source: LMDZ.3.3/branches/rel-LF/libf/phylmd/FCTTRE.h @ 79

Last change on this file since 79 was 2, checked in by lmdz, 25 years ago

Initial revision

  • 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 
1C     ------------------------------------------------------------------
2C     This COMDECK includes the Thermodynamical functions for the cy39
3C       ECMWF Physics package.
4C       Consistent with YOMCST Basic physics constants, assuming the
5C       partial pressure of water vapour is given by a first order
6C       Taylor expansion of Qs(T) w.r.t. to Temperature, using constants
7C       in YOETHF
8C     ------------------------------------------------------------------
9      REAL PTARG, PDELARG, P5ARG, PQSARG, PCOARG
10      REAL FOEEW, FOEDE, qsats, qsatl, dqsats, dqsatl
11      LOGICAL thermcep
12      PARAMETER (thermcep=.TRUE.)
13C
14      FOEEW ( PTARG,PDELARG ) = EXP (
15     S          (R3LES*(1.-PDELARG)+R3IES*PDELARG) * (PTARG-RTT)
16     S / (PTARG-(R4LES*(1.-PDELARG)+R4IES*PDELARG)) )
17C
18      FOEDE ( PTARG,PDELARG,P5ARG,PQSARG,PCOARG ) = PQSARG*PCOARG*P5ARG
19     S / (PTARG-(R4LES*(1.-PDELARG)+R4IES*PDELARG))**2
20c
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))
29c
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.