Changeset 2166 for trunk/LMDZ.MARS/libf

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

File:

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

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",     &