Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Normal
Revision Log

ini_histday.h @ 3670

Last change on this file since 3670 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: 13.3 KB

Rev	Line
[416]	1	IF (ok_journe) 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	write(,)'zx_lon = ',zx_lon(:,1)
	16	write(,)'zx_lat = ',zx_lat(1,:)
	17	CALL histbeg("histday", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
	18	. 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
	19	. nhori, nid_day)
	20	write(,)'Journee ', itau_phy, zjulian
	21	CALL histvert(nid_day, "presnivs", "Vertical levels", "mb",
	22	. klev, presnivs, nvert)
	23	c call histvert(nid_day, 'sig_s', 'Niveaux sigma','-',
	24	c . klev, znivsig, nvert)
	25	c
	26	zsto = dtime
	27	zout = dtime * FLOAT(ecrit_day)
	28	C Essai writephys
	29	c nom_fichier = 'histday1'
	30	c call writephy_ini(fid_day,nom_fichier,klon,iim,jjmp1,klev,
	31	c . rlon,rlat, presnivs,
	32	c . zjulian, dtime)
	33	c call writephy_def(prof2d_on, fid_day, "once", zsto, zout, 0)
	34	c call writephy_def(prof3d_on, fid_day, "once", zsto, zout,
	35	c . klev)
	36	c call writephy_def(prof2d_av, fid_day, "ave(X)", zsto, zout, 0)
	37	c call writephy_def(prof3d_av, fid_day, "ave(X)", zsto, zout,
	38	c . klev)
	39
	40	c
	41	CALL histdef(nid_day, "phis", "Surface geop. height", "-",
	42	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	43	. "once", zsto,zout)
	44	c
	45	CALL histdef(nid_day, "aire", "Grid area", "-",
	46	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	47	. "once", zsto,zout)
	48	c
	49	c Champs 2D:
	50	c
	51	CALL histdef(nid_day, "tsol", "Surface Temperature", "K",
	52	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	53	. "ave(X)", zsto,zout)
	54	c
	55	CALL histdef(nid_day, "tter", "Surface Temperature", "K",
	56	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	57	. "ave(X)", zsto,zout)
	58	c
	59	CALL histdef(nid_day, "tlic", "Surface Temperature", "K",
	60	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	61	. "ave(X)", zsto,zout)
	62	c
	63	CALL histdef(nid_day, "toce", "Surface Temperature", "K",
	64	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	65	. "ave(X)", zsto,zout)
	66	c
	67	CALL histdef(nid_day, "tsic", "Surface Temperature", "K",
	68	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	69	. "ave(X)", zsto,zout)
	70	c
	71	cccIM
	72	c
	73	CALL histdef(nid_day, "t2m", "Temperature 2m", "K",
	74	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	75	. "ave(X)", zsto,zout)
	76	c
	77	CALL histdef(nid_day, "t2mter", "Temp.terre 2m", "K",
	78	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	79	. "ave(X)", zsto,zout)
	80	c
	81	CALL histdef(nid_day, "t2mlic", "Temp.lic 2m", "K",
	82	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	83	. "ave(X)", zsto,zout)
	84	c
	85	CALL histdef(nid_day, "t2moce", "Temp.oce 2m", "K",
	86	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	87	. "ave(X)", zsto,zout)
	88	c
	89	CALL histdef(nid_day, "t2msic", "Temp.sic 2m", "K",
	90	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	91	. "ave(X)", zsto,zout)
	92	c
	93	CALL histdef(nid_day, "t2m_min", "Temp. 2m min.",
	94	. "K",
	95	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	96	. t2mincels, zsto,zout)
	97	c
	98	CALL histdef(nid_day, "t2m_max", "Temp. 2m max.",
	99	. "K",
	100	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	101	. t2maxcels, zsto,zout)
	102	c
	103	CALL histdef(nid_day, "t2mter_min", "Temp.terre 2m min.",
	104	. "K",
	105	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	106	. t2mincels, zsto,zout)
	107	c
	108	CALL histdef(nid_day, "t2mter_max", "Temp.terre 2m max.",
	109	. "K",
	110	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	111	. t2maxcels, zsto,zout)
	112	c
	113	CALL histdef(nid_day, "q2m", "Specific humidity", "Kg/Kg",
	114	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	115	. "ave(X)", zsto,zout)
	116	c
	117	CALL histdef(nid_day, "u10m", "Vent zonal 10m", "m/s",
	118	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	119	. "ave(X)", zsto,zout)
	120	c
	121	CALL histdef(nid_day, "u10mter", "Vent zonal ter 10m", "m/s",
	122	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	123	. "ave(X)", zsto,zout)
	124	c
	125	CALL histdef(nid_day, "u10mlic", "Vent zonal lic 10m", "m/s",
	126	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	127	. "ave(X)", zsto,zout)
	128	c
	129	CALL histdef(nid_day, "u10moce", "Vent zonal oce 10m", "m/s",
	130	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	131	. "ave(X)", zsto,zout)
	132	c
	133	CALL histdef(nid_day, "u10msic", "Vent zonal sic 10m",
	134	. "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
	135	. "ave(X)", zsto,zout)
	136	c
	137	CALL histdef(nid_day, "v10m", "Vent meridien 10m", "m/s",
	138	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	139	. "ave(X)", zsto,zout)
	140	c
	141	CALL histdef(nid_day, "v10mter", "Vent meridien ter 10m",
	142	. "m/s", iim,jjmp1,nhori, 1,1,1, -99, 32,
	143	. "ave(X)", zsto,zout)
	144	c
	145	CALL histdef(nid_day, "v10mlic", "Vent meridien lic 10m",
	146	. "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
	147	. "ave(X)", zsto,zout)
	148	c
	149	CALL histdef(nid_day, "v10moce", "Vent meridien oce 10m",
	150	. "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
	151	. "ave(X)", zsto,zout)
	152	c
	153	CALL histdef(nid_day, "v10msic", "Vent meridien sic 10m",
	154	. "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
	155	. "ave(X)", zsto,zout)
	156	c
	157	CALL histdef(nid_day, "psol", "Surface Pressure", "Pa",
	158	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	159	. "ave(X)", zsto,zout)
	160	c
	161	CALL histdef(nid_day, "precip","Precipitation Totale liq+sol"
	162	. , "kg/s",
	163	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	164	. "ave(X)", zsto,zout)
	165	c
	166	CALL histdef(nid_day, "snow", "Snow fall", "kg/s",
	167	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	168	. "ave(X)", zsto,zout)
	169	c
	170	CALL histdef(nid_day, "snow_mass", "Snow Mass", "kg/m2",
	171	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	172	. "ave(X)", zsto,zout)
	173	c
	174	CALL histdef(nid_day, "evap", "Evaporation", "kg/s",
	175	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	176	. "ave(X)", zsto,zout)
	177	c
	178	CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
	179	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	180	. "ave(X)", zsto,zout)
	181	c
	182	CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
	183	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	184	. "ave(X)", zsto,zout)
	185	c
	186	CALL histdef(nid_day, "sols", "Net Solar rad. at surf.",
	187	. "W/m2",
	188	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	189	. "ave(X)", zsto,zout)
	190	c
	191	CALL histdef(nid_day, "soll", "Net IR rad. at surface", "W/m2",
	192	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	193	. "ave(X)", zsto,zout)
	194	c
	195	CALL histdef(nid_day, "solldown", "Down. IR rad. at surface",
	196	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
	197	. "ave(X)", zsto,zout)
	198	c
	199	CALL histdef(nid_day, "bils", "Surf. total heat flux", "W/m2",
	200	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	201	. "ave(X)", zsto,zout)
	202	c
	203	CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2",
	204	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	205	. "ave(X)", zsto,zout)
	206	c
	207	CALL histdef(nid_day, "fder", "Heat flux derivation", "W/m2",
	208	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	209	. "ave(X)", zsto,zout)
	210	c
	211	CALL histdef(nid_day, "frtu", "Zonal wind stress", "Pa",
	212	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	213	. "ave(X)", zsto,zout)
	214	c
	215	CALL histdef(nid_day, "frtv", "Meridional wind stress", "Pa",
	216	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	217	. "ave(X)", zsto,zout)
	218	c
	219	CXXX PB flux pour chaque sous surface
	220	C
	221	DO nsrf = 1, nbsrf
	222	C
	223	call histdef(nid_day, "pourc_"//clnsurf(nsrf),
	224	$ "Fraction"//clnsurf(nsrf), "W/m2",
	225	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	226	$ "ave(X)", zsto,zout)
	227	C
	228	call histdef(nid_day, "tsol_"//clnsurf(nsrf),
	229	$ "Fraction"//clnsurf(nsrf), "W/m2",
	230	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	231	$ "ave(X)", zsto,zout)
	232	C
	233	call histdef(nid_day, "sens_"//clnsurf(nsrf),
	234	$ "Sensible heat flux "//clnsurf(nsrf), "W/m2",
	235	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	236	$ "ave(X)", zsto,zout)
	237	c
	238	call histdef(nid_day, "lat_"//clnsurf(nsrf),
	239	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
	240	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	241	$ "ave(X)", zsto,zout)
	242	C
	243	call histdef(nid_day, "taux_"//clnsurf(nsrf),
	244	$ "Zonal wind stress"//clnsurf(nsrf),"Pa",
	245	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	246	$ "ave(X)", zsto,zout)
	247
	248	call histdef(nid_day, "tauy_"//clnsurf(nsrf),
	249	$ "Meridional xind stress "//clnsurf(nsrf), "Pa",
	250	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	251	$ "ave(X)", zsto,zout)
	252	C
	253	call histdef(nid_day, "albe_"//clnsurf(nsrf),
	254	$ "Albedo surf. "//clnsurf(nsrf), "W/m2",
	255	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	256	$ "ave(X)", zsto,zout)
	257	C
	258	call histdef(nid_day, "rugs_"//clnsurf(nsrf),
	259	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
	260	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	261	$ "ave(X)", zsto,zout)
	262
	263	CXXX
	264	END DO
	265
	266	CALL histdef(nid_day, "sicf", "Sea-ice fraction", "-",
	267	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	268	. "ave(X)", zsto,zout)
	269	c
	270	CALL histdef(nid_day, "cldl", "Low-level cloudiness", "-",
	271	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	272	. "ave(X)", zsto,zout)
	273	c
	274	CALL histdef(nid_day, "cldm", "Mid-level cloudiness", "-",
	275	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	276	. "ave(X)", zsto,zout)
	277	c
	278	CALL histdef(nid_day, "cldh", "High-level cloudiness", "-",
	279	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	280	. "ave(X)", zsto,zout)
	281	c
	282	CALL histdef(nid_day, "cldt", "Total cloudiness", "-",
	283	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	284	. "ave(X)", zsto,zout)
	285	c
	286	CALL histdef(nid_day, "cldq", "Cloud liquid water path", "-",
	287	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	288	. "ave(X)", zsto,zout)
	289	c
	290	c Champs 3D:
	291	c
	292	CALL histdef(nid_day, "temp", "Air temperature", "K",
	293	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	294	. "ave(X)", zsto,zout)
	295	c
	296	CALL histdef(nid_day, "ovap", "Specific humidity", "Kg/Kg",
	297	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	298	. "ave(X)", zsto,zout)
	299	c
	300	CALL histdef(nid_day, "geop", "Geopotential height", "m",
	301	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	302	. "ave(X)", zsto,zout)
	303	c
	304	CALL histdef(nid_day, "vitu", "Zonal wind", "m/s",
	305	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	306	. "ave(X)", zsto,zout)
	307	c
	308	CALL histdef(nid_day, "vitv", "Meridional wind", "m/s",
	309	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	310	. "ave(X)", zsto,zout)
	311	c
	312	CALL histdef(nid_day, "vitw", "Vertical wind", "m/s",
	313	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	314	. "ave(X)", zsto,zout)
	315	c
	316	CALL histdef(nid_day, "pres", "Air pressure", "Pa",
	317	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	318	. "ave(X)", zsto,zout)
	319	c
	320	cccIM
	321	CALL histdef(nid_day, "SWupTOA", "SWup at TOA","W/m2",
	322	. iim,jjmp1,nhori, 1,1,1,-99,
	323	. 32, "ave(X)", zsto,zout)
	324	c
	325	CALL histdef(nid_day, "SWupSFC", "SWup at surface","W/m2",
	326	. iim,jjmp1,nhori, 1,1,1,-99,
	327	. 32, "ave(X)", zsto,zout)
	328	c
	329	CALL histdef(nid_day, "SWdnTOA", "SWdn at TOA","W/m2",
	330	. iim,jjmp1,nhori, 1,1,1,-99,
	331	. 32, "ave(X)", zsto,zout)
	332	c
	333	CALL histdef(nid_day, "SWdnSFC", "SWdn at surface","W/m2",
	334	. iim,jjmp1,nhori, 1,1,1,-99,
	335	. 32, "ave(X)", zsto,zout)
	336	c
	337	CALL histend(nid_day)
	338	c
	339	ndex2d = 0
	340	ndex3d = 0
	341	c
	342	ENDIF ! fin de test sur ok_journe

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

Context Navigation

source: LMDZ.3.3/trunk/libf/phylmd/ini_histday.h @ 3670

Download in other formats: