Changeset 1985 for trunk/LMDZ.MARS/libf/phymars/rocketduststorm_mod.F90

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

File:

• trunk/LMDZ.MARS/libf/phymars/rocketduststorm_mod.F90 (modified) (1 diff)

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