Context Navigation

source: LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histday.h @ 493

Last change on this file since 493 was 486, checked in by lmdzadmin, 22 years ago
Phasage avec la version de Ionela IM/LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 17.1 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	CALL histbeg("histday", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
	16	. 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
	17	. nhori, nid_day)
	18	write(,)'Journee ', itau_phy, zjulian
	19	CALL histvert(nid_day, "presnivs", "Vertical levels", "mb",
	20	. klev, presnivs, nvert)
	21	c call histvert(nid_day, 'sig_s', 'Niveaux sigma','-',
	22	c . klev, znivsig, nvert)
	23	c
	24	zsto = dtime
	25	zout = dtime * FLOAT(ecrit_day)
	26	C Essai writephys
	27	c nom_fichier = 'histday1'
	28	c call writephy_ini(fid_day,nom_fichier,klon,iim,jjmp1,klev,
	29	c . rlon,rlat, presnivs,
	30	c . zjulian, dtime)
	31	c call writephy_def(prof2d_on, fid_day, "once", zsto, zout, 0)
	32	c call writephy_def(prof3d_on, fid_day, "once", zsto, zout,
	33	c . klev)
	34	c call writephy_def(prof2d_av, fid_day, "ave(X)", zsto, zout, 0)
	35	c call writephy_def(prof3d_av, fid_day, "ave(X)", zsto, zout,
	36	c . klev)
	37
	38	c
	39	CALL histdef(nid_day, "phis", "Surface geop. height", "-",
	40	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	41	. "once", zsto,zout)
	42	c
	43	CALL histdef(nid_day, "aire", "Grid area", "-",
	44	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	45	. "once", zsto,zout)
	46	c
	47	c Champs 2D:
	48	c
	49	CALL histdef(nid_day, "tsol", "Surface Temperature", "K",
	50	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	51	. "ave(X)", zsto,zout)
	52	c
[433]	53	c
	54	CALL histdef(nid_day, "t2m", "Temperature 2m", "K",
	55	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	56	. "ave(X)", zsto,zout)
	57	c
	58	CALL histdef(nid_day, "q2m", "Specific humidity", "Kg/Kg",
	59	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	60	. "ave(X)", zsto,zout)
	61	c
	62	CALL histdef(nid_day, "u10m", "Vent zonal 10m", "m/s",
	63	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	64	. "ave(X)", zsto,zout)
	65	c
	66	CALL histdef(nid_day, "v10m", "Vent meridien 10m", "m/s",
	67	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	68	. "ave(X)", zsto,zout)
	69	c
	70	CALL histdef(nid_day, "psol", "Surface Pressure", "Pa",
[416]	71	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	72	. "ave(X)", zsto,zout)
	73	c
[433]	74	CALL histdef(nid_day, "precip","Precipitation Totale liq+sol"
[451]	75	. , "kg/(s*m2)",
[416]	76	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	77	. "ave(X)", zsto,zout)
	78	c
[451]	79	CALL histdef(nid_day, "snow", "Snow fall", "kg/(s*m2)",
[416]	80	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	81	. "ave(X)", zsto,zout)
	82	c
[433]	83	CALL histdef(nid_day, "snow_mass", "Snow Mass", "kg/m2",
[416]	84	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	85	. "ave(X)", zsto,zout)
	86	c
[451]	87	CALL histdef(nid_day, "evap", "Evaporation", "kg/(s*m2)",
[433]	88	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	89	. "ave(X)", zsto,zout)
[416]	90	c
[433]	91	CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
	92	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	93	. "ave(X)", zsto,zout)
	94	c
[433]	95	CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
	96	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	97	. "ave(X)", zsto,zout)
	98	c
[433]	99	CALL histdef(nid_day, "sols", "Net Solar rad. at surf.",
	100	. "W/m2",
	101	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	102	. "ave(X)", zsto,zout)
	103	c
[433]	104	CALL histdef(nid_day, "soll", "Net IR rad. at surface", "W/m2",
	105	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	106	. "ave(X)", zsto,zout)
	107	c
[433]	108	CALL histdef(nid_day, "solldown", "Down. IR rad. at surface",
	109	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	110	. "ave(X)", zsto,zout)
	111	c
[433]	112	CALL histdef(nid_day, "bils", "Surf. total heat flux", "W/m2",
	113	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	114	. "ave(X)", zsto,zout)
