Changeset 2643 for trunk/LMDZ.MARS/libf/phymars/dust_param_mod.F90

Timestamp:

Mar 17, 2022, 2:36:31 PM (3 years ago)

Author:

abierjon

Message:

Mars GCM:
Implementation of the IRabs-to-VISext dust scenario conversion, depending on the GCM effective radius :

• new flag 'reff_driven_IRtoVIS_scenario', false by default, must be set to true to use this dynamic conversion (otherwise, the coefficient takes the constant value of 2.6, eqv to reff=1.5um, as before) A specific line is added in deftank/callphys.def.GCM6
• this flag requires the 'dustiropacity' to be set to 'tes' to match the IR scenario's wavelength. 'mcs'-like dso diagnostics can only be produced by the use of the post-processing tool util/aeroptical.F90
• the variable IRtoVIScoef is computed in aeropacity via the GCM CDOD ratio : tau_pref_gcm_VIS/tau_pref_gcm_IR (only at the first call to callradite during each physical timestep)
• change read_dust_scenario into module read_dust_scenario_mod

AB

File:

• trunk/LMDZ.MARS/libf/phymars/dust_param_mod.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/libf/phymars/dust_param_mod.F90

@@ -8 +8 @@
   LOGICAL,SAVE :: freedust ! if true: no rescaling (via tauscaling) of the dust mass and number
   LOGICAL,SAVE :: callddevil ! flag to activate dust devil (dust lifing/injection) parametrization
+  LOGICAL,SAVE :: reff_driven_IRtoVIS_scenario ! use GCM dust size to convert IR scenarios to VIS
 
 !$OMP THREADPRIVATE(active, doubleq,submicron,lifting,freedust,     &
@@ -28 +29 @@
                     ! tau_pref_scenario is deemed exact
 
-!$OMP THREADPRIVATE(tauscaling, dustscaling_mode,dust_rad_adjust)
+!$OMP THREADPRIVATE(tauscaling,dustscaling_mode,dust_rad_adjust)
 
 contains