Context Navigation

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

Last change on this file since 703 was 684, checked in by lmdzadmin, 18 years ago
Ajout nouveaux flags (ip_ebil_phy, ok_slab_sicOBS) ds physiq.def + Lecture freq. ecriture en nombre de jours IM
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 24.5 KB

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

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

Download in other formats: