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

Timestamp:

Nov 5, 2025, 6:17:11 PM (4 weeks ago)

Author:

debatzbr

Message:

Pluto PCM: Add cloud initialization and cloud-related diagnostics.
BBT

File:

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

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

@@ -8 +8 @@
     !
     !     The module also contains three methods:
-    !        - mm_initialize(dt,haze_prod_pCH4,p_prod,tx_prod,rc_prod,rm,df,rho_aer,rplanet,g0,air_rad,air_mmol,cfgpath)
+    !        - mm_initialize(dt,call_CH4hazeprod,p_prod,tx_prod,rc_prod,rm,df,rho_aer,rplanet,g0,air_rad,air_mmol,cfgpath)
     !        - read_aprm(parser,sec,pp)
     !        - abort_program(err)
@@ -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,clouds,cfgpath)
+    SUBROUTINE mm_initialize(dt,call_CH4hazeprod,p_prod,tx_prod,rc_prod,rm,df,rho_aer,rplanet,g0,air_rad,air_mmol,clouds,cfgpath)
         !! Initialize global parameters of the model.
         !!
@@ -46 +46 @@
         REAL(kind=mm_wp), INTENT(in)           :: dt
         ! Enable/Disable production from CH4 photolysis.
-        LOGICAL                                :: haze_prod_pCH4
+        LOGICAL                                :: call_CH4hazeprod
         ! Aerosol production pressure level (Pa).
         REAL(kind=mm_wp), INTENT(in)           :: p_prod
@@ -84 +84 @@
         ! Thresholds related parameters.
         REAL(kind=mm_wp)       :: m0as_min,rcs_min,m0af_min,rcf_min
+        REAL(kind=mm_wp)       :: m0ccn_min,drad_min
         ! Debug mode control flag (may print lot of stuff if enabled).
         LOGICAL                :: wdebug
@@ -105 +106 @@
         m0af_min    = 1e-8_mm_wp
         rcf_min     = 1e-9_mm_wp
+        m0ccn_min   = 1e-8_mm_wp
+        drad_min    = 1e-9_mm_wp
         wdebug      = .false.
     
@@ -132 +135 @@
         err = mm_check_opt(cfg_get_value(cparser,"m0af_min",m0af_min)                 ,m0af_min   ,1e-8_mm_wp,mm_log)
         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,"m0ccn_min",m0ccn_min)               ,m0ccn_min  ,1e-8_mm_wp,mm_log)
+        err = mm_check_opt(cfg_get_value(cparser,"drad_min",drad_min)                 ,drad_min   ,1e-9_mm_wp,mm_log)
         err = mm_check_opt(cfg_get_value(cparser,"debug",wdebug)                      ,wdebug     ,.false.   ,mm_log)
 
@@ -154 +159 @@
         ! Transfert probabilities (S --> F):
         ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-        err = mm_check_opt(cfg_get_value(cparser, "transfert_probability", mm_w_ps2s), mm_w_ps2s, wlog=mm_log)
+        err = mm_check_opt(cfg_get_value(cparser, "transfert_probability", mm_call_ps2s), mm_call_ps2s, wlog=mm_log)
         IF (err/=0) call abort_program(err)
 
-        IF (mm_w_haze_coag .AND. mm_w_ps2s) THEN
+        IF (mm_call_hazecoag .AND. mm_call_ps2s) THEN
             err = mm_check_opt(cfg_get_value(cparser, "ps2s_file", pssfile), pssfile)
             IF (err /= 0) call abort_program(err)
@@ -177 +182 @@
         ! Mean electric correction:
         ! ~~~~~~~~~~~~~~~~~~~~~~~~~
-        err = mm_check_opt(cfg_get_value(cparser, "electric_charging", mm_w_qe), mm_w_qe, wlog=mm_log)
+        err = mm_check_opt(cfg_get_value(cparser, "electric_charging", mm_call_qe), mm_call_qe, wlog=mm_log)
         IF (err/=0) call abort_program(err)
 
-        IF (mm_w_haze_coag .AND. mm_w_qe) THEN
+        IF (mm_call_hazecoag .AND. mm_call_qe) THEN
             err = mm_check_opt(cfg_get_value(cparser, "mq_file", mqfile), mqfile)
             IF (err /= 0) call abort_program(err)
@@ -224 +229 @@
         ! 2. YAMMS initialization
         !------------------------
-        err = mm_global_init_0(dt,df,rm,rho_aer,haze_prod_pCH4,p_prod,tx_prod,rc_prod, &
-                               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,                      &
-                               clouds,spcpath,wdebug)
+        err = mm_global_init_0(dt,df,rm,rho_aer,call_CH4hazeprod,p_prod,tx_prod,rc_prod, &
+                               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,                        &
+                               clouds,spcpath,m0ccn_min,drad_min                         &
+                               wdebug)
         IF (err /= 0) call abort_program(err)
         
         ! Dump parameters.
         WRITE(*,'(a)')        "========= MUPHYS PARAMETERS ==========="
-        WRITE(*,'(a,L2)')     "transfert_probability: ", mm_w_ps2s
-        WRITE(*,'(a,L2)')     "electric_charging    : ", mm_w_qe
+        WRITE(*,'(a,L2)')     "transfert_probability: ", mm_call_ps2s
+        WRITE(*,'(a,L2)')     "electric_charging    : ", mm_call_qe
         call mm_dump_parameters()