Context Navigation

ini_histmth.h

Last change on this file was 416, checked in by lmdzadmin, 22 years ago
Inclusion initiale
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 19.0 KB

Rev	Line
[416]	1	IF (ok_mensuel) THEN
	2	c
	3	idayref = day_ref
	4	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
	5	c
	6	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
	7	DO i = 1, iim
	8	zx_lon(i,1) = rlon(i+1)
	9	zx_lon(i,jjmp1) = rlon(i+1)
	10	ENDDO
	11	DO ll=1,klev
	12	znivsig(ll)=float(ll)
	13	ENDDO
	14	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
	15	CALL histbeg("histmth.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
	16	. 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
	17	. nhori, nid_mth)
	18	write(,)'Mensuel ', itau_phy, zjulian
	19	CALL histvert(nid_mth, "presnivs", "Vertical levels", "mb",
	20	. klev, presnivs, nvert)
	21	c call histvert(nid_mth, 'sig_s', 'Niveaux sigma','-',
	22	c . klev, znivsig, nvert)
	23	c
	24	zsto = dtime
	25	zout = dtime * ecrit_mth
	26	c
	27	CALL histdef(nid_mth, "phis", "Surface geop. height", "-",
	28	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	29	. "once", zsto,zout)
	30	c
	31	CALL histdef(nid_mth, "aire", "Grid area", "-",
	32	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	33	. "once", zsto,zout)
	34	c
	35	c Champs 2D:
	36	c
	37	CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
	38	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	39	. "ave(X)", zsto,zout)
	40	c
	41	CALL histdef(nid_mth, "t2m", "Temperature 2m", "K",
	42	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	43	. "ave(X)", zsto,zout)
	44	c
	45	CALL histdef(nid_mth, "q2m", "Specific humidity 2m", "Kg/Kg",
	46	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	47	. "ave(X)", zsto,zout)
	48	c
	49	CALL histdef(nid_mth, "u10m", "Vent zonal 10m", "m/s",
	50	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	51	. "ave(X)", zsto,zout)
	52	c
	53	CALL histdef(nid_mth, "v10m", "Vent meridien 10m", "m/s",
	54	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	55	. "ave(X)", zsto,zout)
	56	c
	57	c
	58	CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
	59	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	60	. "ave(X)", zsto,zout)
	61	c
	62	CALL histdef(nid_mth, "qsol", "Surface humidity", "mm",
	63	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	64	. "ave(X)", zsto,zout)
	65	c
	66	CALL histdef(nid_mth, "precip", "Precipitation Totale liq+sol",
	67	. "kg/s",
	68	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	69	. "ave(X)", zsto,zout)
	70	c
	71	CALL histdef(nid_mth, "plul", "Large-scale Precip.", "kg/s",
	72	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	73	. "ave(X)", zsto,zout)
	74	c
	75	CALL histdef(nid_mth, "pluc", "Convective Precip.", "kg/s",
	76	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	77	. "ave(X)", zsto,zout)
	78	c
	79	CALL histdef(nid_mth, "snow", "Snow fall", "kg/s",
	80	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	81	. "ave(X)", zsto,zout)
	82	c
	83	CALL histdef(nid_mth, "snow_mass", "Snow Mass", "kg/m2",
	84	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	85	. "ave(X)", zsto,zout)
	86	c
	87	CALL histdef(nid_mth, "evap", "Evaporation", "kg/s",
	88	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	89	. "ave(X)", zsto,zout)
	90	c
	91	CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
	92	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	93	. "ave(X)", zsto,zout)
	94	c
	95	CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
	96	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	97	. "ave(X)", zsto,zout)
	98	c
	99	CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
	100	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	101	. "ave(X)", zsto,zout)
	102	c
	103	CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
	104	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	105	. "ave(X)", zsto,zout)
	106	c
	107	CALL histdef(nid_mth, "solldown", "Down. IR rad. at surface",
	108	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
	109	. "ave(X)", zsto,zout)
	110	c
	111	CALL histdef(nid_mth, "tops0", "Solar rad. at TOA", "W/m2",
	112	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	113	. "ave(X)", zsto,zout)
	114	c
	115	CALL histdef(nid_mth, "topl0", "IR rad. at TOA", "W/m2",
	116	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	117	. "ave(X)", zsto,zout)
	118	c
	119	CALL histdef(nid_mth, "sols0", "Solar rad. at surf.", "W/m2",
	120	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	121	. "ave(X)", zsto,zout)
	122	c
	123	CALL histdef(nid_mth, "soll0", "IR rad. at surface", "W/m2",
	124	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	125	. "ave(X)", zsto,zout)
	126	c
	127	CALL histdef(nid_mth, "bils", "Surf. total heat flux", "W/m2",
	128	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	129	. "ave(X)", zsto,zout)
	130	c
	131	CALL histdef(nid_mth, "sens", "Sensible heat flux", "W/m2",
	132	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	133	. "ave(X)", zsto,zout)
	134	c
	135	CALL histdef(nid_mth, "fder", "Heat flux derivation", "W/m2",
	136	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	137	. "ave(X)", zsto,zout)
	138	c
	139	CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
	140	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	141	. "ave(X)", zsto,zout)
	142	c
	143	CALL histdef(nid_mth, "frtv", "Meridional wind stress", "Pa",
	144	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	145	. "ave(X)", zsto,zout)
	146	c
	147	DO nsrf = 1, nbsrf
	148	C
	149	call histdef(nid_mth, "pourc_"//clnsurf(nsrf),
	150	$ "Fraction "//clnsurf(nsrf), "W/m2",
	151	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	152	$ "ave(X)", zsto,zout)
	153	C
	154	call histdef(nid_mth, "tsol_"//clnsurf(nsrf),
	155	$ "Fraction "//clnsurf(nsrf), "W/m2",
	156	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	157	$ "ave(X)", zsto,zout)
	158	C
	159	call histdef(nid_mth, "sens_"//clnsurf(nsrf),
	160	$ "Sensible heat flux "//clnsurf(nsrf), "W/m2",
	161	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	162	$ "ave(X)", zsto,zout)
	163	c
	164	call histdef(nid_mth, "lat_"//clnsurf(nsrf),
	165	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
	166	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	167	$ "ave(X)", zsto,zout)
	168	C
	169	call histdef(nid_mth, "taux_"//clnsurf(nsrf),
	170	$ "Zonal wind stress"//clnsurf(nsrf), "Pa",
	171	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	172	$ "ave(X)", zsto,zout)
	173
	174	call histdef(nid_mth, "tauy_"//clnsurf(nsrf),
	175	$ "Meridional xind stress "//clnsurf(nsrf), "Pa",
	176	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	177	$ "ave(X)", zsto,zout)
	178	c
	179	call histdef(nid_mth, "albe_"//clnsurf(nsrf),
	180	$ "Albedo surf. "//clnsurf(nsrf), "W/m2",
	181	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	182	$ "ave(X)", zsto,zout)
	183	c
	184	call histdef(nid_mth, "rugs_"//clnsurf(nsrf),
	185	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
	186	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	187	$ "ave(X)", zsto,zout)
	188	c
	189	CALL histdef(nid_mth, "ages_"//clnsurf(nsrf), "Snow age","day",
	190	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	191	. "ave(X)", zsto,zout)
	192
	193	END DO
	194	C
	195	CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "-",
	196	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	197	. "ave(X)", zsto,zout)
	198	c
	199	CALL histdef(nid_mth, "albs", "Surface albedo", "-",
	200	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	201	. "ave(X)", zsto,zout)
	202	CALL histdef(nid_mth, "albslw", "Surface albedo LW", "-",
	203	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	204	. "ave(X)", zsto,zout)
	205	c
	206	CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
	207	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	208	. "ave(X)", zsto,zout)
	209	c
	210	CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
	211	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	212	. "ave(X)", zsto,zout)
	213	c
	214	CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-",
	215	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	216	. "ave(X)", zsto,zout)
	217	c
	218	CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-",
	219	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	220	. "ave(X)", zsto,zout)
	221	c
	222	CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-",
	223	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	224	. "ave(X)", zsto,zout)
	225	c
	226	CALL histdef(nid_mth, "cldt", "Total cloudiness", "-",
	227	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	228	. "ave(X)", zsto,zout)
	229	c
	230	CALL histdef(nid_mth, "cldq", "Cloud liquid water path", "-",
	231	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	232	. "ave(X)", zsto,zout)
	233	c
	234	CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
	235	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	236	. "ave(X)", zsto,zout)
	237	c
	238	CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
	239	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	240	. "ave(X)", zsto,zout)
	241	c
	242	CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-",
	243	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	244	. "ave(X)", zsto,zout)
	245	c
	246	CALL histdef(nid_mth, "vq", "Merid humidity transport", "-",
	247	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	248	. "ave(X)", zsto,zout)
	249	cKE43
	250	IF (iflag_con .GE. 3) THEN ! sb
	251	c
	252	CALL histdef(nid_mth, "cape", "Conv avlbl pot ener", "J/Kg",
	253	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	254	. "ave(X)", zsto,zout)
	255	c
	256	CALL histdef(nid_mth, "pbase", "Cld base pressure", "hPa",
	257	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	258	. "ave(X)", zsto,zout)
	259	c
	260	CALL histdef(nid_mth, "ptop", "Cld top pressure", "hPa",
	261	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	262	. "ave(X)", zsto,zout)
	263	c
	264	CALL histdef(nid_mth, "fbase", "Cld base mass flux", "Kg/m2/s",
	265	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	266	. "ave(X)", zsto,zout)
	267	c
	268	c
	269	ENDIF
	270	c34EK
	271	c
	272	c Champs 3D:
	273	c
	274	CALL histdef(nid_mth, "temp", "Air temperature", "K",
	275	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	276	. "ave(X)", zsto,zout)
	277	c
	278	CALL histdef(nid_mth, "ovap", "Specific humidity", "Kg/Kg",
	279	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	280	. "ave(X)", zsto,zout)
	281	c
	282	CALL histdef(nid_mth, "geop", "Geopotential height", "m",
	283	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	284	. "ave(X)", zsto,zout)
	285	c
	286	CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
	287	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	288	. "ave(X)", zsto,zout)
	289	c
	290	CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
	291	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	292	. "ave(X)", zsto,zout)
	293	c
	294	CALL histdef(nid_mth, "vitw", "Vertical wind", "m/s",
	295	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	296	. "ave(X)", zsto,zout)
	297	c
	298	CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
	299	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	300	. "ave(X)", zsto,zout)
	301	c
	302	CALL histdef(nid_mth, "rneb", "Cloud fraction", "-",
	303	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	304	. "ave(X)", zsto,zout)
	305	c
	306	CALL histdef(nid_mth, "rhum", "Relative humidity", "-",
	307	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	308	. "ave(X)", zsto,zout)
	309	c
	310	CALL histdef(nid_mth, "clwcon", "Cloud Liquid water content"
	311	. , "kg/kg",
	312	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	313	. "ave(X)", zsto,zout)
	314	c
	315	CALL histdef(nid_mth, "oliq", "Liquid water content", "kg/kg",
	316	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	317	. "ave(X)", zsto,zout)
	318	c
	319	CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
	320	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	321	. "ave(X)", zsto,zout)
	322	c
	323	CALL histdef(nid_mth, "dqdyn", "Dynamics dQ", "Kg/Kg/s",
	324	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	325	. "ave(X)", zsto,zout)
	326	c
	327	CALL histdef(nid_mth, "dtcon", "Convection dT", "K/s",
	328	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	329	. "ave(X)", zsto,zout)
	330	c
	331	CALL histdef(nid_mth, "ducon", "Convection du", "m/s2",
	332	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	333	. "ave(X)", zsto,zout)
	334	c
	335	CALL histdef(nid_mth, "dqcon", "Convection dQ", "Kg/Kg/s",
	336	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	337	. "ave(X)", zsto,zout)
	338	c
	339	CALL histdef(nid_mth, "dtlsc", "Condensation dT", "K/s",
	340	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	341	. "ave(X)", zsto,zout)
	342	c
	343	CALL histdef(nid_mth, "dqlsc", "Condensation dQ", "Kg/Kg/s",
	344	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	345	. "ave(X)", zsto,zout)
	346	c
	347	CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
	348	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	349	. "ave(X)", zsto,zout)
	350	c
	351	CALL histdef(nid_mth, "dqvdf", "Boundary-layer dQ", "Kg/Kg/s",
	352	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	353	. "ave(X)", zsto,zout)
	354	c
	355	CALL histdef(nid_mth, "dteva", "Reevaporation dT", "K/s",
	356	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	357	. "ave(X)", zsto,zout)
	358	c
	359	CALL histdef(nid_mth, "dqeva", "Reevaporation dQ", "Kg/Kg/s",
	360	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	361	. "ave(X)", zsto,zout)
	362
	363	CALL histdef(nid_mth, "ptconv", "POINTS CONVECTIFS"," ",
	364	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	365	. "ave(X)", zsto,zout)
	366
	367	CALL histdef(nid_mth, "ratqs", "RATQS"," ",
	368	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	369	. "ave(X)", zsto,zout)
	370
	371	c
	372	CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s",
	373	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	374	. "ave(X)", zsto,zout)
	375
	376	CALL histdef(nid_mth, "dqajs", "Dry adjust. dQ", "Kg/Kg/s",
	377	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	378	. "ave(X)", zsto,zout)
	379	c
	380	CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
	381	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	382	. "ave(X)", zsto,zout)
	383	c
	384	CALL histdef(nid_mth, "dtsw0", "SW radiation dT", "K/s",
	385	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	386	. "ave(X)", zsto,zout)
	387	c
	388	CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
	389	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	390	. "ave(X)", zsto,zout)
	391	c
	392	CALL histdef(nid_mth, "dtlw0", "LW radiation dT", "K/s",
	393	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	394	. "ave(X)", zsto,zout)
	395	c
	396	CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
	397	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	398	. "ave(X)", zsto,zout)
	399	c
	400	CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
	401	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	402	. "ave(X)", zsto,zout)
	403	c
	404	CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
	405	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	406	. "ave(X)", zsto,zout)
	407	c
	408	IF (ok_orodr) THEN
	409	CALL histdef(nid_mth, "duoro", "Orography dU", "m/s2",
	410	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	411	. "ave(X)", zsto,zout)
	412	c
	413	CALL histdef(nid_mth, "dvoro", "Orography dV", "m/s2",
	414	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	415	. "ave(X)", zsto,zout)
	416	c
	417	ENDIF
	418	C
	419	IF (ok_orolf) THEN
	420	CALL histdef(nid_mth, "dulif", "Orography dU", "m/s2",
	421	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	422	. "ave(X)", zsto,zout)
	423	c
	424	CALL histdef(nid_mth, "dvlif", "Orography dV", "m/s2",
	425	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	426	. "ave(X)", zsto,zout)
	427	ENDIF
	428	C
	429	CALL histdef(nid_mth, "ozone", "Ozone concentration", "-",
	430	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	431	. "ave(X)", zsto,zout)
	432	c
	433	if (nqmax.GE.3) THEN
	434	DO iq=1,nqmax-2
	435	IF (iq.LE.99) THEN
	436	WRITE(str2,'(i2.2)') iq
	437	CALL histdef(nid_mth, "trac"//str2, "Tracer No."//str2, "-",
	438	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	439	. "ave(X)", zsto,zout)
	440	ELSE
	441	PRINT*, "Trop de traceurs"
	442	CALL abort
	443	ENDIF
	444	ENDDO
	445	ENDIF
	446	c
	447	cKE43
	448	IF (iflag_con.GE.3) THEN ! (sb)
	449	c
	450	CALL histdef(nid_mth, "upwd", "saturated updraft", "Kg/m2/s",
	451	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	452	. "ave(X)", zsto,zout)
	453	c
	454	CALL histdef(nid_mth, "dnwd", "saturated downdraft","Kg/m2/s",
	455	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	456	. "ave(X)", zsto,zout)
	457	c
	458	CALL histdef(nid_mth, "dnwd0", "unsat. downdraft", "Kg/m2/s",
	459	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	460	. "ave(X)", zsto,zout)
	461	c
	462	CALL histdef(nid_mth,"Ma","undilute adiab updraft","Kg/m2/s",
	463	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	464	. "ave(X)", zsto,zout)
	465	cccIM
	466	CALL histdef(nid_mth, "SWupTOA", "SWup at TOA","W/m2",
	467	. iim,jjmp1,nhori, 1,1,1,-99,
	468	. 32, "ave(X)", zsto,zout)
	469	c
	470	CALL histdef(nid_mth, "SWupSFC", "SWup at surface","W/m2",
	471	. iim,jjmp1,nhori, 1,1,1,-99,
	472	. 32, "ave(X)", zsto,zout)
	473	c
	474	CALL histdef(nid_mth, "SWdnTOA", "SWdn at TOA","W/m2",
	475	. iim,jjmp1,nhori, 1,1,1,-99,
	476	. 32, "ave(X)", zsto,zout)
	477	c
	478	CALL histdef(nid_mth, "SWdnSFC", "SWdn at surface","W/m2",
	479	. iim,jjmp1,nhori, 1,1,1,-99,
	480	. 32, "ave(X)", zsto,zout)
	481	c
	482	ENDIF
	483	c34EK
	484	CALL histend(nid_mth)
	485	c
	486	ndex2d = 0
	487	ndex3d = 0
	488	c
	489	ENDIF ! fin de test sur ok_mensuel

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

Context Navigation

source: LMDZ.3.3/trunk/libf/phylmd/ini_histmth.h

Download in other formats: