Context Navigation

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

Last change on this file since 505 was 486, checked in by lmdzadmin, 22 years ago
Phasage avec la version de Ionela IM/LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 24.2 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)
[486]	118	cIM: 071003
	119	CALL histdef(nid_mth,"LWdnSFC","Down. IR rad. at surface",
[416]	120	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
	121	. "ave(X)", zsto,zout)
[486]	122	cIM: 071003
	123	CALL histdef(nid_mth,"LWupSFC","Upwd. IR rad. at surface",
	124	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
	125	. "ave(X)", zsto,zout)
[416]	126	c
	127	CALL histdef(nid_mth, "tops0", "Solar 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, "topl0", "IR rad. at TOA", "W/m2",
	132	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	133	. "ave(X)", zsto,zout)
	134	c
	135	CALL histdef(nid_mth, "sols0", "Solar rad. at surf.", "W/m2",
	136	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	137	. "ave(X)", zsto,zout)
	138	c
	139	CALL histdef(nid_mth, "soll0", "IR rad. at surface", "W/m2",
	140	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	141	. "ave(X)", zsto,zout)
	142	c
	143	CALL histdef(nid_mth, "bils", "Surf. total 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, "sens", "Sensible heat flux", "W/m2",
	148	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	149	. "ave(X)", zsto,zout)
	150	c
	151	CALL histdef(nid_mth, "fder", "Heat flux derivation", "W/m2",
	152	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	153	. "ave(X)", zsto,zout)
	154	c
[470]	155	c CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
	156	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	157	c . "ave(X)", zsto,zout)
[416]	158	c
[470]	159	c CALL histdef(nid_mth, "frtv", "Meridional wind stress", "Pa",
	160	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	161	c . "ave(X)", zsto,zout)
[471]	162	c
[455]	163	CALL histdef(nid_mth, "ffonte","Thermal flux for snow melting",
	164	. "W/m2",iim,jjmp1,nhori, 1,1,1, -99, 32,
	165	. "ave(X)", zsto,zout)
[416]	166	c
[455]	167	CALL histdef(nid_mth, "fqcalving","Ice Calving",
	168	. "kg/m2/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
	169	. "ave(X)", zsto,zout)
	170	c
[416]	171	DO nsrf = 1, nbsrf
	172	C
	173	call histdef(nid_mth, "pourc_"//clnsurf(nsrf),
	174	$ "Fraction "//clnsurf(nsrf), "W/m2",
	175	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	176	$ "ave(X)", zsto,zout)
	177	C
	178	call histdef(nid_mth, "tsol_"//clnsurf(nsrf),
	179	$ "Fraction "//clnsurf(nsrf), "W/m2",
	180	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	181	$ "ave(X)", zsto,zout)
	182	C
	183	call histdef(nid_mth, "sens_"//clnsurf(nsrf),
	184	$ "Sensible heat flux "//clnsurf(nsrf), "W/m2",
	185	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	186	$ "ave(X)", zsto,zout)
	187	c
	188	call histdef(nid_mth, "lat_"//clnsurf(nsrf),
	189	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
	190	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	191	$ "ave(X)", zsto,zout)
	192	C
[433]	193	call histdef(nid_mth, "flw_"//clnsurf(nsrf),
	194	$ "LW "//clnsurf(nsrf), "W/m2",
	195	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	196	$ "ave(X)", zsto,zout)
[471]	197	c
[433]	198	call histdef(nid_mth, "fsw_"//clnsurf(nsrf),
	199	$ "SW "//clnsurf(nsrf), "W/m2",
	200	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	201	$ "ave(X)", zsto,zout)
	202	C
	203	call histdef(nid_mth, "wbils_"//clnsurf(nsrf),
	204	$ "Bilan sol "//clnsurf(nsrf), "W/m2",
	205	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	206	$ "ave(X)", zsto,zout)
	207	C
[416]	208	call histdef(nid_mth, "taux_"//clnsurf(nsrf),
	209	$ "Zonal wind stress"//clnsurf(nsrf), "Pa",
	210	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	211	$ "ave(X)", zsto,zout)
	212
	213	call histdef(nid_mth, "tauy_"//clnsurf(nsrf),
	214	$ "Meridional xind stress "//clnsurf(nsrf), "Pa",
	215	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	216	$ "ave(X)", zsto,zout)
	217	c
	218	call histdef(nid_mth, "albe_"//clnsurf(nsrf),
[440]	219	$ "Albedo surf. "//clnsurf(nsrf), "-",
[416]	220	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	221	$ "ave(X)", zsto,zout)
	222	c
	223	call histdef(nid_mth, "rugs_"//clnsurf(nsrf),
	224	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
	225	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	226	$ "ave(X)", zsto,zout)
	227	c
	228	CALL histdef(nid_mth, "ages_"//clnsurf(nsrf), "Snow age","day",
	229	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	230	. "ave(X)", zsto,zout)
	231
	232	END DO
	233	C
[470]	234	c CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "-",
	235	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	236	c . "ave(X)", zsto,zout)
[416]	237	c
	238	CALL histdef(nid_mth, "albs", "Surface albedo", "-",
	239	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	240	. "ave(X)", zsto,zout)
	241	CALL histdef(nid_mth, "albslw", "Surface albedo LW", "-",
	242	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	243	. "ave(X)", zsto,zout)
	244	c
	245	CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
	246	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	247	. "ave(X)", zsto,zout)
	248	c
	249	CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
	250	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	251	. "ave(X)", zsto,zout)
	252	c
	253	CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-",
	254	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	255	. "ave(X)", zsto,zout)
	256	c
	257	CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-",
	258	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	259	. "ave(X)", zsto,zout)
	260	c
	261	CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-",
	262	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	263	. "ave(X)", zsto,zout)
	264	c
	265	CALL histdef(nid_mth, "cldt", "Total cloudiness", "-",
	266	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	267	. "ave(X)", zsto,zout)
	268	c
[486]	269	CALL histdef(nid_mth,"cldq","Cloud liquid water path","Kg/m2",
[416]	270	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	271	. "ave(X)", zsto,zout)
[486]	272	cIM: 071003
	273	CALL histdef(nid_mth,"lwp","Cloud water path","Kg/m2",
	274	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	275	. "ave(X)", zsto,zout)
[416]	276	c
[486]	277	CALL histdef(nid_mth,"iwp","Cloud ice water path","Kg/m2",
	278	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	279	. "ave(X)", zsto,zout)
	280	cIM: 071003
	281	CALL histdef(nid_mth,"lwcon","Cloud water content","Kg/Kg",
	282	. iim,jjmp1,nhori, klev,1,klev, nvert, 32,
	283	. "ave(X)", zsto,zout)
	284	c
	285	CALL histdef(nid_mth,"iwcon","Cloud ice water content","Kg/Kg",
	286	. iim,jjmp1,nhori, klev,1,klev, nvert, 32,
	287	. "ave(X)", zsto,zout)
	288	c
[416]	289	CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
	290	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	291	. "ave(X)", zsto,zout)
	292	c
	293	CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
	294	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	295	. "ave(X)", zsto,zout)
	296	c
	297	CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-",
	298	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	299	. "ave(X)", zsto,zout)
	300	c
	301	CALL histdef(nid_mth, "vq", "Merid humidity transport", "-",
	302	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	303	. "ave(X)", zsto,zout)
	304	cKE43
	305	IF (iflag_con .GE. 3) THEN ! sb
	306	c
	307	CALL histdef(nid_mth, "cape", "Conv avlbl pot ener", "J/Kg",
	308	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	309	. "ave(X)", zsto,zout)
	310	c
	311	CALL histdef(nid_mth, "pbase", "Cld base pressure", "hPa",
	312	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	313	. "ave(X)", zsto,zout)
	314	c
	315	CALL histdef(nid_mth, "ptop", "Cld top pressure", "hPa",
	316	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	317	. "ave(X)", zsto,zout)
	318	c
	319	CALL histdef(nid_mth, "fbase", "Cld base mass flux", "Kg/m2/s",
	320	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	321	. "ave(X)", zsto,zout)
	322	c
	323	c
