A simulation model of the human body is developed in frequency dependent Crank Nicolson finite difference time domain (FD-CN-FDTD) method. Numerical simulation of electromagnetic wave propagation inside the human head is presented. Advantages of using time discretization beyond the Courant Friedrich-Lewy (CFL) limit in FD-CN-FDTD method are shown. Parallelization using Open Multi-Processing (OpenMP) in a shared memory architecture is performed and the achieved computational efficiencies are shown. © 2011 VSP.