swave.f

xy时域有限差分的fortran数值模拟 xy时域有限差分的fortran数值模拟

*--------------------------------------------------------------------*
      real function swave(i,j)
*
*  function to generate a sine wave source
*
      include 'xfdtmcom.f'
      if ( (i.eq.imid).and.(j.eq.jmid) ) then
         freq=1./(50.*dt)
         swave = sin(2.*pi*freq*t)
      else
         swave=0.
      endif
      return
      end