Context Navigation

ini_histmth.h @ 5456

Last change on this file since 5456 was 574, checked in by Laurent Fairhead, 20 years ago
Erreur sur l'unite de la vitesse verticale LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 29.1 KB

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

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

Context Navigation

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

Download in other formats: