My Project
 All Classes Files Functions Variables Macros
FCTTRE.h
Go to the documentation of this file.
1 !
2 ! $Header$
3 !
4 !
5 ! ATTENTION!!!!: ce fichier include est compatible format fixe/format libre
6 ! veillez n'utiliser que des ! pour les commentaires
7 ! et bien positionner les & des lignes de continuation
8 ! (les placer en colonne 6 et en colonne 73)
9 !
10 ! ------------------------------------------------------------------
11 ! This COMDECK includes the Thermodynamical functions for the cy39
12 ! ECMWF Physics package.
13 ! Consistent with YOMCST Basic physics constants, assuming the
14 ! partial pressure of water vapour is given by a first order
15 ! Taylor expansion of Qs(T) w.r.t. to Temperature, using constants
16 ! in YOETHF
17 ! ------------------------------------------------------------------
18  REAL PTARG, PDELARG, P5ARG, PQSARG, PCOARG
19  REAL FOEEW, FOEDE, qsats, qsatl, dqsats, dqsatl
20  LOGICAL thermcep
21  PARAMETER (thermcep=.TRUE.)
22 !
23  FOEEW ( PTARG,PDELARG ) = EXP ( &
24  & (R3LES*(1.-PDELARG)+R3IES*PDELARG) * (PTARG-RTT) &
25  & / (PTARG-(R4LES*(1.-PDELARG)+R4IES*PDELARG)) )
26 !
27  FOEDE ( PTARG,PDELARG,P5ARG,PQSARG,PCOARG ) = PQSARG*PCOARG*P5ARG &
28  & / (PTARG-(R4LES*(1.-PDELARG)+R4IES*PDELARG))**2
29 !
30  qsats(ptarg) = 100.0 * 0.622 * 10.0 &
31  & ** (2.07023 - 0.00320991 * ptarg &
32  & - 2484.896 / ptarg + 3.56654 * LOG10(ptarg))
33  qsatl(ptarg) = 100.0 * 0.622 * 10.0 &
34  & ** (23.8319 - 2948.964 / ptarg &
35  & - 5.028 * LOG10(ptarg) &
36  & - 29810.16 * EXP( - 0.0699382 * ptarg) &
37  & + 25.21935 * EXP( - 2999.924 / ptarg))
38 !
39  dqsats(ptarg,pqsarg) = RLVTT/RCPD*pqsarg * (3.56654/ptarg &
40  & +2484.896*LOG(10.)/ptarg**2 &
41  & -0.00320991*LOG(10.))
42  dqsatl(ptarg,pqsarg) = RLVTT/RCPD*pqsarg*LOG(10.)* &
43  & (2948.964/ptarg**2-5.028/LOG(10.)/ptarg &
44  & +25.21935*2999.924/ptarg**2*EXP(-2999.924/ptarg) &
45  & +29810.16*0.0699382*EXP(-0.0699382*ptarg))