Changeset 3889

Timestamp:

Aug 12, 2025, 2:32:04 PM (3 days ago)

Author:

debatzbr

Message:

Pluto PCM: Pass haze production parameters as keys in callphys.def
BBT

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• callkeys_mod.F90 (modified) (1 diff)
• hazecloud.F90 (modified) (2 diffs)
• inifis_mod.F90 (modified) (1 diff)
• optci.F90 (modified) (2 diffs)
• optcv.F90 (modified) (2 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/callkeys_mod.F90

@@ -208 +208 @@
 !$OMP THREADPRIVATE(haze_p_prod,haze_tx_prod,haze_rc_prod,haze_rm,haze_df,haze_rho,air_rad)
 
+!! Parameter of conversion precursors to haze (hazecloud.F90)
+      real,save :: tcon_ch4, k_ch4, ncratio_ch4
+!$OMP THREADPRIVATE(tcon_ch4, k_ch4, ncratio_ch4)
+
+!! Variables for aerosol absorption
+      real,save :: Fabs_aers_VI, Fabs_aerf_VI, Fabs_aers_IR, Fabs_aerf_IR
+!$OMP THREADPRIVATE(Fabs_aers_VI, Fabs_aerf_VI, Fabs_aers_IR, Fabs_aerf_IR)
+
       integer,save :: iddist
       integer,save :: iaervar

trunk/LMDZ.PLUTO/libf/phypluto/hazecloud.F90

@@ -8 +8 @@
       use tracer_h, only: igcm_haze, igcm_ch4_gas, igcm_prec_haze, noms, mmol
       use geometry_mod, only: longitude, latitude ! in radians
-      use callkeys_mod, only: hazeconservch4, diurnal
+      use callkeys_mod, only: hazeconservch4, diurnal, tcon_ch4, k_ch4, ncratio_ch4
 
       implicit none
@@ -159 +159 @@
 !! Time constant of conversion in aerosol [second]
 !! To be explore: 1.E3; 1.E5; 1.E7; 1.E9
-      tcon= 1.e5  
+      tcon= tcon_ch4
 !! Parameter of conversion precurseur to haze
-      kch4=0.3
-      ncratio=0.5     ! boost for haze considering nitrogen contribution
-                      ! ratio n/c : =0.25 if there is 1N for 3C
+      kch4=k_ch4
+      ncratio=ncratio_ch4 ! boost for haze considering nitrogen contribution
+                          ! ratio n/c : =0.25 if there is 1N for 3C
       IF (firstcall) then
         write(*,*) 'hazecloud: haze parameters:'

trunk/LMDZ.PLUTO/libf/phypluto/inifis_mod.F90

@@ -776 +776 @@
      if (is_master) write(*,*)" air_rad = ",air_rad
 
+     ! Conversion precursor to haze
+     ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+     if (is_master) write(*,*) "Time constant of conversion in aerosol (s)?"
+     tcon_ch4=1.e7 ! default value
+     call getin_p("tcon_ch4",tcon_ch4)
+     if (is_master) write(*,*)" tcon_ch4 = ",tcon_ch4
+
+     if (is_master) write(*,*) "Ratio CH4 photolysis?"
+     k_ch4=1. ! default value
+     call getin_p("k_ch4",k_ch4)
+     if (is_master) write(*,*)" k_ch4 = ",k_ch4
+
+     if (is_master) write(*,*) "Nitrogen contribution ratio N/C?"
+     ncratio_ch4=0.5 ! default value
+     call getin_p("ncratio_ch4",ncratio_ch4)
+     if (is_master) write(*,*)" ncratio_ch4 = ",ncratio_ch4
+
+     ! Variables for aerosol absorption
+     ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+     if (is_master) write(*,*) "Sph. aer. absorption correction in VI?"
+     Fabs_aers_VI=1. ! default value
+     call getin_p("Fabs_aers_VI",Fabs_aers_VI)
+     if (is_master) write(*,*)" Fabs_aers_VI = ",Fabs_aers_VI
+
+     if (is_master) write(*,*) "Fra. aer. absorption correction in VI?"
+     Fabs_aerf_VI=1. ! default value
+     call getin_p("Fabs_aerf_VI",Fabs_aerf_VI)
+     if (is_master) write(*,*)" Fabs_aerf_VI = ",Fabs_aerf_VI
+
+     if (is_master) write(*,*) "Sph. aer. absorption correction in IR?"
+     Fabs_aers_IR=1. ! default value
+     call getin_p("Fabs_aers_IR",Fabs_aers_IR)
+     if (is_master) write(*,*)" Fabs_aers_IR = ",Fabs_aers_IR
+
+     if (is_master) write(*,*) "Fra. aer. absorption correction in IR?"
+     Fabs_aerf_IR=1. ! default value
+     call getin_p("Fabs_aerf_IR",Fabs_aerf_IR)
+     if (is_master) write(*,*)" Fabs_aerf_IR = ",Fabs_aerf_IR
+
      ! Pluto haze model
      ! ~~~~~~~~~~~~~~~~

trunk/LMDZ.PLUTO/libf/phypluto/optci.F90

@@ -15 +15 @@
                      igas_CH4, igas_N2
   use comcstfi_mod, only: g, r, mugaz
-  use callkeys_mod, only: kastprof,continuum,graybody,callmufi
+  use callkeys_mod, only: kastprof,continuum,graybody,callmufi,Fabs_aers_IR,Fabs_aerf_IR
   use recombin_corrk_mod, only: corrk_recombin, gasi_recomb
   use tpindex_mod, only: tpindex
@@ -144 +144 @@
 
   if(callmufi) then
-   Fabs_aer(1) = 1.2
-   Fabs_aer(2) = 1.3
+   Fabs_aer(1) = Fabs_aers_IR
+   Fabs_aer(2) = Fabs_aerf_IR
   else
    Fabs_aer(:) = 1.0

trunk/LMDZ.PLUTO/libf/phypluto/optcv.F90

@@ -14 +14 @@
                      igas_CH4, igas_N2
   use comcstfi_mod, only: g, r, mugaz
-  use callkeys_mod, only: kastprof,continuum,graybody,callgasvis,callmufi
+  use callkeys_mod, only: kastprof,continuum,graybody,callgasvis,callmufi,Fabs_aers_VI,Fabs_aerf_VI
   use recombin_corrk_mod, only: corrk_recombin, gasv_recomb
   use tpindex_mod, only: tpindex
@@ -141 +141 @@
 
   if(callmufi) then
-   Fabs_aer(1) = 0.7
-   Fabs_aer(2) = 2.5
+   Fabs_aer(1) = Fabs_aers_VI
+   Fabs_aer(2) = Fabs_aerf_VI
   else
    Fabs_aer(:) = 1.0