source: LMDZ5/branches/LMDZ6_rc0/libf/phylmd/clift.F90 @ 3069

Last change on this file since 3069 was 1999, checked in by Laurent Fairhead, 11 years ago

Merged trunk changes r1920:1997 into testing branch

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