Changeset 5144 for LMDZ6/branches/Amaury_dev/libf/phylmdiso/clift.F90

Timestamp:

Jul 29, 2024, 11:01:04 PM (8 weeks ago)

Author:

abarral

Message:

Put YOMCST.h into modules

File:

• LMDZ6/branches/Amaury_dev/libf/phylmdiso/clift.F90 (modified) (4 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmdiso/clift.F90

@@ -1 @@
-
 ! $Header$
 
 SUBROUTINE clift(p, t, rr, rs, plcl, dplcldt, dplcldq)
-IMPLICIT NONE
+  USE lmdz_yomcst
+
+  IMPLICIT NONE
   ! ***************************************************************
   ! *                                                             *
@@ -41 +42 @@
   ! on utilise les constantes thermo du Centre Europeen: (sb)
 
-  include "YOMCST.h"
-  REAL :: p,t,rr,rs,plcl,dplcldt,dplcldq,cpd,cpv,cl,cpvmcl,eps,alv0,a,b
-  REAL :: rh,chi,alv
+  REAL :: p, t, rr, rs, plcl, dplcldt, dplcldq, cpd, cpv, cl, cpvmcl, eps, alv0, a, b
+  REAL :: rh, chi, alv
 
   cpd = rcpd
@@ -49 +49 @@
   cl = rcw
   cpvmcl = cl - cpv
-  eps = rd/rv
+  eps = rd / rv
   alv0 = rlvtt
 
@@ -57 +57 @@
   b = 122.0
 
-  rh = rr/rs
-  chi = t/(a-b*rh-t)
-  plcl = p*(rh**chi)
+  rh = rr / rs
+  chi = t / (a - b * rh - t)
+  plcl = p * (rh**chi)
 
-  alv = alv0 - cpvmcl*(t-273.15)
+  alv = alv0 - cpvmcl * (t - 273.15)
 
   ! -- sb: correction:
   ! DPLCLDQ = PLCL*CHI*( 1./RR - B*CHI/T/RS*ALOG(RH) )
-  dplcldq = plcl*chi*(1./rr+b*chi/t/rs*alog(rh))
+  dplcldq = plcl * chi * (1. / rr + b * chi / t / rs * alog(rh))
   ! sb --
 
-  dplcldt = plcl*chi*((a-b*rh*(1.+alv/rv/t))/t**2*chi*alog(rh)-alv/rv/t**2)
-
-
+  dplcldt = plcl * chi * ((a - b * rh * (1. + alv / rv / t)) / t**2 * chi * alog(rh) - alv / rv / t**2)
 
 END SUBROUTINE clift