Context Navigation

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

Last change on this file since 665 was 653, checked in by Laurent Fairhead, 20 years ago
On avait perdu des modifs de Anne et Marie-Angele LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 68.2 KB

Rev	Line
[644]	1	c
	2	c $Header$
	3	c
[524]	4	IF (ok_mensuel) THEN
	5	c
[644]	6	zstophy = dtime
	7	zstomth = dtime * FLOAT(ecrit_mth) !valeurs mensuelles "instantannees"
	8	zstorad = dtime * FLOAT(radpas) !variables stockees toutes les 2h
	9	zstohf = dtime *FLOAT(ecrit_hf) !variables stockees toutes les 6h
	10	zout = dtime * FLOAT(ecrit_mth) !frequence des sorties : 1 val. par mois
[524]	11	c
[644]	12	cIM PRINT*,' La zstophy,zstorad,zstohf,zstomth,zout=',
	13	cIM . zstophy,zstorad,zstohf,zstomth,zout
[524]	14	c
	15	idayref = day_ref
	16	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
	17	c
	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]	23	c
	24	cIM cf. AM 081204 BEG region
	25	c imin_histmth=1
	26	c imax_histmth=iim
	27	c jmin_histmth=1
	28	c jmax_histmth=jjmp1
	29	c do i=1,iim-1
	30	c if(zx_lon(i,1).lt.lonmin_histmth) imin_histmth=i
	31	c if(zx_lon(i,1).le.lonmax_histmth) imax_histmth=i+1
	32	c enddo
	33	c do j=1,jjmp1-1
	34	c if(zx_lat(1,j).ge.latmin_histmth) jmax_histmth=j
	35	c if(zx_lat(1,j).gt.latmax_histmth) jmin_histmth=j
	36	c enddo
	37	c
	38	c print*,'On stoke le fichier instantanne sur, ',
	39	c s imin_histmth,imax_histmth,jmin_histmth,jmax_histmth
	40	c print*,'On stoke le fichier instantanne sur, ',
	41	c s zx_lon(imin_histmth,1),zx_lon(imax_histmth,1),
	42	c s zx_lat(1,jmin_histmth),zx_lat(1,jmax_histmth)
	43	cIM cf. AM 081204 END region
	44	c
[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]	52	c
	53	cIM cf. AM 081204 BEG region
	54	c CALL histbeg("histmth.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
	55	c . imin_histmth,imax_histmth-imin_histmth+1
	56	c . ,jmin_histmth,jmax_histmth-jmin_histmth+1
	57	c . , itau_phy, zjulian, dtime,
	58	c . nhori, nid_mth)
	59	cIM 081204 END
	60	c
[524]	61	CALL histvert(nid_mth, "presnivs", "Vertical levels", "mb",
[644]	62	. klev, presnivs/100., nvert)
[524]	63	c
[644]	64	IF(type_run.EQ."CLIM".OR.type_run.EQ."ENSP") THEN
[524]	65	c
	66	IF(lev_histmth.GE.1) THEN
[644]	67	c
	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]	72	c
	73	CALL histdef(nid_mth, "aire", "Grid area", "-",
	74	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	75	. "once", zstophy,zout)
[524]	76	c
[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]	81	c
[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)
	86	c
[524]	87	c Champs 2D:
	88	c
	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]	92	c
	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]	96	c
	97	CALL histdef(nid_mth, "t2m", "Temperature 2m", "K",
	98	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	99	. "ave(X)", zstophy,zout)
	100	c
[524]	101	c 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]	106	c
	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]	111	c
	112	c CALL histdef(nid_mth, "tsoil", "Sfce soil Temperature",
	113	c . "K",
	114	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	115	c . "ave(X)", zstophy,zout)
[524]	116	c
	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]	120	c
	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]	124	c
	125	c 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]	129	c
	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]	133	c
	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]	137	c
	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]	141	c
	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]	145	c
	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
	151	c
	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]	156	c
	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]	161	c
	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]	166	c
	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]	171	c
	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]	175	c
[644]	176	c CALL histdef(nid_mth, "snow_mass", "Snow Mass", "kg/m2",
	177	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	178	c . "ave(X)", zstophy,zout)
[524]	179	c
	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]	183	c
	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]	187	c
[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]	191	c
	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]	195	c
	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]	199	c
	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]	203	c
	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]	208	c
	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]	212	c
	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]	217	c
[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]	221	c
	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]	226	c
[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]	230	c
	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]	235	c
[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]	239	c
	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]	244	c
[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]	248	c
	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]	253	c
	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]	257	c
	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]	261	c
	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]	265	c
	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]	269	c
	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]	273	c
	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]	278	c
	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]	282	c
	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]	287	c
	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]	291	c
	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]	295	c
	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]	300	c
	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]	305	c
	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]	309	c
	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]	313	c
	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]	317	c
	318	c CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
	319	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	320	c . "ave(X)", zstophy,zout)
[524]	321	c
	322	c CALL histdef(nid_mth, "frtv", "Meridional wind stress", "Pa",
	323	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	324	c . "ave(X)", zstophy,zout)
[524]	325	c
	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]	329	c
	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)
	333	c
[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
	345	c
	346	DO nsrf = 1, nbsrf
	347	C
[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)
	353	c
	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]	359	C
	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]	364	C
	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]	369	c
	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]	374	C
	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]	379	c
	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]	384	C
	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]	389	C
[644]	390	cIM 250105 BEG
	391	IF(1.EQ.0) THEN
	392	C
	393	call histdef(nid_mth, "evapour_"//clnsurf(nsrf),
	394	$ "Evappourc "//clnsurf(nsrf), "kg/(sm2)",
	395	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	396	$ "ave(X)", zstophy,zout)
	397	C
	398	call histdef(nid_mth, "prepour_"//clnsurf(nsrf),
	399	$ "Precippourc "//clnsurf(nsrf), "kg/(sm2)",
	400	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	401	$ "ave(X)", zstophy,zout)
	402	C
	403	ENDIF !(1.EQ.0) THEN
	404	cIM 250105 END
	405	C
[524]	406	END DO
	407	c
	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]	411	c
	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]	415	c
	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]	419	c
	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]	423	c
	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]	427	c
	428	CALL histdef(nid_mth, "cldt", "Total cloudiness", "-",
	429	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	430	. "ave(X)", zstophy,zout)
[524]	431	c
	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)
	435	c
[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]	439	c
	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]	443	c
	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]	447	c
	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]	451	c
	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]	455	c
	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]	459	cKE43
	460	IF(iflag_con.GE.3) THEN ! sb
	461	c
	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]	465	c
[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]	469	c
[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]	473	c
	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]	477	c
	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]	481	c
[644]	482	ENDIF !iflag_con .GE. 3
[524]	483	c34EK
[644]	484	cIM cf. AM 081204 BEG
	485	c34EK
	486	c 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]	520	c
[644]	521	cIM cf. AM 081204 END
	522	c
[524]	523	c Champs interpolles sur des niveaux de pression
[644]	524	DO k=1, nlevSTD
	525	c
	526	c bb=clevSTD(k)
	527	c
	528	c 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)
	531	c aa=clevSTD(k)
	532	c bb=aa(1:lnblnk1(aa))
	533	c ENDIF
	534	c
	535	IF(bb2.EQ."850".OR.bb2.EQ."700".OR.
	536	$ bb2.EQ."500".OR.bb2.EQ."200") THEN
	537	c
	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]	542	c
[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]	547	c
[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]	552	c
[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]	557	c
[644]	558	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.
	559	ENDDO
	560	c
	561	cIM diagnostiques flux ocean-atm ou ocean-glace de mer
	562	cIM pour utilisation dans un modele de "slab" ocean
	563	c
	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)
	568	c
	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)
	573	c
	574	IF (OCEAN.EQ.'force ') THEN
	575	c
	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)
	580	c
	581	ELSE IF (OCEAN.EQ.'slab ') THEN
	582	c
	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)
	587	c
	588	CALL histdef(nid_mth, "tslab", "Slab SST ", "K",
	589	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	590	. "ave(X)", zstophy,zout)
	591	c
	592	CALL histdef(nid_mth, "seaice", "Slab seaice", "kg/m2",
	593	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	594	. "ave(X)", zstophy,zout)
	595	c
	596	CALL histdef(nid_mth, "siceh", "Slab seaice height", "m",
	597	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	598	. "ave(X)", zstophy,zout)
	599	c
	600	ENDIF
	601	c
