Changeset 4032 for trunk/LMDZ.PLUTO/libf/muphypluto/mp2m_globals.F90

Timestamp:

Jan 28, 2026, 11:32:11 AM (11 days ago)

Author:

debatzbr

Message:

Pluto PCM: Improve haze microphysics - Updated tracers.
BBT

File:

• trunk/LMDZ.PLUTO/libf/muphypluto/mp2m_globals.F90 (modified) (4 diffs)

trunk/LMDZ.PLUTO/libf/muphypluto/mp2m_globals.F90

@@ -30 +30 @@
     !        - Vertical structure part
     !
-    !     The module also contains twenty methods:
+    !     The module also contains twenty-one methods:
     !        - mm_global_init_0
     !        - mm_global_init_1
@@ -36 +36 @@
     !        - mm_aerosols_init
     !        - mm_clouds_init
-    !        - mm_alpha_s, mm_alpha_f
+    !        - mm_alpha_s, mm_alpha_f, mm_alpha_ccn
     !        - mm_set_moments_thresholds
     !        - mm_get_rcs, mm_get_rcf
@@ -80 +80 @@
     ! Moments parameters (mm_aerosols_init)
     PROTECTED :: mm_m0aer_s, mm_m3aer_s, mm_m0aer_f, mm_m3aer_f
-    ! Moments parameters (derived, are updated with moments parameters)
-    PROTECTED :: mm_rcs, mm_rcf
     ! Thresholds parameters
     PROTECTED :: mm_m0as_min, mm_m3as_min, mm_rcs_min, mm_m0af_min, mm_m3af_min, mm_rcf_min
@@ -1170 +1168 @@
     END FUNCTION mm_alpha_f
 
+    PURE FUNCTION mm_alpha_ccn(k) RESULT (res)
+      !! Inter-moment relation for fractal aerosols size distribution law.
+      !! Mk / M0 = rc^k . alpha(k)
+      !!
+      REAL(kind=mm_wp), INTENT(in) :: k ! k order of the moment.
+      REAL(kind=mm_wp) :: sigma         ! Standard deviation.
+      REAL(kind=mm_wp) :: res           ! Alpha value.
+
+      ! Pluto's case
+      !~~~~~~~~~~~~~
+      sigma = 0.2_mm_wp
+      res = exp(k**2 * sigma**2 / 2._mm_wp)
+      
+      RETURN
+    END FUNCTION mm_alpha_ccn
+
 
     !============================================================================