fcte1.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\artemis\fcte1.f
00002 !
00039                   DOUBLE PRECISION FUNCTION FCTE1
00040 !                 *******************************
00041 !
00042      &(XX)
00043 !
00044 !***********************************************************************
00045 ! ARTEMIS   V6P1                                   31/05/2011
00046 !***********************************************************************
00047 !
00048 !
00049 !
00050 !
00051 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00052 !| XX             |-->| K*H
00053 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00054 !
00055       USE INTERFACE_ARTEMIS, EX_FCTE1 => FCTE1
00056       IMPLICIT NONE
00057 !
00058 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00059 !
00060       DOUBLE PRECISION, INTENT(IN) :: XX
00061 !
00062 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00063 !
00064       DOUBLE PRECISION SHHX , SHH2X , CHHX ,XXX
00065       INTRINSIC          SINH, COSH
00066 !
00067       XXX=2.D0*XX
00068 !
00069       IF(XXX.LT.0.001D0) THEN
00070         FCTE1=-1.D0/6.D0
00071       ELSE
00072         SHHX  = SINH(XXX)
00073         SHH2X = SINH(2.D0*XXX)
00074         CHHX  = COSH(XXX)
00075         FCTE1=3.D0*( (CHHX-1.D0)**2 + 2.D0 )*(XXX + 2.D0*SHHX)*XXX
00076         FCTE1=FCTE1 +  SHHX*(XXX**4/SHHX +4.D0*XXX**3 -9.D0*SHH2X)
00077         FCTE1=FCTE1/(3.D0*(XXX+SHHX)**4)
00078       ENDIF
00079 !
00080 !-----------------------------------------------------------------------
00081 !
00082       RETURN
00083       END

Generated on Fri Aug 31 2013 18:12:58 by S.E.Bourban (HRW) using doxygen 1.7.0