Esope source files

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