Changeset 2643 for trunk/LMDZ.MARS/libf/phymars/topmons_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/topmons_mod.F90 (modified) (3 diffs)

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

@@ -22 +22 @@
                                  icount,zday,zls,tsurf,igout,aerosol,  &
                                  tauscaling,dust_rad_adjust,           &
+                                 IRtoVIScoef,                          &
 !             input sub-grid scale rocket dust storm
                                  totstormfract,clearatm,               &
@@ -75 +76 @@
       REAL, INTENT(INOUT) :: tauscaling(ngrid)
       REAL,INTENT(INOUT) :: dust_rad_adjust(ngrid)
+      REAL,INTENT(INOUT) :: IRtoVIScoef(ngrid) ! NB: not modified by this call to callradite,
+                                               ! the OUT is just here because callradite needs it
+                                               
 !     input sub-grid scale rocket dust storm
       LOGICAL, INTENT(IN) :: clearatm
@@ -284 +288 @@
                  zdtlw1,zdtsw1,fluxsurf_lw1,fluxsurf_sw1,fluxtop_lw1,  &
                  fluxtop_sw1,tau_pref_scenario,tau_pref_gcm, &
-                 tau,aerosol,dsodust,tauscaling,dust_rad_adjust,    &
+                 tau,aerosol,dsodust,tauscaling,dust_rad_adjust,IRtoVIScoef, &
                  taucloudtes,rdust,rice,nuice,riceco2,nuiceco2,co2ice,rstormdust,rtopdust, &
                  totstormfract,clearatm,dsords,dsotop,nohmons,&