Changeset 3657 for trunk/LMDZ.TITAN/libf/phytitan/inifis_mod.F90

Timestamp:

Feb 27, 2025, 2:53:07 PM (6 months ago)

Author:

emoisan

Message:

Titan physics (under MESOSCALE precompiling flag):
Change initialization for microphysics model, to be consistent with
initialization for physics.
The parameters are now defined in callphys.def instead of datagcm.
The new callphys.def with all the parameters is in deftank
(callphys.def.allmuparams).
EMo

File:

• trunk/LMDZ.TITAN/libf/phytitan/inifis_mod.F90 (modified) (6 diffs)

trunk/LMDZ.TITAN/libf/phytitan/inifis_mod.F90

@@ -20 +20 @@
   use planete_mod, only: nres
   use planetwide_mod, only: planetwide_sumval
+  use mm_globals, only: mm_global_init_0 
+  use mmp_globals, only: mmp_w_ps2s, mmp_w_qe
   use callkeys_mod
   use mod_phys_lmdz_para, only : is_parallel
+  use errors, only: error
 
 !=======================================================================
@@ -73 +76 @@
   integer,intent(in) :: day_ini
   INTEGER ig,ierr
+  TYPE(error) :: err
  
   EXTERNAL iniorbit,orbite
   EXTERNAL SSUM
   REAL SSUM
- 
+
   ! Initialize flags lunout, prt_level, debug (in print_control_mod)
   CALL init_print_control
@@ -93 +97 @@
 
   ! Initialize some "temporal and calendar" related variables
-#ifndef MESOSCALE
+  ! -- MESOSCALE : we do it for the right print (and avoid crash in debug)
+  ! -- MESOSCALE : and sporadic use in physiq, e.g. daysec used in soil
   CALL init_time(day_ini,pdaysec,nday,ptimestep)
-#endif
 
   ! read in some parameters from "run.def" for physics,
@@ -503 +507 @@
 
      write(*,*) "Equivalent radius production (m) ?"
-     rc_prod=2.0e-8 ! default value
+     rc_prod=1.0e-8 ! default value
      call getin_p("rc_prod",rc_prod)
      write(*,*)" rhc_prod = ",rc_prod
@@ -513 +517 @@
 
      write(*,*) "Path to microphys. config file ?"
-     config_mufi='datagcm/microphysics/config.cfg' ! default value
+     ! default path is set in datadir
+     write(*,*) "microphysics config file data base folder:",trim(datadir),"/microphysics"
+     config_mufi=trim(datadir)//'/microphysics/config.ne15.cfg' ! default value
      call getin_p("config_mufi",config_mufi)
+     config_mufi = trim(datadir)//'/microphysics/'//config_mufi
      write(*,*)" config_mufi = ",config_mufi
