Changeset 5626 for LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_condensation.f90

Timestamp:

Apr 23, 2025, 10:51:22 AM (7 weeks ago)

Author:

aborella

Message:

Corrections coupling with convective clouds

File:

• LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_condensation.f90 (modified) (3 diffs)

LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_condensation.f90

@@ -238 +238 @@
 ! for unadjusted clouds
 REAL :: qiceincld, qvapincld, qvapincld_new
+REAL :: qice_ratio
 !
 ! for deposition / sublimation
@@ -538 +539 @@
               !--Exact explicit integration (qvc exact, qice explicit)
               tauinv_depsub = depo_coef_cirrus * qiceincld**(1./3.) * kappa_depsub
+              qvapincld_new = qsat(i) + ( qvapincld - qsat(i) ) * EXP( - dtime * tauinv_depsub )
             ELSE
               !--If the cloud is initially subsaturated
-              !--Exact explicit integration (qvc exact, qice explicit)
-              !--Same but depo_coef_cirrus = 1
+              !!--Exact explicit integration (qvc exact, qice explicit)
+              !!--Same but depo_coef_cirrus = 1
+              !tauinv_depsub = qiceincld**(1./3.) * kappa_depsub
+              !qvapincld_new = qsat(i) + ( qvapincld - qsat(i) ) * EXP( - dtime * tauinv_depsub )
+              !--Exact explicit integration (qice exact, qvc explicit)
+              !--The barrier is set so that the resulting vapor in cloud
+              !--cannot be greater than qsat
+              !--qice_ratio is the ratio between the new ice content and
+              !--the old one, it is comprised between 0 and 1
               tauinv_depsub = qiceincld**(1./3.) * kappa_depsub
+              qice_ratio = tauinv_depsub * dtime / 1.5 / qiceinmix * ( qsat(i) - qvapincld )
+              !--The new vapor in the cloud is increased with the
+              !--sublimated ice
+              qvapincld_new = qvapincld + qiceincld * ( 1. - MAX(0., 1. - qice_ratio)**1.5 )
+              !--The new vapor in the cloud cannot be greater than qsat
+              qvapincld_new = MIN(qvapincld_new, qsat(i))
+              !--If all the ice is sublimated
+              IF ( qvapincld_new .GE. ( qvapincld + qiceincld ) ) qvapincld_new = 0.
             ENDIF ! qvapincld .GT. qsat
-            qvapincld_new = qsat(i) + ( qvapincld - qsat(i) ) * EXP( - dtime * tauinv_depsub )
-            !--If all the ice is sublimated
-            IF ( qvapincld_new .GE. ( qvapincld + qiceincld ) ) qvapincld_new = 0.
           ELSE
             !--We keep the saturation adjustment hypothesis, and the vapor in the
@@ -836 +850 @@
           qiceinmix = ( qcld(i) - qvc(i) ) / cldfra(i) / ( 1. + clrfra_mix / cldfra_mix )
           tauinv_depsub = qiceinmix**(1./3.) * kappa_depsub
-          qvapinmix_lim = qsat(i) - qiceinmix / ( 1. - EXP( - dtime * tauinv_depsub ) )
+          !qvapinmix_lim = qsat(i) - qiceinmix / ( 1. - EXP( - dtime * tauinv_depsub ) )
+          qvapinmix_lim = qsat(i) - qiceinmix * MAX(1., 1.5 / ( dtime * tauinv_depsub ))
           qvapinclr_lim = qvapinmix_lim * ( 1. + cldfra_mix / clrfra_mix ) &
                         - qvc(i) / cldfra(i) * cldfra_mix / clrfra_mix