source: LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histmth.h @ 2063

Last change on this file since 2063 was 740, checked in by lmdzadmin, 18 years ago

Correction bogues: les ecrit_ sont des REALs lus dans conf_phys.F90 en
nombre de jours sauf pour ecrit_ins et ecrit_tra en secondes!
Les ecrit_ sont initialises dans conf_phys.F90 et peuvent etre modifies dans
physiq.def.
IM

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 68.2 KB
RevLine 
[644]1c
2c $Header$
3c
[524]4      IF (ok_mensuel) THEN
5c
[644]6       zstophy = dtime
[740]7       zstomth = ecrit_mth                !valeurs mensuelles "instantannees"
[644]8       zstorad = dtime * FLOAT(radpas)    !variables stockees toutes les 2h
[740]9       zstohf = ecrit_hf                  !variables stockees toutes les 6h
10       zout = ecrit_mth                   !frequence des sorties : 1 val. par mois
[524]11c
[740]12c     PRINT*,' La zstophy,zstorad,zstohf,zstomth,zout=',
13c   . zstophy,zstorad,zstohf,zstomth,zout
[524]14c
15         idayref = day_ref
16         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
17c
18         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
19         DO i = 1, iim
20            zx_lon(i,1) = rlon(i+1)
21            zx_lon(i,jjmp1) = rlon(i+1)
22         ENDDO
[644]23c
24cIM cf. AM 081204 BEG region
25c         imin_histmth=1
26c         imax_histmth=iim
27c         jmin_histmth=1
28c         jmax_histmth=jjmp1
29c         do i=1,iim-1
30c            if(zx_lon(i,1).lt.lonmin_histmth) imin_histmth=i
31c            if(zx_lon(i,1).le.lonmax_histmth) imax_histmth=i+1
32c         enddo
33c         do j=1,jjmp1-1
34c            if(zx_lat(1,j).ge.latmin_histmth) jmax_histmth=j
35c            if(zx_lat(1,j).gt.latmax_histmth) jmin_histmth=j
36c         enddo
37c
38c         print*,'On stoke le fichier instantanne sur, ',
39c     s   imin_histmth,imax_histmth,jmin_histmth,jmax_histmth
40c         print*,'On stoke le fichier instantanne sur, ',
41c     s   zx_lon(imin_histmth,1),zx_lon(imax_histmth,1),
42c     s   zx_lat(1,jmin_histmth),zx_lat(1,jmax_histmth)
43cIM cf. AM 081204 END region
44c
[524]45         DO ll=1,klev
46            znivsig(ll)=float(ll)
47         ENDDO
48         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
49         CALL histbeg("histmth.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
50     .                 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 
51     .                 nhori, nid_mth)
[644]52c
53cIM cf. AM 081204 BEG region
54c         CALL histbeg("histmth.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
55c     .                 imin_histmth,imax_histmth-imin_histmth+1
56c     .                 ,jmin_histmth,jmax_histmth-jmin_histmth+1
57c     .                 , itau_phy, zjulian, dtime,
58c     .                 nhori, nid_mth)
59cIM 081204 END
60c
[524]61         CALL histvert(nid_mth, "presnivs", "Vertical levels", "mb",
[644]62     .                 klev, presnivs/100., nvert)
[524]63c
[644]64      IF(type_run.EQ."CLIM".OR.type_run.EQ."ENSP") THEN
[524]65c
66      IF(lev_histmth.GE.1) THEN
[644]67c
68         CALL histdef(nid_mth, "phis", 
69     .                "Surface geop. height", "m2/s2",
[524]70     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]71     .                "once",  zstophy,zout)
[524]72c
73         CALL histdef(nid_mth, "aire", "Grid area", "-",
74     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]75     .                "once",  zstophy,zout)
[524]76c
[644]77         CALL histdef(nid_mth, "pourc_"//clnsurf(is_ter),
78     .         "% "//clnsurf(is_ter), "%", 
79     .         iim,jjmp1,nhori, 1,1,1, -99, 32,
80     .         "once", zstophy,zout)
[524]81c
[644]82         CALL histdef(nid_mth, "fract_"//clnsurf(is_ter),
83     .         "Fraction "//clnsurf(is_ter), "1", 
84     .         iim,jjmp1,nhori, 1,1,1, -99, 32,
85     .         "once", zstophy,zout)
86c
[524]87c Champs 2D:
88c
89         CALL histdef(nid_mth, "slp", "Sea Level Pressure", "Pa",
90     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]91     .                "ave(X)", zstophy,zout)
[524]92c
93         CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
94     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]95     .                "ave(X)", zstophy,zout)
[524]96c
97         CALL histdef(nid_mth, "t2m", "Temperature 2m", "K",
98     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]99     .                "ave(X)", zstophy,zout)
100c
[524]101c ENSEMBLES BEG
102         CALL histdef(nid_mth, "t2m_min", "Temp. 2m min.",
103     .                "K",
104     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]105     .                t2mincels, zstophy,zout)
[524]106c
107         CALL histdef(nid_mth, "t2m_max", "Temp. 2m max.",
108     .                "K",
109     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]110     .                t2maxcels, zstophy,zout)
[524]111c
112c        CALL histdef(nid_mth, "tsoil", "Sfce soil Temperature",
113c    .                "K",
114c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]115c    .                "ave(X)", zstophy,zout)
[524]116c
117         CALL histdef(nid_mth, "wind10m","10-m wind speed","m/s",
118     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]119     .                "ave(X)", zstophy,zout)
[524]120c
121         CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "-",
122     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]123     .                "ave(X)", zstophy,zout)
[524]124c
125c ENSEMBLES END
126         CALL histdef(nid_mth, "q2m", "Specific humidity 2m", "kg/kg",
127     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]128     .                "ave(X)", zstophy,zout)
[524]129c
130         CALL histdef(nid_mth, "u10m", "Vent zonal 10m", "m/s",
131     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]132     .                "ave(X)", zstophy,zout)
[524]133c
134         CALL histdef(nid_mth, "v10m", "Vent meridien 10m", "m/s",
135     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]136     .                "ave(X)", zstophy,zout)
[524]137c
138         CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
139     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]140     .                "ave(X)", zstophy,zout)
[524]141c
142         CALL histdef(nid_mth, "qsurf", "Surface Air humidity", "kg/kg",
143     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]144     .                "ave(X)", zstophy,zout)
[524]145c
146      if (.not. ok_veget) then
147         CALL histdef(nid_mth, "qsol", "Soil watter content", "mm",
148     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]149     .                "ave(X)", zstophy,zout)
[524]150      endif
151c
152         CALL histdef(nid_mth, "ndayrain", 
153     .                "Number of day with rain (liq+sol)", "-",
154     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]155     .                "inst(X)", zstomth,zout)
[524]156c
157         CALL histdef(nid_mth, "precip", "Precipitation Totale liq+sol", 
158     .                "kg/(s*m2)",
159     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]160     .                "ave(X)", zstophy,zout)
[524]161c
162         CALL histdef(nid_mth, "plul", "Large-scale Precip.", 
163     .   "kg/(s*m2)",
164     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]165     .                "ave(X)", zstophy,zout)
[524]166c
167         CALL histdef(nid_mth, "pluc", "Convective Precip.", 
168     .   "kg/(s*m2)",
169     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]170     .                "ave(X)", zstophy,zout)
[524]171c
172         CALL histdef(nid_mth, "snow", "Snow fall", "kg/(s*m2)",
173     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]174     .                "ave(X)", zstophy,zout)
[524]175c
[644]176c        CALL histdef(nid_mth, "snow_mass", "Snow Mass", "kg/m2",
177c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
178c    .                "ave(X)", zstophy,zout)
[524]179c
180         CALL histdef(nid_mth, "evap", "Evaporation", "kg/(s*m2)",
181     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]182     .                "ave(X)", zstophy,zout)
[524]183c
184         CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
185     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]186     .                "ave(X)", zstophy,zout)
[524]187c
[644]188         CALL histdef(nid_mth, "tops0", "CS Solar rad. at TOA", "W/m2",
[524]189     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]190     .                "ave(X)", zstophy,zout)
[524]191c
192         CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
193     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]194     .                "ave(X)", zstophy,zout)
[524]195c
196         CALL histdef(nid_mth, "topl0", "IR rad. at TOA", "W/m2",
197     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]198     .                "ave(X)", zstophy,zout)
[524]199c
200         CALL histdef(nid_mth, "SWupTOA", "SWup at TOA","W/m2",
201     .                iim,jjmp1,nhori, 1,1,1,-99,
[644]202     .                32, "ave(X)", zstophy,zout)
[524]203c
204         CALL histdef(nid_mth, "SWupTOAclr", 
205     .                "SWup clear sky at TOA","W/m2",
206     .                iim,jjmp1,nhori, 1,1,1,-99,
[644]207     .                32, "ave(X)", zstophy,zout)
[524]208c
209         CALL histdef(nid_mth, "SWdnTOA", "SWdn at TOA","W/m2",
210     .                iim,jjmp1,nhori, 1,1,1,-99,
[644]211     .                32, "ave(X)", zstophy,zout)
[524]212c
213         CALL histdef(nid_mth, "SWdnTOAclr", 
214     .                "SWdn clear sky at TOA","W/m2",
215     .                iim,jjmp1,nhori, 1,1,1,-99,
[644]216     .                32, "ave(X)", zstophy,zout)
[524]217c
[644]218         CALL histdef(nid_mth, "SWup200", "SWup at 200mb","W/m2",
[524]219     .                iim,jjmp1,nhori, 1,1,1,-99,
[644]220     .                32, "ave(X)", zstophy,zout)
[524]221c
222         CALL histdef(nid_mth, "SWup200clr", 
[644]223     .                "SWup clear sky at 200mb","W/m2",
[524]224     .                iim,jjmp1,nhori, 1,1,1,-99,
[644]225     .                32, "ave(X)", zstophy,zout)
[524]226c
[644]227         CALL histdef(nid_mth, "SWdn200", "SWdn at 200mb","W/m2",
[524]228     .                iim,jjmp1,nhori, 1,1,1,-99,
[644]229     .                32, "ave(X)", zstophy,zout)
[524]230c
231         CALL histdef(nid_mth, "SWdn200clr", 
[644]232     .                "SWdn clear sky at 200mb","W/m2",
[524]233     .                iim,jjmp1,nhori, 1,1,1,-99,
[644]234     .                32, "ave(X)", zstophy,zout)
[524]235c
[644]236         CALL histdef(nid_mth, "LWup200", "LWup at 200mb","W/m2",
[524]237     .                iim,jjmp1,nhori, 1,1,1,-99,
[644]238     .                32, "ave(X)", zstophy,zout)
[524]239c
240         CALL histdef(nid_mth, "LWup200clr", 
[644]241     .                "LWup clear sky at 200mb","W/m2",
[524]242     .                iim,jjmp1,nhori, 1,1,1,-99,
[644]243     .                32, "ave(X)", zstophy,zout)
[524]244c
[644]245         CALL histdef(nid_mth, "LWdn200", "LWdn at 200mb","W/m2",
[524]246     .                iim,jjmp1,nhori, 1,1,1,-99,
[644]247     .                32, "ave(X)", zstophy,zout)
[524]248c
249         CALL histdef(nid_mth, "LWdn200clr", 
[644]250     .                "LWdn clear sky at 200mb","W/m2",
[524]251     .                iim,jjmp1,nhori, 1,1,1,-99,
[644]252     .                32, "ave(X)", zstophy,zout)
[524]253c
254         CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
255     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]256     .                "ave(X)", zstophy,zout)
[524]257c
258         CALL histdef(nid_mth, "sols0", "Solar rad. at surf.", "W/m2",
259     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]260     .                "ave(X)", zstophy,zout)
[524]261c
262         CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
263     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]264     .                "ave(X)", zstophy,zout)
[524]265c
266         CALL histdef(nid_mth, "soll0", "IR rad. at surface", "W/m2",
267     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]268     .                "ave(X)", zstophy,zout)
[524]269c
270         CALL histdef(nid_mth, "SWupSFC", "SWup at surface","W/m2",
271     .                iim,jjmp1,nhori, 1,1,1,-99,
[644]272     .                32, "ave(X)", zstophy,zout)
[524]273c
274         CALL histdef(nid_mth, "SWupSFCclr", 
275     .                "SWup clear sky at surface","W/m2",
276     .                iim,jjmp1,nhori, 1,1,1,-99,
[644]277     .                32, "ave(X)", zstophy,zout)
[524]278c
279         CALL histdef(nid_mth, "SWdnSFC", "SWdn at surface","W/m2",
280     .                iim,jjmp1,nhori, 1,1,1,-99,
[644]281     .                32, "ave(X)", zstophy,zout)
[524]282c
283         CALL histdef(nid_mth, "SWdnSFCclr", 
284     .                "SWdn clear sky at surface","W/m2",
285     .                iim,jjmp1,nhori, 1,1,1,-99,
[644]286     .                32, "ave(X)", zstophy,zout)
[524]287c
288         CALL histdef(nid_mth,"LWupSFC","Upwd. IR rad. at surface", 
289     .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]290     .                "ave(X)", zstophy,zout)
[524]291c
292         CALL histdef(nid_mth,"LWdnSFC","Down. IR rad. at surface", 
293     .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]294     .                "ave(X)", zstophy,zout)
[524]295c
296         CALL histdef(nid_mth,"LWupSFCclr",
297     .                "CS Upwd. IR rad. at surface", 
298     .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]299     .                "ave(X)", zstophy,zout)
[524]300c
301         CALL histdef(nid_mth,"LWdnSFCclr",
302     .                "Down. CS IR rad. at surface", 
303     .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]304     .                "ave(X)", zstophy,zout)
[524]305c
306         CALL histdef(nid_mth, "bils", "Surf. total heat flux", "W/m2",
307     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]308     .                "ave(X)", zstophy,zout)
[524]309c
310         CALL histdef(nid_mth, "sens", "Sensible heat flux", "W/m2",
311     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]312     .                "ave(X)", zstophy,zout)
[524]313c
314         CALL histdef(nid_mth, "fder", "Heat flux derivation", "W/m2",
315     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]316     .                "ave(X)", zstophy,zout)
[524]317c
318c          CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
319c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]320c    .                "ave(X)", zstophy,zout)
[524]321c
322c        CALL histdef(nid_mth, "frtv", "Meridional wind stress", "Pa",
323c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]324c    .                "ave(X)", zstophy,zout)
[524]325c
326         CALL histdef(nid_mth, "ffonte","Thermal flux for snow melting",
327     .                "W/m2",iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]328     .                "ave(X)", zstophy,zout)
[524]329c
330         CALL histdef(nid_mth, "fqcalving","Ice Calving",
331     .                "kg/m2/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]332     .                "ave(X)", zstophy,zout)
333c
[524]334         DO nsrf = 1, nbsrf
335           call histdef(nid_mth, "taux_"//clnsurf(nsrf),
336     $         "Zonal wind stress"//clnsurf(nsrf), "Pa", 
337     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]338     $         "ave(X)", zstophy,zout)
[524]339
340           call histdef(nid_mth, "tauy_"//clnsurf(nsrf),
[644]341     $         "Meridional wind stress "//clnsurf(nsrf), "Pa", 
[524]342     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]343     $         "ave(X)", zstophy,zout)
[524]344         ENDDO
345c
346         DO nsrf = 1, nbsrf
347C
[644]348          IF(nsrf.GT.2) THEN
[524]349           call histdef(nid_mth, "pourc_"//clnsurf(nsrf),
[644]350     $         "% "//clnsurf(nsrf), "%", 
[524]351     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]352     $         "ave(X)", zstophy,zout)
353c
354           call histdef(nid_mth, "fract_"//clnsurf(nsrf),
355     $         "Fraction "//clnsurf(nsrf), "1", 
356     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
357     $         "ave(X)", zstophy,zout)
358          ENDIF !nsrf.GT.2
[524]359C
360           call histdef(nid_mth, "tsol_"//clnsurf(nsrf),
[644]361     $         "Temperature "//clnsurf(nsrf), "K", 
[524]362     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]363     $         "ave(X)", zstophy,zout)
[524]364C
365           call histdef(nid_mth, "sens_"//clnsurf(nsrf),
366     $         "Sensible heat flux "//clnsurf(nsrf), "W/m2", 
367     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]368     $         "ave(X)", zstophy,zout)
[524]369c
370           call histdef(nid_mth, "lat_"//clnsurf(nsrf),
371     $         "Latent heat flux "//clnsurf(nsrf), "W/m2", 
372     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]373     $         "ave(X)", zstophy,zout)
[524]374C
375           call histdef(nid_mth, "flw_"//clnsurf(nsrf),
376     $         "LW "//clnsurf(nsrf), "W/m2",
377     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]378     $         "ave(X)", zstophy,zout)
[524]379c
380           call histdef(nid_mth, "fsw_"//clnsurf(nsrf),
381     $         "SW "//clnsurf(nsrf), "W/m2",
382     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]383     $         "ave(X)", zstophy,zout)
[524]384C
385           call histdef(nid_mth, "wbils_"//clnsurf(nsrf),
386     $         "Bilan sol "//clnsurf(nsrf), "W/m2",
387     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]388     $         "ave(X)", zstophy,zout)
[524]389C
[644]390cIM 250105 BEG
391          IF(1.EQ.0) THEN
392C
393           call histdef(nid_mth, "evapour_"//clnsurf(nsrf),
394     $         "Evap*pourc "//clnsurf(nsrf), "kg/(s*m2)",
395     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
396     $         "ave(X)", zstophy,zout)
397C
398           call histdef(nid_mth, "prepour_"//clnsurf(nsrf),
399     $         "Precip*pourc "//clnsurf(nsrf), "kg/(s*m2)",
400     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
401     $         "ave(X)", zstophy,zout)
402C
403          ENDIF !(1.EQ.0) THEN
404cIM 250105 END
405C
[524]406         END DO
407c
408         CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
409     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]410     .                "ave(X)", zstophy,zout)
[524]411c
412         CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
413     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]414     .                "ave(X)", zstophy,zout)
[524]415c
416         CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-",
417     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]418     .                "ave(X)", zstophy,zout)
[524]419c
420         CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-",
421     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]422     .                "ave(X)", zstophy,zout)
[524]423c
424         CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-",
425     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]426     .                "ave(X)", zstophy,zout)
[524]427c
428         CALL histdef(nid_mth, "cldt", "Total cloudiness", "-",
429     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]430     .                "ave(X)", zstophy,zout)
[524]431c
432         CALL histdef(nid_mth,"cldq","Cloud liquid water path","kg/m2",
433     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]434     .                "ave(X)", zstophy,zout)
435c
[524]436         CALL histdef(nid_mth,"lwp","Cloud water path","kg/m2",
437     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]438     .                "ave(X)", zstophy,zout)
[524]439c
440         CALL histdef(nid_mth,"iwp","Cloud ice water path","kg/m2",
441     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]442     .                "ave(X)", zstophy,zout)
[524]443c
444         CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
445     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]446     .                "ave(X)", zstophy,zout)
[524]447c
448         CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
449     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]450     .                "ave(X)", zstophy,zout)
[524]451c
452         CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-",
453     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]454     .                "ave(X)", zstophy,zout)
[524]455c
456         CALL histdef(nid_mth, "vq", "Merid humidity transport", "-",
457     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]458     .                "ave(X)", zstophy,zout)
[524]459cKE43
460      IF(iflag_con.GE.3) THEN ! sb
461c
462         CALL histdef(nid_mth, "cape", "Conv avlbl pot ener", "J/kg",
463     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]464     .                "ave(X)", zstophy,zout)
[524]465c
[644]466         CALL histdef(nid_mth, "pbase", "Cld base pressure", "mb",
[524]467     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]468     .                "ave(X)", zstophy,zout)
[524]469c
[644]470         CALL histdef(nid_mth, "ptop", "Cld top pressure", "mb",
[524]471     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]472     .                "ave(X)", zstophy,zout)
[524]473c
474         CALL histdef(nid_mth, "fbase", "Cld base mass flux", "kg/m2/s",
475     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]476     .                "ave(X)", zstophy,zout)
[524]477c
478         CALL histdef(nid_mth, "prw", "Precipitable water", "kg/m2",
479     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]480     .                "ave(X)", zstophy,zout)
[524]481c
[644]482      ENDIF !iflag_con .GE. 3
[524]483c34EK
[644]484cIM cf. AM 081204 BEG
485c34EK
486c HBTM2
487         CALL histdef(nid_mth, "s_pblh", "Boundary Layer Height", "m",
488     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
489     .                "ave(X)", zstophy,zout)
490         CALL histdef(nid_mth, "s_pblt", "t at Boundary Layer Height",
491     .                "K",
492     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
493     .                "ave(X)", zstophy,zout)
494         CALL histdef(nid_mth, "s_lcl", "Condensation level", "m",
495     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
496     .                "ave(X)", zstophy,zout)
497         CALL histdef(nid_mth, "s_capCL", "Conv avlbl pot ener
498     .                for ABL", "J/m2",
499     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
500     .                "ave(X)", zstophy,zout)
501         CALL histdef(nid_mth, "s_oliqCL", "Liq Water in BL", "kg/m2",
502     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
503     .                "ave(X)", zstophy,zout)
504         CALL histdef(nid_mth, "s_cteiCL", "Instability criteria
505     .                (ABL)", "K",
506     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
507     .                "ave(X)", zstophy,zout)
508         CALL histdef(nid_mth, "s_therm", "Exces du thermique", "K",
509     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
510     .                "ave(X)", zstophy,zout)
511         CALL histdef(nid_mth, "s_trmb1", "deep_cape(HBTM2)", "J/m2",
512     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
513     .                "ave(X)", zstophy,zout)
514         CALL histdef(nid_mth, "s_trmb2", "inhibition (HBTM2)", "J/m2",
515     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
516     .                "ave(X)", zstophy,zout)
517         CALL histdef(nid_mth, "s_trmb3", "Point Omega (HBTM2)", "m",
518     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
519     .                "ave(X)", zstophy,zout)
[524]520c
[644]521cIM cf. AM 081204 END
522c
[524]523c Champs interpolles sur des niveaux de pression
[644]524         DO k=1, nlevSTD
525c
526c        bb=clevSTD(k)
527c
528c        IF(k.GE.2) THEN
529         IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
530         IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
531c         aa=clevSTD(k)
532c         bb=aa(1:lnblnk1(aa))
533c        ENDIF
534c
535      IF(bb2.EQ."850".OR.bb2.EQ."700".OR.
536     $   bb2.EQ."500".OR.bb2.EQ."200") THEN
537c
538          CALL histdef(nid_mth, "u"//bb2,
539     .                 "Zonal wind "//bb2//"mb","m/s",
[524]540     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]541     .                "inst(X)", zout,zout)
[524]542c
[644]543          CALL histdef(nid_mth, "v"//bb2,
544     .                 "Meridional wind "//bb2//"mb","m/s",
[524]545     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]546     .                "inst(X)", zout,zout)
[524]547c
[644]548          CALL histdef(nid_mth, "w"//bb2,
549     .                 "Vertical wind "//bb2//"mb","m/s",
[524]550     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]551     .                "inst(X)", zout,zout)
[524]552c
[644]553          CALL histdef(nid_mth, "phi"//bb2,
554     .                 "Geopotential "//bb2//"mb","m",
[524]555     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]556     .                "inst(X)", zout,zout)
[524]557c
[644]558        ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.
559       ENDDO
560c
561cIM diagnostiques flux ocean-atm ou ocean-glace de mer
562cIM pour utilisation dans un modele de "slab" ocean
563c
564       CALL histdef(nid_mth, "fluxo",
565     .              "Flux turbulents ocean-atmosphere", "W/m2",
566     .              iim,jjmp1,nhori, 1,1,1, -99, 32,
567     .              "ave(X)", zstophy,zout)
568c
569       CALL histdef(nid_mth, "fluxg",
570     .              "Flux turbulents ocean-glace de mer","W/m2",
571     .              iim,jjmp1,nhori, 1,1,1, -99, 32,
572     .              "ave(X)", zstophy,zout)
573c
574       IF (OCEAN.EQ.'force ') THEN
575c
576        CALL histdef(nid_mth, "lmt_bils",
577     .       "Bilan au sol atmosphere forcee", "W/m2",
578     .              iim,jjmp1,nhori, 1,1,1, -99, 32,
579     .              "ave(X)", zstophy,zout)
580c
581       ELSE IF (OCEAN.EQ.'slab  ') THEN
582c
583        CALL histdef(nid_mth, "slab_bils",
584     .       "Bilan au sol Slab", "W/m2",
585     .              iim,jjmp1,nhori, 1,1,1, -99, 32,
586     .              "ave(X)", zstophy,zout)
587c
588        CALL histdef(nid_mth, "tslab", "Slab SST ", "K",
589     .              iim,jjmp1,nhori, 1,1,1, -99, 32,
590     .              "ave(X)", zstophy,zout)
591c
592        CALL histdef(nid_mth, "seaice", "Slab seaice", "kg/m2",
593     .              iim,jjmp1,nhori, 1,1,1, -99, 32,
594     .              "ave(X)", zstophy,zout)
595c
596        CALL histdef(nid_mth, "siceh", "Slab seaice height", "m",
597     .              iim,jjmp1,nhori, 1,1,1, -99, 32,
598     .              "ave(X)", zstophy,zout)
599c
600       ENDIF
601c
[524]602      ENDIF !lev_histmth.GE.1
603c
604      IF(lev_histmth.GE.2) THEN
605c
606c Champs 3D:
607c
608c
[644]609         CALL histdef(nid_mth,"lwcon","Cloud liquid water content",
610     .                "kg/kg",
[524]611     .                iim,jjmp1,nhori, klev,1,klev, nvert, 32, 
[644]612     .                "ave(X)", zstophy,zout)
[524]613c
614         CALL histdef(nid_mth,"iwcon","Cloud ice water content","kg/kg",
615     .                iim,jjmp1,nhori, klev,1,klev, nvert, 32, 
[644]616     .                "ave(X)", zstophy,zout)
[524]617c
618         CALL histdef(nid_mth, "temp", "Air temperature", "K",
619     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]620     .                "ave(X)", zstophy,zout)
[524]621c
622         CALL histdef(nid_mth, "ovap", "Specific humidity", "kg/kg",
623     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]624     .                "ave(X)", zstophy,zout)
[524]625c
626c        CALL histdef(nid_mth,"wvap","Water vapor mixing ratio","kg/kg",
627c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]628c    .                "ave(X)", zstophy,zout)
[524]629c
[644]630         CALL histdef(nid_mth, "geop", "Geopotential height", "m2/s2",
[524]631     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]632     .                "ave(X)", zstophy,zout)
[524]633c
634         CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
635     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]636     .                "ave(X)", zstophy,zout)
[524]637c
638         CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
639     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]640     .                "ave(X)", zstophy,zout)
[524]641c
[644]642         CALL histdef(nid_mth, "vitw", "Vertical wind", "m/s",
[524]643     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]644     .                "ave(X)", zstophy,zout)
[524]645c
646         CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
647     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]648     .                "ave(X)", zstophy,zout)
[524]649c
650        CALL histdef(nid_mth, "rneb", "Cloud fraction", "-",
651     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]652     .                "ave(X)", zstophy,zout)
[524]653c
654         CALL histdef(nid_mth, "rnebcon", "Convective Cloud Fraction"
655     .                , "-",
656     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]657     .                "ave(X)", zstophy,zout)
[524]658c
659         CALL histdef(nid_mth, "rhum", "Relative humidity", "-",
660     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]661     .                "ave(X)", zstophy,zout)
[524]662c
[644]663         CALL histdef(nid_mth, "ozone", "Ozone concentration", "ppmv",
[524]664     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]665     .                "ave(X)", zstophy,zout)
[524]666c
[644]667         CALL histdef(nid_mth, "upwd", "saturated updraft", "kg/m2/s",
668     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
669     .                "ave(X)", zstophy,zout)
670c
[524]671         CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
672     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]673     .                "ave(X)", zstophy,zout)
[524]674c
[644]675         CALL histdef(nid_mth, "dqphy", "Physics dQ", "(kg/kg)/s",
[524]676     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]677     .                "ave(X)", zstophy,zout)
[524]678c
679      ENDIF !lev_histmth.GE.2
680c
681      IF(lev_histmth.GE.3) THEN
682c
683        DO nsrf=1, nbsrf
684c
685           call histdef(nid_mth, "albe_"//clnsurf(nsrf),
686     $         "Albedo surf. "//clnsurf(nsrf), "-", 
687     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]688     $         "ave(X)", zstophy,zout)
[524]689c
690           call histdef(nid_mth, "rugs_"//clnsurf(nsrf),
691     $         "Latent heat flux "//clnsurf(nsrf), "W/m2", 
692     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]693     $         "ave(X)", zstophy,zout)
[524]694c
695         CALL histdef(nid_mth, "ages_"//clnsurf(nsrf), "Snow age","day",
696     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]697     .                "ave(X)", zstophy,zout)
[524]698c
699         ENDDO !nsrf=1, nbsrf
700c
701         CALL histdef(nid_mth, "albs", "Surface albedo", "-",
702     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]703     .                "ave(X)", zstophy,zout)
[524]704         CALL histdef(nid_mth, "albslw", "Surface albedo LW", "-",
705     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
[644]706     .                "ave(X)", zstophy,zout)
[524]707c
708      ENDIF !lev_histmth.GE.3
709c
710      IF(lev_histmth.GE.4) THEN
711c
712         CALL histdef(nid_mth, "clwcon", 
713     .                "Convective Cloud Liquid water content"
714     .                , "kg/kg",
715     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]716     .                "ave(X)", zstophy,zout)
[524]717c
718         CALL histdef(nid_mth,"Ma","undilute adiab updraft","kg/m2/s",
719     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]720     .                "ave(X)", zstophy,zout)
[524]721c
722         CALL histdef(nid_mth, "dnwd", "saturated downdraft","kg/m2/s",
723     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]724     .                "ave(X)", zstophy,zout)
[524]725c
726         CALL histdef(nid_mth, "dnwd0", "unsat. downdraft", "kg/m2/s",
727     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]728     .                "ave(X)", zstophy,zout)
[524]729c
730         CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
731     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]732     .                "ave(X)", zstophy,zout)
[524]733c
[644]734         CALL histdef(nid_mth, "dqdyn", "Dynamics dQ", "(kg/kg)/s",
[524]735     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]736     .                "ave(X)", zstophy,zout)
[524]737c
738         CALL histdef(nid_mth, "dtcon", "Convection dT", "K/s",
739     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]740     .                "ave(X)", zstophy,zout)
[524]741c
742c        CALL histdef(nid_mth, "ducon", "Convection du", "m/s2",
743c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]744c    .                "ave(X)", zstophy,zout)
[524]745c
[644]746         CALL histdef(nid_mth, "dqcon", "Convection dQ", "(kg/kg)/s",
[524]747     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]748     .                "ave(X)", zstophy,zout)
[524]749c
750         CALL histdef(nid_mth, "dtlsc", "Condensation dT", "K/s",
751     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]752     .                "ave(X)", zstophy,zout)
753c
[524]754         CALL histdef(nid_mth, "dtlschr",
755     $       "Large-scale condensational heating rate", "K/s",iim,jjmp1
[644]756     $       ,nhori, klev,1,klev,nvert, 32,"ave(X)", zstophy,zout)
[524]757c
[644]758         CALL histdef(nid_mth, "dqlsc", "Condensation dQ", "(kg/kg)/s",
[524]759     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]760     .                "ave(X)", zstophy,zout)
[524]761c
762         CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
763     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]764     .                "ave(X)", zstophy,zout)
[524]765c
[644]766         CALL histdef(nid_mth, "dqvdf","Boundary-layer dQ","(kg/kg)/s",
[524]767     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]768     .                "ave(X)", zstophy,zout)
[524]769c
770         CALL histdef(nid_mth, "dteva", "Reevaporation dT", "K/s",
771     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]772     .                "ave(X)", zstophy,zout)
[524]773c
[644]774         CALL histdef(nid_mth, "dqeva","Reevaporation dQ","(kg/kg)/s",
[524]775     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]776     .                "ave(X)", zstophy,zout)
[524]777
778         CALL histdef(nid_mth, "ptconv", "POINTS CONVECTIFS"," ",
779     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]780     .                "ave(X)", zstophy,zout)
[524]781
782         CALL histdef(nid_mth, "ratqs", "RATQS"," ",
783     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]784     .                "ave(X)", zstophy,zout)
[524]785
786c
787         CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s",
788     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]789     .                "ave(X)", zstophy,zout)
[524]790
[644]791         CALL histdef(nid_mth,"dqajs","Dry adjust. dQ","(kg/kg)/s",
[524]792     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]793     .                "ave(X)", zstophy,zout)
[524]794c
795         CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
796     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]797     .                "ave(X)", zstophy,zout)
[524]798c
799         CALL histdef(nid_mth, "dtsw0", "CS SW radiation dT", "K/s",
800     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]801     .                "ave(X)", zstophy,zout)
[524]802c
803         CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
804     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]805     .                "ave(X)", zstophy,zout)
[524]806c
807         CALL histdef(nid_mth,"dtlw0","CS LW radiation dT","K/s",
808     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]809     .                "ave(X)", zstophy,zout)
[524]810c
811         CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
812     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]813     .                "ave(X)", zstophy,zout)
[524]814c
815         CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
816     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]817     .                "ave(X)", zstophy,zout)
[524]818c
819         CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
820     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]821     .                "ave(X)", zstophy,zout)
[524]822c
823         IF (ok_orodr) THEN
824         CALL histdef(nid_mth, "duoro", "Orography dU", "m/s2",
825     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]826     .                "ave(X)", zstophy,zout)
[524]827c
828         CALL histdef(nid_mth, "dvoro", "Orography dV", "m/s2",
829     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]830     .                "ave(X)", zstophy,zout)
[524]831c
832         ENDIF
833C
834         IF (ok_orolf) THEN
835         CALL histdef(nid_mth, "dulif", "Orography dU", "m/s2",
836     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]837     .                "ave(X)", zstophy,zout)
[524]838c
839         CALL histdef(nid_mth, "dvlif", "Orography dV", "m/s2",
840     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]841     .                "ave(X)", zstophy,zout)
[524]842         ENDIF
[644]843C
844         if (nqmax.GE.3) THEN
845         DO iq=1,nqmax-2
846         IF (iq.LE.99) THEN
847         WRITE(str2,'(i2.2)') iq
848         CALL histdef(nid_mth, "trac"//str2, "Tracer No."//str2, "-",
849     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
850     .                "ave(X)", zstophy,zout)
851         ELSE
852         PRINT*, "Trop de traceurs"
853         CALL abort
854         ENDIF
855         ENDDO
856         ENDIF
[524]857c
[644]858      ENDIF !lev_histmth.GE.4
[524]859c
[644]860      ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
861c
862      IF(lev_histmth.GE.1) THEN
863c
864         CALL histdef(nid_mth, "phis", "Surface geop. height", "m",
865     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
866     .                "once",  zstophy,zout)
867c
868         CALL histdef(nid_mth, "aire", "Grid area", "-",
869     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
870     .                "once",  zstophy,zout)
871c
872         CALL histdef(nid_mth, "pourc_"//clnsurf(is_ter),
873     .         "% "//clnsurf(is_ter), "%", 
874     .         iim,jjmp1,nhori, 1,1,1, -99, 32,
875     .         "once", zstophy,zout)
876c
877         CALL histdef(nid_mth, "fract_"//clnsurf(is_ter),
878     .         "Fraction "//clnsurf(is_ter), "1", 
879     .         iim,jjmp1,nhori, 1,1,1, -99, 32,
880     .         "once", zstophy,zout)
881c
882         CALL histdef(nid_mth, "pourc_"//clnsurf(is_lic),
883     .         "% "//clnsurf(is_lic), "%", 
884     .         iim,jjmp1,nhori, 1,1,1, -99, 32,
885     .         "once", zstophy,zout)
886c
887         CALL histdef(nid_mth, "fract_"//clnsurf(is_lic),
888     .         "Fraction "//clnsurf(is_lic), "1", 
889     .         iim,jjmp1,nhori, 1,1,1, -99, 32,
890     .         "once", zstophy,zout)
891c
892         CALL histdef(nid_mth, "pourc_"//clnsurf(is_sic),
893     .         "% "//clnsurf(is_sic), "%",
894     .         iim,jjmp1,nhori, 1,1,1, -99, 32,
895     .                "ave(X)", zstophy,zout)
896c
897         CALL histdef(nid_mth, "fract_"//clnsurf(is_sic),
898     .         "Fraction "//clnsurf(is_sic), "1",
899     .         iim,jjmp1,nhori, 1,1,1, -99, 32,
900     .                "ave(X)", zstophy,zout)
901c
902         CALL histdef(nid_mth, "flat", "Latent heat flux", "W/m2",
903     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
904     .                "ave(X)", zstophy,zout)
905c
906         CALL histdef(nid_mth, "bilTOA",
907     .                "Net radiation at TOA", "W/m2",
908     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
909     .                "ave(X)", zstophy,zout)
910c
911c
912c        CALL histdef(nid_mth, "tops01",
913c    .                "TOA reflected clear-sky shortwave
914c    .                radiation (method 1)", "W/m2",
915c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
916c    .                "ave(X)", zstophy,zout)
917c
918          CALL histdef(nid_mth, "dtdia","dtdia","K/s",
[524]919     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]920     .                "ave(X)", zstophy,zout)
[524]921c
[644]922          CALL histdef(nid_mth, "dtsw","dtsw","K/s",
[524]923     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]924     .                "ave(X)", zstophy,zout)
[524]925c
[644]926          CALL histdef(nid_mth, "dtlw","dtlw","K/s",
[524]927     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]928     .                "ave(X)", zstophy,zout)
[524]929c
[644]930          CALL histdef(nid_mth, "dtcon","dtcon","K/s",
[524]931     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]932     .                "ave(X)", zstophy,zout)
[524]933c
[644]934          CALL histdef(nid_mth, "dtajs","dtajs","K/s",
935     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
936     .                "ave(X)", zstophy,zout)
937c
938          CALL histdef(nid_mth, "dtlsc","dtlsc","K/s",
939     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
940     .                "ave(X)", zstophy,zout)
941c
942          CALL histdef(nid_mth, "dqdia","dqdia","(kg/kg)/s",
943     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
944     .                "ave(X)", zstophy,zout)
945c
946          CALL histdef(nid_mth, "duoli",
947     .         "duoli","m/s2",
948     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
949     .                "ave(X)", zstophy,zout)
950c
951          CALL histdef(nid_mth, "dvoli",
952     .         "dvoli","m/s2",
953     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
954     .                "ave(X)", zstophy,zout)
955c
956          CALL histdef(nid_mth, "ducon",
957     .         "ducon","m/s2",
958     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
959     .                "ave(X)", zstophy,zout)
960c
961          CALL histdef(nid_mth, "dvcon",
962     .         "dvcon","m/s2",
963     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
964     .                "ave(X)", zstophy,zout)
965c
966          CALL histdef(nid_mth, "dutot",
967     .         "dutot","m/s2",
968     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
969     .                "ave(X)", zstophy,zout)
970c
971          CALL histdef(nid_mth, "dvtot",
972     .         "dvtot","m/s2",
973     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
974     .                "ave(X)", zstophy,zout)
975c
976c         CALL histdef(nid_mth, "cldtrue",
977c    .         "TRUE Cloud fraction ","%",
978c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
979c    .         "inst(X)", zout,zout)
980c
981c         CALL histdef(nid_mth, "cldTOA",
982c    .         "Cloud amount (satellite view)","%",
983c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
984c    .         "ave(X)", zstophy,zout)
985c
986c         CALL histdef(nid_mth, "cldSFC",
987c    .         "Cloud amount (surface view)","%",
988c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
989c    .         "ave(X)", zstophy,zout)
990c
991          CALL histdef(nid_mth, "GWDtaux",
992     .         "GWD induced Ewd wind stress","Pa",
993     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
994     .         "ave(X)", zstophy,zout)
995c
996          CALL histdef(nid_mth, "GWDtauy",
997     .         "GWD induced Nwd wind stress","Pa",
998     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
999     .         "ave(X)", zstophy,zout)
1000c
1001c Champs 2D:
1002c
1003         CALL histdef(nid_mth, "slp", "Sea Level Pressure", "Pa",
1004     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1005     .                "ave(X)", zstorad,zout)
1006c
1007         CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
1008     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1009     .                "ave(X)", zstophy,zout)
1010c
1011         CALL histdef(nid_mth, "t2m", "Temperature 2m", "K",
1012     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1013     .                "ave(X)", zstorad,zout)
1014c
1015         CALL histdef(nid_mth, "t2m_min", "Temp. 2m min.",
1016     .                "K",
1017     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1018     .                t2mincels, zstophy,zout)
1019c
1020         CALL histdef(nid_mth, "t2m_max", "Temp. 2m max.",
1021     .                "K",
1022     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1023     .                t2maxcels, zstophy,zout)
1024c
1025c        CALL histdef(nid_mth, "tsoil", "Sfce soil Temperature",
1026c    .                "K",
1027c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1028c    .                "ave(X)", zstophy,zout)
1029c
1030         CALL histdef(nid_mth, "wind10m","10-m wind speed","m/s",
1031     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1032     .                "ave(X)", zstorad,zout)
1033c
1034         CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "-",
1035     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1036     .                "ave(X)", zstophy,zout)
1037
1038c
1039         CALL histdef(nid_mth, "q2m","Specific humidity 2m","kg/kg",
1040     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1041     .                "ave(X)", zstorad,zout)
1042c
1043         CALL histdef(nid_mth, "u10m", "Vent zonal 10m", "m/s",
1044     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1045     .                "ave(X)", zstorad,zout)
1046c
1047         CALL histdef(nid_mth, "v10m", "Vent meridien 10m", "m/s",
1048     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1049     .                "ave(X)", zstorad,zout)
1050c
1051         CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
1052     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1053     .                "ave(X)", zstorad,zout)
1054c
1055         CALL histdef(nid_mth, "qsurf","Surface Air humidity","kg/kg",
1056     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1057     .                "ave(X)", zstophy,zout)
1058c
1059      if (.not. ok_veget) then
1060         CALL histdef(nid_mth, "qsol", "Soil watter content", "mm",
1061     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1062     .                "ave(X)", zstophy,zout)
1063      endif
1064c
1065         CALL histdef(nid_mth, "ndayrain", 
1066     .                "Number of day with rain (liq+sol)", "-",
1067     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1068     .                "inst(X)", zstomth,zout)
1069c
1070         CALL histdef(nid_mth, "precip", "Precipitation Totale liq+sol", 
1071     .                "kg/(s*m2)",
1072     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1073     .                "ave(X)", zstophy,zout)
1074c
1075         CALL histdef(nid_mth, "plul", "Large-scale Precip.", 
1076     .   "kg/(s*m2)",
1077     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1078     .                "ave(X)", zstophy,zout)
1079c
1080         CALL histdef(nid_mth, "pluc", "Convective Precip.", 
1081     .   "kg/(s*m2)",
1082     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1083     .                "ave(X)", zstophy,zout)
1084c
1085         CALL histdef(nid_mth, "snowf", "Snow fall", "kg/(s*m2)",
1086     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1087     .                "ave(X)", zstophy,zout)
1088c
1089c        CALL histdef(nid_mth, "snow_mass", "Snow Mass", "kg/m2",
1090c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1091c    .                "ave(X)", zstophy,zout)
1092c
1093         CALL histdef(nid_mth, "evap", 
1094     .                "Evaporation plus sublimation rate", "kg/(s*m2)",
1095     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1096     .                "ave(X)", zstophy,zout)
1097c
1098         CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
1099     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1100     .                "ave(X)", zstophy,zout)
1101c
1102         CALL histdef(nid_mth, "tops0", "CS Solar rad. at TOA", "W/m2",
1103     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1104     .                "ave(X)", zstophy,zout)
1105c
1106         CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
1107     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1108     .                "ave(X)", zstophy,zout)
1109c
1110         CALL histdef(nid_mth, "topl0", "IR rad. at TOA", "W/m2",
1111     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1112     .                "ave(X)", zstophy,zout)
1113c
1114         CALL histdef(nid_mth, "SWupTOA", "SWup at TOA","W/m2",
1115     .                iim,jjmp1,nhori, 1,1,1,-99,
1116     .                32, "ave(X)", zstophy,zout)
1117c
1118         CALL histdef(nid_mth, "SWupTOAclr", 
1119     .                "SWup clear sky at TOA","W/m2",
1120     .                iim,jjmp1,nhori, 1,1,1,-99,
1121     .                32, "ave(X)", zstophy,zout)
1122c
1123         CALL histdef(nid_mth, "SWdnTOA", "SWdn at TOA","W/m2",
1124     .                iim,jjmp1,nhori, 1,1,1,-99,
1125     .                32, "ave(X)", zstophy,zout)
1126c
1127         CALL histdef(nid_mth, "SWdnTOAclr", 
1128     .                "SWdn clear sky at TOA","W/m2",
1129     .                iim,jjmp1,nhori, 1,1,1,-99,
1130     .                32, "ave(X)", zstophy,zout)
1131c
1132         CALL histdef(nid_mth, "SWup200", "SWup at 200mb","W/m2",
1133     .                iim,jjmp1,nhori, 1,1,1,-99,
1134     .                32, "ave(X)", zstophy,zout)
1135c
1136         CALL histdef(nid_mth, "SWup200clr", 
1137     .                "SWup clear sky at 200mb","W/m2",
1138     .                iim,jjmp1,nhori, 1,1,1,-99,
1139     .                32, "ave(X)", zstophy,zout)
1140c
1141         CALL histdef(nid_mth, "SWdn200", "SWdn at 200mb","W/m2",
1142     .                iim,jjmp1,nhori, 1,1,1,-99,
1143     .                32, "ave(X)", zstophy,zout)
1144c
1145         CALL histdef(nid_mth, "SWdn200clr", 
1146     .                "SWdn clear sky at 200mb","W/m2",
1147     .                iim,jjmp1,nhori, 1,1,1,-99,
1148     .                32, "ave(X)", zstophy,zout)
1149c
1150         CALL histdef(nid_mth, "LWup200", "LWup at 200mb","W/m2",
1151     .                iim,jjmp1,nhori, 1,1,1,-99,
1152     .                32, "ave(X)", zstophy,zout)
1153c
1154         CALL histdef(nid_mth, "LWup200clr", 
1155     .                "LWup clear sky at 200mb","W/m2",
1156     .                iim,jjmp1,nhori, 1,1,1,-99,
1157     .                32, "ave(X)", zstophy,zout)
1158c
1159         CALL histdef(nid_mth, "LWdn200", "LWdn at 200mb","W/m2",
1160     .                iim,jjmp1,nhori, 1,1,1,-99,
1161     .                32, "ave(X)", zstophy,zout)
1162c
1163         CALL histdef(nid_mth, "LWdn200clr", 
1164     .                "LWdn clear sky at 200mb","W/m2",
1165     .                iim,jjmp1,nhori, 1,1,1,-99,
1166     .                32, "ave(X)", zstophy,zout)
1167c
1168         CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
1169     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1170     .                "ave(X)", zstophy,zout)
1171c
1172         CALL histdef(nid_mth, "sols0", "Solar rad. at surf.", "W/m2",
1173     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1174     .                "ave(X)", zstophy,zout)
1175c
1176         CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
1177     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1178     .                "ave(X)", zstophy,zout)
1179c
1180         CALL histdef(nid_mth, "soll0", "IR rad. at surface", "W/m2",
1181     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1182     .                "ave(X)", zstophy,zout)
1183c
1184         CALL histdef(nid_mth, "SWupSFC", "SWup at surface","W/m2",
1185     .                iim,jjmp1,nhori, 1,1,1,-99,
1186     .                32, "ave(X)", zstophy,zout)
1187c
1188         CALL histdef(nid_mth, "SWupSFCclr", 
1189     .                "SWup clear sky at surface","W/m2",
1190     .                iim,jjmp1,nhori, 1,1,1,-99,
1191     .                32, "ave(X)", zstophy,zout)
1192c
1193         CALL histdef(nid_mth, "SWdnSFC", "SWdn at surface","W/m2",
1194     .                iim,jjmp1,nhori, 1,1,1,-99,
1195     .                32, "ave(X)", zstophy,zout)
1196c
1197         CALL histdef(nid_mth, "SWdnSFCclr", 
1198     .                "SWdn clear sky at surface","W/m2",
1199     .                iim,jjmp1,nhori, 1,1,1,-99,
1200     .                32, "ave(X)", zstophy,zout)
1201c
1202         CALL histdef(nid_mth,"LWupSFC","Upwd. IR rad. at surface", 
1203     .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32, 
1204     .                "ave(X)", zstophy,zout)
1205c
1206         CALL histdef(nid_mth,"LWdnSFC","Down. IR rad. at surface", 
1207     .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32, 
1208     .                "ave(X)", zstophy,zout)
1209c
1210         CALL histdef(nid_mth,"LWupSFCclr",
1211     .                "CS Upwd. IR rad. at surface", 
1212     .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32, 
1213     .                "ave(X)", zstophy,zout)
1214c
1215         CALL histdef(nid_mth,"LWdnSFCclr",
1216     .                "Down. CS IR rad. at surface", 
1217     .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32, 
1218     .                "ave(X)", zstophy,zout)
1219c
1220         CALL histdef(nid_mth, "bils", "Surf. total heat flux", "W/m2",
1221     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1222     .                "ave(X)", zstophy,zout)
1223c
1224         CALL histdef(nid_mth, "sens", "Sensible heat flux", "W/m2",
1225     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1226     .                "ave(X)", zstophy,zout)
1227c
1228         CALL histdef(nid_mth, "fder", "Heat flux derivation", "W/m2",
1229     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1230     .                "ave(X)", zstophy,zout)
1231c
1232c          CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
1233c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1234c    .                "ave(X)", zstophy,zout)
1235c
1236c        CALL histdef(nid_mth, "frtv", "Meridional wind stress", "Pa",
1237c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1238c    .                "ave(X)", zstophy,zout)
1239c
1240         CALL histdef(nid_mth, "ffonte","Thermal flux for snow melting",
1241     .                "W/m2",iim,jjmp1,nhori, 1,1,1, -99, 32,
1242     .                "ave(X)", zstophy,zout)
1243c
1244         CALL histdef(nid_mth, "fqcalving","Ice Calving",
1245     .                "kg/m2/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
1246     .                "ave(X)", zstophy,zout)
1247c
1248         DO nsrf = 1, nbsrf
1249           call histdef(nid_mth, "taux_"//clnsurf(nsrf),
1250     $         "Zonal wind stress"//clnsurf(nsrf), "Pa", 
1251     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
1252     $         "ave(X)", zstophy,zout)
1253
1254           call histdef(nid_mth, "tauy_"//clnsurf(nsrf),
1255     $         "Meridional wind stress "//clnsurf(nsrf), "Pa", 
1256     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
1257     $         "ave(X)", zstophy,zout)
1258         ENDDO
1259c
1260           call histdef(nid_mth, "pourc_"//clnsurf(is_oce),
1261     $         "% "//clnsurf(is_oce), "%", 
1262     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
1263     $         "ave(X)", zstophy,zout)
1264c
1265           call histdef(nid_mth, "fract_"//clnsurf(is_oce),
1266     $         "Fraction "//clnsurf(is_oce), "1", 
1267     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
1268     $         "ave(X)", zstophy,zout)
1269C
1270         DO nsrf = 1, nbsrf
1271C
1272           call histdef(nid_mth, "tsol_"//clnsurf(nsrf),
1273     $         "Temperature "//clnsurf(nsrf), "K", 
1274     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
1275     $         "ave(X)", zstophy,zout)
1276C
1277           call histdef(nid_mth, "sens_"//clnsurf(nsrf),
1278     $         "Sensible heat flux "//clnsurf(nsrf), "W/m2", 
1279     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
1280     $         "ave(X)", zstophy,zout)
1281c
1282           call histdef(nid_mth, "lat_"//clnsurf(nsrf),
1283     $         "Latent heat flux "//clnsurf(nsrf), "W/m2", 
1284     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
1285     $         "ave(X)", zstophy,zout)
1286C
1287           call histdef(nid_mth, "flw_"//clnsurf(nsrf),
1288     $         "LW "//clnsurf(nsrf), "W/m2",
1289     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
1290     $         "ave(X)", zstophy,zout)
1291c
1292           call histdef(nid_mth, "fsw_"//clnsurf(nsrf),
1293     $         "SW "//clnsurf(nsrf), "W/m2",
1294     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
1295     $         "ave(X)", zstophy,zout)
1296C
1297           call histdef(nid_mth, "wbils_"//clnsurf(nsrf),
1298     $         "Bilan sol "//clnsurf(nsrf), "W/m2",
1299     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
1300     $         "ave(X)", zstophy,zout)
1301C
1302cIM 250105 BEG
1303          IF(1.EQ.0) THEN
1304C
1305           call histdef(nid_mth, "evapour_"//clnsurf(nsrf),
1306     $         "Evap*pourc "//clnsurf(nsrf), "kg/(s*m2)",
1307     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
1308     $         "ave(X)", zstophy,zout)
1309C
1310           call histdef(nid_mth, "prepour_"//clnsurf(nsrf),
1311     $         "Precip*pourc "//clnsurf(nsrf), "kg/(s*m2)",
1312     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
1313     $         "ave(X)", zstophy,zout)
1314C
1315          ENDIF !(1.EQ.0) THEN
1316cIM 250105 END
1317C
1318         END DO
1319c
1320         CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
1321     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1322     .                "ave(X)", zstophy,zout)
1323c
1324         CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
1325     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1326     .                "ave(X)", zstophy,zout)
1327c
1328         CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-",
1329     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1330     .                "ave(X)", zstophy,zout)
1331c
1332         CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-",
1333     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1334     .                "ave(X)", zstophy,zout)
1335c
1336         CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-",
1337     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1338     .                "ave(X)", zstophy,zout)
1339c
1340         CALL histdef(nid_mth, "cldt", "Total cloudiness", "%",
1341     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1342     .                "ave(X)", zstophy,zout)
1343c
1344         CALL histdef(nid_mth,"cldq","Cloud liquid water path","kg/m2",
1345     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1346     .                "ave(X)", zstophy,zout)
1347c
1348         CALL histdef(nid_mth,"lwp","Cloud water path","kg/m2",
1349     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1350     .                "ave(X)", zstophy,zout)
1351c
1352         CALL histdef(nid_mth,"iwp","Cloud ice water path","kg/m2",
1353     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1354     .                "ave(X)", zstophy,zout)
1355c
1356         CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
1357     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1358     .                "ave(X)", zstophy,zout)
1359c
1360         CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
1361     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1362     .                "ave(X)", zstophy,zout)
1363c
1364         CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-",
1365     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1366     .                "ave(X)", zstophy,zout)
1367c
1368         CALL histdef(nid_mth, "vq", "Merid humidity transport", "-",
1369     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1370     .                "ave(X)", zstophy,zout)
1371cKE43
1372      IF(iflag_con.GE.3) THEN ! sb
1373c
1374         CALL histdef(nid_mth, "cape", "Conv avlbl pot ener", "J/kg",
1375     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1376     .                "ave(X)", zstophy,zout)
1377c
1378         CALL histdef(nid_mth, "pbase", "Cld base pressure", "mb",
1379     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1380     .                "ave(X)", zstophy,zout)
1381c
1382         CALL histdef(nid_mth, "ptop", "Cld top pressure", "mb",
1383     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1384     .                "ave(X)", zstophy,zout)
1385c
1386         CALL histdef(nid_mth, "fbase", "Cld base mass flux", "kg/m2/s",
1387     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1388     .                "ave(X)", zstophy,zout)
1389c
1390         CALL histdef(nid_mth, "prw", "Precipitable water", "kg/m2",
1391     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1392     .                "ave(X)", zstophy,zout)
1393c
[653]1394         CALL histdef(nid_mth, "upwd", "saturated updraft", "kg/m2/s",
1395     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1396     .                "ave(X)", zstophy,zout)
1397c
[644]1398      ENDIF !iflag_con .GE. 3
1399c34EK
1400cIM cf. AM 081204 BEG
1401c34EK
1402c HBTM2
1403         CALL histdef(nid_mth, "s_pblh", "Boundary Layer Height", "m",
1404     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1405     .                "ave(X)", zstophy,zout)
1406         CALL histdef(nid_mth, "s_pblt", "T at Boundary Layer Height",
1407     .                "K",
1408     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1409     .                "ave(X)", zstophy,zout)
1410         CALL histdef(nid_mth, "s_lcl", "Condensation level", "m",
1411     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1412     .                "ave(X)", zstophy,zout)
1413         CALL histdef(nid_mth,"s_capCL","Conv avlbl pot ener for ABL",
1414     .                "J/m2",
1415     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1416     .                "ave(X)", zstophy,zout)
1417         CALL histdef(nid_mth, "s_oliqCL", "Liq Water in BL", "kg/m2",
1418     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1419     .                "ave(X)", zstophy,zout)
1420         CALL histdef(nid_mth,"s_cteiCL","Instability criteria (ABL)",
1421     .                "K",
1422     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1423     .                "ave(X)", zstophy,zout)
1424         CALL histdef(nid_mth, "s_therm", "Exces du thermique", "K",
1425     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1426     .                "ave(X)", zstophy,zout)
1427         CALL histdef(nid_mth, "s_trmb1", "deep_cape(HBTM2)", "J/m2",
1428     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1429     .                "ave(X)", zstophy,zout)
1430         CALL histdef(nid_mth, "s_trmb2", "inhibition (HBTM2)", "J/m2",
1431     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1432     .                "ave(X)", zstophy,zout)
1433         CALL histdef(nid_mth, "s_trmb3", "Point Omega (HBTM2)", "m",
1434     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1435     .                "ave(X)", zstophy,zout)
1436c
1437cIM cf. AM 081204 END
1438c
1439c Champs interpolles sur des niveaux de pression
1440         DO k=1, nlevSTD
1441c
1442c        bb=clevSTD(k)
1443c
1444c        IF(k.GE.2) THEN
1445         IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
1446         IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
1447c         aa=clevSTD(k)
1448c         bb=aa(1:lnblnk1(aa))
1449c        ENDIF
1450c
1451       IF(bb2.EQ."850".OR.bb2.EQ."700"
1452     $   .OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
1453c
1454          CALL histdef(nid_mth, "u"//bb2,
1455     .                 "Zonal wind "//bb2//"mb","m/s",
1456     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1457     .                "inst(X)", zout,zout)
1458c
1459          CALL histdef(nid_mth, "v"//bb2,
1460     .                 "Meridional wind "//bb2//"mb","m/s",
1461     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1462     .                "inst(X)", zout,zout)
1463c
1464          CALL histdef(nid_mth, "w"//bb2,
1465     .                 "Vertical wind "//bb2//"mb","m/s",
1466     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1467     .                "inst(X)", zout,zout)
1468c
1469          CALL histdef(nid_mth, "phi"//bb2,
1470     .                 "Geopotential "//bb2//"mb","m",
1471     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
1472     .                "inst(X)", zout,zout)
1473c
1474         ENDIF !(bb2.EQ."850".OR.bb2.EQ."700"
1475c
1476       ENDDO !nlevSTD
1477c
1478cIM diagnostiques flux ocean-atm ou ocean-glace de mer
1479cIM pour utilisation dans un modele de "slab" ocean
1480c
1481       CALL histdef(nid_mth, "fluxo",
1482     .              "Flux turbulents ocean-atmosphere", "W/m2",
1483     .              iim,jjmp1,nhori, 1,1,1, -99, 32,
1484     .              "ave(X)", zstophy,zout)
1485c
1486       CALL histdef(nid_mth, "fluxg",
1487     .              "Flux turbulents ocean-glace de mer","W/m2",
1488     .              iim,jjmp1,nhori, 1,1,1, -99, 32,
1489     .              "ave(X)", zstophy,zout)
1490c
1491       IF (OCEAN.EQ.'force ') THEN
1492c
1493        CALL histdef(nid_mth, "lmt_bils",
1494     .       "Bilan au sol atmosphere forcee", "W/m2",
1495     .       iim,jjmp1,nhori, 1,1,1, -99, 32,
1496     .       "ave(X)", zstophy,zout)
1497c
1498       ELSE IF (OCEAN.EQ.'slab  ') THEN
1499c
1500        CALL histdef(nid_mth, "slab_bils",
1501     .              "Bilan au sol Slab", "W/m2",
1502     .              iim,jjmp1,nhori, 1,1,1, -99, 32,
1503     .              "ave(X)", zstophy,zout)
1504c
1505        CALL histdef(nid_mth, "tslab", "Slab SST ", "K",
1506     .              iim,jjmp1,nhori, 1,1,1, -99, 32,
1507     .              "ave(X)", zstophy,zout)
1508c
1509        CALL histdef(nid_mth, "seaice", "Slab seaice", "kg/m2",
1510     .              iim,jjmp1,nhori, 1,1,1, -99, 32,
1511     .              "ave(X)", zstophy,zout)
1512c
1513        CALL histdef(nid_mth, "siceh", "Slab seaice height", "m",
1514     .              iim,jjmp1,nhori, 1,1,1, -99, 32,
1515     .              "ave(X)", zstophy,zout)
1516c
1517       ENDIF !OCEAN.EQ.slab/force
1518c
1519      ENDIF !lev_histmth.GE.1
1520c
1521      IF(lev_histmth.GE.2) THEN
1522c
1523c Champs 3D:
1524c
1525c
1526         CALL histdef(nid_mth,"lwcon","Cloud liquid water content",
1527     .                "kg/kg",
1528     .                iim,jjmp1,nhori, klev,1,klev, nvert, 32, 
1529     .                "ave(X)", zstophy,zout)
1530c
1531         CALL histdef(nid_mth,"iwcon","Cloud ice water content","kg/kg",
1532     .                iim,jjmp1,nhori, klev,1,klev, nvert, 32, 
1533     .                "ave(X)", zstophy,zout)
1534c
1535         CALL histdef(nid_mth, "temp", "Air temperature", "K",
1536     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1537     .                "ave(X)", zstophy,zout)
1538c
1539         CALL histdef(nid_mth, "ovap", "Specific humidity", "kg/kg",
1540     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1541     .                "ave(X)", zstophy,zout)
1542c
1543c        CALL histdef(nid_mth,"wvap","Water vapor mixing ratio","kg/kg",
1544c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1545c    .                "ave(X)", zstophy,zout)
1546c
1547         CALL histdef(nid_mth, "geop", "Geopotential height", "m2/s2",
1548     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1549     .                "ave(X)", zstophy,zout)
1550c
1551         CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
1552     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1553     .                "ave(X)", zstophy,zout)
1554c
1555         CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
1556     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1557     .                "ave(X)", zstophy,zout)
1558c
1559         CALL histdef(nid_mth, "vitw", "Vertical wind", "m/s",
1560     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1561     .                "ave(X)", zstophy,zout)
1562c
1563         CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
1564     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1565     .                "ave(X)", zstophy,zout)
1566c
1567        CALL histdef(nid_mth, "rneb", "Cloud fraction", "%",
1568     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1569     .                "ave(X)", zstophy,zout)
1570c
1571         CALL histdef(nid_mth, "rnebcon", "Convective Cloud Fraction"
1572     .                , "-",
1573     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1574     .                "ave(X)", zstophy,zout)
1575c
1576         CALL histdef(nid_mth, "rhum", "Relative humidity", "-",
1577     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1578     .                "ave(X)", zstophy,zout)
1579c
1580         CALL histdef(nid_mth, "ozone", "Ozone concentration", "ppmv",
1581     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1582     .                "ave(X)", zstophy,zout)
1583c
1584         CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
1585     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1586     .                "ave(X)", zstophy,zout)
1587c
1588         CALL histdef(nid_mth, "dqphy","Physics dQ","(kg/kg)/s",
1589     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1590     .                "ave(X)", zstophy,zout)
1591c
1592      ENDIF !lev_histmth.GE.2
1593c
1594      IF(lev_histmth.GE.3) THEN
1595c
1596        DO nsrf=1, nbsrf
1597c
1598           call histdef(nid_mth, "albe_"//clnsurf(nsrf),
1599     $         "Albedo surf. "//clnsurf(nsrf), "-", 
1600     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
1601     $         "ave(X)", zstophy,zout)
1602c
1603           call histdef(nid_mth, "rugs_"//clnsurf(nsrf),
1604     $         "Rugosite "//clnsurf(nsrf), "m", 
1605     $         iim,jjmp1,nhori, 1,1,1, -99, 32,
1606     $         "ave(X)", zstophy,zout)
1607c
1608         CALL histdef(nid_mth, "ages_"//clnsurf(nsrf), "Snow age","day",
1609     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1610     .                "ave(X)", zstophy,zout)
1611c
1612         ENDDO !nsrf=1, nbsrf
1613c
1614         CALL histdef(nid_mth, "albs", "Surface albedo", "-",
1615     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1616     .                "ave(X)", zstophy,zout)
1617         CALL histdef(nid_mth, "albslw", "Surface albedo LW", "-",
1618     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
1619     .                "ave(X)", zstophy,zout)
1620c
1621      ENDIF !lev_histmth.GE.3
1622c
1623      IF(lev_histmth.GE.4) THEN
1624c
1625         CALL histdef(nid_mth, "clwcon", 
1626     .                "Convective Cloud Liquid water content"
1627     .                , "kg/kg",
1628     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1629     .                "ave(X)", zstophy,zout)
1630c
1631         CALL histdef(nid_mth,"Ma","undilute adiab updraft","kg/m2/s",
1632     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1633     .                "ave(X)", zstophy,zout)
1634c
1635         CALL histdef(nid_mth, "dnwd", "saturated downdraft","kg/m2/s",
1636     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1637     .                "ave(X)", zstophy,zout)
1638c
1639         CALL histdef(nid_mth, "dnwd0", "unsat. downdraft", "kg/m2/s",
1640     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1641     .                "ave(X)", zstophy,zout)
1642c
1643         CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
1644     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1645     .                "ave(X)", zstophy,zout)
1646c
1647         CALL histdef(nid_mth, "dqdyn","Dynamics dQ","(kg/kg)/s",
1648     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1649     .                "ave(X)", zstophy,zout)
1650c
1651         CALL histdef(nid_mth, "dqcon","Convection dQ","(kg/kg)/s",
1652     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1653     .                "ave(X)", zstophy,zout)
1654c
1655         CALL histdef(nid_mth, "dtlschr",
1656     $       "Large-scale condensational heating rate", "K/s",iim,jjmp1
1657     $       ,nhori, klev,1,klev,nvert, 32,"ave(X)", zstophy,zout)
1658c
1659         CALL histdef(nid_mth, "dqlsc","Condensation dQ","(kg/kg)/s",
1660     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1661     .                "ave(X)", zstophy,zout)
1662c
1663         CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
1664     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1665     .                "ave(X)", zstophy,zout)
1666c
1667         CALL histdef(nid_mth, "dqvdf","Boundary-layer dQ","(kg/kg)/s",
1668     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1669     .                "ave(X)", zstophy,zout)
1670c
1671         CALL histdef(nid_mth, "dteva", "Reevaporation dT", "K/s",
1672     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1673     .                "ave(X)", zstophy,zout)
1674c
1675         CALL histdef(nid_mth,"dqeva","Reevaporation dQ","(kg/kg)/s",
1676     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1677     .                "ave(X)", zstophy,zout)
1678
1679         CALL histdef(nid_mth, "ptconv", "POINTS CONVECTIFS"," ",
1680     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1681     .                "ave(X)", zstophy,zout)
1682
1683         CALL histdef(nid_mth, "ratqs", "RATQS"," ",
1684     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1685     .                "ave(X)", zstophy,zout)
1686
1687c
1688         CALL histdef(nid_mth, "dqajs","Dry adjust. dQ","(kg/kg)/s",
1689     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1690     .                "ave(X)", zstophy,zout)
1691c
1692         CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
1693     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1694     .                "ave(X)", zstophy,zout)
1695c
1696         CALL histdef(nid_mth, "dtsw0", "CS SW radiation dT", "K/s",
1697     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1698     .                "ave(X)", zstophy,zout)
1699c
1700         CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
1701     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1702     .                "ave(X)", zstophy,zout)
1703c
1704         CALL histdef(nid_mth,"dtlw0","CS LW radiation dT","K/s",
1705     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1706     .                "ave(X)", zstophy,zout)
1707c
1708         CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
1709     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1710     .                "ave(X)", zstophy,zout)
1711c
1712         CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
1713     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1714     .                "ave(X)", zstophy,zout)
1715c
1716         CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
1717     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1718     .                "ave(X)", zstophy,zout)
1719c
1720         IF (ok_orodr) THEN
1721         CALL histdef(nid_mth, "duoro", "Orography dU", "m/s2",
1722     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1723     .                "ave(X)", zstophy,zout)
1724c
1725         CALL histdef(nid_mth, "dvoro", "Orography dV", "m/s2",
1726     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1727     .                "ave(X)", zstophy,zout)
1728c
1729         ENDIF
1730C
1731         IF (ok_orolf) THEN
1732         CALL histdef(nid_mth, "dulif", "Orography dU", "m/s2",
1733     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1734     .                "ave(X)", zstophy,zout)
1735c
1736         CALL histdef(nid_mth, "dvlif", "Orography dV", "m/s2",
1737     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
1738     .                "ave(X)", zstophy,zout)
1739         ENDIF
1740C
[524]1741         if (nqmax.GE.3) THEN
[653]1742         DO iq=3,nqmax
1743         iiq=niadv(iq)
1744         CALL histdef(nid_mth, tnom(iq), ttext(iiq), "-",
[524]1745     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]1746     .                "ave(X)", zstophy,zout)
[524]1747         ENDDO
1748         ENDIF
[653]1749
[524]1750      ENDIF !lev_histmth.GE.4
1751c
[644]1752      ENDIF !type_run.EQ."CLIM".OR.type_run.EQ."ENSP")
1753c
[524]1754         CALL histend(nid_mth)
1755c
1756         ndex2d = 0
1757         ndex3d = 0
1758c
1759      ENDIF ! fin de test sur ok_mensuel
Note: See TracBrowser for help on using the repository browser.