Context Navigation

source: LMDZ4/trunk/libf/phylmd/ini_histmth.h @ 748

Last change on this file since 748 was 684, checked in by lmdzadmin, 19 years ago
Ajout nouveaux flags (ip_ebil_phy, ok_slab_sicOBS) ds physiq.def + Lecture freq. ecriture en nombre de jours IM
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 69.2 KB

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

Note: See TracBrowser for help on using the repository browser.

Download in other formats: