Changeset 673 for LMDZ4/trunk/libf/phylmd/ini_histrac.h

Timestamp:

Nov 16, 2005, 3:46:34 PM (19 years ago)

Author:

lmdzadmin

Message:

suppression de prt_flag_ts (INCA) AC
MAF

File:

• LMDZ4/trunk/libf/phylmd/ini_histrac.h (modified) (10 diffs)

LMDZ4/trunk/libf/phylmd/ini_histrac.h

@@ -21 +21 @@
      .                 nhori, nid_tra2)
 
-         CALL histbeg("histrac_inst", iim,zx_lon(:,1), jjm+1,zx_lat(1,:),
+         CALL histbeg("histrac_inst", iim,zx_lon, jjm+1,zx_lat,
      .                 1,iim,1,jjm+1,  itau_phy, zjulian, pdtphys,
      .                 nhori, nid_tra3)
@@ -27 +27 @@
          call histvert(nid_tra2, "presnivs", "Vertical levels", "mb",
      .                 klev, presnivs, nvert)
-
-         call histvert(nid_tra3, "presnivs", "Vertical levels", "mb",
+         call histvert(nid_tra3, "presnivs", "presnivs", "mb",
      .                 klev, presnivs, nvert)
 #endif
+
 
          zsto = pdtphys
@@ -84 +84 @@
 #ifdef INCA
 #ifdef INCAINFO
+
+#ifdef INCA_CH4
          DO it=1, phtcnt
          WRITE(str2,'(i2.2)') it
@@ -96 +98 @@
      .                "ave(X)", zsto,zout)
          ENDDO
-
          DO it=1, extcnt
          WRITE(str2,'(i2.2)') it
@@ -103 +104 @@
      .                "ave(X)", zsto,zout)
          ENDDO
-
+#endif
+#ifdef INCA_NMHC
+         DO it=1, phtcnt
+         WRITE(str2,'(i2.2)') it
+         CALL histdef(nid_tra, "j"//str2,"j"//str2, "CM-3 S-1",
+     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+         ENDDO
+         DO it=1, hetcnt
+         WRITE(str2,'(i2.2)') it
+         CALL histdef(nid_tra, "w"//str2,"w"//str2, "S-1",
+     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+         ENDDO
+         DO it=1, extcnt
+         WRITE(str2,'(i2.2)') it
+         CALL histdef(nid_tra, "ext"//str2,"ext"//str2, "CM-3 S-1",
+     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+         ENDDO
+#endif
          DO it=1, nfs
          WRITE(str2,'(i2.2)') it
@@ -112 +133 @@
 
 #else
+
+#ifdef INCA_NMHC
          CALL histdef(nid_tra, "jO3","jO3", "CM-3 S-1",
      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
@@ -131 +154 @@
      .                "ave(X)", zsto,zout)
 #endif
-
-         DO it=1, grpcnt
-         CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR",
-     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
-     .                "ave(X)", zsto,zout)
-         ENDDO
+#ifdef INCA_CH4
+         CALL histdef(nid_tra, "jO3","jO3", "CM-3 S-1",
+     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+         CALL histdef(nid_tra, "jNO2","jNO2", "CM-3 S-1",
+     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+         CALL histdef(nid_tra, "jH2O2","jH2O2", "CM-3 S-1",
+     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+         CALL histdef(nid_tra, "wHNO3","wHNO3", "S-1",
+     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+         CALL histdef(nid_tra, "kN2O5", "kN2O5","CM-3 S-1",
+     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+         CALL histdef(nid_tra, "LghtNO","LghtNO", "CM-3 S-1",
+     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+#endif
+#endif
+
+!         DO it=1, grpcnt
+!         CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR",
+!     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
+!     .                "ave(X)", zsto,zout)
+!         ENDDO
 #endif
 
@@ -159 +203 @@
 C champ 2D
 #ifdef INCA
-         IF ( prt_flag_ts(it) == 0 ) CYCLE
 
          CALL histdef(nid_tra, "Emi_"//solsym(it), "Emi_"//solsym(it),
@@ -341 +384 @@
 #endif
 
+#ifdef INCA_CH4
+         DO it=1, grpcnt
+         CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR",
+     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+         ENDDO
+
+#endif
+
+
 #ifdef INCA_NMHC
          CALL histdef(nid_tra, "CO2_basprod", "CO2_basprod", "CM-3 S-1",
@@ -363 +416 @@
      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
      .                "ave(X)", zsto,zout)
+
+         DO it=1, grpcnt
+         CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR",
+     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+         ENDDO
+
+
 #endif
 ! Fin ajout Anne