Changeset 3496

Timestamp:

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

Author:

debatzbr

Message:

Final clean and debug for the microphysical model

Location:

trunk/LMDZ.TITAN/libf/muphytitan

Files:

• mm_clouds.f90 (modified) (3 diffs)
• mm_globals.f90 (modified) (3 diffs)
• mm_methods.f90 (modified) (2 diffs)
• mm_microphysic.f90 (modified) (2 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
 

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

@@ -429 +429 @@
   REAL(kind=mm_wp), DIMENSION(:), ALLOCATABLE, SAVE :: mm_drho
 
-  !> Aerosols precipitations (kg.m-2.s-1).
+  !> Aerosols precipitations (kg.m-2).
   !!
   !! Aerosols precipitations take into account both spherical and fractal modes.
@@ -435 +435 @@
   REAL(kind=mm_wp), SAVE :: mm_aer_prec = 0._mm_wp
 
-  !> CCN precipitations (kg.m-2.s-1).
+  !> CCN precipitations (kg.m-2).
   !! It is updated in [[mm_clouds(module):mm_cloud_microphysics(subroutine)]].
   REAL(kind=mm_wp), SAVE :: mm_ccn_prec = 0._mm_wp
@@ -505 +505 @@
   REAL(kind=mm_wp), DIMENSION(:), ALLOCATABLE, SAVE :: mm_ccn_flux
 
-  !> Ice components precipitations (kg.m-2.s-1).
+  !> Ice components precipitations (kg.m-2).
   !!
   !! It is a vector of [[mm_globals(module):mm_nesp(variable)]] values which share the same indexing

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

@@ -370 +370 @@
       res = (1.0e5 / pres) * exp(1.511e1 - 2.207e3/temp - 2.411e4/temp**2 + 7.744e5/temp**3 - 1.161e7/temp**4 + 6.763e7/temp**5)
     
+    ELSE IF(xESP%name == "AC6H6") THEN
+      ! Fray and Schmidt (2009)
+      res = (1.0e5 / pres) * exp(1.735e1 - 5.663e3/temp)
+
     ELSE IF(xESP%name == "HCN") THEN
       ! Fray and Schmidt (2009)
@@ -406 +410 @@
       res = (1.0e5 / pres) * exp(1.511e1 - 2.207e3/temp - 2.411e4/temp**2 + 7.744e5/temp**3 - 1.161e7/temp**4 + 6.763e7/temp**5)
     
+    ELSE IF(xESP%name == "AC6H6") THEN
+      ! Fray and Schmidt (2009)
+      res = (1.0e5 / pres) * exp(1.735e1 - 5.663e3/temp)
+
     ELSE IF(xESP%name == "HCN") THEN
       ! Fray and Schmidt (2009)

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

@@ -186 +186 @@
     !! microphysics is disabled (see [[mm_globals(module):mm_w_clouds(variable)]] documentation).
     REAL(kind=8), INTENT(IN)                                :: dt         !! Physics timestep (s).
-    REAL(kind=mm_wp), INTENT(out), OPTIONAL                 :: aer_prec   !! Aerosols precipitations (both modes) (kg.m-2.s-1).
-    REAL(kind=mm_wp), INTENT(out), OPTIONAL                 :: ccn_prec   !! CCN precipitations (kg.m-2.s-1).
+    REAL(kind=mm_wp), INTENT(out), OPTIONAL                 :: aer_prec   !! Aerosols precipitations (both modes) (\(kg.m^{-2}.s^{-1}\)).
+    REAL(kind=mm_wp), INTENT(out), OPTIONAL                 :: ccn_prec   !! CCN precipitations (\(kg.m^{-2}.s^{-1}\)).
     REAL(kind=mm_wp), INTENT(out), OPTIONAL, DIMENSION(:)   :: aer_s_w    !! Spherical aerosol settling velocity (\(m.s^{-1}\)).
     REAL(kind=mm_wp), INTENT(out), OPTIONAL, DIMENSION(:)   :: aer_f_w    !! Fractal aerosol settling velocity (\(m.s^{-1}\)).
@@ -196 +196 @@
     REAL(kind=mm_wp), INTENT(out), OPTIONAL, DIMENSION(:,:) :: ice_fluxes !! Ice sedimentation fluxes (\(kg.m^{-2}.s^{-1}\)).
     REAL(kind=mm_wp), INTENT(out), OPTIONAL, DIMENSION(:,:) :: gazs_sat   !! Condensible gaz saturation ratios (--).
-    REAL(kind=mm_wp), INTENT(out), OPTIONAL, DIMENSION(:)   :: ice_prec   !! Ice precipitations (kg.m-2.s-1).
+    REAL(kind=mm_wp), INTENT(out), OPTIONAL, DIMENSION(:)   :: ice_prec   !! Ice precipitations (\(kg.m^{-2}.s^{-1}\)).
 
     IF (PRESENT(aer_prec))   aer_prec   = ABS(mm_aer_prec) / dt