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

Timestamp:

Oct 16, 2020, 9:40:05 AM (4 years ago)

Author:

emillour

Message:

Mars GCM:
Add a new scheme to handle the dust and its radiative impact. Triggered using
a new flag dustscaling_mode=2 (dustscaling_mod=0: no rescaling at all, and
dustscaling_mode=1: full rescaling using tauscaling, GCMv5.3 style). Rescaling
is then only done on the radiative impact (see dust_scaling_mod.F90) of dust.
Moreover the scaling factor "dust_rad_adjust" is evaluated using the target dust
scenario opacity for the next sol and left to evolve linearly until then to not
impose the diurnal evolution of dust.
In practice, main changes or additions in the code are:

• renamed flag "tauscaling_mode" as "dustscaling_mode"
• moved parameter "t_scenario_sol" to "dust_param_mod"
• adapted "compute_dustscaling" routine in "dust_scaling_mod"
• added module "dust_rad_adjust_mod"
• 2D fields "dust_rad_adjust_prev" and "dust_rad_adjust_next" required to compute coefficient "dust_rad_adjust" need to be stored in (re)startfi files

EM

File:

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

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

@@ -15 +15 @@
   
   REAL,SAVE,ALLOCATABLE :: tauscaling(:)   ! Convertion factor for qdust and Ndust
-  INTEGER,SAVE :: tauscaling_mode ! =0, no rescaling (freedust)
-                                  ! =1, prescribed scaling GCM5.3 style
+  INTEGER,SAVE :: dustscaling_mode ! dust scaling modes
+                  ! =0, no rescaling (freedust)
+                  ! =1, prescribed scaling GCM5.3 style (using tauscaling)
+                  ! =2, only radiative scaling (using dust_rad_adjust)
+  REAL,SAVE,ALLOCATABLE :: dust_rad_adjust(:) ! radiative scaling for dust
+  REAL,PARAMETER :: t_scenario_sol=14/24. ! time of day (sol) at which
+                    ! tau_pref_scenario is deemed exact
 
 contains
@@ -25 +30 @@
   
     allocate(tauscaling(ngrid))
+    allocate(dust_rad_adjust(ngrid))
   
   end subroutine ini_dust_param_mod
@@ -32 +38 @@
     
     if (allocated(tauscaling))  deallocate(tauscaling)
+    if (allocated(dust_rad_adjust)) deallocate(dust_rad_adjust)
 
   end subroutine end_dust_param_mod