Changeset 3682 for trunk/LMDZ.TITAN/libf/muphytitan/mm_globals.f90

Timestamp:

Mar 14, 2025, 2:27:06 PM (3 months ago)

Author:

lrosset

Message:

Titan microphysics and physics : Outputs the nucleation rates and growth rates for the condensible species. (LR)

File:

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

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

@@ -522 +522 @@
   REAL(kind=mm_wp), DIMENSION(:,:), ALLOCATABLE, SAVE :: mm_ice_fluxes
 
+  !> Condensible components nucleation rates.
+  !!
+  !! It is a 2D-array with the vertical layers in first dimension and the number of condensible
+  !! species in the second.
+  !! It is updated in [[mm_clouds(module):mm_cloud_microphysics(subroutine)]].
+  REAL(kind=mm_wp), DIMENSION(:,:), ALLOCATABLE, SAVE :: mm_nrate
+
+  !> Condensible components growth rates.
+  !!
+  !! It is a 2D-array with the vertical layers in first dimension and the number of condensible
+  !! species in the second.
+  !! It is updated in [[mm_clouds(module):mm_cloud_microphysics(subroutine)]].
+  REAL(kind=mm_wp), DIMENSION(:,:), ALLOCATABLE, SAVE :: mm_grate
+
   !> Condensible species saturation ratio (--).
   !!
@@ -583 +597 @@
   !$OMP THREADPRIVATE(mm_rcs,mm_rcf,mm_drad,mm_drho)
   !$OMP THREADPRIVATE(mm_m0as_vsed,mm_m3as_vsed,mm_m0af_vsed,mm_m3af_vsed)
-  !$OMP THREADPRIVATE(mm_aer_s_flux,mm_aer_f_flux,mm_ccn_vsed,mm_ccn_flux,mm_ice_prec,mm_ice_fluxes,mm_gazs_sat)
+  !$OMP THREADPRIVATE(mm_aer_s_flux,mm_aer_f_flux,mm_ccn_vsed,mm_ccn_flux,mm_ice_prec,mm_ice_fluxes,mm_gazs_sat, mm_nrate, mm_grate)
   !$OMP THREADPRIVATE(mm_m0as_min,mm_m3as_min,mm_rcs_min,mm_m0af_min,mm_m3af_min,mm_rcf_min,mm_m0n_min,mm_m3cld_min)
   !$OMP THREADPRIVATE(mm_nla,mm_nle)
@@ -1251 +1265 @@
       ALLOCATE(mm_gazs_sat(mm_nla,mm_nesp)) ; mm_gazs_sat(:,:) = 0._mm_wp
     ENDIF
+    ! ADDED : Allocate memory for nucleation and growth rate
+    IF (.NOT.ALLOCATED(mm_nrate)) THEN
+      ALLOCATE(mm_nrate(mm_nla,mm_nesp)) ; mm_nrate(:,:) = 0._mm_wp
+    ENDIF
+      IF (.NOT.ALLOCATED(mm_grate)) THEN
+      ALLOCATE(mm_grate(mm_nla,mm_nesp)) ; mm_grate(:,:) = 0._mm_wp
+    ENDIF
 
     ! note mm_dzlev already from top to ground