Changeset 2823 for trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F

Timestamp:

Nov 18, 2022, 10:43:43 AM (2 years ago)

Author:

emillour

Message:

Mars GCM:
Remove the "tracer" (logical) flag as we now always run with at least
one tracer.
EM

File:

• trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F (modified) (6 diffs)

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

@@ -211 +211 @@
          ico2=0
 
-         if (tracer) then
 c          Prepare Special treatment if one of the tracer is CO2 gas
            do iq=1,nq
@@ -226 +225 @@
              endif
            enddo
-         end if
 
         firstcall=.false.
@@ -351 +349 @@
          ENDDO
       ENDDO
-      if(tracer) then
-        DO iq =1, nq
-         DO ilev=1,nlay
-           DO ig=1,ngrid
-              zq(ig,ilev,iq)=pq(ig,ilev,iq)+pdqfi(ig,ilev,iq)*ptimestep
-           ENDDO
-         ENDDO
-        ENDDO
-      end if
+      zq(1:ngrid,1:nlay,1:nq)=pq(1:ngrid,1:nlay,1:nq)+
+     &                        pdqfi(1:ngrid,1:nlay,1:nq)*ptimestep
 
 c-----------------------------------------------------------------------
@@ -705 +696 @@
 c   TRACERS
 c   -------
-
-      if(tracer) then
 
 c     Using the wind modified by friction for lifting and  sublimation
@@ -1056 +1045 @@
 
         enddo ! of do iq=1,nq
-      end if ! of if(tracer)
 
 c --------- end of tracers  ----------------------------
@@ -1087 +1075 @@
       ENDDO
 
-      if (tracer) then 
-        DO iq = 1, nq
-          DO ilev = 1, nlay
-            DO ig=1,ngrid
-              pdqdif(ig,ilev,iq)=(zq(ig,ilev,iq)-
-     $      (pq(ig,ilev,iq) + pdqfi(ig,ilev,iq)*ptimestep))/ptimestep
-            ENDDO
-          ENDDO
-        ENDDO
-      end if
+      pdqdif(1:ngrid,1:nlay,1:nq)=(zq(1:ngrid,1:nlay,1:nq)-
+     &                             (pq(1:ngrid,1:nlay,1:nq)
+     &                              +pdqfi(1:ngrid,1:nlay,1:nq)
+     &                                    *ptimestep))/ptimestep
 
 c    ** diagnostique final