source: LMDZ4/branches/V3_test/libf/phylmd/FCTTRE.h @ 735

Last change on this file since 735 was 703, checked in by (none), 19 years ago

This commit was manufactured by cvs2svn to create branch 'V3_test'.

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