[416]	115	c
[433]	116	CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2",
	117	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	118	. "ave(X)", zsto,zout)
[416]	119	c
[433]	120	CALL histdef(nid_day, "fder", "Heat flux derivation", "W/m2",
	121	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	122	. "ave(X)", zsto,zout)
[416]	123	c
[433]	124	c CALL histdef(nid_day, "frtu", "Zonal wind stress", "Pa",
	125	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	126	c . "ave(X)", zsto,zout)
[416]	127	c
[433]	128	c CALL histdef(nid_day, "frtv", "Meridional wind stress", "Pa",
	129	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	130	c . "ave(X)", zsto,zout)
[416]	131	c
[433]	132	c CALL histdef(nid_day, "sicf", "Sea-ice fraction", "-",
	133	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	134	c . "ave(X)", zsto,zout)
	135	c
	136	CALL histdef(nid_day, "cldl", "Low-level cloudiness", "-",
	137	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	138	. "ave(X)", zsto,zout)
	139	c
[433]	140	CALL histdef(nid_day, "cldm", "Mid-level cloudiness", "-",
	141	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	142	. "ave(X)", zsto,zout)
	143	c
[433]	144	CALL histdef(nid_day, "cldh", "High-level cloudiness", "-",
	145	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	146	. "ave(X)", zsto,zout)
	147	c
[433]	148	CALL histdef(nid_day, "cldt", "Total cloudiness", "-",
	149	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	150	. "ave(X)", zsto,zout)
	151	c
[433]	152	CALL histdef(nid_day, "cldq", "Cloud liquid water path", "-",
	153	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	154	. "ave(X)", zsto,zout)
	155	c
[433]	156	CALL histdef(nid_day, "SWupTOA", "SWup at TOA","W/m2",
	157	. iim,jjmp1,nhori, 1,1,1,-99,
	158	. 32, "ave(X)", zsto,zout)
	159	c
	160	CALL histdef(nid_day, "SWupSFC", "SWup at surface","W/m2",
	161	. iim,jjmp1,nhori, 1,1,1,-99,
	162	. 32, "ave(X)", zsto,zout)
	163	c
	164	CALL histdef(nid_day, "SWdnTOA", "SWdn at TOA","W/m2",
	165	. iim,jjmp1,nhori, 1,1,1,-99,
	166	. 32, "ave(X)", zsto,zout)
	167	c
	168	CALL histdef(nid_day, "SWdnSFC", "SWdn at surface","W/m2",
	169	. iim,jjmp1,nhori, 1,1,1,-99,
	170	. 32, "ave(X)", zsto,zout)
	171	c
[467]	172	CALL histdef(nid_day, "SWupTOAclr",
	173	. "SWup clear sky at TOA","W/m2",
	174	. iim,jjmp1,nhori, 1,1,1,-99,
	175	. 32, "ave(X)", zsto,zout)
	176
	177	CALL histdef(nid_day, "SWupSFCclr",
	178	. "SWup clear sky at surface","W/m2",
	179	. iim,jjmp1,nhori, 1,1,1,-99,
	180	. 32, "ave(X)", zsto,zout)
	181
	182	CALL histdef(nid_day, "SWdnTOAclr",
	183	. "SWdn clear sky at TOA","W/m2",
	184	. iim,jjmp1,nhori, 1,1,1,-99,
	185	. 32, "ave(X)", zsto,zout)
	186
	187	CALL histdef(nid_day, "SWdnSFCclr",
	188	. "SWdn clear sky at surface","W/m2",
	189	. iim,jjmp1,nhori, 1,1,1,-99,
	190	. 32, "ave(X)", zsto,zout)
	191	cccIM
	192	CALL histdef(nid_day, "prw", "Precipitable water", "kg/m2",
	193	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	194	. "ave(X)", zsto,zout)
	195	c
[433]	196	c Champs dynamiques sur niveaux de pression
	197
	198	CALL histdef(nid_day, "u850", "Zonal wind 850mb", "m/s",
[416]	199	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	200	. "ave(X)", zsto,zout)
[433]	201
	202	CALL histdef(nid_day, "v850", "Meridional wind 850mb", "m/s",
	203	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	204	. "ave(X)", zsto,zout)
	205	c
[433]	206	CALL histdef(nid_day, "u500", "Zonal wind 500mb", "m/s",
	207	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	208	. "ave(X)", zsto,zout)
