Context Navigation

ini_histmth.h @ 694

Last change on this file since 694 was 634, checked in by Laurent Fairhead, 19 years ago
Modifications faites à la physique pour la rendre parallele YM Une branche de travail LMDZ4_par_0 a été créée provisoirement afin de tester les modifs pleinement avant leurs inclusions dans le tronc principal LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 29.5 KB

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

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

Context Navigation

source: LMDZ4/branches/LMDZ4_par_0/libf/phylmd/ini_histmth.h @ 694

Download in other formats: