Changeset 2081

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

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/physiq_mod.F (modified) (2 diffs)

trunk/LMDZ.MARS/README

@@ -2685 +2685 @@
 - typo in vdifc.F (visible in debug mode) in the case of rdstorm=false, condition for defining pdqsdif(stormdust_mass/number)
 - typo in physiq_mod.F (visible in debug mode): wrong indices for tracers in aerosol(:,:,:)
+
+== 23/01/2019 == MV
+- follow-up of the last change in physiq_mod.F: put iaer indices (defined in dimradmars_mod.F) for aerosol(:,:,:) in writediagfi calls.

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))