Changeset 2166

Timestamp:

Oct 4, 2019, 3:07:15 PM (5 years ago)

Author:

abierjon

Message:

MARS GCM:
Handle case when injecting without dust scenario (iaervar=1) in compute_dtau_mod
AB

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/compute_dtau_mod.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/README

@@ -2757 +2757 @@
 be computed every "ichemistry" physics step (but chemistry tendencies are
 added every physics step). Default is ichemistry=1.
+
+== 04/10/2019 == AB+EM
+Handle case when injecting without dust scenario (iaervar=1) in compute_dtau_mod

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

@@ -19 +19 @@
         USE comcstfi_h, only: g
         USE tracer_mod, only: alpha_lift,igcm_dust_mass,igcm_dust_number
+        USE dimradmars_mod, only: tauvis
         
         IMPLICIT NONE
@@ -61 +62 @@
         
         ! 1. Obtain tauref_scenario from dust scenario at zday+1
-        zday_scenario=zday-modulo(zday,1.) ! integer value of the day: the scenario opacity is measured at 14:00
-        zday_scenario=zday_scenario+1      ! opacity of the dust scenario is read the day after
-        call read_dust_scenario(ngrid,nlayer,zday_scenario,pplev,       &
-       &                        tauref_scenario)
+        if (iaervar.eq.1) then
+          tauref_scenario = tauvis
+        else
+          zday_scenario=zday-modulo(zday,1.) ! integer value of the day: the scenario opacity is measured at 14:00
+          zday_scenario=zday_scenario+1      ! opacity of the dust scenario is read the day after
+          call read_dust_scenario(ngrid,nlayer,zday_scenario,pplev,     &
+       &                          tauref_scenario)
+        endif
        ! for diagnostics
         call WRITEDIAGFI(ngrid,"tauref_scenario","tauref_scenario",     &