Changeset 2189 for trunk/LMDZ.MARS/libf

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

File:

• trunk/LMDZ.MARS/libf/phymars/co2condens_mod.F (modified) (1 diff)

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