Context Navigation

source: LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histmth.h @ 456

Last change on this file since 456 was 455, checked in by lmdzadmin, 22 years ago
Ajout slp, ffonte, fqcalving IM/FH/JLD
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 22.7 KB

Rev	Line
[416]	1	IF (ok_mensuel) THEN
	2	c
	3	idayref = day_ref
	4	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
	5	c
	6	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
	7	DO i = 1, iim
	8	zx_lon(i,1) = rlon(i+1)
	9	zx_lon(i,jjmp1) = rlon(i+1)
	10	ENDDO
	11	DO ll=1,klev
	12	znivsig(ll)=float(ll)
	13	ENDDO
	14	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
	15	CALL histbeg("histmth.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
	16	. 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
	17	. nhori, nid_mth)
	18	write(,)'Mensuel ', itau_phy, zjulian
	19	CALL histvert(nid_mth, "presnivs", "Vertical levels", "mb",
	20	. klev, presnivs, nvert)
	21	c call histvert(nid_mth, 'sig_s', 'Niveaux sigma','-',
	22	c . klev, znivsig, nvert)
	23	c
	24	zsto = dtime
	25	zout = dtime * ecrit_mth
	26	c
	27	CALL histdef(nid_mth, "phis", "Surface geop. height", "-",
	28	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	29	. "once", zsto,zout)
	30	c
	31	CALL histdef(nid_mth, "aire", "Grid area", "-",
	32	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	33	. "once", zsto,zout)
	34	c
	35	c Champs 2D:
	36	c
[455]	37	CALL histdef(nid_mth, "slp", "Sea Level Pressure", "Pa",
	38	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	39	. "ave(X)", zsto,zout)
	40	c
[416]	41	CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
	42	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	43	. "ave(X)", zsto,zout)
	44	c
	45	CALL histdef(nid_mth, "t2m", "Temperature 2m", "K",
	46	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	47	. "ave(X)", zsto,zout)
	48	c
	49	CALL histdef(nid_mth, "q2m", "Specific humidity 2m", "Kg/Kg",
	50	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	51	. "ave(X)", zsto,zout)
	52	c
	53	CALL histdef(nid_mth, "u10m", "Vent zonal 10m", "m/s",
	54	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	55	. "ave(X)", zsto,zout)
	56	c
	57	CALL histdef(nid_mth, "v10m", "Vent meridien 10m", "m/s",
	58	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	59	. "ave(X)", zsto,zout)
	60	c
	61	c
	62	CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
	63	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	64	. "ave(X)", zsto,zout)
	65	c
[440]	66	CALL histdef(nid_mth, "qsurf", "Surface Air humidity", "Kg/Kg",
[416]	67	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	68	. "ave(X)", zsto,zout)
	69	c
[440]	70	if (.not. ok_veget) then
	71	CALL histdef(nid_mth, "qsol", "Soil watter content", "mm",
	72	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	73	. "ave(X)", zsto,zout)
	74	endif
	75	c
[416]	76	CALL histdef(nid_mth, "precip", "Precipitation Totale liq+sol",
[451]	77	. "kg/(s*m2)",
[416]	78	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	79	. "ave(X)", zsto,zout)
	80	c
[451]	81	CALL histdef(nid_mth, "plul", "Large-scale Precip.",
	82	. "kg/(s*m2)",
[416]	83	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	84	. "ave(X)", zsto,zout)
	85	c
[451]	86	CALL histdef(nid_mth, "pluc", "Convective Precip.",
	87	. "kg/(s*m2)",
[416]	88	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	89	. "ave(X)", zsto,zout)
	90	c
[451]	91	CALL histdef(nid_mth, "snow", "Snow fall", "kg/(s*m2)",
[416]	92	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	93	. "ave(X)", zsto,zout)
	94	c
	95	CALL histdef(nid_mth, "snow_mass", "Snow Mass", "kg/m2",
	96	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	97	. "ave(X)", zsto,zout)
	98	c
[451]	99	CALL histdef(nid_mth, "evap", "Evaporation", "kg/(s*m2)",
[416]	100	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	101	. "ave(X)", zsto,zout)
	102	c
	103	CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
	104	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	105	. "ave(X)", zsto,zout)
	106	c
	107	CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
	108	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	109	. "ave(X)", zsto,zout)
	110	c
	111	CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
	112	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	113	. "ave(X)", zsto,zout)
	114	c
	115	CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
	116	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	117	. "ave(X)", zsto,zout)
	118	c
	119	CALL histdef(nid_mth, "solldown", "Down. IR rad. at surface",
	120	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
	121	. "ave(X)", zsto,zout)
	122	c
	123	CALL histdef(nid_mth, "tops0", "Solar rad. at TOA", "W/m2",
	124	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	125	. "ave(X)", zsto,zout)
	126	c
	127	CALL histdef(nid_mth, "topl0", "IR rad. at TOA", "W/m2",
	128	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	129	. "ave(X)", zsto,zout)
	130	c
	131	CALL histdef(nid_mth, "sols0", "Solar rad. at surf.", "W/m2",
	132	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	133	. "ave(X)", zsto,zout)
	134	c
	135	CALL histdef(nid_mth, "soll0", "IR rad. at surface", "W/m2",
	136	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	137	. "ave(X)", zsto,zout)
	138	c
	139	CALL histdef(nid_mth, "bils", "Surf. total heat flux", "W/m2",
	140	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	141	. "ave(X)", zsto,zout)
	142	c
	143	CALL histdef(nid_mth, "sens", "Sensible heat flux", "W/m2",
	144	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	145	. "ave(X)", zsto,zout)
	146	c
	147	CALL histdef(nid_mth, "fder", "Heat flux derivation", "W/m2",
	148	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	149	. "ave(X)", zsto,zout)
	150	c
	151	CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
	152	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	153	. "ave(X)", zsto,zout)
	154	c
	155	CALL histdef(nid_mth, "frtv", "Meridional wind stress", "Pa",
	156	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	157	. "ave(X)", zsto,zout)
[455]	158	cIM cf JLD
	159	CALL histdef(nid_mth, "ffonte","Thermal flux for snow melting",
	160	. "W/m2",iim,jjmp1,nhori, 1,1,1, -99, 32,
	161	. "ave(X)", zsto,zout)
[416]	162	c
[455]	163	CALL histdef(nid_mth, "fqcalving","Ice Calving",
	164	. "kg/m2/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
	165	. "ave(X)", zsto,zout)
	166	c
[416]	167	DO nsrf = 1, nbsrf
	168	C
	169	call histdef(nid_mth, "pourc_"//clnsurf(nsrf),
	170	$ "Fraction "//clnsurf(nsrf), "W/m2",
	171	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	172	$ "ave(X)", zsto,zout)
	173	C
	174	call histdef(nid_mth, "tsol_"//clnsurf(nsrf),
	175	$ "Fraction "//clnsurf(nsrf), "W/m2",
	176	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	177	$ "ave(X)", zsto,zout)
	178	C
	179	call histdef(nid_mth, "sens_"//clnsurf(nsrf),
	180	$ "Sensible heat flux "//clnsurf(nsrf), "W/m2",
	181	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	182	$ "ave(X)", zsto,zout)
	183	c
	184	call histdef(nid_mth, "lat_"//clnsurf(nsrf),
	185	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
	186	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	187	$ "ave(X)", zsto,zout)
	188	C
[433]	189	call histdef(nid_mth, "flw_"//clnsurf(nsrf),
	190	$ "LW "//clnsurf(nsrf), "W/m2",
	191	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	192	$ "ave(X)", zsto,zout)
	193	cIM cf. JLD
	194	call histdef(nid_mth, "fsw_"//clnsurf(nsrf),
	195	$ "SW "//clnsurf(nsrf), "W/m2",
	196	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	197	$ "ave(X)", zsto,zout)
	198	C
	199	call histdef(nid_mth, "wbils_"//clnsurf(nsrf),
	200	$ "Bilan sol "//clnsurf(nsrf), "W/m2",
	201	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	202	$ "ave(X)", zsto,zout)
	203	C
