pulse.f

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

*--------------------------------------------------------------------*
      real function pulse(i,j)
      include 'xfdtmcom.f'
      if ( (i.eq.imid).and.(j.eq.jmid) ) then
         arg=t/(20.*dt)
         if(arg.le.1.)then
           pulse = (sin(2.*pi*arg - pi/2) + 1.)/2.
         else
           pulse = 0.
         endif
      else
         pulse=0.
      endif
      return
      end