[524]	602	ENDIF !lev_histmth.GE.1
	603	c
	604	IF(lev_histmth.GE.2) THEN
	605	c
	606	c Champs 3D:
	607	c
	608	c
[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]	613	c
	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]	617	c
	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]	621	c
	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]	625	c
	626	c CALL histdef(nid_mth,"wvap","Water vapor mixing ratio","kg/kg",
	627	c . iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]	628	c . "ave(X)", zstophy,zout)
[524]	629	c
[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]	633	c
	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]	637	c
	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]	641	c
[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]	645	c
	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]	649	c
	650	CALL histdef(nid_mth, "rneb", "Cloud fraction", "-",
	651	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]	652	. "ave(X)", zstophy,zout)
[524]	653	c
	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]	658	c
	659	CALL histdef(nid_mth, "rhum", "Relative humidity", "-",
	660	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]	661	. "ave(X)", zstophy,zout)
[524]	662	c
[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]	666	c
[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)
	670	c
[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]	674	c
[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]	678	c
	679	ENDIF !lev_histmth.GE.2
	680	c
	681	IF(lev_histmth.GE.3) THEN
	682	c
	683	DO nsrf=1, nbsrf
	684	c
	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]	689	c
	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]	694	c
	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]	698	c
	699	ENDDO !nsrf=1, nbsrf
	700	c
	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]	707	c
	708	ENDIF !lev_histmth.GE.3
	709	c
	710	IF(lev_histmth.GE.4) THEN
	711	c
	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]	717	c
	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]	721	c
	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]	725	c
	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]	729	c
	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]	733	c
[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]	737	c
	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]	741	c
	742	c CALL histdef(nid_mth, "ducon", "Convection du", "m/s2",
	743	c . iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]	744	c . "ave(X)", zstophy,zout)
[524]	745	c
[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]	749	c
	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)
	753	c
[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]	757	c
[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]	761	c
	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]	765	c
[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]	769	c
	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]	773	c
[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
	786	c
	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]	794	c
	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]	798	c
	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]	802	c
	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]	806	c
	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]	810	c
	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]	814	c
	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]	818	c
	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]	822	c
	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]	827	c
	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]	831	c
	832	ENDIF
	833	C
	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]	838	c
	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]	843	C
	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]	857	c
[644]	858	ENDIF !lev_histmth.GE.4
[524]	859	c
[644]	860	ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
	861	c
	862	IF(lev_histmth.GE.1) THEN
	863	c
	864	CALL histdef(nid_mth, "phis", "Surface geop. height", "m",
	865	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	866	. "once", zstophy,zout)
	867	c
	868	CALL histdef(nid_mth, "aire", "Grid area", "-",
	869	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	870	. "once", zstophy,zout)
	871	c
	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)
	876	c
	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)
	881	c
	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)
	886	c
	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)
	891	c
	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)
	896	c
	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)
	901	c
	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)
	905	c
	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)
	910	c
	911	c
	912	c CALL histdef(nid_mth, "tops01",
	913	c . "TOA reflected clear-sky shortwave
	914	c . radiation (method 1)", "W/m2",
	915	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	916	c . "ave(X)", zstophy,zout)
	917	c
	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]	921	c
