Context Navigation

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

Last change on this file since 510 was 506, checked in by lmdzadmin, 21 years ago

IM: ajout lev_hist(hf,day,mth) niveaux output definis dans physiq.def et

variables ENSEMBLES + rearrangement

Property svn:eol-style set to native
Property svn:keywords set to Author Date Id Revision

File size: 27.9 KB

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

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

Download in other formats: