total1.f90

fortran source code for fdtd parall computing

!*************************************************
!建模
!*************************************************


!c=====================================================
Program main
real ob(-700:700,-700:700)              !存放各网络介质的编号
real Media(1:4,1:15)                    !各介质的电磁参数
character *10 FileName                  !目标信息输出到FileName文件中
real WaveLength,WL                      !入射波长和每波长网格数
real Ismin,Ismax,Jsmin,Jsmax            !目标区域界限
real MediaNo                            !目标区域所有介质数目
integer Record,i,j

!c======================Initiate===============================
FileName='rect.sct'
WaveLength=1
WL=40.
Ismin=-60
Ismax=60
Jsmin=-60
Jsmax=60

MediaNo=1                                 !本例中只有金属一种介质，因此为1

Media(1,1)=1.                             ! epsilon_r
Media(2,1)=1.                             ! mu_r
Media(3,1)=3.72e+7                        ! sigma
Media(4,1)=0.                             ! sigma_m

!====================建立散射体================================
 
 do i=Ismin,Ismax                         !空间所有网格初始值设为背景
   do j=Jsmin,Jsmax
     Ob(i,j)=0.0
   end do
 end do
 
 do i=Ismin,Ismax
    do j=Jsmin,Jsmax
	if((i.ge.-WL).and.(i.le.WL-1).and.(j.ge.-WL).and.(j.le.WL-1)) then    !根据(F-2-1)式判断网格是否在目标上
      Ob(i,j)=1.0
	 end if
   end do
end do 

!c=====================输出======================================
open(1,File=FileName,access='direct',recl=4)
write(1,rec=1)WaveLength
write(1,rec=2)WL
write(1,rec=3)Ismin
write(1,rec=4)Ismax
write(1,rec=5)Jsmin
write(1,rec=6)Jsmax
write(1,rec=7)MediaNo
Record=7
do j=1,MediaNo
  do i=1,4
     Record=Record+1
	 write(1,rec=Record)Media(i,j)
  end do
end do

do i=Ismin,Ismax
  do j=Jsmin,Jsmax
     Record=Record+1
	 write(1,rec=Record) ob(i,j)
  end do 
end do
close(1)

write(*,*) Ismin,Ismax,Jsmin,Jsmax

end

!==================END OF MAIN PROGRAM===============================