/[lmdze]/trunk/libf/phylmd/clift.f
ViewVC logotype

Annotation of /trunk/libf/phylmd/clift.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 38 - (hide annotations)
Thu Jan 6 17:52:19 2011 UTC (13 years, 4 months ago) by guez
File size: 2180 byte(s) 
Extracted ASCII art from "inigeom" into a separate text file in the
documentation.

"test_disvert" now creates a separate file for layer thicknesses.

Moved variables from module "yomcst" to module "suphec_m" because this
is where those variables are defined. Kept in "yomcst" only parameters
of Earth orbit. Gave the attribute "parameter" to some variables of
module "suphec_m".

Variables of module "yoethf" were defined in procedure "suphec". Moved
these definitions to a new procedure "yoethf" in module "yoethf_m".
1 guez 3 !
2     ! $Header: /home/cvsroot/LMDZ4/libf/phylmd/clift.F,v 1.1.1.1 2004/05/19 12:53:08 lmdzadmin Exp $
3     !
4     SUBROUTINE CLIFT (P,T,RR,RS,PLCL,DPLCLDT,DPLCLDQ)
5     C***************************************************************
6     C* *
7     C* CLIFT : COMPUTE LIFTING CONDENSATION LEVEL AND ITS *
8     C* DERIVATIVES RELATIVE TO T AND R *
9     C* (WITHIN 0.2% OF FORMULA OF BOLTON, MON. WEA. REV.,1980) *
10     C* *
11     C* written by : GRANDPEIX Jean-Yves, 17/11/98, 12.39.01 *
12     C* modified by : *
13     C***************************************************************
14     C*
15     C*Arguments :
16     C*
17     C* Input : P = pressure of level from wich lifting is performed
18     C* T = temperature of level P
19     C* RR = vapour mixing ratio at level P
20     C* RS = vapour saturation mixing ratio at level P
21     C*
22     C* Output : PLCL = lifting condensation level
23     C* DPLCLDT = derivative of PLCL relative to T
24     C* DPLCLDQ = derivative of PLCL relative to R
25     C*
26     ccccccccccccccccccccccc
27     c constantes coherentes avec le modele du Centre Europeen
28     c RD = 1000.0 * 1.380658E-23 * 6.0221367E+23 / 28.9644
29     c RV = 1000.0 * 1.380658E-23 * 6.0221367E+23 / 18.0153
30     c CPD = 3.5 * RD
31     c CPV = 4.0 * RV
32     c CL = 4218.0
33     c CI=2090.0
34     c CPVMCL=CL-CPV
35     c CLMCI=CL-CI
36     c EPS=RD/RV
37     c ALV0=2.5008E+06
38     c ALF0=3.34E+05
39     c
40     c on utilise les constantes thermo du Centre Europeen: (sb)
41     c
42 guez 38 use SUPHEC_M
43 guez 3 c
44     CPD = RCPD
45     CPV = RCPV
46     CL = RCW
47     CPVMCL = CL-CPV
48     EPS = RD/RV
49     ALV0 = RLVTT
50     c
51     c
52     c Bolton formula coefficients :
53     A = 1669.0
54     B = 122.0
55     c
56     RH=RR/RS
57     CHI=T/(A-B*RH-T)
58     PLCL=P*(RH**CHI)
59     c
60     ALV = ALV0 - CPVMCL*(T-273.15)
61     c
62     c -- sb: correction:
63     c DPLCLDQ = PLCL*CHI*( 1./RR - B*CHI/T/RS*ALOG(RH) )
64     DPLCLDQ = PLCL*CHI*( 1./RR + B*CHI/T/RS*ALOG(RH) )
65     c sb --
66     c
67     DPLCLDT = PLCL*CHI*((A-B*RH*(1.+ALV/RV/T))/T**2*CHI*ALOG(RH)
68     $ - ALV/RV/T**2 )
69     c
70     c
71     RETURN
72     END
  ViewVC Help
Powered by ViewVC 1.1.21