Changeset 1985 for trunk

Timestamp:

Aug 9, 2018, 4:27:52 PM (6 years ago)

Author:

mvals

Message:

Mars GCM:
Correction of the dust storm scheme:

• in rocketduststorm_mod.F90 the case of not daytime and no storm was not taken into account (it was actually in comment), which lead to NaNs? in the calculation of pdqrds (scheme=0, division by alpha=0.)

Useful:

• in compute_dtau_mod.F90, the writediagfi of tauref_scenario has been added

MV

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/compute_dtau_mod.F90 (modified) (1 diff)
• libf/phymars/rocketduststorm_mod.F90 (modified) (1 diff)

trunk/LMDZ.MARS/README

@@ -2632 +2632 @@
 Correction:
 - in callsedim_mod, declaration of variable tau(ngrid,nlay) corrected to tau(ngrid,naerkind)
+
+== 09/08/2018 == MV
+Correction of the dust storm scheme:
+- in rocketduststorm_mod.F90 the case of not daytime and no storm was not taken into account (it was actually in comment), which lead to NaNs in the calculation of pdqrds (scheme=0, division by alpha=0.)
+Useful:
+- in compute_dtau_mod.F90, the writediagfi of tauref_scenario has been added

trunk/LMDZ.MARS/libf/phymars/compute_dtau_mod.F90

@@ -59 +59 @@
         call read_dust_scenario(ngrid,nlayer,zday_scenario,pplev,       &
        &                        tauref_scenario)
+       ! for diagnostics
+        call WRITEDIAGFI(ngrid,"tauref_scenario","tauref_scenario",     &
+       &                  "",2,tauref_scenario)
 
         ! 2. Compute dtau

trunk/LMDZ.MARS/libf/phymars/rocketduststorm_mod.F90

@@ -331 +331 @@
       DO ig=1,ngrid
         !! **********************************************************************
-        !! 2.1 case 1: Nothing to do when no storm and no slope
-!        IF ((mu0(ig) .LE. mu0lim) .AND. .NOT.(storm(ig)) ) then
-!          scheme(ig)=1
-!          cycle
-!        endif
-        IF ((alpha_hmons(ig) .EQ. 0.) .AND. .NOT.(storm(ig))) then
+        !! 2.1 case 1: Nothing to do when no storm and no slope, or
+        !!             no storm and not daytime
+        IF (((alpha_hmons(ig) .EQ. 0.) .AND. .NOT.(storm(ig))) &
+             .OR. ((mu0(ig) .LE. mu0lim) .AND. .NOT.(storm(ig))) ) then
           scheme(ig)=1
-          cycle !!no slope 
+          cycle
         endif