[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]	925	c
[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]	929	c
[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]	933	c
[644]	934	CALL histdef(nid_mth, "dtajs","dtajs","K/s",
	935	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	936	. "ave(X)", zstophy,zout)
	937	c
	938	CALL histdef(nid_mth, "dtlsc","dtlsc","K/s",
	939	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	940	. "ave(X)", zstophy,zout)
	941	c
	942	CALL histdef(nid_mth, "dqdia","dqdia","(kg/kg)/s",
	943	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	944	. "ave(X)", zstophy,zout)
	945	c
	946	CALL histdef(nid_mth, "duoli",
	947	. "duoli","m/s2",
	948	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	949	. "ave(X)", zstophy,zout)
	950	c
	951	CALL histdef(nid_mth, "dvoli",
	952	. "dvoli","m/s2",
	953	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	954	. "ave(X)", zstophy,zout)
	955	c
	956	CALL histdef(nid_mth, "ducon",
	957	. "ducon","m/s2",
	958	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	959	. "ave(X)", zstophy,zout)
	960	c
	961	CALL histdef(nid_mth, "dvcon",
	962	. "dvcon","m/s2",
	963	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	964	. "ave(X)", zstophy,zout)
	965	c
	966	CALL histdef(nid_mth, "dutot",
	967	. "dutot","m/s2",
	968	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	969	. "ave(X)", zstophy,zout)
	970	c
	971	CALL histdef(nid_mth, "dvtot",
	972	. "dvtot","m/s2",
	973	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	974	. "ave(X)", zstophy,zout)
	975	c
	976	c CALL histdef(nid_mth, "cldtrue",
	977	c . "TRUE Cloud fraction ","%",
	978	c . iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	979	c . "inst(X)", zout,zout)
	980	c
	981	c CALL histdef(nid_mth, "cldTOA",
	982	c . "Cloud amount (satellite view)","%",
	983	c . iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	984	c . "ave(X)", zstophy,zout)
	985	c
	986	c CALL histdef(nid_mth, "cldSFC",
	987	c . "Cloud amount (surface view)","%",
	988	c . iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	989	c . "ave(X)", zstophy,zout)
	990	c
	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)
	995	c
	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)
	1000	c
	1001	c Champs 2D:
	1002	c
	1003	CALL histdef(nid_mth, "slp", "Sea Level Pressure", "Pa",
	1004	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	1005	. "ave(X)", zstorad,zout)
	1006	c
	1007	CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
	1008	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	1009	. "ave(X)", zstophy,zout)
	1010	c
	1011	CALL histdef(nid_mth, "t2m", "Temperature 2m", "K",
	1012	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	1013	. "ave(X)", zstorad,zout)
	1014	c
	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)
	1019	c
	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)
	1024	c
	1025	c CALL histdef(nid_mth, "tsoil", "Sfce soil Temperature",
	1026	c . "K",
	1027	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	1028	c . "ave(X)", zstophy,zout)
	1029	c
	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)
	1033	c
	1034	CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "-",
	1035	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	1036	. "ave(X)", zstophy,zout)
	1037
	1038	c
	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)
	1042	c
	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)
	1046	c
	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)
	1050	c
	1051	CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
	1052	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	1053	. "ave(X)", zstorad,zout)
	1054	c
	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)
	1058	c
	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
	1064	c
	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)
	1069	c
	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)
	1074	c
	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)
	1079	c
	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)
	1084	c
	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)
	1088	c
	1089	c CALL histdef(nid_mth, "snow_mass", "Snow Mass", "kg/m2",
	1090	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	1091	c . "ave(X)", zstophy,zout)
	1092	c
	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)
	1097	c
	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)
	1101	c
	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)
	1105	c
	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)
	1109	c
	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)
	1113	c
	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)
	1117	c
	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)
	1122	c
	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)
	1126	c
	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)
	1131	c
	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)
	1135	c
	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)
	1140	c
	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)
	1144	c
	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)
	1149	c
	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)
	1153	c
	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)
	1158	c
	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)
	1162	c
	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)
	1167	c
	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)
	1171	c
	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)
	1175	c
	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)
	1179	c
	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)
	1183	c
	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)
	1187	c
	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)
	1192	c
	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)
	1196	c
	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)
	1201	c
	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)
	1205	c
	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)
	1209	c
	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)
	1214	c
	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)
	1219	c
	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)
	1223	c
	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)
	1227	c
	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)
	1231	c
	1232	c CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
	1233	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	1234	c . "ave(X)", zstophy,zout)
	1235	c
	1236	c CALL histdef(nid_mth, "frtv", "Meridional wind stress", "Pa",
	1237	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	1238	c . "ave(X)", zstophy,zout)
	1239	c
	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)
	1243	c
	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)
	1247	c
	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
	1259	c
	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)
	1264	c
	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)
	1269	C
	1270	DO nsrf = 1, nbsrf
	1271	C
	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)
	1276	C
	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)
	1281	c
	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)
	1286	C
	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)
	1291	c
	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)
	1296	C
	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)
	1301	C
	1302	cIM 250105 BEG
	1303	IF(1.EQ.0) THEN
	1304	C
	1305	call histdef(nid_mth, "evapour_"//clnsurf(nsrf),
	1306	$ "Evappourc "//clnsurf(nsrf), "kg/(sm2)",
	1307	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	1308	$ "ave(X)", zstophy,zout)
	1309	C
	1310	call histdef(nid_mth, "prepour_"//clnsurf(nsrf),
	1311	$ "Precippourc "//clnsurf(nsrf), "kg/(sm2)",
	1312	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	1313	$ "ave(X)", zstophy,zout)
	1314	C
	1315	ENDIF !(1.EQ.0) THEN
	1316	cIM 250105 END
	1317	C
	1318	END DO
	1319	c
	1320	CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
	1321	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	1322	. "ave(X)", zstophy,zout)
	1323	c
	1324	CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
	1325	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	1326	. "ave(X)", zstophy,zout)
	1327	c
	1328	CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-",
	1329	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	1330	. "ave(X)", zstophy,zout)
	1331	c
	1332	CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-",
	1333	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	1334	. "ave(X)", zstophy,zout)
	1335	c
	1336	CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-",
	1337	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	1338	. "ave(X)", zstophy,zout)
	1339	c
	1340	CALL histdef(nid_mth, "cldt", "Total cloudiness", "%",
	1341	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	1342	. "ave(X)", zstophy,zout)
	1343	c
	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)
	1347	c
	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)
	1351	c
	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)
	1355	c
	1356	CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
	1357	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	1358	. "ave(X)", zstophy,zout)
	1359	c
	1360	CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
	1361	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	1362	. "ave(X)", zstophy,zout)
	1363	c
	1364	CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-",
	1365	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	1366	. "ave(X)", zstophy,zout)
	1367	c
	1368	CALL histdef(nid_mth, "vq", "Merid humidity transport", "-",
	1369	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	1370	. "ave(X)", zstophy,zout)
	1371	cKE43
	1372	IF(iflag_con.GE.3) THEN ! sb
	1373	c
	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)
	1377	c
	1378	CALL histdef(nid_mth, "pbase", "Cld base pressure", "mb",
	1379	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	1380	. "ave(X)", zstophy,zout)
	1381	c
	1382	CALL histdef(nid_mth, "ptop", "Cld top pressure", "mb",
	1383	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	1384	. "ave(X)", zstophy,zout)
	1385	c
	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)
	1389	c
	1390	CALL histdef(nid_mth, "prw", "Precipitable water", "kg/m2",
	1391	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	1392	. "ave(X)", zstophy,zout)
	1393	c
