Esope source files

   1 : C FAU_M4    SOURCE    BECC      10/05/05    21:15:05     6674           
   2 :       function fau_m4(M,coef)
   3 : C     
   4 :       real*8 M, coef
   5 :      &     , fau_m1, fau_m2, fau_m4, beta
   6 :      &     , f1, f2, f2m
   7 :       parameter(beta = 1.0D0/8.0D0)
   8 : C
   9 :       f1 = fau_m1(M,coef)
  10 :       f2 = fau_m2(M,coef)
  11 :       f2m = fau_m2(M,(-1.0D0 * coef))
  12 :       if (abs(M) .gt. 1.0D0) then
  13 :          fau_m4 = f1
  14 :       else
  15 :          fau_m4 = 1.0D0 - (coef * 16.0D0 * beta * f2m)
  16 :          fau_m4 = f2 * fau_m4
  17 :       endif
  18 : C
  19 :       return
  20 :       end
  21 : 
  22 :          
  23 :