source: LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/clift.F90 @ 2404

Last change on this file since 2404 was 1992, checked in by lguez, 11 years ago

Converted to free source form files in libf/phylmd which were still in
fixed source form. The conversion was done using the polish mode of
the NAG Fortran Compiler.

In addition to converting to free source form, the processing of the
files also:

-- indented the code (including comments);

-- set Fortran keywords to uppercase, and set all other identifiers
to lower case;

-- added qualifiers to end statements (for example "end subroutine
conflx", instead of "end");

-- changed the terminating statements of all DO loops so that each
loop ends with an ENDDO statement (instead of a labeled continue).

-- replaced #include by include.
  • 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.