source: LMDZ4/trunk/libf/phylmd/clift.F @ 3193

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