Changeset 4146 for trunk/LMDZ.GENERIC/libf/phygeneric/rain.F90

Timestamp:

Mar 19, 2026, 2:35:46 PM (3 weeks ago)

Author:

gmilcareck

Message:

Thermodynamics update on LMDZ.GENERIC

File:

• trunk/LMDZ.GENERIC/libf/phygeneric/rain.F90 (modified) (3 diffs)

trunk/LMDZ.GENERIC/libf/phygeneric/rain.F90

@@ -6 +6 @@
   use aerosol_radius, only: aerosol_radius_h2o_liquid_ice_separate
   USE tracer_h, only: igcm_h2o_vap, igcm_h2o_ice
-  use comcstfi_mod, only: g, r
+  use comcstfi_mod, only: g, rd_ref
   implicit none
 
@@ -215 +215 @@
                      zqev = MAX (0.0, (zqs(i,k)-q(i,k)))*dmass(i,k)/ptimestep !there was a bug here
                      zqevt= 2.0e-5*(1.0-q(i,k)/zqs(i,k))    & !default was 2.e-5
-                        *sqrt(precip_rate(i))*dmass(i,k)/pplay(i,k)*zt(i,k)*R ! BC modif here
+                        *sqrt(precip_rate(i))*dmass(i,k)/pplay(i,k)*zt(i,k)*rd_ref ! BC modif here
                      zqevt = MAX (zqevt, 0.0)
                      zqev  = MIN (zqev, zqevt)
@@ -268 +268 @@
                IF (rneb(i,k).GT.0.0) THEN
                   zoliq(i) = ql(i,k)
-                  zrho(i)  = pplay(i,k) / ( zt(i,k) * R )
+                  zrho(i)  = pplay(i,k) / ( zt(i,k) * rd_ref )
                   zdz(i)   = (pplev(i,k)-pplev(i,k+1)) / (zrho(i)*g)
                   zfrac(i) = (zt(i,k)-T_h2O_ice_clouds) / (T_h2O_ice_liq-T_h2O_ice_clouds)