[资料共享]看看我编的二维FDTD程序

以下是我编写的二维TM波的FDTD 的FORTRAN程序 ;cKH1  
是自由空间中的电磁波的传播模拟，没有其他的限制条件，只是想实现模拟电磁波在二维真空中的传播情况，但这个程序的结果却是不稳定的，就是电场随时间无限增大，我看了下这个程序是符合时间稳定条件的，现在怎么也找不到毛病出在哪里了，请高手帮忙看看问题出在哪里： Oz{%k#X-  
rbyY8 bX  
没有加边界的，但我想只要时间步小点应该还不会发散，但不知道为什么还是不收敛 &[_ZXVva~  
P~RhUKfd  
      program fdtd2D m'x;,xfY&F  
real*8,dimension(:,:),allocatable::hx,hy,ez >w.'KR0L  
real eps0,mu0,dt,dx,dy,c0,pi,f #?Mj$ZB  
integer pmin,pmax k4{:9zL1#?  
pi=3.14159265 y+U83a[L*  
c0=3.0e8  !光速 ,l HLH  
f=3.e12    !源频率 {)@D`{$  
dx=3.e-7  !300纳米 Ei5wel6!  
dy=dx          mS%4gx~~_n  
dt=dx/(2*c0)        lb~E0U`\E`  
eps0=8.854e-12 2z_2.0/3  
mu0=4.0e-7 pu4,0bw  
nn=100  !网格数 p|qyTeg  
pmin=0 >(3\kiYS  
pmax=nn nY_?Jq  
  open(1,file='D:\MSDEV\ex.txt') ^Ga&
}-  
allocate (hx(pmin:pmax,pmin:pmax)) ?6P.b6m}0  
allocate (hy(pmin:pmax,pmin:pmax)) *(QH{!-$s  
allocate (ez(pmin:pmax,pmin:pmax)) sJZ!sznn  
do i=pmin,pmax #NNewzC<*  
  do j=pmin,pmax OBOwz4<  
    hx(i,j)=0. s_ bR]G  
  end do >jx.R  
end do rzeLx Wt  
do i=pmin,pmax /ty?<24ko  
  do j=pmin,pmax @r;wobt  
    hy(i,j)=0. t6lE#<xZV;  
  end do 8=;k"  
end do          WE6\dhJ<  
do i=pmin,pmax 2.lgT|p  
  do j=pmin,pmax  ?=Db@97  
    ez(i,j)=0. #0"~G][#  
  end do Yr+23Ro  
end do qOy=O [+9  
SVWIEH0?  
nt=20 !时间步数 bw4oLu?  
7J*N_8?2  
do n=1,nt 1n=_y o  
do j=pmin+1,pmax Ne6}oQy(S`  
  if(j>=pmin+1) then `bG7"o`  
    do i=pmin+1,pmax z[nS$]u  
          if (i.eq.50.and.j.eq.50) then h7yqk4'Lq  
                ez(i,j)=sin(2*pi*f*n*dt) !激励源        ! Bh]!WMAw.  
          else aFj.i8+  
                ez(i,j)=ez(i,j)+(dt/eps0)*((hy(i,j)-hy(i-1,j))/dx-(hx(i,j)-hx(i,j-1))/dy) 6'C2SihYp  
          end if h|;qG)f^  
    end do m:+8J,jW  
  end if \7v)iG|#G&  
  if(j<=pmax-1) then xJwG=$o  
      do i=pmin,pmax-1 T:iP="?{  
              hx(i,j)=hx(i,j)-(dt/mu0)*(ez(i,j+1)-ez(i,j))/dy $m| V :/  
  &nbs .. aM=D84@  
9X}I>  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册

感谢“憨企鹅”，已经解决了 ;#/@+4@a&  
,e722wz  
真空磁导系数少乘个Pi

有没有c的源程序啊