minmod.for

2D有限差分法程序,未加边界条件,适合一般情况下的波动方程模拟

!     copyright:wangaifeng 
!     DATE     :2007-07-25
!     Company  :GUCAS
      real function minmod(x,y)
      implicit real(a-h,o-z)
      if(x*y.le.0.) then
        minmod = 0.
      elseif(abs(x).ge.abs(y)) then
        minmod = y
      else
        minmod = x
      endif
       
      return
      end