zj.f

a 2-d FDTD program using PML

*--------------------------------------------------------------------*
      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