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

Timestamp:

May 25, 2005, 3:10:09 PM (19 years ago)

Author:

Laurent Fairhead

Message:

Synchronisation avec tous les diagnostiques de Ionela IM
Inclusion du slab ocean IM
LF

File:

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

LMDZ4/trunk/libf/phylmd/ini_histISCCP.h

@@ -7 +7 @@
 c zout et zsto.
 c dtime est passe par ailleurs a histbeg
-c zsto = frequence de stockage des champs
+c zstophy = frequence de stockage des champs tous les pdt physiques
 c zout = frequence d'ecriture des champs
-         zsto = dtime
+        zstophy = dtime
 c
 c ecriture 8 fois par jour
@@ -15 +15 @@
 c ecriture toutes les 2h (12 fois par jour)
 c       zout = dtime * 4.
-c ecriture toutes les 1/h (48 fois par jour)
+c ecriture toutes les 1/2 h (48 fois par jour)
 c       zout = dtime
 c
 c ecriture mensuelle
         zout = dtime * ecrit_mth
-c
-        print*,'ISCCP zout,zsto=',zout,zsto
+c       zout = dtime * ecrit_day
 c
 c       PRINT*, 'La frequence de sortie ISCCP est de ', ecrit_isccp
@@ -48 +47 @@
      .                 lmaxm1, cldtopres, nvert,'down')
 c
+        IF(type_run.EQ."ENSP".OR.type_run.EQ."CLIM") THEN
+c
 c variables a ecrire
 c
@@ -54 +55 @@
      .                "LMDZ ISCCP cld", "%",
      .                iim, jjmp1,nhori,lmaxm1,1,lmaxm1,nvert,32,
-     .                "ave(X)", zsto,zout)
+     .                "ave(X)", zstophy,zout)
          ENDDO
 c
@@ -60 +61 @@
      .                "Nb of calls with sunlit ", "%",
      .                iim, jjmp1,nhori,1,1,1,-99,32,
-     .                "ave(X)", zsto,zout)
+     .                "ave(X)", zstophy,zout)
+c
+        ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
+c
+         DO k=1, kmaxm1
+          DO l=1, lmaxm1
+c
+          CALL histdef(nid_isccp, pclev(l)//taulev(k),
+     .                "LMDZ ISCCP cld "//cnameisccp(l,k), "%",
+     .                iim, jjmp1,nhori,1,1,1,-99,32,
+     .                "ave(X)", zstophy,zout)
+c
+          ENDDO
+         ENDDO
+c
+         CALL histdef(nid_isccp, "nsunlit",
+     .                "Nb of calls with sunlit ", "%",
+     .                iim, jjmp1,nhori,1,1,1,-99,32,
+     .                "ave(X)", zstophy,zout)
+c
+        ENDIF 
 c
         CALL histend(nid_isccp)