[416]	204	call histdef(nid_mth, "taux_"//clnsurf(nsrf),
	205	$ "Zonal wind stress"//clnsurf(nsrf), "Pa",
	206	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	207	$ "ave(X)", zsto,zout)
	208
	209	call histdef(nid_mth, "tauy_"//clnsurf(nsrf),
	210	$ "Meridional xind stress "//clnsurf(nsrf), "Pa",
	211	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	212	$ "ave(X)", zsto,zout)
	213	c
	214	call histdef(nid_mth, "albe_"//clnsurf(nsrf),
[440]	215	$ "Albedo surf. "//clnsurf(nsrf), "-",
[416]	216	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	217	$ "ave(X)", zsto,zout)
	218	c
	219	call histdef(nid_mth, "rugs_"//clnsurf(nsrf),
	220	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
	221	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	222	$ "ave(X)", zsto,zout)
	223	c
	224	CALL histdef(nid_mth, "ages_"//clnsurf(nsrf), "Snow age","day",
	225	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	226	. "ave(X)", zsto,zout)
	227
	228	END DO
	229	C
	230	CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "-",
	231	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	232	. "ave(X)", zsto,zout)
	233	c
	234	CALL histdef(nid_mth, "albs", "Surface albedo", "-",
	235	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	236	. "ave(X)", zsto,zout)
	237	CALL histdef(nid_mth, "albslw", "Surface albedo LW", "-",
	238	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	239	. "ave(X)", zsto,zout)
	240	c
	241	CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
	242	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	243	. "ave(X)", zsto,zout)
	244	c
	245	CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
	246	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	247	. "ave(X)", zsto,zout)
	248	c
	249	CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-",
	250	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	251	. "ave(X)", zsto,zout)
	252	c
	253	CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-",
	254	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	255	. "ave(X)", zsto,zout)
	256	c
	257	CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-",
	258	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	259	. "ave(X)", zsto,zout)
	260	c
	261	CALL histdef(nid_mth, "cldt", "Total cloudiness", "-",
	262	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	263	. "ave(X)", zsto,zout)
	264	c
	265	CALL histdef(nid_mth, "cldq", "Cloud liquid water path", "-",
	266	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	267	. "ave(X)", zsto,zout)
	268	c
	269	CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
	270	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	271	. "ave(X)", zsto,zout)
	272	c
	273	CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
	274	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	275	. "ave(X)", zsto,zout)
	276	c
	277	CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-",
	278	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	279	. "ave(X)", zsto,zout)
	280	c
	281	CALL histdef(nid_mth, "vq", "Merid humidity transport", "-",
	282	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	283	. "ave(X)", zsto,zout)
	284	cKE43
	285	IF (iflag_con .GE. 3) THEN ! sb
	286	c
	287	CALL histdef(nid_mth, "cape", "Conv avlbl pot ener", "J/Kg",
	288	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	289	. "ave(X)", zsto,zout)
	290	c
	291	CALL histdef(nid_mth, "pbase", "Cld base pressure", "hPa",
	292	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	293	. "ave(X)", zsto,zout)
	294	c
	295	CALL histdef(nid_mth, "ptop", "Cld top pressure", "hPa",
	296	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	297	. "ave(X)", zsto,zout)
	298	c
	299	CALL histdef(nid_mth, "fbase", "Cld base mass flux", "Kg/m2/s",
	300	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	301	. "ave(X)", zsto,zout)
	302	c
	303	c
	304	ENDIF
	305	c34EK
	306	c
	307	c Champs 3D:
	308	c
	309	CALL histdef(nid_mth, "temp", "Air temperature", "K",
	310	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	311	. "ave(X)", zsto,zout)
	312	c
	313	CALL histdef(nid_mth, "ovap", "Specific humidity", "Kg/Kg",
	314	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	315	. "ave(X)", zsto,zout)
	316	c
	317	CALL histdef(nid_mth, "geop", "Geopotential height", "m",
	318	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	319	. "ave(X)", zsto,zout)
	320	c
	321	CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
	322	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	323	. "ave(X)", zsto,zout)
	324	c
	325	CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
	326	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	327	. "ave(X)", zsto,zout)
	328	c
	329	CALL histdef(nid_mth, "vitw", "Vertical wind", "m/s",
	330	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	331	. "ave(X)", zsto,zout)
	332	c
	333	CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
	334	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	335	. "ave(X)", zsto,zout)
	336	c
	337	CALL histdef(nid_mth, "rneb", "Cloud fraction", "-",
	338	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	339	. "ave(X)", zsto,zout)
[433]	340	cIM cf. FH
	341	CALL histdef(nid_mth, "rnebcon", "Convective Cloud Fraction"
	342	. , "-",
	343	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	344	. "ave(X)", zsto,zout)
[416]	345	c
	346	CALL histdef(nid_mth, "rhum", "Relative humidity", "-",
	347	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	348	. "ave(X)", zsto,zout)
	349	c
	350	CALL histdef(nid_mth, "clwcon", "Cloud Liquid water content"
	351	. , "kg/kg",
	352	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	353	. "ave(X)", zsto,zout)
	354	c
	355	CALL histdef(nid_mth, "oliq", "Liquid water content", "kg/kg",
	356	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	357	. "ave(X)", zsto,zout)
	358	c
	359	CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
	360	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	361	. "ave(X)", zsto,zout)
	362	c
	363	CALL histdef(nid_mth, "dqdyn", "Dynamics dQ", "Kg/Kg/s",
	364	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	365	. "ave(X)", zsto,zout)
	366	c
	367	CALL histdef(nid_mth, "dtcon", "Convection dT", "K/s",
	368	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	369	. "ave(X)", zsto,zout)
	370	c
	371	CALL histdef(nid_mth, "ducon", "Convection du", "m/s2",
	372	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	373	. "ave(X)", zsto,zout)
	374	c
	375	CALL histdef(nid_mth, "dqcon", "Convection dQ", "Kg/Kg/s",
	376	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	377	. "ave(X)", zsto,zout)
	378	c
	379	CALL histdef(nid_mth, "dtlsc", "Condensation dT", "K/s",
	380	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	381	. "ave(X)", zsto,zout)
	382	c
	383	CALL histdef(nid_mth, "dqlsc", "Condensation dQ", "Kg/Kg/s",
	384	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	385	. "ave(X)", zsto,zout)
	386	c
	387	CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
	388	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	389	. "ave(X)", zsto,zout)
	390	c
	391	CALL histdef(nid_mth, "dqvdf", "Boundary-layer dQ", "Kg/Kg/s",
	392	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	393	. "ave(X)", zsto,zout)
	394	c
	395	CALL histdef(nid_mth, "dteva", "Reevaporation dT", "K/s",
	396	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	397	. "ave(X)", zsto,zout)
	398	c
	399	CALL histdef(nid_mth, "dqeva", "Reevaporation dQ", "Kg/Kg/s",
	400	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	401	. "ave(X)", zsto,zout)
	402
	403	CALL histdef(nid_mth, "ptconv", "POINTS CONVECTIFS"," ",
	404	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	405	. "ave(X)", zsto,zout)
	406
	407	CALL histdef(nid_mth, "ratqs", "RATQS"," ",
	408	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	409	. "ave(X)", zsto,zout)
	410
	411	c
	412	CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s",
	413	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	414	. "ave(X)", zsto,zout)
	415
	416	CALL histdef(nid_mth, "dqajs", "Dry adjust. dQ", "Kg/Kg/s",
	417	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	418	. "ave(X)", zsto,zout)
	419	c
	420	CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
	421	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	422	. "ave(X)", zsto,zout)
	423	c
	424	CALL histdef(nid_mth, "dtsw0", "SW radiation dT", "K/s",
	425	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	426	. "ave(X)", zsto,zout)
	427	c
	428	CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
	429	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	430	. "ave(X)", zsto,zout)
	431	c
	432	CALL histdef(nid_mth, "dtlw0", "LW radiation dT", "K/s",
	433	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	434	. "ave(X)", zsto,zout)
	435	c
	436	CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
	437	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	438	. "ave(X)", zsto,zout)
	439	c
	440	CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
	441	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	442	. "ave(X)", zsto,zout)
	443	c
	444	CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
	445	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	446	. "ave(X)", zsto,zout)
	447	c
	448	IF (ok_orodr) THEN
	449	CALL histdef(nid_mth, "duoro", "Orography dU", "m/s2",
	450	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	451	. "ave(X)", zsto,zout)
	452	c
	453	CALL histdef(nid_mth, "dvoro", "Orography dV", "m/s2",
	454	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	455	. "ave(X)", zsto,zout)
	456	c
	457	ENDIF
	458	C
	459	IF (ok_orolf) THEN
	460	CALL histdef(nid_mth, "dulif", "Orography dU", "m/s2",
	461	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	462	. "ave(X)", zsto,zout)
	463	c
	464	CALL histdef(nid_mth, "dvlif", "Orography dV", "m/s2",
	465	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	466	. "ave(X)", zsto,zout)
	467	ENDIF
	468	C
	469	CALL histdef(nid_mth, "ozone", "Ozone concentration", "-",
	470	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	471	. "ave(X)", zsto,zout)
	472	c
	473	if (nqmax.GE.3) THEN
	474	DO iq=1,nqmax-2
	475	IF (iq.LE.99) THEN
	476	WRITE(str2,'(i2.2)') iq
	477	CALL histdef(nid_mth, "trac"//str2, "Tracer No."//str2, "-",
	478	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	479	. "ave(X)", zsto,zout)
	480	ELSE
	481	PRINT*, "Trop de traceurs"
	482	CALL abort
	483	ENDIF
	484	ENDDO
	485	ENDIF
	486	c
	487	cKE43
	488	IF (iflag_con.GE.3) THEN ! (sb)
	489	c
	490	CALL histdef(nid_mth, "upwd", "saturated updraft", "Kg/m2/s",
	491	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	492	. "ave(X)", zsto,zout)
	493	c
	494	CALL histdef(nid_mth, "dnwd", "saturated downdraft","Kg/m2/s",
	495	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	496	. "ave(X)", zsto,zout)
	497	c
	498	CALL histdef(nid_mth, "dnwd0", "unsat. downdraft", "Kg/m2/s",
	499	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	500	. "ave(X)", zsto,zout)
	501	c
	502	CALL histdef(nid_mth,"Ma","undilute adiab updraft","Kg/m2/s",
	503	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	504	. "ave(X)", zsto,zout)
	505	cccIM
	506	CALL histdef(nid_mth, "SWupTOA", "SWup at TOA","W/m2",
	507	. iim,jjmp1,nhori, 1,1,1,-99,
	508	. 32, "ave(X)", zsto,zout)
	509	c
	510	CALL histdef(nid_mth, "SWupSFC", "SWup at surface","W/m2",
	511	. iim,jjmp1,nhori, 1,1,1,-99,
	512	. 32, "ave(X)", zsto,zout)
	513	c
	514	CALL histdef(nid_mth, "SWdnTOA", "SWdn at TOA","W/m2",
	515	. iim,jjmp1,nhori, 1,1,1,-99,
	516	. 32, "ave(X)", zsto,zout)
	517	c
	518	CALL histdef(nid_mth, "SWdnSFC", "SWdn at surface","W/m2",
	519	. iim,jjmp1,nhori, 1,1,1,-99,
	520	. 32, "ave(X)", zsto,zout)
[433]	521	cccIM clear sky
	522	CALL histdef(nid_mth, "SWupTOAclr",
	523	. "SWup clear sky at TOA","W/m2",
	524	. iim,jjmp1,nhori, 1,1,1,-99,
	525	. 32, "ave(X)", zsto,zout)
[416]	526	c
[433]	527	CALL histdef(nid_mth, "SWupSFCclr",
	528	. "SWup clear sky at surface","W/m2",
	529	. iim,jjmp1,nhori, 1,1,1,-99,
	530	. 32, "ave(X)", zsto,zout)
	531	c
	532	CALL histdef(nid_mth, "SWdnTOAclr",
	533	. "SWdn clear sky at TOA","W/m2",
	534	. iim,jjmp1,nhori, 1,1,1,-99,
	535	. 32, "ave(X)", zsto,zout)
[440]	536	c
	537	CALL histdef(nid_mth, "SWdnSFCclr",
	538	. "SWdn clear sky at surface","W/m2",
	539	. iim,jjmp1,nhori, 1,1,1,-99,
	540	. 32, "ave(X)", zsto,zout)
[433]	541	c Champs interpolles sur des niveaux de pression
	542
	543	CALL histdef(nid_mth, "u850", "Zonal wind 850mb", "m/s",
	544	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	545	. "ave(X)", zsto,zout)
	546
	547	CALL histdef(nid_mth, "v850", "Meridional wind 850mb", "m/s",
	548	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	549	. "ave(X)", zsto,zout)
	550	c
	551	CALL histdef(nid_mth, "u500", "Zonal wind 500mb", "m/s",
	552	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	553	. "ave(X)", zsto,zout)
	554
	555	CALL histdef(nid_mth, "v500", "Meridional wind 500mb", "m/s",
	556	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	557	. "ave(X)", zsto,zout)
	558
	559	CALL histdef(nid_mth, "u200", "Zonal wind 200mb", "m/s",
	560	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	561	. "ave(X)", zsto,zout)
	562
	563	CALL histdef(nid_mth, "v200", "Meridional wind 200mb", "m/s",
	564	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	565	. "ave(X)", zsto,zout)
	566
	567	CALL histdef(nid_mth, "phi500","Geopotentiel à 500mb","m2/s2",
	568	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	569	. "ave(X)", zsto,zout)
	570
	571	CALL histdef(nid_mth, "w500", "Vertical wind 500mb", "Pa/s",
	572	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	573	. "ave(X)", zsto,zout)
	574
	575
[416]	576	ENDIF
	577	c34EK
	578	CALL histend(nid_mth)
	579	c
	580	ndex2d = 0
	581	ndex3d = 0
	582	c
	583	ENDIF ! fin de test sur ok_mensuel

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

Download in other formats: