source: LMDZ4/branches/unlabeled-1.1.1/libf/phylmd/FCTTRE.inc @ 1329

Last change on this file since 1329 was 524, checked in by lmdzadmin, 20 years ago

Initial revision

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