Sources Esope | Cast3M

fau_m4
C FAU_M4    SOURCE    BECC      10/05/05    21:15:05     6674
      function fau_m4(M,coef)
C
      real*8 M, coef
     &     , fau_m1, fau_m2, fau_m4, beta
     &     , f1, f2, f2m
      parameter(beta = 1.0D0/8.0D0)
C
      f1 = fau_m1(M,coef)
      f2 = fau_m2(M,coef)
      f2m = fau_m2(M,(-1.0D0 * coef))
      if (abs(M) .gt. 1.0D0) then
         fau_m4 = f1
      else
         fau_m4 = 1.0D0 - (coef * 16.0D0 * beta * f2m)
         fau_m4 = f2 * fau_m4
      endif
C
      return
      end