Context Navigation

ini_histday.h @ 1006

Last change on this file since 1006 was 574, checked in by Laurent Fairhead, 20 years ago
Erreur sur l'unite de la vitesse verticale LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 17.7 KB

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

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

Context Navigation

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

Download in other formats: