source: LMDZ6/branches/SETHET_DECOUPLE/libf/phylmd/FCTTRE.h @ 5182

Last change on this file since 5182 was 1907, checked in by lguez, 11 years ago

Added a copyright property to every file of the distribution, except
for the fcm files (which have their own copyright). Use svn propget on
a file to see the copyright. For instance:

$ svn propget copyright libf/phylmd/physiq.F90
Name of program: LMDZ
Creation date: 1984
Version: LMDZ5
License: CeCILL version 2
Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
See the license file in the root directory

Also added the files defining the CeCILL version 2 license, in French
and English, at the top of the LMDZ tree.
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 2.3 KB
Line 
1!
2! $Header$
3!
4!
5!  ATTENTION!!!!: ce fichier include est compatible format fixe/format libre
6!                 veillez  n'utiliser que des ! pour les commentaires
7!                 et  bien positionner les & des lignes de continuation
8!                 (les placer en colonne 6 et en colonne 73)
9!
10!     ------------------------------------------------------------------
11!     This COMDECK includes the Thermodynamical functions for the cy39
12!       ECMWF Physics package.
13!       Consistent with YOMCST Basic physics constants, assuming the
14!       partial pressure of water vapour is given by a first order
15!       Taylor expansion of Qs(T) w.r.t. to Temperature, using constants
16!       in YOETHF
17!     ------------------------------------------------------------------
18      REAL PTARG, PDELARG, P5ARG, PQSARG, PCOARG
19      REAL FOEEW, FOEDE, qsats, qsatl, dqsats, dqsatl
20      LOGICAL thermcep
21      PARAMETER (thermcep=.TRUE.)
22!
23      FOEEW ( PTARG,PDELARG ) = EXP (                                   &
24     &          (R3LES*(1.-PDELARG)+R3IES*PDELARG) * (PTARG-RTT)        &
25     & / (PTARG-(R4LES*(1.-PDELARG)+R4IES*PDELARG)) )
26!
27      FOEDE ( PTARG,PDELARG,P5ARG,PQSARG,PCOARG ) = PQSARG*PCOARG*P5ARG &
28     & / (PTARG-(R4LES*(1.-PDELARG)+R4IES*PDELARG))**2
29!
30      qsats(ptarg) = 100.0 * 0.622 * 10.0                               &
31     &           ** (2.07023 - 0.00320991 * ptarg                       &
32     &           - 2484.896 / ptarg + 3.56654 * LOG10(ptarg))
33      qsatl(ptarg) = 100.0 * 0.622 * 10.0                               &
34     &           ** (23.8319 - 2948.964 / ptarg                         &
35     &           - 5.028 * LOG10(ptarg)                                 &
36     &           - 29810.16 * EXP( - 0.0699382 * ptarg)                 &
37     &           + 25.21935 * EXP( - 2999.924 / ptarg))
38!
39      dqsats(ptarg,pqsarg) = RLVTT/RCPD*pqsarg * (3.56654/ptarg         &
40     &                     +2484.896*LOG(10.)/ptarg**2                  &
41     &                     -0.00320991*LOG(10.))
42      dqsatl(ptarg,pqsarg) = RLVTT/RCPD*pqsarg*LOG(10.)*                &
43     &                (2948.964/ptarg**2-5.028/LOG(10.)/ptarg           &
44     &                +25.21935*2999.924/ptarg**2*EXP(-2999.924/ptarg)  &
45     &                +29810.16*0.0699382*EXP(-0.0699382*ptarg))
Note: See TracBrowser for help on using the repository browser.