Context Navigation

source: LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/ini_histmth.h @ 892

Last change on this file since 892 was 892, checked in by Laurent Fairhead, 16 years ago
Mise sous flag ok_ISCCP d'une sortie d'une variable dependant du simulateur ISCCP Correction coquille sur des sorties sous-surfaces SD/LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 72.7 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",
[850]	287	. iim,jj_nb,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,
[850]	312	. "ave(X)", zstophy,zout)
[766]	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
[892]	328	IF(nsrf.GE.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)
[892]	338	ENDIF !nsrf.GE.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
[869]	632	CALL histdef(nid_mth, "vitw", "Vertical wind", "Pa/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
[862]	834	#if defined(INCA_AER) && defined(CPP_COUPLE)
	835	CALL histdef(nid_mth, "tau", "cloud opt thickness", "-",
	836	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
	837	. "ave(X)", zstophy,zout)
	838	CALL histdef(nid_mth, "CCM1", "CCM1", "",
	839	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
	840	. "ave(X)", zstophy,zout)
	841	CALL histdef(nid_mth, "CCM2", "CCM2", "",
	842	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
	843	. "ave(X)", zstophy,zout)
	844	CALL histdef(nid_mth, "taupi", "cloud opt thickn. (pi)", "-",
	845	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
	846	. "ave(X)", zstophy,zout)
	847	CALL histdef(nid_mth, "cldemi", "emissivite IR", "-",
	848	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
	849	. "ave(X)", zstophy,zout)
	850	CALL histdef(nid_mth, "SWTOAAS_ZERO", "?", "?",
	851	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	852	. "ave(X)", zstophy,zout)
	853	CALL histdef(nid_mth, "SWTOAAS_AER", "?", "?",
	854	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	855	. "ave(X)", zstophy,zout)
	856	CALL histdef(nid_mth, "SWTOAAS_NAT", "?", "?",
	857	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	858	. "ave(X)", zstophy,zout)
	859	CALL histdef(nid_mth, "SWTOAAS_AD", "?", "?",
	860	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	861	. "ave(X)", zstophy,zout)
	862	CALL histdef(nid_mth, "SWTOAAS_AI", "?", "?",
	863	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	864	. "ave(X)", zstophy,zout)
	865	CALL histdef(nid_mth, "SWTOACS_ZERO", "?", "?",
	866	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	867	. "ave(X)", zstophy,zout)
	868	CALL histdef(nid_mth, "SWTOACS_AER", "?", "?",
	869	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	870	. "ave(X)", zstophy,zout)
	871	CALL histdef(nid_mth, "SWTOACS_AD", "?", "?",
	872	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	873	. "ave(X)", zstophy,zout)
	874	CALL histdef(nid_mth, "SWTOACS_NAT", "?", "?",
	875	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	876	. "ave(X)", zstophy,zout)
	877	#endif
	878	#ifndef INCA
[644]	879	if (nqmax.GE.3) THEN
[845]	880	DO iq=3,nqmax
	881	iiq=niadv(iq)
	882	CALL histdef(nid_mth, tnom(iq), ttext(iiq), "-",
[776]	883	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	884	. "ave(X)", zstophy,zout)
[845]	885	ENDDO
[644]	886	ENDIF
[862]	887	#endif
[524]	888	c
[644]	889	ENDIF !lev_histmth.GE.4
[524]	890	c
[644]	891	ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
	892	c
	893	IF(lev_histmth.GE.1) THEN
	894	c
	895	CALL histdef(nid_mth, "phis", "Surface geop. height", "m",
[776]	896	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	897	. "once", zstophy,zout)
	898	c
	899	CALL histdef(nid_mth, "aire", "Grid area", "-",
[776]	900	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	901	. "once", zstophy,zout)
	902	c
	903	CALL histdef(nid_mth, "pourc_"//clnsurf(is_ter),
	904	. "% "//clnsurf(is_ter), "%",
[776]	905	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	906	. "once", zstophy,zout)
	907	c
	908	CALL histdef(nid_mth, "fract_"//clnsurf(is_ter),
	909	. "Fraction "//clnsurf(is_ter), "1",
[776]	910	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	911	. "once", zstophy,zout)
	912	c
	913	CALL histdef(nid_mth, "pourc_"//clnsurf(is_lic),
	914	. "% "//clnsurf(is_lic), "%",
[776]	915	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	916	. "once", zstophy,zout)
	917	c
	918	CALL histdef(nid_mth, "fract_"//clnsurf(is_lic),
	919	. "Fraction "//clnsurf(is_lic), "1",
[776]	920	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	921	. "once", zstophy,zout)
	922	c
	923	CALL histdef(nid_mth, "pourc_"//clnsurf(is_sic),
	924	. "% "//clnsurf(is_sic), "%",
[776]	925	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	926	. "ave(X)", zstophy,zout)
	927	c
	928	CALL histdef(nid_mth, "fract_"//clnsurf(is_sic),
	929	. "Fraction "//clnsurf(is_sic), "1",
[776]	930	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	931	. "ave(X)", zstophy,zout)
	932	c
	933	CALL histdef(nid_mth, "flat", "Latent heat flux", "W/m2",
[776]	934	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	935	. "ave(X)", zstophy,zout)
	936	c
	937	CALL histdef(nid_mth, "bilTOA",
	938	. "Net radiation at TOA", "W/m2",
[776]	939	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	940	. "ave(X)", zstophy,zout)
	941	c
	942	c
	943	c CALL histdef(nid_mth, "tops01",
	944	c . "TOA reflected clear-sky shortwave
	945	c . radiation (method 1)", "W/m2",
[776]	946	c . iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	947	c . "ave(X)", zstophy,zout)
	948	c
	949	CALL histdef(nid_mth, "dtdia","dtdia","K/s",
[776]	950	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	951	. "ave(X)", zstophy,zout)
[524]	952	c
[644]	953	CALL histdef(nid_mth, "dtsw","dtsw","K/s",
[776]	954	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	955	. "ave(X)", zstophy,zout)
[524]	956	c
[644]	957	CALL histdef(nid_mth, "dtlw","dtlw","K/s",
[776]	958	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	959	. "ave(X)", zstophy,zout)
[524]	960	c
[644]	961	CALL histdef(nid_mth, "dtcon","dtcon","K/s",
[776]	962	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	963	. "ave(X)", zstophy,zout)
[524]	964	c
[644]	965	CALL histdef(nid_mth, "dtajs","dtajs","K/s",
[776]	966	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	967	. "ave(X)", zstophy,zout)
	968	c
	969	CALL histdef(nid_mth, "dtlsc","dtlsc","K/s",
[776]	970	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	971	. "ave(X)", zstophy,zout)
	972	c
	973	CALL histdef(nid_mth, "dqdia","dqdia","(kg/kg)/s",
[776]	974	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	975	. "ave(X)", zstophy,zout)
	976	c
	977	CALL histdef(nid_mth, "duoli",
	978	. "duoli","m/s2",
[776]	979	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	980	. "ave(X)", zstophy,zout)
	981	c
	982	CALL histdef(nid_mth, "dvoli",
	983	. "dvoli","m/s2",
[776]	984	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	985	. "ave(X)", zstophy,zout)
	986	c
	987	CALL histdef(nid_mth, "ducon",
	988	. "ducon","m/s2",
[776]	989	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	990	. "ave(X)", zstophy,zout)
	991	c
	992	CALL histdef(nid_mth, "dvcon",
	993	. "dvcon","m/s2",
[776]	994	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	995	. "ave(X)", zstophy,zout)
	996	c
	997	CALL histdef(nid_mth, "dutot",
	998	. "dutot","m/s2",
[776]	999	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1000	. "ave(X)", zstophy,zout)
	1001	c
	1002	CALL histdef(nid_mth, "dvtot",
	1003	. "dvtot","m/s2",
[776]	1004	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1005	. "ave(X)", zstophy,zout)
	1006	c
	1007	c CALL histdef(nid_mth, "cldtrue",
	1008	c . "TRUE Cloud fraction ","%",
[776]	1009	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1010	c . "inst(X)", zout,zout)
	1011	c
	1012	c CALL histdef(nid_mth, "cldTOA",
	1013	c . "Cloud amount (satellite view)","%",
[776]	1014	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1015	c . "ave(X)", zstophy,zout)
	1016	c
	1017	c CALL histdef(nid_mth, "cldSFC",
	1018	c . "Cloud amount (surface view)","%",
[776]	1019	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1020	c . "ave(X)", zstophy,zout)
	1021	c
	1022	CALL histdef(nid_mth, "GWDtaux",
	1023	. "GWD induced Ewd wind stress","Pa",
[776]	1024	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1025	. "ave(X)", zstophy,zout)
	1026	c
	1027	CALL histdef(nid_mth, "GWDtauy",
	1028	. "GWD induced Nwd wind stress","Pa",
[776]	1029	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1030	. "ave(X)", zstophy,zout)
	1031	c
	1032	c Champs 2D:
	1033	c
	1034	CALL histdef(nid_mth, "slp", "Sea Level Pressure", "Pa",
[776]	1035	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1036	. "ave(X)", zstorad,zout)
	1037	c
	1038	CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
[776]	1039	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1040	. "ave(X)", zstophy,zout)
	1041	c
	1042	CALL histdef(nid_mth, "t2m", "Temperature 2m", "K",
[776]	1043	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1044	. "ave(X)", zstorad,zout)
	1045	c
	1046	CALL histdef(nid_mth, "t2m_min", "Temp. 2m min.",
	1047	. "K",
[776]	1048	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1049	. t2mincels, zstophy,zout)
	1050	c
	1051	CALL histdef(nid_mth, "t2m_max", "Temp. 2m max.",
	1052	. "K",
[776]	1053	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1054	. t2maxcels, zstophy,zout)
	1055	c
	1056	c CALL histdef(nid_mth, "tsoil", "Sfce soil Temperature",
	1057	c . "K",
[776]	1058	c . iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1059	c . "ave(X)", zstophy,zout)
	1060	c
	1061	CALL histdef(nid_mth, "wind10m","10-m wind speed","m/s",
[776]	1062	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1063	. "ave(X)", zstorad,zout)
	1064	c
	1065	CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "-",
[776]	1066	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1067	. "ave(X)", zstophy,zout)
	1068
	1069	c
	1070	CALL histdef(nid_mth, "q2m","Specific humidity 2m","kg/kg",
[776]	1071	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1072	. "ave(X)", zstorad,zout)
	1073	c
	1074	CALL histdef(nid_mth, "u10m", "Vent zonal 10m", "m/s",
[776]	1075	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1076	. "ave(X)", zstorad,zout)
	1077	c
	1078	CALL histdef(nid_mth, "v10m", "Vent meridien 10m", "m/s",
[776]	1079	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1080	. "ave(X)", zstorad,zout)
	1081	c
	1082	CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
[776]	1083	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1084	. "ave(X)", zstorad,zout)
	1085	c
	1086	CALL histdef(nid_mth, "qsurf","Surface Air humidity","kg/kg",
[776]	1087	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1088	. "ave(X)", zstophy,zout)
	1089	c
	1090	if (.not. ok_veget) then
	1091	CALL histdef(nid_mth, "qsol", "Soil watter content", "mm",
[776]	1092	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1093	. "ave(X)", zstophy,zout)
	1094	endif
	1095	c
	1096	CALL histdef(nid_mth, "ndayrain",
	1097	. "Number of day with rain (liq+sol)", "-",
[776]	1098	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1099	. "inst(X)", zstomth,zout)
	1100	c
	1101	CALL histdef(nid_mth, "precip", "Precipitation Totale liq+sol",
	1102	. "kg/(s*m2)",
[776]	1103	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1104	. "ave(X)", zstophy,zout)
	1105	c
	1106	CALL histdef(nid_mth, "plul", "Large-scale Precip.",
	1107	. "kg/(s*m2)",
[776]	1108	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1109	. "ave(X)", zstophy,zout)
	1110	c
	1111	CALL histdef(nid_mth, "pluc", "Convective Precip.",
	1112	. "kg/(s*m2)",
[776]	1113	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1114	. "ave(X)", zstophy,zout)
	1115	c
	1116	CALL histdef(nid_mth, "snowf", "Snow fall", "kg/(s*m2)",
[776]	1117	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1118	. "ave(X)", zstophy,zout)
	1119	c
	1120	c CALL histdef(nid_mth, "snow_mass", "Snow Mass", "kg/m2",
[776]	1121	c . iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1122	c . "ave(X)", zstophy,zout)
	1123	c
	1124	CALL histdef(nid_mth, "evap",
	1125	. "Evaporation plus sublimation rate", "kg/(s*m2)",
[776]	1126	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1127	. "ave(X)", zstophy,zout)
	1128	c
	1129	CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
[776]	1130	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1131	. "ave(X)", zstophy,zout)
	1132	c
	1133	CALL histdef(nid_mth, "tops0", "CS Solar rad. at TOA", "W/m2",
[776]	1134	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1135	. "ave(X)", zstophy,zout)
	1136	c
	1137	CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
[776]	1138	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1139	. "ave(X)", zstophy,zout)
	1140	c
	1141	CALL histdef(nid_mth, "topl0", "IR rad. at TOA", "W/m2",
[776]	1142	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1143	. "ave(X)", zstophy,zout)
	1144	c
	1145	CALL histdef(nid_mth, "SWupTOA", "SWup at TOA","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, "SWupTOAclr",
	1150	. "SWup clear sky at TOA","W/m2",
[776]	1151	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	1152	. 32, "ave(X)", zstophy,zout)
	1153	c
	1154	CALL histdef(nid_mth, "SWdnTOA", "SWdn at TOA","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, "SWdnTOAclr",
	1159	. "SWdn clear sky at TOA","W/m2",
[776]	1160	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	1161	. 32, "ave(X)", zstophy,zout)
	1162	c
	1163	CALL histdef(nid_mth, "SWup200", "SWup at 200mb","W/m2",
[776]	1164	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	1165	. 32, "ave(X)", zstophy,zout)
	1166	c
	1167	CALL histdef(nid_mth, "SWup200clr",
	1168	. "SWup clear sky at 200mb","W/m2",
[776]	1169	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	1170	. 32, "ave(X)", zstophy,zout)
	1171	c
	1172	CALL histdef(nid_mth, "SWdn200", "SWdn at 200mb","W/m2",
[776]	1173	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	1174	. 32, "ave(X)", zstophy,zout)
	1175	c
	1176	CALL histdef(nid_mth, "SWdn200clr",
	1177	. "SWdn clear sky at 200mb","W/m2",
[776]	1178	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	1179	. 32, "ave(X)", zstophy,zout)
	1180	c
	1181	CALL histdef(nid_mth, "LWup200", "LWup at 200mb","W/m2",
[776]	1182	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	1183	. 32, "ave(X)", zstophy,zout)
	1184	c
	1185	CALL histdef(nid_mth, "LWup200clr",
	1186	. "LWup clear sky at 200mb","W/m2",
[776]	1187	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	1188	. 32, "ave(X)", zstophy,zout)
	1189	c
	1190	CALL histdef(nid_mth, "LWdn200", "LWdn at 200mb","W/m2",
[776]	1191	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	1192	. 32, "ave(X)", zstophy,zout)
	1193	c
	1194	CALL histdef(nid_mth, "LWdn200clr",
	1195	. "LWdn clear sky at 200mb","W/m2",
[776]	1196	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	1197	. 32, "ave(X)", zstophy,zout)
	1198	c
	1199	CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
[776]	1200	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1201	. "ave(X)", zstophy,zout)
	1202	c
	1203	CALL histdef(nid_mth, "sols0", "Solar rad. at surf.", "W/m2",
[776]	1204	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1205	. "ave(X)", zstophy,zout)
	1206	c
	1207	CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
[776]	1208	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1209	. "ave(X)", zstophy,zout)
	1210	c
[684]	1211	CALL histdef(nid_mth, "radsol", "Rayonnement au sol", "W/m2",
[776]	1212	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[684]	1213	. "ave(X)", zstophy,zout)
	1214	c
[644]	1215	CALL histdef(nid_mth, "soll0", "IR rad. at surface", "W/m2",
[776]	1216	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1217	. "ave(X)", zstophy,zout)
	1218	c
	1219	CALL histdef(nid_mth, "SWupSFC", "SWup at surface","W/m2",
[776]	1220	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	1221	. 32, "ave(X)", zstophy,zout)
	1222	c
	1223	CALL histdef(nid_mth, "SWupSFCclr",
	1224	. "SWup clear sky at surface","W/m2",
[776]	1225	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	1226	. 32, "ave(X)", zstophy,zout)
	1227	c
	1228	CALL histdef(nid_mth, "SWdnSFC", "SWdn at surface","W/m2",
[776]	1229	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	1230	. 32, "ave(X)", zstophy,zout)
	1231	c
	1232	CALL histdef(nid_mth, "SWdnSFCclr",
	1233	. "SWdn clear sky at surface","W/m2",
[776]	1234	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	1235	. 32, "ave(X)", zstophy,zout)
	1236	c
	1237	CALL histdef(nid_mth,"LWupSFC","Upwd. IR rad. at surface",
[776]	1238	. "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1239	. "ave(X)", zstophy,zout)
	1240	c
	1241	CALL histdef(nid_mth,"LWdnSFC","Down. IR rad. at surface",
[776]	1242	. "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1243	. "ave(X)", zstophy,zout)
	1244	c
	1245	CALL histdef(nid_mth,"LWupSFCclr",
	1246	. "CS Upwd. IR rad. at surface",
[776]	1247	. "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1248	. "ave(X)", zstophy,zout)
	1249	c
	1250	CALL histdef(nid_mth,"LWdnSFCclr",
	1251	. "Down. CS IR rad. at surface",
[776]	1252	. "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1253	. "ave(X)", zstophy,zout)
	1254	c
	1255	CALL histdef(nid_mth, "bils", "Surf. total heat flux", "W/m2",
[776]	1256	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1257	. "ave(X)", zstophy,zout)
	1258	c
	1259	CALL histdef(nid_mth, "sens", "Sensible heat flux", "W/m2",
[776]	1260	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1261	. "ave(X)", zstophy,zout)
	1262	c
	1263	CALL histdef(nid_mth, "fder", "Heat flux derivation", "W/m2",
[776]	1264	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1265	. "ave(X)", zstophy,zout)
	1266	c
	1267	c CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
[776]	1268	c . iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1269	c . "ave(X)", zstophy,zout)
	1270	c
	1271	c CALL histdef(nid_mth, "frtv", "Meridional wind stress", "Pa",
[776]	1272	c . iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1273	c . "ave(X)", zstophy,zout)
	1274	c
	1275	CALL histdef(nid_mth, "ffonte","Thermal flux for snow melting",
[776]	1276	. "W/m2",iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1277	. "ave(X)", zstophy,zout)
	1278	c
	1279	CALL histdef(nid_mth, "fqcalving","Ice Calving",
[776]	1280	. "kg/m2/s",iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1281	. "ave(X)", zstophy,zout)
	1282	c
	1283	DO nsrf = 1, nbsrf
	1284	call histdef(nid_mth, "taux_"//clnsurf(nsrf),
	1285	$ "Zonal wind stress"//clnsurf(nsrf), "Pa",
[776]	1286	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1287	$ "ave(X)", zstophy,zout)
	1288
	1289	call histdef(nid_mth, "tauy_"//clnsurf(nsrf),
	1290	$ "Meridional wind stress "//clnsurf(nsrf), "Pa",
[776]	1291	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1292	$ "ave(X)", zstophy,zout)
	1293	ENDDO
	1294	c
	1295	call histdef(nid_mth, "pourc_"//clnsurf(is_oce),
	1296	$ "% "//clnsurf(is_oce), "%",
[776]	1297	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1298	$ "ave(X)", zstophy,zout)
	1299	c
	1300	call histdef(nid_mth, "fract_"//clnsurf(is_oce),
	1301	$ "Fraction "//clnsurf(is_oce), "1",
[776]	1302	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1303	$ "ave(X)", zstophy,zout)
	1304	C
	1305	DO nsrf = 1, nbsrf
	1306	C
	1307	call histdef(nid_mth, "tsol_"//clnsurf(nsrf),
	1308	$ "Temperature "//clnsurf(nsrf), "K",
[776]	1309	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1310	$ "ave(X)", zstophy,zout)
	1311	C
	1312	call histdef(nid_mth, "sens_"//clnsurf(nsrf),
	1313	$ "Sensible heat flux "//clnsurf(nsrf), "W/m2",
[776]	1314	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1315	$ "ave(X)", zstophy,zout)
	1316	c
	1317	call histdef(nid_mth, "lat_"//clnsurf(nsrf),
	1318	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
[776]	1319	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1320	$ "ave(X)", zstophy,zout)
	1321	C
	1322	call histdef(nid_mth, "flw_"//clnsurf(nsrf),
	1323	$ "LW "//clnsurf(nsrf), "W/m2",
[776]	1324	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1325	$ "ave(X)", zstophy,zout)
	1326	c
	1327	call histdef(nid_mth, "fsw_"//clnsurf(nsrf),
	1328	$ "SW "//clnsurf(nsrf), "W/m2",
[776]	1329	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1330	$ "ave(X)", zstophy,zout)
	1331	C
	1332	call histdef(nid_mth, "wbils_"//clnsurf(nsrf),
	1333	$ "Bilan sol "//clnsurf(nsrf), "W/m2",
[776]	1334	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1335	$ "ave(X)", zstophy,zout)
	1336	C
[684]	1337	call histdef(nid_mth, "wbilo_"//clnsurf(nsrf),
	1338	$ "Bilan eau "//clnsurf(nsrf), "kg/(m2*s)",
[776]	1339	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[684]	1340	$ "ave(X)", zstophy,zout)
	1341	C
[644]	1342	cIM 250105 BEG
	1343	IF(1.EQ.0) THEN
	1344	C
	1345	call histdef(nid_mth, "evapour_"//clnsurf(nsrf),
	1346	$ "Evappourc "//clnsurf(nsrf), "kg/(sm2)",
[776]	1347	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1348	$ "ave(X)", zstophy,zout)
	1349	C
	1350	call histdef(nid_mth, "prepour_"//clnsurf(nsrf),
	1351	$ "Precippourc "//clnsurf(nsrf), "kg/(sm2)",
[776]	1352	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1353	$ "ave(X)", zstophy,zout)
	1354	C
	1355	ENDIF !(1.EQ.0) THEN
	1356	cIM 250105 END
	1357	C
	1358	END DO
	1359	c
	1360	CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
[776]	1361	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1362	. "ave(X)", zstophy,zout)
	1363	c
	1364	CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
[776]	1365	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1366	. "ave(X)", zstophy,zout)
	1367	c
	1368	CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-",
[776]	1369	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1370	. "ave(X)", zstophy,zout)
	1371	c
	1372	CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-",
[776]	1373	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1374	. "ave(X)", zstophy,zout)
	1375	c
	1376	CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-",
[776]	1377	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1378	. "ave(X)", zstophy,zout)
	1379	c
	1380	CALL histdef(nid_mth, "cldt", "Total cloudiness", "%",
[776]	1381	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1382	. "ave(X)", zstophy,zout)
	1383	c
	1384	CALL histdef(nid_mth,"cldq","Cloud liquid water path","kg/m2",
[776]	1385	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1386	. "ave(X)", zstophy,zout)
	1387	c
	1388	CALL histdef(nid_mth,"lwp","Cloud water path","kg/m2",
[776]	1389	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1390	. "ave(X)", zstophy,zout)
	1391	c
	1392	CALL histdef(nid_mth,"iwp","Cloud ice water path","kg/m2",
[776]	1393	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1394	. "ave(X)", zstophy,zout)
	1395	c
	1396	CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
[776]	1397	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1398	. "ave(X)", zstophy,zout)
	1399	c
	1400	CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
[776]	1401	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1402	. "ave(X)", zstophy,zout)
	1403	c
	1404	CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-",
[776]	1405	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1406	. "ave(X)", zstophy,zout)
	1407	c
	1408	CALL histdef(nid_mth, "vq", "Merid humidity transport", "-",
[776]	1409	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1410	. "ave(X)", zstophy,zout)
	1411	cKE43
	1412	IF(iflag_con.GE.3) THEN ! sb
	1413	c
	1414	CALL histdef(nid_mth, "cape", "Conv avlbl pot ener", "J/kg",
[776]	1415	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1416	. "ave(X)", zstophy,zout)
	1417	c
	1418	CALL histdef(nid_mth, "pbase", "Cld base pressure", "mb",
[776]	1419	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1420	. "ave(X)", zstophy,zout)
	1421	c
	1422	CALL histdef(nid_mth, "ptop", "Cld top pressure", "mb",
[776]	1423	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1424	. "ave(X)", zstophy,zout)
	1425	c
	1426	CALL histdef(nid_mth, "fbase", "Cld base mass flux", "kg/m2/s",
[776]	1427	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1428	. "ave(X)", zstophy,zout)
	1429	c
	1430	CALL histdef(nid_mth, "prw", "Precipitable water", "kg/m2",
[776]	1431	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1432	. "ave(X)", zstophy,zout)
	1433	c
[653]	1434	CALL histdef(nid_mth, "upwd", "saturated updraft", "kg/m2/s",
[776]	1435	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[653]	1436	. "ave(X)", zstophy,zout)
	1437	c
[644]	1438	ENDIF !iflag_con .GE. 3
	1439	c34EK
	1440	cIM cf. AM 081204 BEG
	1441	c34EK
	1442	c HBTM2
	1443	CALL histdef(nid_mth, "s_pblh", "Boundary Layer Height", "m",
[776]	1444	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1445	. "ave(X)", zstophy,zout)
	1446	CALL histdef(nid_mth, "s_pblt", "T at Boundary Layer Height",
	1447	. "K",
[776]	1448	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1449	. "ave(X)", zstophy,zout)
	1450	CALL histdef(nid_mth, "s_lcl", "Condensation level", "m",
[776]	1451	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1452	. "ave(X)", zstophy,zout)
	1453	CALL histdef(nid_mth,"s_capCL","Conv avlbl pot ener for ABL",
	1454	. "J/m2",
[776]	1455	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1456	. "ave(X)", zstophy,zout)
	1457	CALL histdef(nid_mth, "s_oliqCL", "Liq Water in BL", "kg/m2",
[776]	1458	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1459	. "ave(X)", zstophy,zout)
	1460	CALL histdef(nid_mth,"s_cteiCL","Instability criteria (ABL)",
	1461	. "K",
[776]	1462	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1463	. "ave(X)", zstophy,zout)
	1464	CALL histdef(nid_mth, "s_therm", "Exces du thermique", "K",
[776]	1465	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1466	. "ave(X)", zstophy,zout)
	1467	CALL histdef(nid_mth, "s_trmb1", "deep_cape(HBTM2)", "J/m2",
[776]	1468	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1469	. "ave(X)", zstophy,zout)
	1470	CALL histdef(nid_mth, "s_trmb2", "inhibition (HBTM2)", "J/m2",
[776]	1471	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1472	. "ave(X)", zstophy,zout)
	1473	CALL histdef(nid_mth, "s_trmb3", "Point Omega (HBTM2)", "m",
[776]	1474	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1475	. "ave(X)", zstophy,zout)
	1476	c
	1477	cIM cf. AM 081204 END
	1478	c
	1479	c Champs interpolles sur des niveaux de pression
	1480	DO k=1, nlevSTD
	1481	c
	1482	c bb=clevSTD(k)
	1483	c
	1484	c IF(k.GE.2) THEN
	1485	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
	1486	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
	1487	c aa=clevSTD(k)
	1488	c bb=aa(1:lnblnk1(aa))
	1489	c ENDIF
	1490	c
	1491	IF(bb2.EQ."850".OR.bb2.EQ."700"
	1492	$ .OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
	1493	c
	1494	CALL histdef(nid_mth, "u"//bb2,
	1495	. "Zonal wind "//bb2//"mb","m/s",
[776]	1496	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1497	. "inst(X)", zout,zout)
	1498	c
	1499	CALL histdef(nid_mth, "v"//bb2,
	1500	. "Meridional wind "//bb2//"mb","m/s",
[776]	1501	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1502	. "inst(X)", zout,zout)
	1503	c
	1504	CALL histdef(nid_mth, "w"//bb2,
	1505	. "Vertical wind "//bb2//"mb","m/s",
[776]	1506	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1507	. "inst(X)", zout,zout)
	1508	c
	1509	CALL histdef(nid_mth, "phi"//bb2,
	1510	. "Geopotential "//bb2//"mb","m",
[776]	1511	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1512	. "inst(X)", zout,zout)
	1513	c
	1514	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700"
	1515	c
	1516	ENDDO !nlevSTD
	1517	c
	1518	cIM diagnostiques flux ocean-atm ou ocean-glace de mer
	1519	cIM pour utilisation dans un modele de "slab" ocean
	1520	c
	1521	CALL histdef(nid_mth, "fluxo",
	1522	. "Flux turbulents ocean-atmosphere", "W/m2",
[776]	1523	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1524	. "ave(X)", zstophy,zout)
	1525	c
	1526	CALL histdef(nid_mth, "fluxg",
	1527	. "Flux turbulents ocean-glace de mer","W/m2",
[776]	1528	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1529	. "ave(X)", zstophy,zout)
	1530	c
[684]	1531	CALL histdef(nid_mth, "t_oce_sic",
	1532	. "Temp. mixte oce-sic","K",
[776]	1533	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[684]	1534	. "ave(X)", zstophy,zout)
	1535	c
[644]	1536	IF (OCEAN.EQ.'force ') THEN
	1537	c
	1538	CALL histdef(nid_mth, "lmt_bils",
	1539	. "Bilan au sol atmosphere forcee", "W/m2",
[776]	1540	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1541	. "ave(X)", zstophy,zout)
	1542	c
	1543	ELSE IF (OCEAN.EQ.'slab ') THEN
	1544	c
	1545	CALL histdef(nid_mth, "slab_bils",
	1546	. "Bilan au sol Slab", "W/m2",
[776]	1547	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1548	. "ave(X)", zstophy,zout)
	1549	c
	1550	CALL histdef(nid_mth, "tslab", "Slab SST ", "K",
[776]	1551	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1552	. "ave(X)", zstophy,zout)
	1553	c
	1554	CALL histdef(nid_mth, "seaice", "Slab seaice", "kg/m2",
[776]	1555	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1556	. "ave(X)", zstophy,zout)
	1557	c
	1558	CALL histdef(nid_mth, "siceh", "Slab seaice height", "m",
[776]	1559	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1560	. "ave(X)", zstophy,zout)
	1561	c
	1562	ENDIF !OCEAN.EQ.slab/force
	1563	c
	1564	ENDIF !lev_histmth.GE.1
	1565	c
	1566	IF(lev_histmth.GE.2) THEN
	1567	c
	1568	c Champs 3D:
	1569	c
	1570	c
	1571	CALL histdef(nid_mth,"lwcon","Cloud liquid water content",
	1572	. "kg/kg",
[776]	1573	. iim,jj_nb,nhori, klev,1,klev, nvert, 32,
[644]	1574	. "ave(X)", zstophy,zout)
	1575	c
	1576	CALL histdef(nid_mth,"iwcon","Cloud ice water content","kg/kg",
[776]	1577	. iim,jj_nb,nhori, klev,1,klev, nvert, 32,
[644]	1578	. "ave(X)", zstophy,zout)
	1579	c
	1580	CALL histdef(nid_mth, "temp", "Air temperature", "K",
[776]	1581	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1582	. "ave(X)", zstophy,zout)
	1583	c
	1584	CALL histdef(nid_mth, "ovap", "Specific humidity", "kg/kg",
[776]	1585	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1586	. "ave(X)", zstophy,zout)
	1587	c
	1588	c CALL histdef(nid_mth,"wvap","Water vapor mixing ratio","kg/kg",
[776]	1589	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1590	c . "ave(X)", zstophy,zout)
	1591	c
	1592	CALL histdef(nid_mth, "geop", "Geopotential height", "m2/s2",
[776]	1593	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1594	. "ave(X)", zstophy,zout)
	1595	c
	1596	CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
[776]	1597	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1598	. "ave(X)", zstophy,zout)
	1599	c
	1600	CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
[776]	1601	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1602	. "ave(X)", zstophy,zout)
	1603	c
[869]	1604	CALL histdef(nid_mth, "vitw", "Vertical wind", "Pa/s",
[776]	1605	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1606	. "ave(X)", zstophy,zout)
	1607	c
	1608	CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
[776]	1609	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1610	. "ave(X)", zstophy,zout)
	1611	c
	1612	CALL histdef(nid_mth, "rneb", "Cloud fraction", "%",
[776]	1613	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1614	. "ave(X)", zstophy,zout)
	1615	c
	1616	CALL histdef(nid_mth, "rnebcon", "Convective Cloud Fraction"
	1617	. , "-",
[776]	1618	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1619	. "ave(X)", zstophy,zout)
	1620	c
	1621	CALL histdef(nid_mth, "rhum", "Relative humidity", "-",
[776]	1622	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1623	. "ave(X)", zstophy,zout)
	1624	c
	1625	CALL histdef(nid_mth, "ozone", "Ozone concentration", "ppmv",
[776]	1626	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1627	. "ave(X)", zstophy,zout)
	1628	c
	1629	CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
[776]	1630	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1631	. "ave(X)", zstophy,zout)
	1632	c
	1633	CALL histdef(nid_mth, "dqphy","Physics dQ","(kg/kg)/s",
[776]	1634	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1635	. "ave(X)", zstophy,zout)
	1636	c
	1637	ENDIF !lev_histmth.GE.2
	1638	c
	1639	IF(lev_histmth.GE.3) THEN
	1640	c
	1641	DO nsrf=1, nbsrf
	1642	c
	1643	call histdef(nid_mth, "albe_"//clnsurf(nsrf),
	1644	$ "Albedo surf. "//clnsurf(nsrf), "-",
[776]	1645	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1646	$ "ave(X)", zstophy,zout)
	1647	c
	1648	call histdef(nid_mth, "rugs_"//clnsurf(nsrf),
	1649	$ "Rugosite "//clnsurf(nsrf), "m",
[776]	1650	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1651	$ "ave(X)", zstophy,zout)
	1652	c
	1653	CALL histdef(nid_mth, "ages_"//clnsurf(nsrf), "Snow age","day",
[776]	1654	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1655	. "ave(X)", zstophy,zout)
	1656	c
	1657	ENDDO !nsrf=1, nbsrf
	1658	c
	1659	CALL histdef(nid_mth, "albs", "Surface albedo", "-",
[776]	1660	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1661	. "ave(X)", zstophy,zout)
	1662	CALL histdef(nid_mth, "albslw", "Surface albedo LW", "-",
[776]	1663	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	1664	. "ave(X)", zstophy,zout)
	1665	c
	1666	ENDIF !lev_histmth.GE.3
	1667	c
	1668	IF(lev_histmth.GE.4) THEN
	1669	c
	1670	CALL histdef(nid_mth, "clwcon",
	1671	. "Convective Cloud Liquid water content"
	1672	. , "kg/kg",
[776]	1673	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1674	. "ave(X)", zstophy,zout)
	1675	c
	1676	CALL histdef(nid_mth,"Ma","undilute adiab updraft","kg/m2/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, "dnwd", "saturated downdraft","kg/m2/s",
[776]	1681	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1682	. "ave(X)", zstophy,zout)
	1683	c
	1684	CALL histdef(nid_mth, "dnwd0", "unsat. downdraft", "kg/m2/s",
[776]	1685	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1686	. "ave(X)", zstophy,zout)
	1687	c
	1688	CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
[776]	1689	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1690	. "ave(X)", zstophy,zout)
	1691	c
	1692	CALL histdef(nid_mth, "dqdyn","Dynamics dQ","(kg/kg)/s",
[776]	1693	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1694	. "ave(X)", zstophy,zout)
	1695	c
	1696	CALL histdef(nid_mth, "dqcon","Convection dQ","(kg/kg)/s",
[776]	1697	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1698	. "ave(X)", zstophy,zout)
	1699	c
	1700	CALL histdef(nid_mth, "dtlschr",
[766]	1701	$ "Large-scale condensational heating rate", "K/s",
[776]	1702	$ iim,jj_nb
[644]	1703	$ ,nhori, klev,1,klev,nvert, 32,"ave(X)", zstophy,zout)
	1704	c
	1705	CALL histdef(nid_mth, "dqlsc","Condensation dQ","(kg/kg)/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, "dtvdf", "Boundary-layer 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, "dqvdf","Boundary-layer dQ","(kg/kg)/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, "dteva", "Reevaporation dT", "K/s",
[776]	1718	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1719	. "ave(X)", zstophy,zout)
	1720	c
	1721	CALL histdef(nid_mth,"dqeva","Reevaporation dQ","(kg/kg)/s",
[776]	1722	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1723	. "ave(X)", zstophy,zout)
	1724
	1725	CALL histdef(nid_mth, "ptconv", "POINTS CONVECTIFS"," ",
[776]	1726	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1727	. "ave(X)", zstophy,zout)
	1728
	1729	CALL histdef(nid_mth, "ratqs", "RATQS"," ",
[776]	1730	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1731	. "ave(X)", zstophy,zout)
	1732
	1733	c
	1734	CALL histdef(nid_mth, "dqajs","Dry adjust. dQ","(kg/kg)/s",
[776]	1735	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1736	. "ave(X)", zstophy,zout)
	1737	c
	1738	CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
[776]	1739	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1740	. "ave(X)", zstophy,zout)
	1741	c
	1742	CALL histdef(nid_mth, "dtsw0", "CS SW radiation dT", "K/s",
[776]	1743	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1744	. "ave(X)", zstophy,zout)
	1745	c
	1746	CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
[776]	1747	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1748	. "ave(X)", zstophy,zout)
	1749	c
	1750	CALL histdef(nid_mth,"dtlw0","CS LW radiation dT","K/s",
[776]	1751	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1752	. "ave(X)", zstophy,zout)
	1753	c
	1754	CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
[776]	1755	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1756	. "ave(X)", zstophy,zout)
	1757	c
	1758	CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
[776]	1759	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1760	. "ave(X)", zstophy,zout)
	1761	c
	1762	CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
[776]	1763	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1764	. "ave(X)", zstophy,zout)
	1765	c
	1766	IF (ok_orodr) THEN
	1767	CALL histdef(nid_mth, "duoro", "Orography dU", "m/s2",
[776]	1768	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1769	. "ave(X)", zstophy,zout)
	1770	c
	1771	CALL histdef(nid_mth, "dvoro", "Orography dV", "m/s2",
[776]	1772	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1773	. "ave(X)", zstophy,zout)
	1774	c
	1775	ENDIF
	1776	C
	1777	IF (ok_orolf) THEN
	1778	CALL histdef(nid_mth, "dulif", "Orography dU", "m/s2",
[776]	1779	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1780	. "ave(X)", zstophy,zout)
	1781	c
	1782	CALL histdef(nid_mth, "dvlif", "Orography dV", "m/s2",
[776]	1783	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1784	. "ave(X)", zstophy,zout)
	1785	ENDIF
	1786	C
[862]	1787	#if defined(INCA_AER) && defined(CPP_COUPLE)
	1788	CALL histdef(nid_mth, "tau", "cloud opt thickness", "-",
	1789	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
	1790	. "ave(X)", zstophy,zout)
	1791	CALL histdef(nid_mth, "CCM1", "CCM1", "",
	1792	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
	1793	. "ave(X)", zstophy,zout)
	1794	CALL histdef(nid_mth, "CCM2", "CCM2", "",
	1795	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
	1796	. "ave(X)", zstophy,zout)
	1797	CALL histdef(nid_mth, "taupi", "cloud opt thickn. (pi)", "-",
	1798	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
	1799	. "ave(X)", zstophy,zout)
	1800	CALL histdef(nid_mth, "cldemi", "emissivite IR", "-",
	1801	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
	1802	. "ave(X)", zstophy,zout)
	1803	CALL histdef(nid_mth, "SWTOAAS_ZERO", "?", "?",
	1804	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	1805	. "ave(X)", zstophy,zout)
	1806	CALL histdef(nid_mth, "SWTOAAS_AER", "?", "?",
	1807	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	1808	. "ave(X)", zstophy,zout)
	1809	CALL histdef(nid_mth, "SWTOAAS_NAT", "?", "?",
	1810	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	1811	. "ave(X)", zstophy,zout)
	1812	CALL histdef(nid_mth, "SWTOAAS_AD", "?", "?",
	1813	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	1814	. "ave(X)", zstophy,zout)
	1815	CALL histdef(nid_mth, "SWTOAAS_AI", "?", "?",
	1816	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	1817	. "ave(X)", zstophy,zout)
	1818	CALL histdef(nid_mth, "SWTOACS_ZERO", "?", "?",
	1819	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	1820	. "ave(X)", zstophy,zout)
	1821	CALL histdef(nid_mth, "SWTOACS_AER", "?", "?",
	1822	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	1823	. "ave(X)", zstophy,zout)
	1824	CALL histdef(nid_mth, "SWTOACS_AD", "?", "?",
	1825	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	1826	. "ave(X)", zstophy,zout)
	1827	CALL histdef(nid_mth, "SWTOACS_NAT", "?", "?",
	1828	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	1829	. "ave(X)", zstophy,zout)
	1830	#endif
	1831	#ifndef INCA
[524]	1832	if (nqmax.GE.3) THEN
[845]	1833	DO iq=3,nqmax
	1834	iiq=niadv(iq)
	1835	CALL histdef(nid_mth, tnom(iq), ttext(iiq), "-",
[776]	1836	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	1837	. "ave(X)", zstophy,zout)
[845]	1838	ENDDO
	1839	ENDIF
[862]	1840	#endif
[524]	1841	ENDIF !lev_histmth.GE.4
	1842	c
[644]	1843	ENDIF !type_run.EQ."CLIM".OR.type_run.EQ."ENSP")
	1844	c
[524]	1845	CALL histend(nid_mth)
	1846	c
	1847	ndex2d = 0
	1848	ndex3d = 0
	1849	c
[766]	1850	c$OMP END MASTER
[524]	1851	ENDIF ! fin de test sur ok_mensuel

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

Download in other formats: