libf/phylmd/clift.f

!
! $Header: /home/cvsroot/LMDZ4/libf/phylmd/clift.F,v 1.1.1.1 2004/05/19 12:53:08 lmdzadmin Exp $
!
        SUBROUTINE CLIFT (P,T,RR,RS,PLCL,DPLCLDT,DPLCLDQ)
C***************************************************************
C*                                                             *
C* CLIFT : COMPUTE LIFTING CONDENSATION LEVEL AND ITS          *
C*         DERIVATIVES RELATIVE TO T AND R                     *
C*   (WITHIN 0.2% OF FORMULA OF BOLTON, MON. WEA. REV.,1980)   *
C*                                                             *
C* written by   : GRANDPEIX Jean-Yves, 17/11/98, 12.39.01      *
C* modified by :                                               *
C***************************************************************
C*
C*Arguments :
C*
C* Input :  P = pressure of level from wich lifting is performed
C*          T = temperature of level P
C*          RR = vapour mixing ratio at level P
C*          RS = vapour saturation mixing ratio at level P
C*
C* Output : PLCL = lifting condensation level
C*          DPLCLDT = derivative of PLCL relative to T
C*          DPLCLDQ = derivative of PLCL relative to R
C*
ccccccccccccccccccccccc
c constantes coherentes avec le modele du Centre Europeen
c      RD = 1000.0 * 1.380658E-23 * 6.0221367E+23 / 28.9644
c      RV = 1000.0 * 1.380658E-23 * 6.0221367E+23 / 18.0153
c      CPD = 3.5 * RD
c      CPV = 4.0 * RV
c      CL = 4218.0
c      CI=2090.0
c      CPVMCL=CL-CPV
c      CLMCI=CL-CI
c      EPS=RD/RV
c      ALV0=2.5008E+06
c      ALF0=3.34E+05
c
c on utilise les constantes thermo du Centre Europeen: (sb)
c
      use YOMCST
c
       CPD = RCPD
       CPV = RCPV
       CL = RCW
       CPVMCL = CL-CPV
       EPS = RD/RV
       ALV0 = RLVTT
c
c
c      Bolton formula coefficients :
      A = 1669.0
      B = 122.0
c
      RH=RR/RS
      CHI=T/(A-B*RH-T)
      PLCL=P*(RH**CHI)
c
      ALV = ALV0 - CPVMCL*(T-273.15)
c
c -- sb: correction:
c       DPLCLDQ = PLCL*CHI*( 1./RR - B*CHI/T/RS*ALOG(RH) )
      DPLCLDQ = PLCL*CHI*( 1./RR + B*CHI/T/RS*ALOG(RH) )
c sb --
c
      DPLCLDT = PLCL*CHI*((A-B*RH*(1.+ALV/RV/T))/T**2*CHI*ALOG(RH)
     $                    - ALV/RV/T**2 )
c
c
      RETURN
      END
1	!
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	use YOMCST
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 = PLCLCHI( 1./RR - BCHI/T/RSALOG(RH) )
64	DPLCLDQ = PLCLCHI( 1./RR + BCHI/T/RSALOG(RH) )
65	c sb --
66	c
67	DPLCLDT = PLCLCHI((A-BRH(1.+ALV/RV/T))/T*2CHI*ALOG(RH)
68	$ - ALV/RV/T**2 )
69	c
70	c
71	RETURN
72	END