fcte1.f

Go to the documentation of this file.

!                 *******************************
                  DOUBLE PRECISION FUNCTION fcte1
!                 *******************************
!
     &(xx)
!
!***********************************************************************
! ARTEMIS   V6P1                                   31/05/2011
!***********************************************************************
!
!brief    EVALUATE FUNCTION E1(KH) FOR SECOND
!+        ORDER BOTTOM EFFECTS (GRADIENT)
!
!
!history  C.PEYRARD
!+        31/05/2011
!+        V6P1
!+
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| XX             |-->| K*H
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE interface_artemis, ex_fcte1 => fcte1
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      DOUBLE PRECISION, INTENT(IN) :: XX
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      DOUBLE PRECISION SHHX , SHH2X , CHHX ,XXX
      INTRINSIC          sinh, cosh
!
      xxx=2.d0*xx
!
      IF(xxx.LT.0.001d0) THEN
        fcte1=-1.d0/6.d0
      ELSE
        shhx  = sinh(xxx)
        shh2x = sinh(2.d0*xxx)
        chhx  = cosh(xxx)
        fcte1=3.d0*( (chhx-1.d0)**2 + 2.d0 )*(xxx + 2.d0*shhx)*xxx
        fcte1=fcte1 +  shhx*(xxx**4/shhx +4.d0*xxx**3 -9.d0*shh2x)
        fcte1=fcte1/(3.d0*(xxx+shhx)**4)
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END