Context Navigation

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

Last change on this file since 777 was 776, checked in by Laurent Fairhead, 17 years ago
Suite du merge entre la version et la HEAD: quelques modifications de Yann sur le LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 68.5 KB

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

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

Download in other formats: