Changeset 3496 for trunk/LMDZ.TITAN/libf/muphytitan/mm_clouds.f90

Timestamp:

Nov 7, 2024, 11:09:51 AM (2 weeks ago)

Author:

debatzbr

Message:

Final clean and debug for the microphysical model

File:

• trunk/LMDZ.TITAN/libf/muphytitan/mm_clouds.f90 (modified) (3 diffs)

trunk/LMDZ.TITAN/libf/muphytitan/mm_clouds.f90

@@ -112 +112 @@
       mm_ccn_vsed(:) = wsettle(mm_play,mm_temp,mm_zlay,mm_drho,mm_drad)
 
-      ! Computes flux [kg.m-2.s-1] and precipitation [kg.m-2.iphysiq] of ccn
+      ! Computes flux [kg.m-2.s-1] and precipitation [kg.m-2] of ccn
       mm_ccn_flux(:) = get_mass_flux(mm_rhoaer,mm_m3ccn(:))
       mm_ccn_prec = SUM(zdm3n*mm_dzlev*mm_rhoaer)
 
-      ! Computes flux [kg.m-2.s-1] and precipitation [kg.m-2.iphysiq] of ices
+      ! Computes flux [kg.m-2.s-1] and precipitation [kg.m-2] of ices
       DO i = 1, mm_nesp
         mm_ice_fluxes(:,i) = get_mass_flux(mm_xESPS(i)%rho,(3._mm_wp*mm_m3ice(:,i))/(4._mm_wp*mm_pi))
@@ -729 +729 @@
     REAL(kind=mm_wp), INTENT(in) :: rad !! Radius of the particle (m).
     REAL(kind=mm_wp) :: w               !! Settling velocity (\(m.s^{-1}\)).
-    REAL(kind=mm_wp) :: Us, Fc, kn
+    REAL(kind=mm_wp) :: Us, Fc, kn, wtmp, wmax
     REAL(kind=mm_wp), PARAMETER :: ra = 1.75e-10_mm_wp
     
@@ -743 +743 @@
     ! Computes settling velocity (correction factor : x3.0)
     w = Us * Fc * 3._mm_wp
+
+    ! Imposes a velocity limit
+    wmax = 20._mm_wp ! 20 m/s [Lorenz 1993]
+    wtmp = (1._mm_wp / w) + (1._mm_wp / wmax)
+    w = 1._mm_wp / wtmp
   END FUNCTION wsettle