source: lmdz_wrf/trunk/WRFV3/lmdz/FCTTRE.h @ 1425

Last change on this file since 1425 was 1, checked in by lfita, 10 years ago
  • -- --- Opening of the WRF+LMDZ coupling repository --- -- -

WRF: version v3.3
LMDZ: version v1818

More details in:

File size: 2.3 KB
Line 
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))
Note: See TracBrowser for help on using the repository browser.