+
+! in mesoscale the following variables are defined in callphys.def
+! before they were defined in the config file, in datagcm (EMoi)
+#ifdef MESOSCALE
+     write(*,*) "Monomer radius (m)"
+     rm=50e-9 ! default value
+     call getin_p("rm",rm)
+     write(*,*)" rm = ",rm
+ 
+     write(*,*) "Fractal dimension"
+     df=2.0 ! default value
+     call getin_p("df",df)
+     write(*,*)" df = ",df
+ 
+     write(*,*) "Aerosol density (kg.m-3)"
+     rho_aer=600. ! default value
+     call getin_p("rho_aer",rho_aer)
+     write(*,*)" rho_aer = ",rho_aer
+ 
+     write(*,*) "Enable/disable Haze production process"
+     w_h_prod=.true. ! default value
+     call getin_p("haze_production",w_h_prod)
+     write(*,*)" haze_production = ",w_h_prod
+ 
+     write(*,*) "Enable/disable Haze coagulation process"
+     w_h_coag=.true. ! default value
+     call getin_p("haze_coagulation",w_h_coag)
+     write(*,*)" haze_coagulation = ",w_h_coag
+ 
+     write(*,*) "Coagulation interactions"
+     coag_choice=7 ! default value
+     call getin_p("haze_coag_interactions",coag_choice)
+     write(*,*)" haze_coag_interactions = ",coag_choice
+ 
+     write(*,*) "Enable/disable Haze sedimentation process"
+     w_h_sed=.true. ! default value
+     call getin_p("haze_sedimentation",w_h_sed)
+     write(*,*)" haze_sedimentation = ",w_h_sed
+ 
+     write(*,*) "Disable Fiadero correction for sedimentation process"
+     no_fiadero=.true. ! default value
+     call getin_p("no_fiadero",no_fiadero)
+     write(*,*)" no_fiadero = ",no_fiadero
+ 
+     write(*,*) "Fiadero correction minimum ratio threshold"
+     fiad_min=0.1 ! default value
+     call getin_p("fiadero_min_ratio",fiad_min)
+     write(*,*)" fiadero_min_ratio = ",fiad_min
+ 
+     write(*,*) "Fiadero correction maximum ratio threshold"
+     fiad_max=10. ! default value
+     call getin_p("fiadero_max_ratio",fiad_max)
+     write(*,*)" fiadero_max_ratio = ",fiad_max
+ 
+     write(*,*) "Force settling velocity to M0"
+     fwsed_m0=.true. ! default value
+     call getin_p("wsed_m0",fwsed_m0)
+     write(*,*)" wsed_m0 = ",fwsed_m0
+ 
+     write(*,*) "Force settling velocity to M3"
+     fwsed_m3=.false. ! default value
+     call getin_p("wsed_m3",fwsed_m3)
+     write(*,*)" wsed_m3 = ",fwsed_m3
+ 
+     write(*,*) "Enable/disable clouds sedimentation process"
+     w_c_sed=.true. ! default value
+     call getin_p("clouds_sedimentation",w_c_sed)
+     write(*,*)" clouds_sedimentation = ",w_c_sed
+ 
+     write(*,*) "Enable/disable clouds nucleation and condensation processes"
+     w_c_nucond=.true. ! default value
+     call getin_p("clouds_nuc_cond",w_c_nucond)
+     write(*,*)" clouds_nuc_cond = ",w_c_nucond
+ 
+     write(*,*) "Path to condensible species configuration file"
+     ! default path is set in datadir
+     write(*,*) "condensible species config file data base folder:",trim(datadir),"/microphysics"
+     spcpath=trim(datadir)//"/microphysics/mp2m_species.cfg" ! default value
+     call getin_p("specie_cfg",spcpath)
+     spcpath = trim(datadir)//'/microphysics/'//spcpath
+     write(*,*)" specie_cfg = ",spcpath
+ 
+     write(*,*) "Enable/disable spherical mode transfert probability"
+     mmp_w_ps2s=.true. ! default value
+     call getin_p("transfert_probability",mmp_w_ps2s)
+     write(*,*)" transfert_probability = ",mmp_w_ps2s
+ 
+     write(*,*) "Path to spherical mode transfert probability look-up tables file"
+     ! default path is set in datadir
+     write(*,*) "spherical mode transfert probability look-up tables file data base folder:",trim(datadir),"/microphysics"
+     pssfile=trim(datadir)//"/microphysics/mmp_ps2s_rm50_ne15.nc" ! default value
+     call getin_p("ps2s_file",pssfile)
+     pssfile = trim(datadir)//'/microphysics/'//pssfile
+     write(*,*)" ps2s_file = ",pssfile
+ 
+     write(*,*) "Electric charging coagulation correction"
+     mmp_w_qe=.true. ! default value
+     call getin_p("electric_charging",mmp_w_qe)
+     write(*,*)" electric_charging = ",mmp_w_qe
+ 
+     write(*,*) "Electric charging correction factor file"
+     ! default path is set in datadir
+     write(*,*) "Electric charging correction factor file data base folder:",trim(datadir),"/microphysics"
+     mqfile=trim(datadir)//"/microphysics/mmp_qmean_rm50_ne15.nc" ! default value
+     call getin_p("mq_file",mqfile)
+     mqfile = trim(datadir)//'/microphysics/'//mqfile
+     write(*,*)" mq_file = ",mqfile
+ 
+     write(*,*) "Aerosol spherical mode total number min. threshold"
+     m0as_min=1e-10 ! default value
+     call getin_p("m0as_min",m0as_min)
+     write(*,*)" m0as_min = ",m0as_min
+ 
+     write(*,*) "Aerosol spherical mode min. caracteristic radius threshold"
+     rcs_min=1e-9 ! default value
+     call getin_p("rcs_min",rcs_min)
+     write(*,*)" rcs_min = ",rcs_min
+ 
+     write(*,*) "Aerosol fractal mode total number min. threshold"
+     m0af_min=1e-10 ! default value
+     call getin_p("m0af_min",m0af_min)
+     write(*,*)" m0af_min = ",m0af_min
+ 
+     write(*,*) "Aerosol fractal mode min. caracteristic radius threshold (updated to monomer radius at runtime if needed)"
+     rcf_min=1e-9 ! default value
+     call getin_p("rcf_min",rcf_min)
+     write(*,*)" rcf_min = ",rcf_min
+ 
+     write(*,*) "cloud drop total number min. threshold"
+     m0n_min=1e-10 ! default value
+     call getin_p("m0n_min",m0n_min)
+     write(*,*)" m0n_min = ",m0n_min
+
+     write(*,*) "microphysics debug flag"
+     wdebug=.false. ! default value
+     call getin_p("wdebug",wdebug)
+     write(*,*)" wdebug = ",wdebug
+#endif
 
 ! Soil model
@@ -669 +814 @@
   ! allocate "comsoil_h" arrays
   call ini_comsoil_h(ngrid)
-    
+
+#ifdef MESOSCALE
+  ! YAMMS initialization.
+  err = mm_global_init_0(ptimestep,df,rm,rho_aer,p_prod,tx_prod,rc_prod,rad,g, &
+                         air_rad,mugaz,coag_choice,callclouds,spcpath,  &
+                         w_h_prod,w_h_sed,w_h_coag,w_c_nucond,  &
+                         w_c_sed,fwsed_m0,fwsed_m3, &
+                         no_fiadero,fiad_min,fiad_max, &
+                         m0as_min,rcs_min,m0af_min,rcf_min,m0n_min,wdebug)
+  IF (ierr /= 0) THEN
+       PRINT *, "PROBLEM WITH YAMMS INITIATILIZATION"
+  ENDIF
+#endif
+
+
   END SUBROUTINE inifis