Changeset 2081 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Jan 23, 2019, 10:49:43 AM (6 years ago)

Author:

mvals

Message:

Mars GCM:

• follow-up of the last change in physiq_mod.F: put iaer indices (defined in dimradmars_mod.F) for aerosol(:,:,:) in writediagfi calls.

MV

File:

• trunk/LMDZ.MARS/libf/phymars/physiq_mod.F (modified) (2 diffs)

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -49 +49 @@
       use dimradmars_mod, only: tauscaling, aerosol, totcloudfrac,
      &                          dtrad, fluxrad_sky, fluxrad, albedo,
-     &                          naerkind
+     &                          naerkind, iaer_dust_doubleq, 
+     &                          iaer_stormdust_doubleq
       use turb_mod, only: q2, wstar, ustar, sensibFlux, 
      &                    zmax_th, hfmax_th, turb_resolved
@@ -3068 +3069 @@
            if (rdstorm) then
              call writediagfi(1,'aerosol_dust','opacity of env. dust',''
-     &                        ,1,aerosol(:,:,1))
+     &                        ,1,aerosol(:,:,iaer_dust_doubleq))
              call writediagfi(1,'aerosol_stormdust',
      &                         'opacity of storm dust',''
-     &                        ,1,aerosol(:,:,2))
+     &                        ,1,aerosol(:,:,iaer_stormdust_doubleq))
              call WRITEDIAGFI(ngrid,'dqsdifdustq','diffusion',
      &                       'kg.m-2.s-1',0,zdqsdif(1,igcm_dust_mass))