Changeset 2189

Timestamp:

Dec 12, 2019, 5:17:03 PM (5 years ago)

Author:

emillour

Message:

Mars GCM:
Bug fix in the scavenging routine (missing initialisation of tendencies of
surface tracers) in co2condens.
EM

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/co2condens_mod.F (modified) (1 diff)

trunk/LMDZ.MARS/README

@@ -2793 +2793 @@
 == 04/12/2019 == AB
 Add the instantaneous scavenging by CO2 of dust, ccn and water ice in co2condens_mod. It can be activated or deactivated with the new flag "scavco2cond" (=false by default). Expected to be replaced by the CO2 clouds microphysics in the future.
+
+== 12/12/2019 == EM
+Bug fix in the scavenging routine (missing initialisation of tendencies of
+surface tracers) in co2condens.
+
+ 

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

@@ -209 +209 @@
 c     Tendencies set to 0
 c     -------------------------------------
-      DO l=1,nlayer
-         DO ig=1,ngrid
-           zcondicea(ig,l) = 0.
-           zfallice(ig,l) = 0.
-           pduc(ig,l)  = 0
-           pdvc(ig,l)  = 0
-           pdtc(ig,l) = 0.
-         END DO
-      END DO
+      zcondicea(1:ngrid,1:nlayer) = 0.
+      zfallice(1:ngrid,1:nlayer+1) = 0.
+      pduc(1:ngrid,1:nlayer)  = 0
+      pdvc(1:ngrid,1:nlayer)  = 0
+      pdtc(1:ngrid,1:nlayer) = 0.
+      pdqsc(1:ngrid,1:nq) = 0
          
-      DO iq=1,nq
-      DO l=1,nlayer
-         DO ig=1,ngrid
-           pdqc(ig,l,iq)  = 0
-         END DO
-      END DO
-      END DO
-
-      DO ig=1,ngrid
-         zfallice(ig,nlayer+1) = 0.
-         zcondices(ig) = 0.
-         pdtsrfc(ig) = 0.
-         pdpsrf(ig) = 0.
-         condsub(ig) = .false.
-         zdiceco2(ig) = 0.
-      ENDDO
+      pdqc(1:ngrid,1:nlayer,1:nq)  = 0
+
+      zcondices(1:ngrid) = 0.
+      pdtsrfc(1:ngrid) = 0.
+      pdpsrf(1:ngrid) = 0.
+      condsub(1:ngrid) = .false.
+      zdiceco2(1:ngrid) = 0.
+
       zfallheat=0