[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)
	1397	c
[644]	1398	ENDIF !iflag_con .GE. 3
	1399	c34EK
	1400	cIM cf. AM 081204 BEG
	1401	c34EK
	1402	c 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)
	1436	c
	1437	cIM cf. AM 081204 END
	1438	c
	1439	c Champs interpolles sur des niveaux de pression
	1440	DO k=1, nlevSTD
	1441	c
	1442	c bb=clevSTD(k)
	1443	c
	1444	c 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)
	1447	c aa=clevSTD(k)
	1448	c bb=aa(1:lnblnk1(aa))
	1449	c ENDIF
	1450	c
	1451	IF(bb2.EQ."850".OR.bb2.EQ."700"
	1452	$ .OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
	1453	c
	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)
	1458	c
	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)
	1463	c
	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)
	1468	c
	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)
	1473	c
	1474	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700"
	1475	c
	1476	ENDDO !nlevSTD
	1477	c
	1478	cIM diagnostiques flux ocean-atm ou ocean-glace de mer
	1479	cIM pour utilisation dans un modele de "slab" ocean
	1480	c
	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)
	1485	c
	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)
	1490	c
	1491	IF (OCEAN.EQ.'force ') THEN
	1492	c
	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)
	1497	c
	1498	ELSE IF (OCEAN.EQ.'slab ') THEN
	1499	c
	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)
	1504	c
	1505	CALL histdef(nid_mth, "tslab", "Slab SST ", "K",
	1506	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	1507	. "ave(X)", zstophy,zout)
	1508	c
	1509	CALL histdef(nid_mth, "seaice", "Slab seaice", "kg/m2",
	1510	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	1511	. "ave(X)", zstophy,zout)
	1512	c
	1513	CALL histdef(nid_mth, "siceh", "Slab seaice height", "m",
	1514	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	1515	. "ave(X)", zstophy,zout)
	1516	c
	1517	ENDIF !OCEAN.EQ.slab/force
	1518	c
	1519	ENDIF !lev_histmth.GE.1
	1520	c
	1521	IF(lev_histmth.GE.2) THEN
	1522	c
	1523	c Champs 3D:
	1524	c
	1525	c
	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)
	1530	c
	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)
	1534	c
	1535	CALL histdef(nid_mth, "temp", "Air temperature", "K",
	1536	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	1537	. "ave(X)", zstophy,zout)
	1538	c
	1539	CALL histdef(nid_mth, "ovap", "Specific humidity", "kg/kg",
	1540	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	1541	. "ave(X)", zstophy,zout)
	1542	c
	1543	c CALL histdef(nid_mth,"wvap","Water vapor mixing ratio","kg/kg",
	1544	c . iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	1545	c . "ave(X)", zstophy,zout)
	1546	c
	1547	CALL histdef(nid_mth, "geop", "Geopotential height", "m2/s2",
	1548	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	1549	. "ave(X)", zstophy,zout)
	1550	c
	1551	CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
	1552	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	1553	. "ave(X)", zstophy,zout)
	1554	c
	1555	CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
	1556	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	1557	. "ave(X)", zstophy,zout)
	1558	c
	1559	CALL histdef(nid_mth, "vitw", "Vertical wind", "m/s",
	1560	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	1561	. "ave(X)", zstophy,zout)
	1562	c
	1563	CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
	1564	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	1565	. "ave(X)", zstophy,zout)
	1566	c
	1567	CALL histdef(nid_mth, "rneb", "Cloud fraction", "%",
	1568	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	1569	. "ave(X)", zstophy,zout)
	1570	c
	1571	CALL histdef(nid_mth, "rnebcon", "Convective Cloud Fraction"
	1572	. , "-",
	1573	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	1574	. "ave(X)", zstophy,zout)
	1575	c
	1576	CALL histdef(nid_mth, "rhum", "Relative humidity", "-",
	1577	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	1578	. "ave(X)", zstophy,zout)
	1579	c
	1580	CALL histdef(nid_mth, "ozone", "Ozone concentration", "ppmv",
	1581	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	1582	. "ave(X)", zstophy,zout)
	1583	c
	1584	CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
	1585	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	1586	. "ave(X)", zstophy,zout)
	1587	c
	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)
	1591	c
	1592	ENDIF !lev_histmth.GE.2
	1593	c
	1594	IF(lev_histmth.GE.3) THEN
	1595	c
	1596	DO nsrf=1, nbsrf
	1597	c
	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)
	1602	c
	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)
	1607	c
	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)
	1611	c
	1612	ENDDO !nsrf=1, nbsrf
	1613	c
	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)
	1620	c
	1621	ENDIF !lev_histmth.GE.3
	1622	c
	1623	IF(lev_histmth.GE.4) THEN
	1624	c
	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)
	1630	c
	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)
	1634	c
	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)
	1638	c
	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)
	1642	c
	1643	CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
	1644	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	1645	. "ave(X)", zstophy,zout)
	1646	c
	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)
	1650	c
	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)
	1654	c
	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)
	1658	c
	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)
	1662	c
	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)
	1666	c
	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)
	1670	c
	1671	CALL histdef(nid_mth, "dteva", "Reevaporation dT", "K/s",
	1672	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	1673	. "ave(X)", zstophy,zout)
	1674	c
	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
	1687	c
	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)
	1691	c
	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)
	1695	c
	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)
	1699	c
	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)
	1703	c
	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)
	1707	c
	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)
	1711	c
	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)
	1715	c
	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)
	1719	c
	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)
	1724	c
	1725	CALL histdef(nid_mth, "dvoro", "Orography dV", "m/s2",
	1726	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	1727	. "ave(X)", zstophy,zout)
	1728	c
	1729	ENDIF
	1730	C
	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)
	1735	c
	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
	1740	C
[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
	1751	c
[644]	1752	ENDIF !type_run.EQ."CLIM".OR.type_run.EQ."ENSP")
	1753	c
[524]	1754	CALL histend(nid_mth)
	1755	c
	1756	ndex2d = 0
	1757	ndex3d = 0
	1758	c
	1759	ENDIF ! fin de test sur ok_mensuel

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

Download in other formats: