Changeset 2080 for trunk/LMDZ.MARS/libf

Timestamp:

Jan 22, 2019, 7:05:11 PM (6 years ago)

Author:

mvals

Message:

Mars GCM:

• 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(:,:,:)

MV

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• physiq_mod.F (modified) (1 diff)
• vdifc_mod.F (modified) (2 diffs)

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

@@ -3068 +3068 @@
            if (rdstorm) then
              call writediagfi(1,'aerosol_dust','opacity of env. dust',''
-     &                        ,1,aerosol(:,:,igcm_dust_mass))
+     &                        ,1,aerosol(:,:,1))
              call writediagfi(1,'aerosol_stormdust',
      &                         'opacity of storm dust',''
-     &                        ,1,aerosol(:,:,igcm_stormdust_mass))
+     &                        ,1,aerosol(:,:,2))
              call WRITEDIAGFI(ngrid,'dqsdifdustq','diffusion',
      &                       'kg.m-2.s-1',0,zdqsdif(1,igcm_dust_mass))

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

@@ -548 +548 @@
 
 ! calcul de zc et zd pour la couche top en prenant en compte le terme
-! de variation de masse (on fait une boucle pour que ça converge)
+! de variation de masse (on fait une boucle pour que \E7a converge)
 
 ! Identification des points de grilles qui ont besoin de la correction
@@ -740 +740 @@
                   endif ! if (submicron)
                 ELSE ! outside dust injection time frame
-                        pdqsdif(ig,igcm_dust_mass)= 0.
-                        pdqsdif(ig,igcm_dust_number)= 0.
+                  pdqsdif(ig,igcm_dust_mass)= 0.
+                  pdqsdif(ig,igcm_dust_number)= 0.
+                  if (rdstorm) then      
                         pdqsdif(ig,igcm_stormdust_mass)= 0.
                         pdqsdif(ig,igcm_stormdust_number)= 0.
+                  end if
                 ENDIF