Changeset 827 for LMDZ4/trunk/libf/phylmd/ini_histISCCP.h

Timestamp:

Aug 28, 2007, 2:17:59 PM (17 years ago)

Author:

lmdzadmin

Message:

Correction bug sorties ISCCP (7pres X 7tau) + ENSEMBLES + ajout freq_ISCCP/ecrit_ISCCP
IM

File:

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

LMDZ4/trunk/libf/phylmd/ini_histISCCP.h

@@ -15 +15 @@
 cIM 300505     zstophy = dtime 
 c appel du simulateur toutes les 3heures
-          zcals(1) = dtime *6.  !toutes les 3h (en s)
+!IM on lit la frequence d'appel dans physiq.def
+!         zcals(1) = dtime *6.  !toutes les 3h (en s)
+          zcals(1) = freq_ISCCP !toutes les freq_ISCCP secondes
         DO n=1, napisccp
           zcalh(n) = zcals(n)/3600. !stoutes les Xh (en heures)
@@ -29 +31 @@
 c       IF(freqout_isccp.EQ.1.) THEN
 c ecriture jounaliere
-          zout_isccp(1) = ecrit_day !(en s)
+!IM on ecrit les resultats du simulateur ISCCP toutes les 
+! ecrit_ISCCP secondes      zout_isccp(1) = ecrit_day !(en s)
+          zout_isccp(1) = ecrit_ISCCP !(en s)
 c ecriture mensuelle
 c         zout = dtime * ecrit_mth !(en s)
@@ -37 +41 @@
 c le nombre de sous-colonnes ncol : ncol=(100.*zcalh)/zoutd
           ncol(n)=NINT((100.*zcalh(n))/zoutj(n))
-          PRINT*,'n ncol',n,ncol(n)
+          IF(ncol(n).GT.ncolmx) THEN
+           PRINT*,'Warning: Augmenter le nombre colonnes du simulateur'
+           PRINT*,'         ISCCP ncol=', ncol,' ncolmx=',ncolmx
+c          PRINT*,'n ncol',n,ncol(n)
+           CALL abort
+          ENDIF
 c
         DO l=1, ncol(n)
@@ -80 +89 @@
 cym          enddo
 c
-          print*,'On stoke le fichier histISCCP instantanne sur ',
+          print*,'On stoke le fichier histISCCP sur ',
      s   imin_ins,imax_ins,jmin_ins,jmax_ins
 cym          print*,'On stoke le fichier histISCCP instantanne sur ',
@@ -94 +103 @@
          CALL histbeg_phy("histISCCP.nc", itau_phy, zjulian, dtime,
      .                 nhori, nid_isccp)
-
         ENDIF !(1.EQ.0) THEN
 c
@@ -106 +114 @@
          CALL histbeg_phy("histISCCP.nc", itau_phy, zjulian, dtime,
      .                 nhori, nid_isccp)
-
 c
         IF(type_run.EQ."ENSP".OR.type_run.EQ."CLIM") THEN
@@ -128 +135 @@
 c
          DO k=1, kmaxm1
-          CALL histdef(nid_isccp, "cldISCCP_"//taulev(k)//typinout(n),
+          CALL histdef(nid_isccp, "cldISCCP_"//taulev(k)//verticaxe(n),
      .                "LMDZ ISCCP cld", "%",
      .                iim, jj_nb,nhori,lmaxm1,1,lmaxm1,nvert,32,
@@ -134 +141 @@
          ENDDO
 c
-         CALL histdef(nid_isccp, "nsunlit"//typinout(n),
+         CALL histdef(nid_isccp, "nsunlit"//verticaxe(n),
      .                "Nb of calls with sunlit ", "%",
      .                iim, jj_nb,nhori,1,1,1,-99,32,
@@ -150 +157 @@
            DO l=1, lmaxm1
 c
-           CALL histdef(nid_isccp, pclev(l)//taulev(k)//typinout(n),
+           CALL histdef(nid_isccp, pclev(l)//taulev(k)//verticaxe(n),
      .                "LMDZ ISCCP cld "//cnameisccp(l,k), "%",
      .                iim, jj_nb,nhori,1,1,1,-99,32,
@@ -159 +166 @@
 c
 c         print*,'n=',n,' avant histdef(..Nb of calls sunlit'
-          CALL histdef(nid_isccp, "nsunlit"//typinout(n),
+          CALL histdef(nid_isccp, "nsunlit"//verticaxe(n),
      .                "Nb of calls with sunlit ", "%",
      .                iim, jj_nb,nhori,1,1,1,-99,32,