Changeset 3951 for trunk/LMDZ.PLUTO/libf/muphypluto/mp2m_intgcm.F90

Timestamp:

Nov 4, 2025, 5:51:18 PM (4 weeks ago)

Author:

debatzbr

Message:

Pluto PCM: Add variables, indices, and flags related to microphysical clouds
BBT

File:

• trunk/LMDZ.PLUTO/libf/muphypluto/mp2m_intgcm.F90 (modified) (7 diffs)

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

@@ -29 +29 @@
     CONTAINS
 
-    SUBROUTINE mm_initialize(dt,haze_prod_pCH4,p_prod,tx_prod,rc_prod,rm,df,rho_aer,rplanet,g0,air_rad,air_mmol,cfgpath)
+    SUBROUTINE mm_initialize(dt,haze_prod_pCH4,p_prod,tx_prod,rc_prod,rm,df,rho_aer,rplanet,g0,air_rad,air_mmol,clouds,cfgpath)
         !! Initialize global parameters of the model.
         !!
@@ -67 +67 @@
         ! Mean molar mass of air molecules (kg.mol-1).
         REAL(kind=mm_wp), INTENT(in)           :: air_mmol
+        ! Clouds microphysics control flag.
+        LOGICAL, INTENT(in)                    :: clouds
         ! Internal microphysics configuration file.
         CHARACTER(len=*), INTENT(in), OPTIONAL :: cfgpath
@@ -72 +74 @@
         ! Microphysical configuration file.
         TYPE(cfgparser)        :: cparser
-        ! Look-up tables (transfert probabilities, mean electric correction).
+        ! Look-up tables - Haze related: transfert probabilities, mean electric correction.
         CHARACTER(len=st_slen) :: pssfile,mqfile
+        ! Look-up tables - Cloud related: species properties.
+        CHARACTER(len=st_slen) :: spcpath
         ! Enable/disable Haze process.
         LOGICAL                :: w_h_prod,w_h_sed,w_h_coag,fwsed_m0,fwsed_m3
@@ -84 +88 @@
 
         ! Local variables.
-        TYPE(error)                                 :: err
+        INTEGER     :: i
+        TYPE(error) :: err
         REAL(kind=mm_wp), DIMENSION(:), ALLOCATABLE :: tmp
 
@@ -128 +133 @@
         err = mm_check_opt(cfg_get_value(cparser,"rcf_min",rcf_min)                   ,rcf_min    ,rm        ,mm_log)
         err = mm_check_opt(cfg_get_value(cparser,"debug",wdebug)                      ,wdebug     ,.false.   ,mm_log)
+
+        ! Retrieve clouds species configuration file
+        spcpath = ''
+        IF (clouds) THEN
+          err = mm_check_opt(cfg_get_value(cparser,"species_cfg",spcpath), spcpath, wlog=mm_log)
+          IF (err/=0) call abort_program(err)
+        ENDIF
     
         ! Alpha function: 
@@ -215 +227 @@
                                rplanet,g0,air_rad,air_mmol,coag_choice,                &
                                w_h_prod,w_h_sed,w_h_coag,fwsed_m0,fwsed_m3,            &
-                               m0as_min,rcs_min,m0af_min,rcf_min,wdebug)
+                               m0as_min,rcs_min,m0af_min,rcf_min,                      &
+                               clouds,spcpath,wdebug)
         IF (err /= 0) call abort_program(err)
         
@@ -223 +236 @@
         WRITE(*,'(a,L2)')     "electric_charging    : ", mm_w_qe
         call mm_dump_parameters()
+
+        IF (clouds) THEN
+            DO i = 1, size(mm_xESPS)
+                print*, TRIM(mm_xESPS(i)%name), " fmol2fmas = ", mm_xESPS(i)%fmol2fmas
+            ENDDO
+        ENDIF
         
     END SUBROUTINE mm_initialize