Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Normal
Revision Log

ini_histday.h @ 1717

Last change on this file since 1717 was 741, checked in by (none), 18 years ago
This commit was manufactured by cvs2svn to create tag 'LMDZ4_V2'.
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 22.5 KB

Rev	Line
[644]	1	c
	2	c $Header$
	3	c
[524]	4	IF (ok_journe) THEN
	5	c
[644]	6	zstophy = dtime
[740]	7	zstoday= ecrit_day
	8	zout = ecrit_day
[524]	9	c
	10	idayref = day_ref
	11	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
	12	c
	13	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
	14	DO i = 1, iim
	15	zx_lon(i,1) = rlon(i+1)
	16	zx_lon(i,jjmp1) = rlon(i+1)
	17	ENDDO
	18	DO ll=1,klev
	19	znivsig(ll)=float(ll)
	20	ENDDO
	21	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
	22	CALL histbeg("histday", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
	23	. 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
	24	. nhori, nid_day)
	25	CALL histvert(nid_day, "presnivs", "Vertical levels", "mb",
[644]	26	. klev, presnivs/100., nvert)
[524]	27	c
	28	IF(lev_histday.GE.1) THEN
	29	c
[644]	30	CALL histdef(nid_day, "phis", "Surface geop. height", "m2/s2",
[524]	31	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	32	. "once", zstophy,zout)
[524]	33	c
	34	CALL histdef(nid_day, "aire", "Grid area", "-",
	35	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	36	. "once", zstophy,zout)
[524]	37	c
	38	CALL histdef(nid_day, "contfracATM","% sfce ter+lic ","-",
	39	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	40	. "once", zstophy,zout)
[524]	41	c
	42	CALL histdef(nid_day, "contfracOR","% sfce terre OR", "-",
	43	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	44	. "once", zstophy,zout)
[524]	45	c
	46	c Champs 2D:
	47	c
	48	CALL histdef(nid_day, "tsol", "Surface Temperature", "K",
	49	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	50	. "ave(X)", zstophy,zout)
[524]	51	c
	52	CALL histdef(nid_day, "t2m", "Temperature 2m", "K",
	53	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	54	. "ave(X)", zstophy,zout)
[524]	55	c
	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",
	61	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	62	. t2mincels, zstophy,zout)
[524]	63	c
	64	CALL histdef(nid_day, "t2m_max", "Temp. 2m max.",
	65	. "K",
	66	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	67	. t2maxcels, zstophy,zout)
[524]	68	c
	69	CALL histdef(nid_day, "plul", "Large-scale Precip.",
	70	. "kg/(s*m2)",
	71	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	72	. "ave(X)", zstophy,zout)
[524]	73	c
	74	CALL histdef(nid_day, "pluc", "Convective Precip.",
	75	. "kg/(s*m2)",
	76	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	77	. "ave(X)", zstophy,zout)
[524]	78	c
	79	CALL histdef(nid_day, "snowl", "Solid Large-scale Precip.",
[644]	80	. "kg/(m2*s)",
[524]	81	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	82	. "ave(X)", zstophy,zout)
[524]	83	c
	84	IF(1.EQ.0) THEN !snowc=0.
	85	CALL histdef(nid_day, "snowc", "Solid Convective Precip.",
[644]	86	. "kg/(m2*s)",
[524]	87	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	88	. "ave(X)", zstophy,zout)
[524]	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,
[644]	93	. "ave(X)", zstophy,zout)
[524]	94	c
	95	CALL histdef(nid_day, "sicf", "Sea-ice fraction", "-",
	96	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	97	. "inst(X)", zstoday,zout)
[524]	98	c
	99	CALL histdef(nid_day, "q2m", "Specific humidity", "kg/kg",
	100	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	101	. "ave(X)", zstophy,zout)
[524]	102	c
	103	CALL histdef(nid_day, "u10m", "Vent zonal 10m", "m/s",
	104	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	105	. "ave(X)", zstophy,zout)
[524]	106	c
	107	CALL histdef(nid_day, "v10m", "Vent meridien 10m", "m/s",
	108	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	109	. "ave(X)", zstophy,zout)
