fcte2.f

Go to the documentation of this file.

!                 *******************************
                  DOUBLE PRECISION FUNCTION fcte2
!                 *******************************
!
     &(xx)
!
!***********************************************************************
! ARTEMIS   V6P1                                   31/05/2011
!***********************************************************************
!
!brief    EVALUATE FUNCTION E2(KH) FOR SECOND
!+        ORDER BOTTOM EFFECTS (CURVATURE)
!
!history  C.PEYRARD
!+        31/05/2011
!+        V6P1
!+
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| XX             |-->| K*H
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE interface_artemis, ex_fcte2 => fcte2
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      DOUBLE PRECISION, INTENT(IN) :: XX
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      DOUBLE PRECISION  CHHX,SHHX , THHX,XXX
      INTRINSIC         sinh, cosh
!
      xxx = 2.d0*xx
!
      IF(xxx.LT.0.001d0) THEN
        fcte2=-1.d0/12.d0
      ELSE
        shhx = sinh(xxx)
        chhx = cosh(xxx)
        thhx = shhx/chhx
        fcte2=(thhx-xxx)*chhx/(xxx*(shhx+xxx)**2)
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END