[471]	324	CALL histdef(nid_mth, "upwd", "saturated updraft", "Kg/m2/s",
	325	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	326	. "ave(X)", zsto,zout)
	327	c
	328	CALL histdef(nid_mth, "dnwd", "saturated downdraft","Kg/m2/s",
	329	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	330	. "ave(X)", zsto,zout)
	331	c
	332	CALL histdef(nid_mth, "dnwd0", "unsat. downdraft", "Kg/m2/s",
	333	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	334	. "ave(X)", zsto,zout)
	335	c
	336	CALL histdef(nid_mth,"Ma","undilute adiab updraft","Kg/m2/s",
	337	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	338	. "ave(X)", zsto,zout)
	339	c
	340	CALL histdef(nid_mth, "prw", "Precipitable water", "kg/m2",
	341	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	342	. "ave(X)", zsto,zout)
	343	c
[416]	344	ENDIF
	345	c34EK
	346	c
	347	c Champs 3D:
	348	c
	349	CALL histdef(nid_mth, "temp", "Air temperature", "K",
	350	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	351	. "ave(X)", zsto,zout)
	352	c
	353	CALL histdef(nid_mth, "ovap", "Specific humidity", "Kg/Kg",
	354	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	355	. "ave(X)", zsto,zout)
[486]	356	cIM: 071003
	357	CALL histdef(nid_mth,"wvap","Water vapor mixing ratio","Kg/Kg",
	358	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	359	. "ave(X)", zsto,zout)
[416]	360	c
	361	CALL histdef(nid_mth, "geop", "Geopotential height", "m",
	362	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	363	. "ave(X)", zsto,zout)
	364	c
	365	CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
	366	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	367	. "ave(X)", zsto,zout)
	368	c
	369	CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
	370	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	371	. "ave(X)", zsto,zout)
	372	c
	373	CALL histdef(nid_mth, "vitw", "Vertical wind", "m/s",
	374	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	375	. "ave(X)", zsto,zout)
	376	c
	377	CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
	378	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	379	. "ave(X)", zsto,zout)
	380	c
	381	CALL histdef(nid_mth, "rneb", "Cloud fraction", "-",
	382	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	383	. "ave(X)", zsto,zout)
[471]	384	c
[433]	385	CALL histdef(nid_mth, "rnebcon", "Convective Cloud Fraction"
	386	. , "-",
	387	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	388	. "ave(X)", zsto,zout)
[416]	389	c
	390	CALL histdef(nid_mth, "rhum", "Relative humidity", "-",
	391	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	392	. "ave(X)", zsto,zout)
	393	c
	394	CALL histdef(nid_mth, "clwcon", "Cloud Liquid water content"
	395	. , "kg/kg",
	396	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	397	. "ave(X)", zsto,zout)
	398	c
	399	CALL histdef(nid_mth, "oliq", "Liquid water content", "kg/kg",
	400	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	401	. "ave(X)", zsto,zout)
	402	c
[486]	403	CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
	404	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	405	. "ave(X)", zsto,zout)
	406	c
[416]	407	CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
	408	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	409	. "ave(X)", zsto,zout)
	410	c
	411	CALL histdef(nid_mth, "dqdyn", "Dynamics dQ", "Kg/Kg/s",
	412	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	413	. "ave(X)", zsto,zout)
	414	c
	415	CALL histdef(nid_mth, "dtcon", "Convection dT", "K/s",
	416	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	417	. "ave(X)", zsto,zout)
	418	c
[470]	419	c CALL histdef(nid_mth, "ducon", "Convection du", "m/s2",
	420	c . iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	421	c . "ave(X)", zsto,zout)
[416]	422	c
	423	CALL histdef(nid_mth, "dqcon", "Convection dQ", "Kg/Kg/s",
	424	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	425	. "ave(X)", zsto,zout)
	426	c
	427	CALL histdef(nid_mth, "dtlsc", "Condensation dT", "K/s",
	428	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	429	. "ave(X)", zsto,zout)
[486]	430	cIM: 071003
	431	CALL histdef(nid_mth, "dtlschr",
	432	$ "Large-scale condensational heating rate", "K/s",iim,jjmp1
	433	$ ,nhori, klev,1,klev,nvert, 32,"ave(X)", zsto,zout)
[416]	434	c
	435	CALL histdef(nid_mth, "dqlsc", "Condensation dQ", "Kg/Kg/s",
	436	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	437	. "ave(X)", zsto,zout)
	438	c
	439	CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
	440	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	441	. "ave(X)", zsto,zout)
	442	c
	443	CALL histdef(nid_mth, "dqvdf", "Boundary-layer dQ", "Kg/Kg/s",
	444	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	445	. "ave(X)", zsto,zout)
	446	c
	447	CALL histdef(nid_mth, "dteva", "Reevaporation dT", "K/s",
	448	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	449	. "ave(X)", zsto,zout)
	450	c
	451	CALL histdef(nid_mth, "dqeva", "Reevaporation dQ", "Kg/Kg/s",
	452	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	453	. "ave(X)", zsto,zout)
	454
	455	CALL histdef(nid_mth, "ptconv", "POINTS CONVECTIFS"," ",
	456	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	457	. "ave(X)", zsto,zout)
	458
	459	CALL histdef(nid_mth, "ratqs", "RATQS"," ",
	460	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	461	. "ave(X)", zsto,zout)
	462
	463	c
	464	CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s",
	465	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	466	. "ave(X)", zsto,zout)
	467
	468	CALL histdef(nid_mth, "dqajs", "Dry adjust. dQ", "Kg/Kg/s",
	469	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	470	. "ave(X)", zsto,zout)
	471	c
	472	CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
	473	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	474	. "ave(X)", zsto,zout)
	475	c
[486]	476	CALL histdef(nid_mth, "dtsw0", "CS SW radiation dT", "K/s",
[416]	477	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	478	. "ave(X)", zsto,zout)
	479	c
	480	CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
	481	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	482	. "ave(X)", zsto,zout)
	483	c
[486]	484	CALL histdef(nid_mth,"dtlw0","CS LW radiation dT","K/s",
[416]	485	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	486	. "ave(X)", zsto,zout)
	487	c
	488	CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
	489	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	490	. "ave(X)", zsto,zout)
	491	c
	492	CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
	493	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	494	. "ave(X)", zsto,zout)
	495	c
	496	CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
	497	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	498	. "ave(X)", zsto,zout)
	499	c
	500	IF (ok_orodr) THEN
	501	CALL histdef(nid_mth, "duoro", "Orography dU", "m/s2",
	502	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	503	. "ave(X)", zsto,zout)
	504	c
	505	CALL histdef(nid_mth, "dvoro", "Orography dV", "m/s2",
	506	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	507	. "ave(X)", zsto,zout)
	508	c
	509	ENDIF
	510	C
	511	IF (ok_orolf) THEN
	512	CALL histdef(nid_mth, "dulif", "Orography dU", "m/s2",
	513	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	514	. "ave(X)", zsto,zout)
	515	c
	516	CALL histdef(nid_mth, "dvlif", "Orography dV", "m/s2",
	517	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	518	. "ave(X)", zsto,zout)
	519	ENDIF
	520	C
	521	CALL histdef(nid_mth, "ozone", "Ozone concentration", "-",
	522	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	523	. "ave(X)", zsto,zout)
	524	c
	525	if (nqmax.GE.3) THEN
	526	DO iq=1,nqmax-2
	527	IF (iq.LE.99) THEN
	528	WRITE(str2,'(i2.2)') iq
	529	CALL histdef(nid_mth, "trac"//str2, "Tracer No."//str2, "-",
	530	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	531	. "ave(X)", zsto,zout)
	532	ELSE
	533	PRINT*, "Trop de traceurs"
	534	CALL abort
	535	ENDIF
	536	ENDDO
	537	ENDIF
	538	c
	539	CALL histdef(nid_mth, "SWupTOA", "SWup at TOA","W/m2",
	540	. iim,jjmp1,nhori, 1,1,1,-99,
	541	. 32, "ave(X)", zsto,zout)
	542	c
	543	CALL histdef(nid_mth, "SWupSFC", "SWup at surface","W/m2",
	544	. iim,jjmp1,nhori, 1,1,1,-99,
	545	. 32, "ave(X)", zsto,zout)
	546	c
	547	CALL histdef(nid_mth, "SWdnTOA", "SWdn at TOA","W/m2",
	548	. iim,jjmp1,nhori, 1,1,1,-99,
	549	. 32, "ave(X)", zsto,zout)
	550	c
	551	CALL histdef(nid_mth, "SWdnSFC", "SWdn at surface","W/m2",
	552	. iim,jjmp1,nhori, 1,1,1,-99,
	553	. 32, "ave(X)", zsto,zout)
[471]	554	c
[433]	555	CALL histdef(nid_mth, "SWupTOAclr",
	556	. "SWup clear sky at TOA","W/m2",
	557	. iim,jjmp1,nhori, 1,1,1,-99,
	558	. 32, "ave(X)", zsto,zout)
[416]	559	c
[433]	560	CALL histdef(nid_mth, "SWupSFCclr",
	561	. "SWup clear sky at surface","W/m2",
	562	. iim,jjmp1,nhori, 1,1,1,-99,
	563	. 32, "ave(X)", zsto,zout)
	564	c
	565	CALL histdef(nid_mth, "SWdnTOAclr",
	566	. "SWdn clear sky at TOA","W/m2",
	567	. iim,jjmp1,nhori, 1,1,1,-99,
	568	. 32, "ave(X)", zsto,zout)
[440]	569	c
	570	CALL histdef(nid_mth, "SWdnSFCclr",
	571	. "SWdn clear sky at surface","W/m2",
	572	. iim,jjmp1,nhori, 1,1,1,-99,
	573	. 32, "ave(X)", zsto,zout)
[467]	574
[433]	575	c Champs interpolles sur des niveaux de pression
	576
	577	CALL histdef(nid_mth, "u850", "Zonal wind 850mb", "m/s",
	578	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	579	. "ave(X)", zsto,zout)
	580
	581	CALL histdef(nid_mth, "v850", "Meridional wind 850mb", "m/s",
	582	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	583	. "ave(X)", zsto,zout)
	584	c
	585	CALL histdef(nid_mth, "u500", "Zonal wind 500mb", "m/s",
	586	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	587	. "ave(X)", zsto,zout)
	588
	589	CALL histdef(nid_mth, "v500", "Meridional wind 500mb", "m/s",
	590	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	591	. "ave(X)", zsto,zout)
	592
	593	CALL histdef(nid_mth, "u200", "Zonal wind 200mb", "m/s",
	594	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	595	. "ave(X)", zsto,zout)
	596
	597	CALL histdef(nid_mth, "v200", "Meridional wind 200mb", "m/s",
	598	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	599	. "ave(X)", zsto,zout)
	600
	601	CALL histdef(nid_mth, "phi500","Geopotentiel à 500mb","m2/s2",
	602	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	603	. "ave(X)", zsto,zout)
	604
	605	CALL histdef(nid_mth, "w500", "Vertical wind 500mb", "Pa/s",
	606	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	607	. "ave(X)", zsto,zout)
[471]	608	c
[416]	609	CALL histend(nid_mth)
	610	c
	611	ndex2d = 0
	612	ndex3d = 0
	613	c
	614	ENDIF ! fin de test sur ok_mensuel

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

Download in other formats: