Changeset 2080 for trunk

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

Files:

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

trunk/LMDZ.MARS/README

@@ -2681 +2681 @@
 We want to separate both parametrizations related to the formation of the detached dust layers. Therefore, rocketduststorm_mod.F90 now only comprises the rocket dust storm scheme, whereas it contained also before the calculation of the vertical velocity induced by the presence of the sub-grid scale topography. This latter part is under development and will be integrated as a fully independant parametrization: the aim is to simulate the entrainment of dust by slope winds, from the boundary layer up to the top of the sub-grid scale topography.
 - Addition of initial surface parameters "summit" and "base" to prepare the previously described slope wind parametrization
+
+== 22/01/2019 == MV
+- 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(:,:,:)

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