constants.f90

来自「三维FDTD,matlal编程,mur边界条件,平面波光源」· F90 代码 · 共 31 行

!===========================================================
!           set common parameters 
!=========================================================== 

module constants
implicit none

	real(8), parameter ::	pi=3.14159265358979d0,&
						nm=1d-9,&           
						e0=(1d-9)/(36.0d0*pi),&  
						u0=pi*4d-7,&         
						g0=0d0,&
						c=(u0*e0)**(-0.5d0)
	integer,parameter :: kind=8

	complex(kind),allocatable,dimension(:,:,:),save ::	Ex,Ey,Ez,&
									Exy_old,Exz_old,Eyx_old,Eyz_old,Ezx_old,Ezy_old,&  ! field value
													Hx,Hy,Hz,&
													CEx,DEx,CEy,DEy,CEz,DEz  ! material distribution	
	real(4),allocatable,dimension(:,:,:),save ::		outdata    ! output data 
	complex(kind),allocatable,dimension(:),save ::	c_p,d_p,& ! material 
												e_p,& ! dielectric constant
												u_p,&
												g_p
	integer(kind),save :: nx,ny,nz,nt,mat_num !numbers of material
	 
	real(kind),save  :: lamda,& ! wavelength
					w,&		! wave frequency
					ds,dx,dy,dz,dt ! space and time step
     
end module