[524]	110	c
	111	CALL histdef(nid_day, "wind10m","10-m wind speed","m/s",
	112	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	113	. "ave(X)", zstophy,zout)
[524]	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,
[644]	118	. t2maxcels, zstophy,zout)
[524]	119	c
	120	CALL histdef(nid_day, "psol", "Surface Pressure", "Pa",
	121	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	122	. "ave(X)", zstophy,zout)
[524]	123	c
	124	CALL histdef(nid_day, "precip","Precipitation Totale liq+sol"
	125	. , "kg/(s*m2)",
	126	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	127	. "ave(X)", zstophy,zout)
[524]	128	c
[644]	129	CALL histdef(nid_day, "snowf", "Snow fall", "kg/(m2*s)",
[524]	130	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	131	. "ave(X)", zstophy,zout)
[524]	132	c
[644]	133	c CALL histdef(nid_day, "snow_mass", "Snow Mass", "kg/m2",
	134	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	135	c . "ave(X)", zstophy,zout)
[524]	136	c
[644]	137	CALL histdef(nid_day, "evap", "Evaporation", "kg/(m2*s)",
[524]	138	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	139	. "ave(X)", zstophy,zout)
[524]	140	c
	141	CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
	142	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	143	. "ave(X)", zstophy,zout)
[524]	144	c
	145	CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
	146	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	147	. "ave(X)", zstophy,zout)
[524]	148	c
	149	CALL histdef(nid_day, "sols", "Net Solar rad. at surf.",
	150	. "W/m2",
	151	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	152	. "ave(X)", zstophy,zout)
[524]	153	c
	154	CALL histdef(nid_day, "soll", "Net IR rad. at surface", "W/m2",
	155	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	156	. "ave(X)", zstophy,zout)
[524]	157	c
	158	CALL histdef(nid_day, "cldl", "Low-level cloudiness", "-",
	159	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	160	. "ave(X)", zstophy,zout)
[524]	161	c
	162	CALL histdef(nid_day, "cldm", "Mid-level cloudiness", "-",
	163	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	164	. "ave(X)", zstophy,zout)
[524]	165	c
	166	CALL histdef(nid_day, "cldh", "High-level cloudiness", "-",
	167	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	168	. "ave(X)", zstophy,zout)
[524]	169	c
	170	CALL histdef(nid_day, "cldt", "Total cloudiness", "-",
	171	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	172	. "ave(X)", zstophy,zout)
[524]	173	c
	174	CALL histdef(nid_day, "cldq", "Cloud liquid water path",
	175	. "kg/m2",
	176	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	177	. "ave(X)", zstophy,zout)
[524]	178	c
	179	CALL histdef(nid_day, "prw", "Precipitable water", "kg/m2",
	180	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	181	. "ave(X)", zstophy,zout)
[524]	182	c
	183	c Champs dynamiques sur niveaux de pression
	184	c
[644]	185	c DO k=1, nlevENS
	186	DO k=1, nlevSTD
	187	c
	188	c bb=clevSTD(k)
	189	c
	190	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
	191	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
	192	c aa=clevSTD(k)
	193	c bb2=clevSTD(k)
	194	c ENDIF
	195	c
	196	IF(bb2.EQ."850".OR.bb2.EQ."700"
	197	. .OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
	198	c
	199	CALL histdef(nid_day, "u"//bb2,
	200	. "Zonal wind "//bb2//"mb","m/s",
[524]	201	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	202	. "inst(X)", zout,zout)
[524]	203	c
[644]	204	CALL histdef(nid_day, "v"//bb2,
	205	. "Meridional wind "//bb2//"mb","m/s",
[524]	206	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	207	. "inst(X)", zout,zout)
[524]	208	c
[644]	209	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
[524]	210	c
[644]	211	IF(bb2.EQ."500") THEN
	212	c
	213	CALL histdef(nid_day, "w"//bb2, "Vertical wind "//bb2//"mb",
	214	. "m/s",
[524]	215	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	216	. "inst(X)", zout,zout)
[524]	217	c
[644]	218	CALL histdef(nid_day,"phi"//bb2, "Geopotential "//bb2//"mb",
	219	. "m",
[524]	220	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	221	. "inst(X)", zout,zout)
[524]	222	c
[644]	223	ENDIF !(bb2.EQ."500")
	224	c
	225	ENDDO !nlevSTD
	226	c
[524]	227	CALL histdef(nid_day, "slp", "Sea Level Pressure", "Pa",
	228	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	229	. "ave(X)", zstophy,zout)
[524]	230	c
	231	CALL histdef(nid_day, "cape_max", "CAPE max.",
	232	. "J/kg",
	233	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	234	. capemaxcels, zstophy,zout)
[524]	235	c
	236	CALL histdef(nid_day, "solldown", "Down. IR rad. at surface",
	237	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	238	. "ave(X)", zstophy,zout)
[524]	239	c
	240	CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2",
	241	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	242	. "ave(X)", zstophy,zout)
[524]	243	c
	244	CALL histdef(nid_day, "SWdnSFC", "SWdn at surface","W/m2",
	245	. iim,jjmp1,nhori, 1,1,1,-99,
[644]	246	. 32, "ave(X)", zstophy,zout)
[524]	247	c
[644]	248	cIM diagnostiques flux ocean-atm ou ocean-glace de mer
	249	cIM pour utilisation dans un modele de "slab" ocean
	250	c
	251	IF (OCEAN.EQ.'force ') THEN
	252	c
	253	CALL histdef(nid_day, "lmt_bils",
	254	. "Bilan au sol atmosphere forcee", "W/m2",
	255	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	256	. "ave(X)", zstophy,zout)
	257	c
	258	ELSE IF (OCEAN.EQ.'slab ') THEN
	259	c
	260	CALL histdef(nid_day, "slab_bils",
	261	. "Bilan au sol slab", "W/m2",
	262	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	263	. "ave(X)", zstophy,zout)
	264	c
	265	CALL histdef(nid_day, "tslab", "Slab SST ", "K",
	266	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	267	. "ave(X)", zstophy,zout)
	268	c
	269	CALL histdef(nid_day, "seaice", "Slab seaice", "kg/m2",
	270	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	271	. "ave(X)", zstophy,zout)
	272	c
	273	CALL histdef(nid_day, "siceh", "Slab seaice height", "m",
	274	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	275	. "ave(X)", zstophy,zout)
	276	c
	277	ENDIF !OCEAN.EQ.force/slab
	278	c
	279	CALL histdef(nid_day, "fluxo",
	280	. "Flux turbulents ocean-atmosphere", "W/m2",
	281	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	282	. "ave(X)", zstophy,zout)
	283	c
	284	CALL histdef(nid_day, "fluxg",
	285	. "Flux ocean-glace de mer", "W/m2",
	286	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	287	. "ave(X)", zstophy,zout)
	288	c
[524]	289	ENDIF !lev_histday.GE.1
	290	c
	291	IF (lev_histday.GE.2) THEN
	292	c
	293	CALL histdef(nid_day, "bils", "Surf. total heat flux", "W/m2",
	294	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	295	. "ave(X)", zstophy,zout)
[524]	296	c
	297	CALL histdef(nid_day, "fder", "Heat flux derivation", "W/m2/K",
	298	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	299	. "ave(X)", zstophy,zout)
[524]	300	c
	301	ENDIF !lev_histday.GE.2
	302	c=================================================================
	303	IF(lev_histday.GE.3) THEN
	304	c=================================================================
	305	c INITIALISATION DES CHAMPS 3D
	306	c=================================================================
	307	c=================================================================
	308	c=================================================================
	309	c Champs 3D:
	310	c
	311	CALL histdef(nid_day, "temp", "Air temperature", "K",
	312	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]	313	. "ave(X)", zstophy,zout)
[524]	314	c
	315	CALL histdef(nid_day, "ovap", "Specific humidity", "kg/kg",
	316	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]	317	. "ave(X)", zstophy,zout)
