zj.f

来自「xy时域有限差分的fortran数值模拟 xy时域有限差分的fortran数值」· F 代码 · 共 15 行

*--------------------------------------------------------------------*
      real function zj(i,j)
*
*  current source at i=imid, j=jmid
*
      include 'xfdtmcom.f'
      integer i,j
      if (wave.eq.'sine') then
         zj = swave(i,j)
      else
         zj = pulse(i,j)
      endif
      return
      end