Esope source files

   1 : C CALFF     SOURCE    CHAT      05/01/12    21:45:56     5004           
   2 :       SUBROUTINE CALFF(A,FE)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8 (A-H,O-Z)
   5 :       IF((1+A).LT.1E-3) THEN
   6 :                       FE = -1.5
   7 :                  ELSE
   8 :          IF ((A-1).GT.-1E-3)  THEN
   9 :                       FE = -1.5 + DLOG(4.D0)
  10 :                  ELSE
  11 :                       B =  DSQRT(1.-A*A)
  12 :                       FE = -1.5 +((1+A)/2)*DLOG(2*(1+A))
  13 :      *                     + B* (DATAN((1+A)/B) - DATAN(A/B) )
  14 :          ENDIF
  15 :       ENDIF
  16 :       RETURN
  17 :       END
  18 :