[524]	318	c
[644]	319	CALL histdef(nid_day, "geop", "Geopotential height", "m2/s2",
[524]	320	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]	321	. "ave(X)", zstophy,zout)
[524]	322	c
	323	CALL histdef(nid_day, "vitu", "Zonal wind", "m/s",
	324	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]	325	. "ave(X)", zstophy,zout)
[524]	326	c
	327	CALL histdef(nid_day, "vitv", "Meridional wind", "m/s",
	328	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]	329	. "ave(X)", zstophy,zout)
[524]	330	c
[644]	331	CALL histdef(nid_day, "vitw", "Vertical wind", "m/s",
[524]	332	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]	333	. "ave(X)", zstophy,zout)
[524]	334	c
	335	CALL histdef(nid_day, "pres", "Air pressure", "Pa",
	336	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
[644]	337	. "ave(X)", zstophy,zout)
[524]	338	c
	339	c=================================================================
	340	c FIN INITIALISATION DES CHAMPS 3D
	341	c=================================================================
	342	ENDIF ! lev_histday.GE.3
	343	c
	344	IF (lev_histday.GE.4) THEN
	345	c=================================================================
	346	c
	347	c INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
	348	c
	349	c=================================================================
	350	c
	351	CALL histdef(nid_day, "SWupTOA", "SWup at TOA","W/m2",
	352	. iim,jjmp1,nhori, 1,1,1,-99,
[644]	353	. 32, "ave(X)", zstophy,zout)
[524]	354	c
	355	CALL histdef(nid_day, "SWupSFC", "SWup at surface","W/m2",
	356	. iim,jjmp1,nhori, 1,1,1,-99,
[644]	357	. 32, "ave(X)", zstophy,zout)
[524]	358	c
	359	CALL histdef(nid_day, "SWdnTOA", "SWdn at TOA","W/m2",
	360	. iim,jjmp1,nhori, 1,1,1,-99,
[644]	361	. 32, "ave(X)", zstophy,zout)
[524]	362	c
	363	CALL histdef(nid_day, "SWupTOAclr",
	364	. "SWup clear sky at TOA","W/m2",
	365	. iim,jjmp1,nhori, 1,1,1,-99,
[644]	366	. 32, "ave(X)", zstophy,zout)
[524]	367
	368	CALL histdef(nid_day, "SWupSFCclr",
	369	. "SWup clear sky at surface","W/m2",
	370	. iim,jjmp1,nhori, 1,1,1,-99,
[644]	371	. 32, "ave(X)", zstophy,zout)
[524]	372
	373	CALL histdef(nid_day, "SWdnTOAclr",
	374	. "SWdn clear sky at TOA","W/m2",
	375	. iim,jjmp1,nhori, 1,1,1,-99,
[644]	376	. 32, "ave(X)", zstophy,zout)
[524]	377
	378	CALL histdef(nid_day, "SWdnSFCclr",
	379	. "SWdn clear sky at surface","W/m2",
	380	. iim,jjmp1,nhori, 1,1,1,-99,
[644]	381	. 32, "ave(X)", zstophy,zout)
[524]	382	c
	383	c=================================================================
	384	c=================================================================
	385	c=================================================================
	386	c INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
	387	c=================================================================
	388	c
[644]	389	CALL histdef(nid_day, "tter", "Surface Temperature ter", "K",
[524]	390	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	391	. "ave(X)", zstophy,zout)
[524]	392	c
[644]	393	CALL histdef(nid_day, "tlic", "Surface Temperature lic", "K",
[524]	394	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	395	. "ave(X)", zstophy,zout)
[524]	396	c
[644]	397	CALL histdef(nid_day, "toce", "Surface Temperature oce", "K",
[524]	398	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	399	. "ave(X)", zstophy,zout)
[524]	400	c
[644]	401	CALL histdef(nid_day, "tsic", "Surface Temperature sic", "K",
[524]	402	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	403	. "ave(X)", zstophy,zout)
[524]	404	c
[644]	405	CALL histdef(nid_day, "t2mter", "Temp.tere 2m", "K",
[524]	406	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	407	. "ave(X)", zstophy,zout)
[524]	408	c
	409	CALL histdef(nid_day, "t2mlic", "Temp.lic 2m", "K",
	410	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	411	. "ave(X)", zstophy,zout)
[524]	412	c
	413	CALL histdef(nid_day, "t2moce", "Temp.oce 2m", "K",
	414	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	415	. "ave(X)", zstophy,zout)
[524]	416	c
	417	CALL histdef(nid_day, "t2msic", "Temp.sic 2m", "K",
	418	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	419	. "ave(X)", zstophy,zout)
[524]	420	c
	421	CALL histdef(nid_day, "t2mter_min", "Temp.terre 2m min.",
	422	. "K",
	423	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	424	. t2mincels, zstophy,zout)
[524]	425	c
	426	CALL histdef(nid_day, "t2mter_max", "Temp.terre 2m max.",
	427	. "K",
	428	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	429	. t2maxcels, zstophy,zout)
[524]	430	c
	431	CALL histdef(nid_day, "u10mter", "Vent zonal ter 10m", "m/s",
	432	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	433	. "ave(X)", zstophy,zout)
[524]	434	c
	435	CALL histdef(nid_day, "u10mlic", "Vent zonal lic 10m", "m/s",
	436	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	437	. "ave(X)", zstophy,zout)
[524]	438	c
	439	CALL histdef(nid_day, "u10moce", "Vent zonal oce 10m", "m/s",
	440	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	441	. "ave(X)", zstophy,zout)
[524]	442	c
	443	CALL histdef(nid_day, "u10msic", "Vent zonal sic 10m",
	444	. "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	445	. "ave(X)", zstophy,zout)
[524]	446	c
	447	CALL histdef(nid_day, "v10mter", "Vent meridien ter 10m",
	448	. "m/s", iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	449	. "ave(X)", zstophy,zout)
[524]	450	c
	451	CALL histdef(nid_day, "v10mlic", "Vent meridien lic 10m",
	452	. "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	453	. "ave(X)", zstophy,zout)
[524]	454	c
	455	CALL histdef(nid_day, "v10moce", "Vent meridien oce 10m",
	456	. "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	457	. "ave(X)", zstophy,zout)
[524]	458	c
	459	CALL histdef(nid_day, "v10msic", "Vent meridien sic 10m",
	460	. "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	461	. "ave(X)", zstophy,zout)
[524]	462	c
	463	DO nsrf = 1, nbsrf
	464	C
	465	call histdef(nid_day, "pourc_"//clnsurf(nsrf),
[644]	466	$ "% "//clnsurf(nsrf), "%",
[524]	467	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	468	$ "ave(X)", zstophy,zout)
[524]	469	C
[644]	470	call histdef(nid_day, "fract_"//clnsurf(nsrf),
	471	$ "Fraction "//clnsurf(nsrf), "1",
	472	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	473	$ "ave(X)", zstophy,zout)
	474	C
[524]	475	call histdef(nid_day, "tsol_"//clnsurf(nsrf),
[644]	476	$ "Temperature "//clnsurf(nsrf), "K",
[524]	477	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	478	$ "ave(X)", zstophy,zout)
[524]	479	C
	480	call histdef(nid_day, "sens_"//clnsurf(nsrf),
	481	$ "Sensible heat flux "//clnsurf(nsrf), "W/m2",
	482	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	483	$ "ave(X)", zstophy,zout)
[524]	484	c
	485	call histdef(nid_day, "lat_"//clnsurf(nsrf),
	486	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
	487	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	488	$ "ave(X)", zstophy,zout)
[524]	489	C
	490	call histdef(nid_day, "taux_"//clnsurf(nsrf),
[644]	491	$ "Zonal wind stress "//clnsurf(nsrf),"Pa",
[524]	492	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	493	$ "ave(X)", zstophy,zout)
[524]	494
	495	call histdef(nid_day, "tauy_"//clnsurf(nsrf),
[644]	496	$ "Meridional wind stress "//clnsurf(nsrf), "Pa",
[524]	497	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	498	$ "ave(X)", zstophy,zout)
[524]	499	C
	500	call histdef(nid_day, "albe_"//clnsurf(nsrf),
[644]	501	$ "Albedo surf. "//clnsurf(nsrf), "-",
[524]	502	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	503	$ "ave(X)", zstophy,zout)
[524]	504	C
	505	call histdef(nid_day, "rugs_"//clnsurf(nsrf),
[644]	506	$ "Rugosite "//clnsurf(nsrf), "m",
[524]	507	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	508	$ "ave(X)", zstophy,zout)
[524]	509	C
	510	END DO
	511	C
	512	c=================================================================
	513	c
	514	c FIN INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
	515	c
	516	c=================================================================
	517	ENDIF !lev_histday.GE.4
[644]	518	c
	519	IF(lev_histday.GE.5) THEN !lev_histday.GE.5
	520	c
	521	c rajout sorties F. Aires
	522	CALL histdef(nid_day, "tops0", "CS Solar rad. at TOA",
	523	. "W/m2",
	524	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	525	. "ave(X)", zstophy,zout)
	526	c
	527	CALL histdef(nid_day, "topl0", "CS IR rad. at TOA",
	528	. "W/m2",
	529	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	530	. "ave(X)", zstophy,zout)
	531	c
	532	CALL histdef(nid_day, "sols0", "CS Net Solar rad. at surf.",
	533	. "W/m2",
	534	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	535	. "ave(X)", zstophy,zout)
	536	c
	537	CALL histdef(nid_day, "soll0", "CS Net IR rad. at surface",
	538	. "W/m2",
	539	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	540	. "ave(X)", zstophy,zout)
	541	c
	542	CALL histdef(nid_day, "rneb", "Cloud fraction", "-",
	543	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	544	. "ave(X)", zstophy,zout)
	545	c
	546	CALL histdef(nid_day, "rnebcon", "Convective Cloud Fraction"
	547	. , "-",
	548	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	549	. "ave(X)", zstophy,zout)
	550	c
	551	CALL histdef(nid_day,"lwcon","Cloud liquid water content",
	552	. "kg/kg",
	553	. iim,jjmp1,nhori, klev,1,klev, nvert, 32,
	554	. "ave(X)", zstophy,zout)
	555	c
	556	CALL histdef(nid_day,"iwcon","Cloud ice water content","kg/kg",
	557	. iim,jjmp1,nhori, klev,1,klev, nvert, 32,
	558	. "ave(X)", zstophy,zout)
	559	c
	560	CALL histdef(nid_day,"lwp","Cloud water path","kg/m2",
	561	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	562	. "ave(X)", zstophy,zout)
	563	c
	564	CALL histdef(nid_day,"iwp","Cloud ice water path","kg/m2",
	565	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	566	. "ave(X)", zstophy,zout)
	567	c
	568	CALL histdef(nid_day,"meantaucld",
	569	. "ISCCP mean cloud optical thickness","1",
	570	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	571	. "ave(X)", zstophy,zout)
	572	c
	573	CALL histdef(nid_day,"cldtau","Cloud optical thickness","1",
	574	. iim,jjmp1,nhori, klev,1,klev, nvert, 32,
	575	. "ave(X)", zstophy,zout)
	576	c
	577	CALL histdef(nid_day,"cldemi","Cloud optical emissivity","1",
	578	. iim,jjmp1,nhori, klev,1,klev, nvert, 32,
	579	. "ave(X)", zstophy,zout)
	580	c
	581	ENDIF !lev_histday.GE.5
[524]	582	c=================================================================
	583	c
	584	CALL histend(nid_day)
	585	c
	586	ndex2d = 0
	587	ndex3d = 0
	588	c
	589	c=================================================================
	590	ENDIF ! fin de test sur ok_journe

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

Context Navigation

source: LMDZ4/tags/LMDZ4_V2/libf/phylmd/ini_histday.h @ 1717

Download in other formats: