Changeset 850 for trunk/LMDZ.GENERIC/libf

Timestamp:

Nov 23, 2012, 3:53:54 PM (12 years ago)

Author:

aslmd

Message:

LMDZ.GENERIC. bug fix for start2archive (there was a problem because initracer was not called). also corrected iniadvtrac so that it could be compiled with -t 0. compiler was complaining about lines with iadv(nqmx). those lines are actually not necessary (this could be included in the loop above). added an if statement in iniadvtrac so that traceur.def does not even need to be here when the model is compiled with -t 0

Location:

trunk/LMDZ.GENERIC/libf/dyn3d

Files:

• iniadvtrac.F (modified) (4 diffs)
• start2archive.F (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/dyn3d/iniadvtrac.F

@@ -21 +21 @@
       INTEGER :: ierr
 
-#ifdef NOTRAC
+
+      if (nqmx > 0) then
 
 ! Look for file traceur.def
@@ -79 +80 @@
 c
 c
-      DO  iq = 1, nqmx-1
+      DO  iq = 1, nqmx
        IF( iadv(iq).EQ.1 ) PRINT *,' Choix du shema humidite specifique'
      * ,' pour le traceur no ', iq
@@ -102 +103 @@
       ENDDO
 
-       IF( iadv(nqmx).EQ.1 ) PRINT *,' Choix du shema humidite '
-     * ,'specifique pour la vapeur d''eau'
-       IF( iadv(nqmx).EQ.2 ) PRINT *,' Choix du shema  amont',' pour la'
-     * ,' vapeur d''eau '
-       IF( iadv(nqmx).EQ.3 ) PRINT *,' Choix du shema  Van-Leer '
-     * ,' pour la vapeur d''eau'
-       IF( iadv(nqmx).EQ.4 ) PRINT *,' Choix du shema  Van-Leer + '
-     * ,' humidite specifique pour la vapeur d''eau'
-c
-!       IF( (iadv(nqmx).LE.0).OR.(iadv(nqmx).GT.4) )   THEN
-!     MODIFICATION TO TEST WITHOUT TRACER ADVECTION BY RDW
-       IF( (iadv(nqmx).LT.0).OR.(iadv(nqmx).GT.4) )   THEN
-        PRINT *,' Erreur dans le choix de iadv (nqmx).Corriger et '
-     * ,' repasser car  iadv(nqmx) = ', iadv(nqmx)
-         CALL ABORT
-       ENDIF
+!!!! AS: compiler complains about iadv(nqmx) when there is nqmx=0
+!!!! AS: so I commented those lines and changed nqmx-1 for nqmx above
+!       IF( iadv(nqmx).EQ.1 ) PRINT *,' Choix du shema humidite '
+!     * ,'specifique pour la vapeur d''eau'
+!       IF( iadv(nqmx).EQ.2 ) PRINT *,' Choix du shema  amont',' pour la'
+!     * ,' vapeur d''eau '
+!       IF( iadv(nqmx).EQ.3 ) PRINT *,' Choix du shema  Van-Leer '
+!     * ,' pour la vapeur d''eau'
+!       IF( iadv(nqmx).EQ.4 ) PRINT *,' Choix du shema  Van-Leer + '
+!     * ,' humidite specifique pour la vapeur d''eau'
+!
+!c
+!!       IF( (iadv(nqmx).LE.0).OR.(iadv(nqmx).GT.4) )   THEN
+!!     MODIFICATION TO TEST WITHOUT TRACER ADVECTION BY RDW
+!       IF( (iadv(nqmx).LT.0).OR.(iadv(nqmx).GT.4) )   THEN
+!        PRINT *,' Erreur dans le choix de iadv (nqmx).Corriger et '
+!     * ,' repasser car  iadv(nqmx) = ', iadv(nqmx)
+!         CALL ABORT
+!       ENDIF
 
       first = .TRUE.
@@ -138 +142 @@
       ENDDO
 c
-
-#endif
+      endif ! of if nqmx > 0
 
       end

trunk/LMDZ.GENERIC/libf/dyn3d/start2archive.F

@@ -137 +137 @@
 ! Load tracer names:
       call iniadvtrac(nq,numvanle)
+! Initialize global tracer indexes (stored in tracer.h)
+! ... this has to be done before phyetat0
+      call initracer(ngridmx,nqmx,tnom)
 
       fichnom = 'start.nc'