My Project
 All Classes Files Functions Variables Macros
clift.F
Go to the documentation of this file.
1 !
2 ! $Header$
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 #include "YOMCST.h"
43 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