Context Navigation

source: LMDZ4/trunk/libf/phylmd/ini_histday.h @ 844

Last change on this file since 844 was 776, checked in by Laurent Fairhead, 18 years ago
Suite du merge entre la version et la HEAD: quelques modifications de Yann sur le LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 24.7 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
	194	c Champs dynamiques sur niveaux de pression
	195	c
[644]	196	c DO k=1, nlevENS
	197	DO k=1, nlevSTD
	198	c
	199	c bb=clevSTD(k)
	200	c
	201	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
	202	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
	203	c aa=clevSTD(k)
	204	c bb2=clevSTD(k)
	205	c ENDIF
	206	c
	207	IF(bb2.EQ."850".OR.bb2.EQ."700"
	208	. .OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
	209	c
	210	CALL histdef(nid_day, "u"//bb2,
	211	. "Zonal wind "//bb2//"mb","m/s",
[776]	212	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	213	. "inst(X)", zout,zout)
[524]	214	c
[644]	215	CALL histdef(nid_day, "v"//bb2,
	216	. "Meridional wind "//bb2//"mb","m/s",
[776]	217	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	218	. "inst(X)", zout,zout)
[524]	219	c
[644]	220	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
[524]	221	c
[644]	222	IF(bb2.EQ."500") THEN
	223	c
	224	CALL histdef(nid_day, "w"//bb2, "Vertical wind "//bb2//"mb",
	225	. "m/s",
[776]	226	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	227	. "inst(X)", zout,zout)
[524]	228	c
[644]	229	CALL histdef(nid_day,"phi"//bb2, "Geopotential "//bb2//"mb",
	230	. "m",
[776]	231	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	232	. "inst(X)", zout,zout)
[524]	233	c
[644]	234	ENDIF !(bb2.EQ."500")
	235	c
	236	ENDDO !nlevSTD
	237	c
[524]	238	CALL histdef(nid_day, "slp", "Sea Level Pressure", "Pa",
[776]	239	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	240	. "ave(X)", zstophy,zout)
[524]	241	c
	242	CALL histdef(nid_day, "cape_max", "CAPE max.",
	243	. "J/kg",
[776]	244	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	245	. capemaxcels, zstophy,zout)
[524]	246	c
	247	CALL histdef(nid_day, "solldown", "Down. IR rad. at surface",
[776]	248	. "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	249	. "ave(X)", zstophy,zout)
[524]	250	c
	251	CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2",
[776]	252	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	253	. "ave(X)", zstophy,zout)
[524]	254	c
	255	CALL histdef(nid_day, "SWdnSFC", "SWdn at surface","W/m2",
[776]	256	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	257	. 32, "ave(X)", zstophy,zout)
[524]	258	c
[644]	259	cIM diagnostiques flux ocean-atm ou ocean-glace de mer
	260	cIM pour utilisation dans un modele de "slab" ocean
	261	c
	262	IF (OCEAN.EQ.'force ') THEN
	263	c
	264	CALL histdef(nid_day, "lmt_bils",
	265	. "Bilan au sol atmosphere forcee", "W/m2",
[776]	266	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	267	. "ave(X)", zstophy,zout)
	268	c
	269	ELSE IF (OCEAN.EQ.'slab ') THEN
	270	c
	271	CALL histdef(nid_day, "slab_bils",
	272	. "Bilan au sol slab", "W/m2",
[776]	273	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	274	. "ave(X)", zstophy,zout)
	275	c
	276	CALL histdef(nid_day, "tslab", "Slab SST ", "K",
[776]	277	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	278	. "ave(X)", zstophy,zout)
	279	c
	280	CALL histdef(nid_day, "seaice", "Slab seaice", "kg/m2",
[776]	281	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	282	. "ave(X)", zstophy,zout)
	283	c
	284	CALL histdef(nid_day, "siceh", "Slab seaice height", "m",
[776]	285	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	286	. "ave(X)", zstophy,zout)
	287	c
	288	ENDIF !OCEAN.EQ.force/slab
	289	c
	290	CALL histdef(nid_day, "fluxo",
	291	. "Flux turbulents ocean-atmosphere", "W/m2",
[776]	292	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	293	. "ave(X)", zstophy,zout)
	294	c
	295	CALL histdef(nid_day, "fluxg",
	296	. "Flux ocean-glace de mer", "W/m2",
[776]	297	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	298	. "ave(X)", zstophy,zout)
	299	c
[524]	300	ENDIF !lev_histday.GE.1
	301	c
	302	IF (lev_histday.GE.2) THEN
	303	c
	304	CALL histdef(nid_day, "bils", "Surf. total heat flux", "W/m2",
[776]	305	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	306	. "ave(X)", zstophy,zout)
[524]	307	c
	308	CALL histdef(nid_day, "fder", "Heat flux derivation", "W/m2/K",
[776]	309	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	310	. "ave(X)", zstophy,zout)
[524]	311	c
	312	ENDIF !lev_histday.GE.2
	313	c=================================================================
	314	IF(lev_histday.GE.3) THEN
	315	c=================================================================
	316	c INITIALISATION DES CHAMPS 3D
	317	c=================================================================
	318	c=================================================================
	319	c=================================================================
	320	c Champs 3D:
	321	c
	322	CALL histdef(nid_day, "temp", "Air temperature", "K",
[776]	323	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	324	. "ave(X)", zstophy,zout)
[524]	325	c
	326	CALL histdef(nid_day, "ovap", "Specific humidity", "kg/kg",
[776]	327	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	328	. "ave(X)", zstophy,zout)
[524]	329	c
[644]	330	CALL histdef(nid_day, "geop", "Geopotential height", "m2/s2",
[776]	331	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	332	. "ave(X)", zstophy,zout)
[524]	333	c
	334	CALL histdef(nid_day, "vitu", "Zonal wind", "m/s",
[776]	335	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	336	. "ave(X)", zstophy,zout)
[524]	337	c
	338	CALL histdef(nid_day, "vitv", "Meridional wind", "m/s",
[776]	339	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	340	. "ave(X)", zstophy,zout)
[524]	341	c
[644]	342	CALL histdef(nid_day, "vitw", "Vertical wind", "m/s",
[776]	343	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	344	. "ave(X)", zstophy,zout)
[524]	345	c
	346	CALL histdef(nid_day, "pres", "Air pressure", "Pa",
[776]	347	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	348	. "ave(X)", zstophy,zout)
[524]	349	c
	350	c=================================================================
	351	c FIN INITIALISATION DES CHAMPS 3D
	352	c=================================================================
	353	ENDIF ! lev_histday.GE.3
	354	c
	355	IF (lev_histday.GE.4) THEN
	356	c=================================================================
	357	c
	358	c INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
	359	c
	360	c=================================================================
	361	c
	362	CALL histdef(nid_day, "SWupTOA", "SWup at TOA","W/m2",
[776]	363	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	364	. 32, "ave(X)", zstophy,zout)
[524]	365	c
	366	CALL histdef(nid_day, "SWupSFC", "SWup at surface","W/m2",
[776]	367	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	368	. 32, "ave(X)", zstophy,zout)
[524]	369	c
	370	CALL histdef(nid_day, "SWdnTOA", "SWdn at TOA","W/m2",
[776]	371	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	372	. 32, "ave(X)", zstophy,zout)
[524]	373	c
	374	CALL histdef(nid_day, "SWupTOAclr",
	375	. "SWup clear sky at TOA","W/m2",
[776]	376	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	377	. 32, "ave(X)", zstophy,zout)
[524]	378
	379	CALL histdef(nid_day, "SWupSFCclr",
	380	. "SWup clear sky at surface","W/m2",
[776]	381	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	382	. 32, "ave(X)", zstophy,zout)
[524]	383
	384	CALL histdef(nid_day, "SWdnTOAclr",
	385	. "SWdn clear sky at TOA","W/m2",
[776]	386	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	387	. 32, "ave(X)", zstophy,zout)
[524]	388
	389	CALL histdef(nid_day, "SWdnSFCclr",
	390	. "SWdn clear sky at surface","W/m2",
[776]	391	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	392	. 32, "ave(X)", zstophy,zout)
[524]	393	c
[684]	394	cIM AMMA-MIP
	395	CALL histdef(nid_day, "LWdnSFC",
	396	. "LW down at surface","W/m2",
[776]	397	. iim,jj_nb,nhori, 1,1,1,-99,
[684]	398	. 32, "ave(X)", zstophy,zout)
	399	c
	400	CALL histdef(nid_day, "LWupSFC",
	401	. "LW down at surface","W/m2",
[776]	402	. iim,jj_nb,nhori, 1,1,1,-99,
[684]	403	. 32, "ave(X)", zstophy,zout)
	404	c
	405	cIM IF (iflag_con.GE.3) THEN
	406	CALL histdef(nid_day, "ptop",
	407	c . "Cld top pressure","mb",
	408	. "Cld top pressure (cf ISCCP simulator)","Pa",
[776]	409	. iim,jj_nb,nhori, 1,1,1,-99,
[684]	410	. 32, "ave(X)", zstophy,zout)
	411	cIM ENDIF !(iflag_con.GE.3) THEN
	412	c
[524]	413	c=================================================================
	414	c=================================================================
	415	c=================================================================
	416	c INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
	417	c=================================================================
	418	c
[644]	419	CALL histdef(nid_day, "tter", "Surface Temperature ter", "K",
[776]	420	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	421	. "ave(X)", zstophy,zout)
[524]	422	c
[644]	423	CALL histdef(nid_day, "tlic", "Surface Temperature lic", "K",
[776]	424	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	425	. "ave(X)", zstophy,zout)
[524]	426	c
[644]	427	CALL histdef(nid_day, "toce", "Surface Temperature oce", "K",
[776]	428	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	429	. "ave(X)", zstophy,zout)
[524]	430	c
[644]	431	CALL histdef(nid_day, "tsic", "Surface Temperature sic", "K",
[776]	432	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	433	. "ave(X)", zstophy,zout)
[524]	434	c
[644]	435	CALL histdef(nid_day, "t2mter", "Temp.tere 2m", "K",
[776]	436	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	437	. "ave(X)", zstophy,zout)
[524]	438	c
	439	CALL histdef(nid_day, "t2mlic", "Temp.lic 2m", "K",
[776]	440	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	441	. "ave(X)", zstophy,zout)
[524]	442	c
	443	CALL histdef(nid_day, "t2moce", "Temp.oce 2m", "K",
[776]	444	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	445	. "ave(X)", zstophy,zout)
[524]	446	c
	447	CALL histdef(nid_day, "t2msic", "Temp.sic 2m", "K",
[776]	448	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	449	. "ave(X)", zstophy,zout)
[524]	450	c
	451	CALL histdef(nid_day, "t2mter_min", "Temp.terre 2m min.",
	452	. "K",
[776]	453	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	454	. t2mincels, zstophy,zout)
[524]	455	c
	456	CALL histdef(nid_day, "t2mter_max", "Temp.terre 2m max.",
	457	. "K",
[776]	458	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	459	. t2maxcels, zstophy,zout)
[524]	460	c
	461	CALL histdef(nid_day, "u10mter", "Vent zonal ter 10m", "m/s",
[776]	462	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	463	. "ave(X)", zstophy,zout)
[524]	464	c
	465	CALL histdef(nid_day, "u10mlic", "Vent zonal lic 10m", "m/s",
[776]	466	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	467	. "ave(X)", zstophy,zout)
[524]	468	c
	469	CALL histdef(nid_day, "u10moce", "Vent zonal oce 10m", "m/s",
[776]	470	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	471	. "ave(X)", zstophy,zout)
[524]	472	c
	473	CALL histdef(nid_day, "u10msic", "Vent zonal sic 10m",
[776]	474	. "m/s",iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	475	. "ave(X)", zstophy,zout)
[524]	476	c
	477	CALL histdef(nid_day, "v10mter", "Vent meridien ter 10m",
[776]	478	. "m/s", iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	479	. "ave(X)", zstophy,zout)
[524]	480	c
	481	CALL histdef(nid_day, "v10mlic", "Vent meridien lic 10m",
[776]	482	. "m/s",iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	483	. "ave(X)", zstophy,zout)
[524]	484	c
	485	CALL histdef(nid_day, "v10moce", "Vent meridien oce 10m",
[776]	486	. "m/s",iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	487	. "ave(X)", zstophy,zout)
[524]	488	c
	489	CALL histdef(nid_day, "v10msic", "Vent meridien 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	DO nsrf = 1, nbsrf
	494	C
	495	call histdef(nid_day, "pourc_"//clnsurf(nsrf),
[644]	496	$ "% "//clnsurf(nsrf), "%",
[776]	497	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	498	$ "ave(X)", zstophy,zout)
[524]	499	C
[644]	500	call histdef(nid_day, "fract_"//clnsurf(nsrf),
	501	$ "Fraction "//clnsurf(nsrf), "1",
[776]	502	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	503	$ "ave(X)", zstophy,zout)
	504	C
[524]	505	call histdef(nid_day, "tsol_"//clnsurf(nsrf),
[644]	506	$ "Temperature "//clnsurf(nsrf), "K",
[776]	507	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	508	$ "ave(X)", zstophy,zout)
[524]	509	C
	510	call histdef(nid_day, "sens_"//clnsurf(nsrf),
	511	$ "Sensible heat flux "//clnsurf(nsrf), "W/m2",
[776]	512	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	513	$ "ave(X)", zstophy,zout)
[524]	514	c
	515	call histdef(nid_day, "lat_"//clnsurf(nsrf),
	516	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
[776]	517	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	518	$ "ave(X)", zstophy,zout)
[524]	519	C
	520	call histdef(nid_day, "taux_"//clnsurf(nsrf),
[644]	521	$ "Zonal wind stress "//clnsurf(nsrf),"Pa",
[776]	522	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	523	$ "ave(X)", zstophy,zout)
[524]	524
	525	call histdef(nid_day, "tauy_"//clnsurf(nsrf),
[644]	526	$ "Meridional wind stress "//clnsurf(nsrf), "Pa",
[776]	527	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	528	$ "ave(X)", zstophy,zout)
[524]	529	C
	530	call histdef(nid_day, "albe_"//clnsurf(nsrf),
[644]	531	$ "Albedo surf. "//clnsurf(nsrf), "-",
[776]	532	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	533	$ "ave(X)", zstophy,zout)
[524]	534	C
	535	call histdef(nid_day, "rugs_"//clnsurf(nsrf),
[644]	536	$ "Rugosite "//clnsurf(nsrf), "m",
[776]	537	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	538	$ "ave(X)", zstophy,zout)
[524]	539	C
	540	END DO
	541	C
	542	c=================================================================
	543	c
	544	c FIN INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
	545	c
	546	c=================================================================
	547	ENDIF !lev_histday.GE.4
[644]	548	c
	549	IF(lev_histday.GE.5) THEN !lev_histday.GE.5
	550	c
	551	c rajout sorties F. Aires
	552	CALL histdef(nid_day, "tops0", "CS Solar rad. at TOA",
	553	. "W/m2",
[776]	554	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	555	. "ave(X)", zstophy,zout)
	556	c
	557	CALL histdef(nid_day, "topl0", "CS IR rad. at TOA",
	558	. "W/m2",
[776]	559	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	560	. "ave(X)", zstophy,zout)
	561	c
	562	CALL histdef(nid_day, "sols0", "CS Net Solar rad. at surf.",
	563	. "W/m2",
[776]	564	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	565	. "ave(X)", zstophy,zout)
	566	c
	567	CALL histdef(nid_day, "soll0", "CS Net IR rad. at surface",
	568	. "W/m2",
[776]	569	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	570	. "ave(X)", zstophy,zout)
	571	c
	572	CALL histdef(nid_day, "rneb", "Cloud fraction", "-",
[776]	573	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	574	. "ave(X)", zstophy,zout)
	575	c
	576	CALL histdef(nid_day, "rnebcon", "Convective Cloud Fraction"
	577	. , "-",
[776]	578	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	579	. "ave(X)", zstophy,zout)
	580	c
	581	CALL histdef(nid_day,"lwcon","Cloud liquid water content",
	582	. "kg/kg",
[776]	583	. iim,jj_nb,nhori, klev,1,klev, nvert, 32,
[644]	584	. "ave(X)", zstophy,zout)
	585	c
	586	CALL histdef(nid_day,"iwcon","Cloud ice water content","kg/kg",
[776]	587	. iim,jj_nb,nhori, klev,1,klev, nvert, 32,
[644]	588	. "ave(X)", zstophy,zout)
	589	c
	590	CALL histdef(nid_day,"lwp","Cloud water path","kg/m2",
[776]	591	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	592	. "ave(X)", zstophy,zout)
	593	c
	594	CALL histdef(nid_day,"iwp","Cloud ice water path","kg/m2",
[776]	595	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	596	. "ave(X)", zstophy,zout)
	597	c
	598	CALL histdef(nid_day,"meantaucld",
	599	. "ISCCP mean cloud optical thickness","1",
[776]	600	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	601	. "ave(X)", zstophy,zout)
	602	c
	603	CALL histdef(nid_day,"cldtau","Cloud optical thickness","1",
[776]	604	. iim,jj_nb,nhori, klev,1,klev, nvert, 32,
[644]	605	. "ave(X)", zstophy,zout)
	606	c
	607	CALL histdef(nid_day,"cldemi","Cloud optical emissivity","1",
[776]	608	. iim,jj_nb,nhori, klev,1,klev, nvert, 32,
[644]	609	. "ave(X)", zstophy,zout)
	610	c
[684]	611	cIM rajout AMMA-MIP
	612	c
	613	CALL histdef(nid_day, "dtcon", "Convection dT", "K/s",
[776]	614	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[684]	615	. "ave(X)", zstophy,zout)
	616	c
	617	CALL histdef(nid_day, "dqcon", "Convection dQ", "(kg/kg)/s",
[776]	618	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[684]	619	. "ave(X)", zstophy,zout)
	620	c
	621	CALL histdef(nid_day, "rh2m", "Relative humidity at 2m", "%",
[776]	622	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[684]	623	. "ave(X)", zstophy,zout)
	624	c
	625	CALL histdef(nid_day, "qsat2m", "Saturant humidity at 2m", "%",
[776]	626	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[684]	627	. "ave(X)", zstophy,zout)
	628	c
	629	CALL histdef(nid_day, "tpot",
	630	. "Surface air potential temperature", "K",
[776]	631	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[684]	632	. "ave(X)", zstophy,zout)
	633	c
	634	CALL histdef(nid_day, "tpote",
	635	. "Surface air equivalent potential temperature", "K",
[776]	636	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[684]	637	. "ave(X)", zstophy,zout)
	638	c
[644]	639	ENDIF !lev_histday.GE.5
[524]	640	c=================================================================
	641	c
	642	CALL histend(nid_day)
	643	c
	644	ndex2d = 0
	645	ndex3d = 0
	646	c
	647	c=================================================================
[766]	648	c$OMP END MASTER
[524]	649	ENDIF ! fin de test sur ok_journe

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

Download in other formats: