Context Navigation

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

Last change on this file since 525 was 524, checked in by lmdzadmin, 21 years ago
Initial revision
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 29.1 KB

Rev	Line
[524]	1	!
	2	! $Header$
	3	!
	4	IF (ok_mensuel) THEN
	5	c
	6	zsto = dtime
	7	zout = dtime * ecrit_mth
	8	c
	9	c zsto1: pour des valeurs "instantannees" mensuelles
	10	zsto1 = dtime * ecrit_mth
	11	c zsto2: pour des flux radiatifs calcules tous les 2 heures
	12	zsto2 = dtime * radpas
	13	PRINT*,' zsto,zsto1,zsto2,zout=',zsto, zsto1, zsto2,zout
	14	c
	15	idayref = day_ref
	16	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
	17	c
	18	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
	19	DO i = 1, iim
	20	zx_lon(i,1) = rlon(i+1)
	21	zx_lon(i,jjmp1) = rlon(i+1)
	22	ENDDO
	23	DO ll=1,klev
	24	znivsig(ll)=float(ll)
	25	ENDDO
	26	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
	27	CALL histbeg("histmth.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
	28	. 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
	29	. nhori, nid_mth)
	30	write(,)'Mensuel ', itau_phy, zjulian
	31	CALL histvert(nid_mth, "presnivs", "Vertical levels", "mb",
	32	. klev, presnivs, nvert)
	33	c call histvert(nid_mth, 'sig_s', 'Niveaux sigma','-',
	34	c . klev, znivsig, nvert)
	35	c
	36	c
	37	IF(lev_histmth.GE.1) THEN
	38	CALL histdef(nid_mth, "phis", "Surface geop. height", "-",
	39	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	40	. "once", zsto,zout)
	41	c
	42	CALL histdef(nid_mth, "aire", "Grid area", "-",
	43	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	44	. "once", zsto,zout)
	45	c
	46	call histdef(nid_mth, "pourc_"//clnsurf(is_ter),
	47	$ "Fraction "//clnsurf(is_ter), "W/m2",
	48	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	49	$ "once", zsto,zout)
	50	c
	51	c Champs 2D:
	52	c
	53	CALL histdef(nid_mth, "slp", "Sea Level Pressure", "Pa",
	54	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	55	. "ave(X)", zsto,zout)
	56	c
	57	CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
	58	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	59	. "ave(X)", zsto,zout)
	60	c
	61	CALL histdef(nid_mth, "t2m", "Temperature 2m", "K",
	62	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	63	. "ave(X)", zsto,zout)
	64	c ENSEMBLES BEG
	65	CALL histdef(nid_mth, "t2m_min", "Temp. 2m min.",
	66	. "K",
	67	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	68	. t2mincels, zsto,zout)
	69	c
	70	CALL histdef(nid_mth, "t2m_max", "Temp. 2m max.",
	71	. "K",
	72	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	73	. t2maxcels, zsto,zout)
	74	c
	75	c CALL histdef(nid_mth, "tsoil", "Sfce soil Temperature",
	76	c . "K",
	77	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	78	c . "ave(X)", zsto,zout)
	79	c
	80	CALL histdef(nid_mth, "wind10m","10-m wind speed","m/s",
	81	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	82	. "ave(X)", zsto,zout)
	83	c
	84	CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "-",
	85	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	86	. "ave(X)", zsto,zout)
	87
	88	c
	89	c ENSEMBLES END
	90	CALL histdef(nid_mth, "q2m", "Specific humidity 2m", "kg/kg",
	91	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	92	. "ave(X)", zsto,zout)
	93	c
	94	CALL histdef(nid_mth, "u10m", "Vent zonal 10m", "m/s",
	95	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	96	. "ave(X)", zsto,zout)
	97	c
	98	CALL histdef(nid_mth, "v10m", "Vent meridien 10m", "m/s",
	99	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	100	. "ave(X)", zsto,zout)
	101	c
	102	CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
	103	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	104	. "ave(X)", zsto,zout)
	105	c
	106	CALL histdef(nid_mth, "qsurf", "Surface Air humidity", "kg/kg",
	107	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	108	. "ave(X)", zsto,zout)
	109	c
	110	if (.not. ok_veget) then
	111	CALL histdef(nid_mth, "qsol", "Soil watter content", "mm",
	112	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	113	. "ave(X)", zsto,zout)
	114	endif
	115	c
	116	CALL histdef(nid_mth, "ndayrain",
	117	. "Number of day with rain (liq+sol)", "-",
	118	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	119	. "inst(X)", zsto1,zout)
	120	c
	121	CALL histdef(nid_mth, "precip", "Precipitation Totale liq+sol",
	122	. "kg/(s*m2)",
	123	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	124	. "ave(X)", zsto,zout)
	125	c
	126	CALL histdef(nid_mth, "plul", "Large-scale Precip.",
	127	. "kg/(s*m2)",
	128	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	129	. "ave(X)", zsto,zout)
	130	c
	131	CALL histdef(nid_mth, "pluc", "Convective Precip.",
	132	. "kg/(s*m2)",
	133	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	134	. "ave(X)", zsto,zout)
	135	c
	136	CALL histdef(nid_mth, "snow", "Snow fall", "kg/(s*m2)",
	137	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	138	. "ave(X)", zsto,zout)
	139	c
	140	CALL histdef(nid_mth, "snow_mass", "Snow Mass", "kg/m2",
	141	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	142	. "ave(X)", zsto,zout)
	143	c
	144	CALL histdef(nid_mth, "evap", "Evaporation", "kg/(s*m2)",
	145	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	146	. "ave(X)", zsto,zout)
	147	c
	148	CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
	149	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	150	. "ave(X)", zsto,zout)
	151	c
	152	CALL histdef(nid_mth, "tops0", "Solar rad. at TOA", "W/m2",
	153	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	154	. "ave(X)", zsto,zout)
	155	c
	156	CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
	157	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	158	. "ave(X)", zsto,zout)
	159	c
	160	CALL histdef(nid_mth, "topl0", "IR rad. at TOA", "W/m2",
	161	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	162	. "ave(X)", zsto,zout)
	163	c
	164	CALL histdef(nid_mth, "SWupTOA", "SWup at TOA","W/m2",
	165	. iim,jjmp1,nhori, 1,1,1,-99,
	166	. 32, "ave(X)", zsto2,zout)
	167	c
	168	CALL histdef(nid_mth, "SWupTOAclr",
	169	. "SWup clear sky at TOA","W/m2",
	170	. iim,jjmp1,nhori, 1,1,1,-99,
	171	. 32, "ave(X)", zsto2,zout)
	172	c
	173	CALL histdef(nid_mth, "SWdnTOA", "SWdn at TOA","W/m2",
	174	. iim,jjmp1,nhori, 1,1,1,-99,
	175	. 32, "ave(X)", zsto2,zout)
	176	c
	177	CALL histdef(nid_mth, "SWdnTOAclr",
	178	. "SWdn clear sky at TOA","W/m2",
	179	. iim,jjmp1,nhori, 1,1,1,-99,
	180	. 32, "ave(X)", zsto2,zout)
	181	c
	182	CALL histdef(nid_mth, "SWup200", "SWup at 200hPa","W/m2",
	183	. iim,jjmp1,nhori, 1,1,1,-99,
	184	. 32, "ave(X)", zsto2,zout)
	185	c
	186	CALL histdef(nid_mth, "SWup200clr",
	187	. "SWup clear sky at 200hPa","W/m2",
	188	. iim,jjmp1,nhori, 1,1,1,-99,
	189	. 32, "ave(X)", zsto2,zout)
	190	c
	191	CALL histdef(nid_mth, "SWdn200", "SWdn at 200hPa","W/m2",
	192	. iim,jjmp1,nhori, 1,1,1,-99,
	193	. 32, "ave(X)", zsto2,zout)
	194	c
	195	CALL histdef(nid_mth, "SWdn200clr",
	196	. "SWdn clear sky at 200hPa","W/m2",
	197	. iim,jjmp1,nhori, 1,1,1,-99,
	198	. 32, "ave(X)", zsto2,zout)
	199	c
	200	CALL histdef(nid_mth, "LWup200", "LWup at 200hPa","W/m2",
	201	. iim,jjmp1,nhori, 1,1,1,-99,
	202	. 32, "ave(X)", zsto2,zout)
	203	c
	204	CALL histdef(nid_mth, "LWup200clr",
	205	. "LWup clear sky at 200hPa","W/m2",
	206	. iim,jjmp1,nhori, 1,1,1,-99,
	207	. 32, "ave(X)", zsto2,zout)
	208	c
	209	CALL histdef(nid_mth, "LWdn200", "LWdn at 200hPa","W/m2",
	210	. iim,jjmp1,nhori, 1,1,1,-99,
	211	. 32, "ave(X)", zsto2,zout)
	212	c
	213	CALL histdef(nid_mth, "LWdn200clr",
	214	. "LWdn clear sky at 200hPa","W/m2",
	215	. iim,jjmp1,nhori, 1,1,1,-99,
	216	. 32, "ave(X)", zsto2,zout)
	217	c
	218	CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
	219	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	220	. "ave(X)", zsto,zout)
	221	c
	222	CALL histdef(nid_mth, "sols0", "Solar rad. at surf.", "W/m2",
	223	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	224	. "ave(X)", zsto,zout)
	225	c
	226	CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
	227	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	228	. "ave(X)", zsto,zout)
	229	c
	230	CALL histdef(nid_mth, "soll0", "IR rad. at surface", "W/m2",
	231	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	232	. "ave(X)", zsto,zout)
	233	c
	234	CALL histdef(nid_mth, "SWupSFC", "SWup at surface","W/m2",
	235	. iim,jjmp1,nhori, 1,1,1,-99,
	236	. 32, "ave(X)", zsto2,zout)
	237	c
	238	CALL histdef(nid_mth, "SWupSFCclr",
	239	. "SWup clear sky at surface","W/m2",
	240	. iim,jjmp1,nhori, 1,1,1,-99,
	241	. 32, "ave(X)", zsto2,zout)
	242	c
	243	CALL histdef(nid_mth, "SWdnSFC", "SWdn at surface","W/m2",
	244	. iim,jjmp1,nhori, 1,1,1,-99,
	245	. 32, "ave(X)", zsto2,zout)
	246	c
	247	CALL histdef(nid_mth, "SWdnSFCclr",
	248	. "SWdn clear sky at surface","W/m2",
	249	. iim,jjmp1,nhori, 1,1,1,-99,
	250	. 32, "ave(X)", zsto2,zout)
	251	c
	252	CALL histdef(nid_mth,"LWupSFC","Upwd. IR rad. at surface",
	253	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
	254	. "ave(X)", zsto2,zout)
	255	c
	256	CALL histdef(nid_mth,"LWdnSFC","Down. IR rad. at surface",
	257	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
	258	. "ave(X)", zsto2,zout)
	259	c
	260	CALL histdef(nid_mth,"LWupSFCclr",
	261	. "CS Upwd. IR rad. at surface",
	262	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
	263	. "ave(X)", zsto2,zout)
	264	c
	265	CALL histdef(nid_mth,"LWdnSFCclr",
	266	. "Down. CS IR rad. at surface",
	267	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
	268	. "ave(X)", zsto2,zout)
	269	c
	270	CALL histdef(nid_mth, "bils", "Surf. total heat flux", "W/m2",
	271	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	272	. "ave(X)", zsto,zout)
	273	c
	274	CALL histdef(nid_mth, "sens", "Sensible heat flux", "W/m2",
	275	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	276	. "ave(X)", zsto,zout)
	277	c
	278	CALL histdef(nid_mth, "fder", "Heat flux derivation", "W/m2",
	279	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	280	. "ave(X)", zsto,zout)
	281	c
	282	c Effets des aerosols
	283	c
	284	c IF (ok_ade.OR.ok_aie) THEN
	285	CALL histdef(nid_mth, "topsad", "ADE at TOA", "W/m2",
	286	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	287	. "ave(X)", zsto,zout)
	288	c
	289	CALL histdef(nid_mth, "solsad", "ADE at sfc", "W/m2",
	290	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	291	. "ave(X)", zsto,zout)
	292	c
	293	CALL histdef(nid_mth, "topsai", "AIE at TOA", "W/m2",
	294	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	295	. "ave(X)", zsto,zout)
	296	c
	297	CALL histdef(nid_mth, "solsai", "AIE at sfc", "W/m2",
	298	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	299	. "ave(X)", zsto,zout)
	300	c endif
	301	c
	302
	303	c
	304	c CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
	305	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	306	c . "ave(X)", zsto,zout)
	307	c
	308	c CALL histdef(nid_mth, "frtv", "Meridional wind stress", "Pa",
	309	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	310	c . "ave(X)", zsto,zout)
	311	c
	312	CALL histdef(nid_mth, "ffonte","Thermal flux for snow melting",
	313	. "W/m2",iim,jjmp1,nhori, 1,1,1, -99, 32,
	314	. "ave(X)", zsto,zout)
	315	c
	316	CALL histdef(nid_mth, "fqcalving","Ice Calving",
	317	. "kg/m2/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
	318	. "ave(X)", zsto,zout)
	319	cIM: 171003
	320	DO nsrf = 1, nbsrf
	321	call histdef(nid_mth, "taux_"//clnsurf(nsrf),
	322	$ "Zonal wind stress"//clnsurf(nsrf), "Pa",
	323	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	324	$ "ave(X)", zsto,zout)
	325
	326	call histdef(nid_mth, "tauy_"//clnsurf(nsrf),
	327	$ "Meridional xind stress "//clnsurf(nsrf), "Pa",
	328	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	329	$ "ave(X)", zsto,zout)
	330	ENDDO
	331	cIM: 171003
	332	c
	333	DO nsrf = 1, nbsrf
	334	C
	335	IF(nsrf.GT.1) THEN
	336	call histdef(nid_mth, "pourc_"//clnsurf(nsrf),
	337	$ "Fraction "//clnsurf(nsrf), "W/m2",
	338	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	339	$ "ave(X)", zsto,zout)
	340	ENDIF !nsrf.GT.1
	341	C
	342	call histdef(nid_mth, "tsol_"//clnsurf(nsrf),
	343	$ "Fraction "//clnsurf(nsrf), "W/m2",
	344	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	345	$ "ave(X)", zsto,zout)
	346	C
	347	call histdef(nid_mth, "sens_"//clnsurf(nsrf),
	348	$ "Sensible heat flux "//clnsurf(nsrf), "W/m2",
	349	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	350	$ "ave(X)", zsto,zout)
	351	c
	352	call histdef(nid_mth, "lat_"//clnsurf(nsrf),
	353	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
	354	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	355	$ "ave(X)", zsto,zout)
	356	C
	357	call histdef(nid_mth, "flw_"//clnsurf(nsrf),
	358	$ "LW "//clnsurf(nsrf), "W/m2",
	359	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	360	$ "ave(X)", zsto2,zout)
	361	c
	362	call histdef(nid_mth, "fsw_"//clnsurf(nsrf),
	363	$ "SW "//clnsurf(nsrf), "W/m2",
	364	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	365	$ "ave(X)", zsto2,zout)
	366	C
	367	call histdef(nid_mth, "wbils_"//clnsurf(nsrf),
	368	$ "Bilan sol "//clnsurf(nsrf), "W/m2",
	369	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	370	$ "ave(X)", zsto,zout)
	371	C
	372	END DO
	373	c
	374	CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
	375	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	376	. "ave(X)", zsto,zout)
	377	c
	378	CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
	379	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	380	. "ave(X)", zsto,zout)
	381	c
	382	CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-",
	383	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	384	. "ave(X)", zsto,zout)
	385	c
	386	CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-",
	387	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	388	. "ave(X)", zsto,zout)
	389	c
	390	CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-",
	391	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	392	. "ave(X)", zsto,zout)
	393	c
	394	CALL histdef(nid_mth, "cldt", "Total cloudiness", "-",
	395	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	396	. "ave(X)", zsto,zout)
	397	c
	398	CALL histdef(nid_mth,"cldq","Cloud liquid water path","kg/m2",
	399	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	400	. "ave(X)", zsto,zout)
	401	cIM: 071003
	402	CALL histdef(nid_mth,"lwp","Cloud water path","kg/m2",
	403	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	404	. "ave(X)", zsto,zout)
	405	c
	406	CALL histdef(nid_mth,"iwp","Cloud ice water path","kg/m2",
	407	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	408	. "ave(X)", zsto,zout)
	409	cIM: 071003
	410	c
	411	CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
	412	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	413	. "ave(X)", zsto,zout)
	414	c
	415	CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
	416	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	417	. "ave(X)", zsto,zout)
	418	c
	419	CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-",
	420	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	421	. "ave(X)", zsto,zout)
	422	c
	423	CALL histdef(nid_mth, "vq", "Merid humidity transport", "-",
	424	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	425	. "ave(X)", zsto,zout)
	426	cKE43
	427	IF(iflag_con.GE.3) THEN ! sb
	428	c
	429	CALL histdef(nid_mth, "cape", "Conv avlbl pot ener", "J/kg",
	430	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	431	. "ave(X)", zsto,zout)
	432	c
	433	CALL histdef(nid_mth, "pbase", "Cld base pressure", "hPa",
	434	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	435	. "ave(X)", zsto,zout)
	436	c
	437	CALL histdef(nid_mth, "ptop", "Cld top pressure", "hPa",
	438	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	439	. "ave(X)", zsto,zout)
	440	c
	441	CALL histdef(nid_mth, "fbase", "Cld base mass flux", "kg/m2/s",
	442	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	443	. "ave(X)", zsto,zout)
	444	c
	445	CALL histdef(nid_mth, "prw", "Precipitable water", "kg/m2",
	446	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	447	. "ave(X)", zsto,zout)
	448	c
	449	ENDIF !iflag_con .GE. 3
	450	c34EK
	451	c
	452	c Champs interpolles sur des niveaux de pression
	453	DO k=1, nlevENS
	454	CALL histdef(nid_mth, "u"//clev(k),
	455	. "Zonal wind"//clev(k)//"mb","m/s",
	456	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	457	. "ave(X)", zsto,zout)
	458	c
	459	CALL histdef(nid_mth, "v"//clev(k),
	460	. "Meridional wind"//clev(k)//"mb","m/s",
	461	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	462	. "ave(X)", zsto,zout)
	463	c
	464	CALL histdef(nid_mth, "w"//clev(k),
	465	. "Vertical wind"//clev(k)//"mb","m/s",
	466	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	467	. "ave(X)", zsto,zout)
	468	c
	469	CALL histdef(nid_mth, "phi"//clev(k),
	470	. "Geopotential"//clev(k)//"mb","m2/s2",
	471	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	472	. "ave(X)", zsto,zout)
	473	c
	474	ENDDO
	475	ENDIF !lev_histmth.GE.1
	476	c
	477	IF(lev_histmth.GE.2) THEN
	478	c
	479	c Champs 3D:
	480	c
	481	c
	482	CALL histdef(nid_mth,"lwcon","Cloud water content","kg/kg",
	483	. iim,jjmp1,nhori, klev,1,klev, nvert, 32,
	484	. "ave(X)", zsto,zout)
	485	c
	486	CALL histdef(nid_mth,"iwcon","Cloud ice water content","kg/kg",
	487	. iim,jjmp1,nhori, klev,1,klev, nvert, 32,
	488	. "ave(X)", zsto,zout)
	489	c
	490	CALL histdef(nid_mth, "temp", "Air temperature", "K",
	491	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	492	. "ave(X)", zsto,zout)
	493	c
	494	CALL histdef(nid_mth, "ovap", "Specific humidity", "kg/kg",
	495	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	496	. "ave(X)", zsto,zout)
	497	c
	498	c CALL histdef(nid_mth,"wvap","Water vapor mixing ratio","kg/kg",
	499	c . iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	500	c . "ave(X)", zsto,zout)
	501	c
	502	CALL histdef(nid_mth, "geop", "Geopotential height", "m",
	503	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	504	. "ave(X)", zsto,zout)
	505	c
	506	CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
	507	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	508	. "ave(X)", zsto,zout)
	509	c
	510	CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
	511	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	512	. "ave(X)", zsto,zout)
	513	c
	514	CALL histdef(nid_mth, "vitw", "Vertical wind", "m/s",
	515	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	516	. "ave(X)", zsto,zout)
	517	c
	518	CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
	519	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	520	. "ave(X)", zsto,zout)
	521	c
	522	CALL histdef(nid_mth, "rneb", "Cloud fraction", "-",
	523	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	524	. "ave(X)", zsto,zout)
	525	c
	526	CALL histdef(nid_mth, "rnebcon", "Convective Cloud Fraction"
	527	. , "-",
	528	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	529	. "ave(X)", zsto,zout)
	530	c
	531	CALL histdef(nid_mth, "rhum", "Relative humidity", "-",
	532	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	533	. "ave(X)", zsto,zout)
	534	c
	535	CALL histdef(nid_mth, "ozone", "Ozone concentration", "-",
	536	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	537	. "ave(X)", zsto,zout)
	538	c
	539	CALL histdef(nid_mth, "upwd", "saturated updraft", "kg/m2/s",
	540	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	541	. "ave(X)", zsto,zout)
	542	c
	543	CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
	544	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	545	. "ave(X)", zsto,zout)
	546	c
	547	CALL histdef(nid_mth, "dqphy", "Physics dQ", "kg/kg/s",
	548	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	549	. "ave(X)", zsto,zout)
	550	c
	551	c#define histmthNMC
	552	c#ifdef histmthNMC
	553	c#include "ini_histmthNMC.h"
	554	c#endif
	555	c
	556	ENDIF !lev_histmth.GE.2
	557	c
	558	IF(lev_histmth.GE.3) THEN
	559	c
	560	DO nsrf=1, nbsrf
	561	c
	562	call histdef(nid_mth, "albe_"//clnsurf(nsrf),
	563	$ "Albedo surf. "//clnsurf(nsrf), "-",
	564	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	565	$ "ave(X)", zsto,zout)
	566	c
	567	call histdef(nid_mth, "rugs_"//clnsurf(nsrf),
	568	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
	569	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	570	$ "ave(X)", zsto,zout)
	571	c
	572	CALL histdef(nid_mth, "ages_"//clnsurf(nsrf), "Snow age","day",
	573	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	574	. "ave(X)", zsto,zout)
	575	c
	576	ENDDO !nsrf=1, nbsrf
	577	c
	578	CALL histdef(nid_mth, "albs", "Surface albedo", "-",
	579	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	580	. "ave(X)", zsto,zout)
	581	CALL histdef(nid_mth, "albslw", "Surface albedo LW", "-",
	582	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	583	. "ave(X)", zsto,zout)
	584	c
	585	ENDIF !lev_histmth.GE.3
	586	c
	587	IF(lev_histmth.GE.4) THEN
	588	c
	589	CALL histdef(nid_mth, "clwcon",
	590	. "Convective Cloud Liquid water content"
	591	. , "kg/kg",
	592	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	593	. "ave(X)", zsto,zout)
	594	c
	595	CALL histdef(nid_mth,"Ma","undilute adiab updraft","kg/m2/s",
	596	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	597	. "ave(X)", zsto,zout)
	598	c
	599	CALL histdef(nid_mth, "dnwd", "saturated downdraft","kg/m2/s",
	600	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	601	. "ave(X)", zsto,zout)
	602	c
	603	CALL histdef(nid_mth, "dnwd0", "unsat. downdraft", "kg/m2/s",
	604	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	605	. "ave(X)", zsto,zout)
	606	c
	607	CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
	608	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	609	. "ave(X)", zsto,zout)
	610	c
	611	CALL histdef(nid_mth, "dqdyn", "Dynamics dQ", "kg/kg/s",
	612	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	613	. "ave(X)", zsto,zout)
	614	c
	615	CALL histdef(nid_mth, "dtcon", "Convection dT", "K/s",
	616	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	617	. "ave(X)", zsto,zout)
	618	c
	619	c CALL histdef(nid_mth, "ducon", "Convection du", "m/s2",
	620	c . iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	621	c . "ave(X)", zsto,zout)
	622	c
	623	CALL histdef(nid_mth, "dqcon", "Convection dQ", "kg/kg/s",
	624	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	625	. "ave(X)", zsto,zout)
	626	c
	627	CALL histdef(nid_mth, "dtlsc", "Condensation dT", "K/s",
	628	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	629	. "ave(X)", zsto,zout)
	630	cIM: 071003
	631	CALL histdef(nid_mth, "dtlschr",
	632	$ "Large-scale condensational heating rate", "K/s",iim,jjmp1
	633	$ ,nhori, klev,1,klev,nvert, 32,"ave(X)", zsto,zout)
	634	c
	635	CALL histdef(nid_mth, "dqlsc", "Condensation dQ", "kg/kg/s",
	636	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	637	. "ave(X)", zsto,zout)
	638	c
	639	CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
	640	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	641	. "ave(X)", zsto,zout)
	642	c
	643	CALL histdef(nid_mth, "dqvdf", "Boundary-layer dQ", "kg/kg/s",
	644	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	645	. "ave(X)", zsto,zout)
	646	c
	647	CALL histdef(nid_mth, "dteva", "Reevaporation dT", "K/s",
	648	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	649	. "ave(X)", zsto,zout)
	650	c
	651	CALL histdef(nid_mth, "dqeva", "Reevaporation dQ", "kg/kg/s",
	652	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	653	. "ave(X)", zsto,zout)
	654
	655	CALL histdef(nid_mth, "ptconv", "POINTS CONVECTIFS"," ",
	656	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	657	. "ave(X)", zsto,zout)
	658
	659	CALL histdef(nid_mth, "ratqs", "RATQS"," ",
	660	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	661	. "ave(X)", zsto,zout)
	662
	663	c
	664	CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s",
	665	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	666	. "ave(X)", zsto,zout)
	667
	668	CALL histdef(nid_mth, "dqajs", "Dry adjust. dQ", "kg/kg/s",
	669	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	670	. "ave(X)", zsto,zout)
	671	c
	672	CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
	673	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	674	. "ave(X)", zsto,zout)
	675	c
	676	CALL histdef(nid_mth, "dtsw0", "CS SW radiation dT", "K/s",
	677	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	678	. "ave(X)", zsto,zout)
	679	c
	680	CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
	681	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	682	. "ave(X)", zsto,zout)
	683	c
	684	CALL histdef(nid_mth,"dtlw0","CS LW radiation dT","K/s",
	685	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	686	. "ave(X)", zsto,zout)
	687	c
	688	CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
	689	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	690	. "ave(X)", zsto,zout)
	691	c
	692	CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
	693	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	694	. "ave(X)", zsto,zout)
	695	c
	696	CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
	697	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	698	. "ave(X)", zsto,zout)
	699	c
	700	IF (ok_orodr) THEN
	701	CALL histdef(nid_mth, "duoro", "Orography dU", "m/s2",
	702	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	703	. "ave(X)", zsto,zout)
	704	c
	705	CALL histdef(nid_mth, "dvoro", "Orography dV", "m/s2",
	706	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	707	. "ave(X)", zsto,zout)
	708	c
	709	ENDIF
	710	C
	711	IF (ok_orolf) THEN
	712	CALL histdef(nid_mth, "dulif", "Orography dU", "m/s2",
	713	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	714	. "ave(X)", zsto,zout)
	715	c
	716	CALL histdef(nid_mth, "dvlif", "Orography dV", "m/s2",
	717	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	718	. "ave(X)", zsto,zout)
	719	ENDIF
	720	c
	721	c Effets des aerosols
	722	c
	723	c IF (ok_ade.OR.ok_aie) THEN
	724	CALL histdef(nid_mth, "re", "CDR", "um",
	725	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	726	. "ave(X)", zsto,zout)
	727	c
	728	CALL histdef(nid_mth, "redenom", "CDR denominator", "-",
	729	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	730	. "ave(X)", zsto,zout)
	731	c
	732	CALL histdef(nid_mth, "tau", "cloud opt thickness", "-",
	733	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	734	. "ave(X)", zsto,zout)
	735	c
	736	CALL histdef(nid_mth, "taupi", "cloud opt thickn. (pi)", "-",
	737	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	738	. "ave(X)", zsto,zout)
	739	c endif
	740	c
	741	CALL histdef(nid_mth, "ozone", "Ozone concentration", "-",
	742	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	743	. "ave(X)", zsto,zout)
	744	c
	745	if (nqmax.GE.3) THEN
	746	DO iq=1,nqmax-2
	747	IF (iq.LE.99) THEN
	748	WRITE(str2,'(i2.2)') iq
	749	CALL histdef(nid_mth, "trac"//str2, "Tracer No."//str2, "-",
	750	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	751	. "ave(X)", zsto,zout)
	752	ELSE
	753	PRINT*, "Trop de traceurs"
	754	CALL abort
	755	ENDIF
	756	ENDDO
	757	ENDIF
	758	c
	759	ENDIF !lev_histmth.GE.4
	760	c
	761	CALL histend(nid_mth)
	762	c
	763	ndex2d = 0
	764	ndex3d = 0
	765	c
	766	ENDIF ! fin de test sur ok_mensuel

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

Download in other formats: