Context Navigation

source: LMDZ4/trunk/libf/phytherm/ini_histday.h @ 945

Last change on this file since 945 was 860, checked in by lmdzadmin, 17 years ago
Unite de vitw etait fausse dans histday.nc et histmth.nc IM
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 25.6 KB

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

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

Download in other formats: