Context Navigation

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

Last change on this file since 544 was 517, checked in by lmdzadmin, 21 years ago
Inclusion des modifications de O. Boucher et de J. Quaas pour le calcul des premiers effets directs et indirects dus aux aerosols 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),
	462	. "Vertical wind"//clev(k)//"mb","m/s",
	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
	511	CALL histdef(nid_mth, "vitw", "Vertical wind", "m/s",
	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.

Download in other formats: