Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Normal
Revision Log

ini_histday.h @ 1309

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

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

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

Context Navigation

source: LMDZ4/tags/LMDZ4_V3_2/libf/phylmd/ini_histday.h @ 1309

Download in other formats: