source: LMDZ.3.3/branches/rel-LF/libf/phylmd/clift.F @ 275

Last change on this file since 275 was 230, checked in by lmdzadmin, 23 years ago

Merge de la physique avec la branche principale
LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 2.0 KB
Line 
1c $Header$
2        SUBROUTINE CLIFT (P,T,RR,RS,PLCL,DPLCLDT,DPLCLDQ)
3C***************************************************************
4C*                                                             *
5C* CLIFT : COMPUTE LIFTING CONDENSATION LEVEL AND ITS          *
6C*         DERIVATIVES RELATIVE TO T AND R                     *
7C*   (WITHIN 0.2% OF FORMULA OF BOLTON, MON. WEA. REV.,1980)   *
8C*                                                             *
9C* written by   : GRANDPEIX Jean-Yves, 17/11/98, 12.39.01      *
10C* modified by :                                               *
11C***************************************************************
12C*
13C*Arguments :
14C*
15C* Input :  P = pressure of level from wich lifting is performed
16C*          T = temperature of level P
17C*          RR = vapour mixing ratio at level P
18C*          RS = vapour saturation mixing ratio at level P
19C*
20C* Output : PLCL = lifting condensation level
21C*          DPLCLDT = derivative of PLCL relative to T
22C*          DPLCLDQ = derivative of PLCL relative to R
23C*
24ccccccccccccccccccccccc
25c constantes coherentes avec le modele du Centre Europeen
26c      RD = 1000.0 * 1.380658E-23 * 6.0221367E+23 / 28.9644
27c      RV = 1000.0 * 1.380658E-23 * 6.0221367E+23 / 18.0153
28c      CPD = 3.5 * RD
29c      CPV = 4.0 * RV
30c      CL = 4218.0
31c      CI=2090.0
32c      CPVMCL=CL-CPV
33c      CLMCI=CL-CI
34c      EPS=RD/RV
35c      ALV0=2.5008E+06
36c      ALF0=3.34E+05
37c
38c on utilise les constantes thermo du Centre Europeen: (sb)
39c
40#include "YOMCST.h"
41c
42       CPD = RCPD
43       CPV = RCPV
44       CL = RCW
45       CPVMCL = CL-CPV
46       EPS = RD/RV
47       ALV0 = RLVTT
48c
49c
50c      Bolton formula coefficients :
51      A = 1669.0
52      B = 122.0
53c
54      RH=RR/RS
55      CHI=T/(A-B*RH-T)
56      PLCL=P*(RH**CHI)
57c
58      ALV = ALV0 - CPVMCL*(T-273.15)
59c
60c -- sb: correction:
61c       DPLCLDQ = PLCL*CHI*( 1./RR - B*CHI/T/RS*ALOG(RH) )
62      DPLCLDQ = PLCL*CHI*( 1./RR + B*CHI/T/RS*ALOG(RH) )
63c sb --
64c
65      DPLCLDT = PLCL*CHI*((A-B*RH*(1.+ALV/RV/T))/T**2*CHI*ALOG(RH)
66     $                    - ALV/RV/T**2 )
67c
68c
69      RETURN
70      END
Note: See TracBrowser for help on using the repository browser.