source: LMDZ4/branches/LMDZ4_par_0/libf/phylmd/FCTTRE.inc @ 2725

Last change on this file since 2725 was 633, checked in by (none), 20 years ago

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

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 2.0 KB
RevLine 
[524]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 ( &
[541]18!    &          (R3LES*(1.-PDELARG)+R3IES*PDELARG) * (PTARG-TEST_RTT) &
19     &          (R3LES*(1.-PDELARG)+R3IES*PDELARG) * (PTARG-273.16) &
[524]20     & / (PTARG-(R4LES*(1.-PDELARG)+R4IES*PDELARG)) )
21!
22      FOEDE ( PTARG,PDELARG,P5ARG,PQSARG,PCOARG ) = PQSARG*PCOARG*P5ARG &
23     & / (PTARG-(R4LES*(1.-PDELARG)+R4IES*PDELARG))**2
24!
25      qsats(ptarg) = 100.0 * 0.622 * 10.0 &
26     &           ** (2.07023 - 0.00320991 * ptarg &
27     &           - 2484.896 / ptarg + 3.56654 * LOG10(ptarg))
28      qsatl(ptarg) = 100.0 * 0.622 * 10.0 &
29     &           ** (23.8319 - 2948.964 / ptarg &
30     &           - 5.028 * LOG10(ptarg) &
31     &           - 29810.16 * EXP( - 0.0699382 * ptarg) &
32     &           + 25.21935 * EXP( - 2999.924 / ptarg))
33!
[541]34!      dqsats(ptarg,pqsarg) = TEST_RLVTT/TEST_RCPD*pqsarg * (3.56654/ptarg &
35      dqsats(ptarg,pqsarg) = 0.2500800E+07/0.1004709E+04*pqsarg * (3.56654/ptarg &
[524]36     &                     +2484.896*LOG(10.)/ptarg**2 &
37     &                     -0.00320991*LOG(10.))
[541]38!     dqsatl(ptarg,pqsarg) = TEST_RLVTT/TEST_RCPD*pqsarg*LOG(10.)* &
39      dqsatl(ptarg,pqsarg) = 0.2500800E+07/0.1004709E+04*pqsarg*LOG(10.)* &
[524]40     &                (2948.964/ptarg**2-5.028/LOG(10.)/ptarg &
41     &                +25.21935*2999.924/ptarg**2*EXP(-2999.924/ptarg) &
42     &                +29810.16*0.0699382*EXP(-0.0699382*ptarg))
Note: See TracBrowser for help on using the repository browser.