source: LMDZ4/branches/V3_test/libf/phylmd/FCTTRE.inc

Last change on this file was 541, checked in by lmdzadmin, 20 years ago

Convergence avec la version d'Olivia Coindreau incluant:

  • le offline
  • les thermiques
  • mellor & yamada dans la couche limite

LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 2.0 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-TEST_RTT) &
19     &          (R3LES*(1.-PDELARG)+R3IES*PDELARG) * (PTARG-273.16) &
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!
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 &
36     &                     +2484.896*LOG(10.)/ptarg**2 &
37     &                     -0.00320991*LOG(10.))
38!     dqsatl(ptarg,pqsarg) = TEST_RLVTT/TEST_RCPD*pqsarg*LOG(10.)* &
39      dqsatl(ptarg,pqsarg) = 0.2500800E+07/0.1004709E+04*pqsarg*LOG(10.)* &
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.