Sources Esope | Cast3M

calff
C CALFF     SOURCE    CB215821  16/04/21    21:15:28     8920
      SUBROUTINE CALFF(A,FE)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
      IF((1+A).LT.1E-3) THEN
                      FE = -1.5
                 ELSE
         IF ((A-1).GT.-1E-3)  THEN
                      FE = -1.5 + LOG(4.D0)
                 ELSE
                      B =  SQRT(1.-A*A)
                      FE = -1.5 +((1+A)/2)*LOG(2*(1+A))
     *                     + B* (ATAN((1+A)/B) - ATAN(A/B) )
         ENDIF
      ENDIF
      RETURN
      END