[433]	209
	210	CALL histdef(nid_day, "v500", "Meridional wind 500mb", "m/s",
	211	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	212	. "ave(X)", zsto,zout)
[433]	213
	214	CALL histdef(nid_day, "u200", "Zonal wind 200mb", "m/s",
	215	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	216	. "ave(X)", zsto,zout)
[433]	217
	218	CALL histdef(nid_day, "v200", "Meridional wind 200mb", "m/s",
	219	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	220	. "ave(X)", zsto,zout)
	221
	222	CALL histdef(nid_day,"phi500", "Geopotentiel à 500mb", "m2/s2",
	223	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	224	. "ave(X)", zsto,zout)
	225
	226	CALL histdef(nid_day, "slp", "Sea Level Pressure", "Pa",
	227	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	228	. "ave(X)", zsto,zout)
	229
[447]	230	CALL histdef(nid_day, "cape_max", "CAPE max.",
	231	. "J/Kg",
	232	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	233	. capemaxcels, zsto,zout)
[433]	234
[447]	235
[433]	236	c=================================================================
	237	c CI DESSOUS, SORTIES A UTILISER EN DEBUG
	238	c=================================================================
	239	if ( lev_histday.gt.1 ) then
	240
	241	c Champs retires momentannéement en attendant un hypothetique
	242	c debugage
	243
	244	CALL histdef(nid_day, "t2m_min", "Temp. 2m min.",
	245	. "K",
	246	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	247	. t2mincels, zsto,zout)
[416]	248	c
[433]	249	CALL histdef(nid_day, "t2m_max", "Temp. 2m max.",
	250	. "K",
	251	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	252	. t2maxcels, zsto,zout)
	253	c
	254	c=================================================================
	255	c=================================================================
	256	c=================================================================
	257	c INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
	258	c=================================================================
	259
	260	CALL histdef(nid_day, "tter", "Surface Temperature", "K",
[416]	261	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	262	. "ave(X)", zsto,zout)
	263	c
[433]	264	CALL histdef(nid_day, "tlic", "Surface Temperature", "K",
[416]	265	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	266	. "ave(X)", zsto,zout)
	267	c
[433]	268	CALL histdef(nid_day, "toce", "Surface Temperature", "K",
[416]	269	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	270	. "ave(X)", zsto,zout)
	271	c
[433]	272	CALL histdef(nid_day, "tsic", "Surface Temperature", "K",
[416]	273	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	274	. "ave(X)", zsto,zout)
	275	c
[433]	276	cccIM
	277	c
	278	CALL histdef(nid_day, "t2mter", "Temp.terre 2m", "K",
	279	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	280	. "ave(X)", zsto,zout)
	281	c
[433]	282	CALL histdef(nid_day, "t2mlic", "Temp.lic 2m", "K",
	283	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	284	. "ave(X)", zsto,zout)
	285	c
[433]	286	CALL histdef(nid_day, "t2moce", "Temp.oce 2m", "K",
	287	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	288	. "ave(X)", zsto,zout)
	289	c
[433]	290	CALL histdef(nid_day, "t2msic", "Temp.sic 2m", "K",
	291	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	292	. "ave(X)", zsto,zout)
	293	c
[433]	294	CALL histdef(nid_day, "t2mter_min", "Temp.terre 2m min.",
	295	. "K",
	296	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	297	. t2mincels, zsto,zout)
	298	c
	299	CALL histdef(nid_day, "t2mter_max", "Temp.terre 2m max.",
	300	. "K",
	301	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	302	. t2maxcels, zsto,zout)
	303	c
	304	CALL histdef(nid_day, "u10mter", "Vent zonal ter 10m", "m/s",
	305	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	306	. "ave(X)", zsto,zout)
	307	c
[433]	308	CALL histdef(nid_day, "u10mlic", "Vent zonal lic 10m", "m/s",
	309	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	310	. "ave(X)", zsto,zout)
	311	c
[433]	312	CALL histdef(nid_day, "u10moce", "Vent zonal oce 10m", "m/s",
	313	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	314	. "ave(X)", zsto,zout)
	315	c
[433]	316	CALL histdef(nid_day, "u10msic", "Vent zonal sic 10m",
	317	. "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	318	. "ave(X)", zsto,zout)
	319	c
[433]	320	CALL histdef(nid_day, "v10mter", "Vent meridien ter 10m",
	321	. "m/s", iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	322	. "ave(X)", zsto,zout)
	323	c
[433]	324	CALL histdef(nid_day, "v10mlic", "Vent meridien lic 10m",
	325	. "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	326	. "ave(X)", zsto,zout)
	327	c
[433]	328	CALL histdef(nid_day, "v10moce", "Vent meridien oce 10m",
	329	. "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	330	. "ave(X)", zsto,zout)
	331	c
[433]	332	CALL histdef(nid_day, "v10msic", "Vent meridien sic 10m",
	333	. "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
	334	. "ave(X)", zsto,zout)
	335	c
	336	c
[416]	337	CXXX PB flux pour chaque sous surface
	338	C
	339	DO nsrf = 1, nbsrf
	340	C
	341	call histdef(nid_day, "pourc_"//clnsurf(nsrf),
	342	$ "Fraction"//clnsurf(nsrf), "W/m2",
	343	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	344	$ "ave(X)", zsto,zout)
	345	C
	346	call histdef(nid_day, "tsol_"//clnsurf(nsrf),
	347	$ "Fraction"//clnsurf(nsrf), "W/m2",
	348	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	349	$ "ave(X)", zsto,zout)
	350	C
	351	call histdef(nid_day, "sens_"//clnsurf(nsrf),
	352	$ "Sensible heat flux "//clnsurf(nsrf), "W/m2",
	353	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	354	$ "ave(X)", zsto,zout)
	355	c
	356	call histdef(nid_day, "lat_"//clnsurf(nsrf),
	357	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
	358	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	359	$ "ave(X)", zsto,zout)
	360	C
	361	call histdef(nid_day, "taux_"//clnsurf(nsrf),
	362	$ "Zonal wind stress"//clnsurf(nsrf),"Pa",
	363	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	364	$ "ave(X)", zsto,zout)
	365
	366	call histdef(nid_day, "tauy_"//clnsurf(nsrf),
	367	$ "Meridional xind stress "//clnsurf(nsrf), "Pa",
	368	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	369	$ "ave(X)", zsto,zout)
	370	C
	371	call histdef(nid_day, "albe_"//clnsurf(nsrf),
	372	$ "Albedo surf. "//clnsurf(nsrf), "W/m2",
	373	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	374	$ "ave(X)", zsto,zout)
	375	C
	376	call histdef(nid_day, "rugs_"//clnsurf(nsrf),
	377	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
	378	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	379	$ "ave(X)", zsto,zout)
	380
	381	CXXX
	382	END DO
	383
[433]	384	c=================================================================
	385	c FIN DES CHAMPS SUR LES SOUS SURFACES
	386	c=================================================================
	387	c=================================================================
	388	c FIN INITIALISATION DES CHAMPS 3D
	389	c=================================================================
	390
[416]	391	c Champs 3D:
	392	c
	393	CALL histdef(nid_day, "temp", "Air temperature", "K",
	394	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	395	. "ave(X)", zsto,zout)
	396	c
	397	CALL histdef(nid_day, "ovap", "Specific humidity", "Kg/Kg",
	398	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	399	. "ave(X)", zsto,zout)
	400	c
	401	CALL histdef(nid_day, "geop", "Geopotential height", "m",
	402	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	403	. "ave(X)", zsto,zout)
	404	c
	405	CALL histdef(nid_day, "vitu", "Zonal wind", "m/s",
	406	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	407	. "ave(X)", zsto,zout)
	408	c
	409	CALL histdef(nid_day, "vitv", "Meridional wind", "m/s",
	410	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	411	. "ave(X)", zsto,zout)
	412	c
	413	CALL histdef(nid_day, "vitw", "Vertical wind", "m/s",
	414	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	415	. "ave(X)", zsto,zout)
	416	c
	417	CALL histdef(nid_day, "pres", "Air pressure", "Pa",
	418	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	419	. "ave(X)", zsto,zout)
	420	c
	421	c
[433]	422	c=================================================================
	423	c FIN INITIALISATION DES CHAMPS 3D
	424	c=================================================================
	425
	426
	427	c=================================================================
	428	c=================================================================
	429	endif ! lev_histday.gt.1
	430	c=================================================================
	431
	432	CALL histend(nid_day)
	433	c
[416]	434	ndex2d = 0
	435	ndex3d = 0
	436	c
[433]	437	c=================================================================
[416]	438	ENDIF ! fin de test sur ok_journe

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

Download in other formats: