source: LMDZ.3.3/branches/rel-LF/libf/phylmd/FCTTRE.inc @ 456

Last change on this file since 456 was 97, checked in by (none), 24 years ago

This commit was manufactured by cvs2svn to create branch 'rel-LF'.

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 1.7 KB
Line 
1!      ------------------------------------------------------------------
2!      This COMDECK includes the Thermodynamical functions for the cy39
3!       ECMWF Physics package.
4!       Consistent with YOMCST Basic physics constants, assuming the
5!       partial pressure of water vapour is given by a first order
6!       Taylor expansion of Qs(T) w.r.t. to Temperature, using constants
7!       in YOETHF
8!     ------------------------------------------------------------------
9      REAL PTARG, PDELARG, P5ARG, PQSARG, PCOARG
10      REAL FOEEW, FOEDE, qsats, qsatl, dqsats, dqsatl
11      LOGICAL thermcep
12      PARAMETER (thermcep=.TRUE.)
13!
14      FOEEW ( PTARG,PDELARG ) = EXP ( &
15     &          (R3LES*(1.-PDELARG)+R3IES*PDELARG) * (PTARG-RTT) &
16     & / (PTARG-(R4LES*(1.-PDELARG)+R4IES*PDELARG)) )
17!
18      FOEDE ( PTARG,PDELARG,P5ARG,PQSARG,PCOARG ) = PQSARG*PCOARG*P5ARG &
19     & / (PTARG-(R4LES*(1.-PDELARG)+R4IES*PDELARG))**2
20!
21      qsats(ptarg) = 100.0 * 0.622 * 10.0 &
22     &           ** (2.07023 - 0.00320991 * ptarg &
23     &           - 2484.896 / ptarg + 3.56654 * LOG10(ptarg))
24      qsatl(ptarg) = 100.0 * 0.622 * 10.0 &
25     &           ** (23.8319 - 2948.964 / ptarg &
26     &           - 5.028 * LOG10(ptarg) &
27     &           - 29810.16 * EXP( - 0.0699382 * ptarg) &
28     &           + 25.21935 * EXP( - 2999.924 / ptarg))
29!
30      dqsats(ptarg,pqsarg) = RLVTT/RCPD*pqsarg * (3.56654/ptarg &
31     &                     +2484.896*LOG(10.)/ptarg**2 &
32     &                     -0.00320991*LOG(10.))
33      dqsatl(ptarg,pqsarg) = RLVTT/RCPD*pqsarg*LOG(10.)* &
34     &                (2948.964/ptarg**2-5.028/LOG(10.)/ptarg &
35     &                +25.21935*2999.924/ptarg**2*EXP(-2999.924/ptarg) &
36     &                +29810.16*0.0699382*EXP(-0.0699382*ptarg))
Note: See TracBrowser for help on using the repository browser.