Context Navigation

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

Last change on this file since 1080 was 895, checked in by Laurent Fairhead, 17 years ago
Un peu de menage sur les fichiers de sortie apres les tests aquaplanete FH LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 28.1 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
[878]	55	CALL histdef(nid_day, "weakinv", "Weak inversion", "",
	56	s iim,jjmp1,nhori, 1,1,1, -99, 32,
	57	s "ave(X)", zstophy,zout)
	58
	59	CALL histdef(nid_day, "dthmin", "dTheta mini", "K/m",
	60	s iim,jjmp1,nhori, 1,1,1, -99, 32,
	61	s "ave(X)", zstophy,zout)
	62
	63
[524]	64	CALL histdef(nid_day, "tsol", "Surface Temperature", "K",
[776]	65	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	66	. "ave(X)", zstophy,zout)
[524]	67	c
	68	CALL histdef(nid_day, "t2m", "Temperature 2m", "K",
[776]	69	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	70	. "ave(X)", zstophy,zout)
[524]	71	c
[766]	72	c Champs retires momentannement en attendant un hypothetique
[524]	73	c debugage
	74
	75	CALL histdef(nid_day, "t2m_min", "Temp. 2m min.",
	76	. "K",
[776]	77	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	78	. t2mincels, zstophy,zout)
[524]	79	c
	80	CALL histdef(nid_day, "t2m_max", "Temp. 2m max.",
	81	. "K",
[776]	82	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	83	. t2maxcels, zstophy,zout)
[524]	84	c
[895]	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)
[524]	89	c
	90	CALL histdef(nid_day, "pluc", "Convective Precip.",
	91	. "kg/(s*m2)",
[776]	92	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	93	. "ave(X)", zstophy,zout)
[524]	94	c
	95	CALL histdef(nid_day, "snowl", "Solid Large-scale Precip.",
[644]	96	. "kg/(m2*s)",
[776]	97	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	98	. "ave(X)", zstophy,zout)
[524]	99	c
	100	IF(1.EQ.0) THEN !snowc=0.
	101	CALL histdef(nid_day, "snowc", "Solid Convective Precip.",
[644]	102	. "kg/(m2*s)",
[776]	103	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	104	. "ave(X)", zstophy,zout)
[524]	105	ENDIF !snowc=0.
	106	c
	107	CALL histdef(nid_day, "flat", "Latent heat flux", "W/m2",
[776]	108	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	109	. "ave(X)", zstophy,zout)
[524]	110	c
	111	CALL histdef(nid_day, "sicf", "Sea-ice fraction", "-",
[776]	112	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	113	. "inst(X)", zstoday,zout)
[524]	114	c
	115	CALL histdef(nid_day, "q2m", "Specific humidity", "kg/kg",
[776]	116	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	117	. "ave(X)", zstophy,zout)
[524]	118	c
	119	CALL histdef(nid_day, "u10m", "Vent zonal 10m", "m/s",
[776]	120	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	121	. "ave(X)", zstophy,zout)
[524]	122	c
	123	CALL histdef(nid_day, "v10m", "Vent meridien 10m", "m/s",
[776]	124	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	125	. "ave(X)", zstophy,zout)
[524]	126	c
	127	CALL histdef(nid_day, "wind10m","10-m wind speed","m/s",
[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, "wind10max", "10-m wind speed max.",
	132	. "m/s",
[776]	133	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	134	. t2maxcels, zstophy,zout)
[524]	135	c
	136	CALL histdef(nid_day, "psol", "Surface Pressure", "Pa",
[776]	137	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	138	. "ave(X)", zstophy,zout)
[524]	139	c
	140	CALL histdef(nid_day, "precip","Precipitation Totale liq+sol"
	141	. , "kg/(s*m2)",
[776]	142	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	143	. "ave(X)", zstophy,zout)
[524]	144	c
[879]	145	CALL histdef(nid_day, "plul","Precipitation ls liq+sol"
	146	. , "kg/(s*m2)",
	147	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	148	. "ave(X)", zstophy,zout)
	149	c
[644]	150	CALL histdef(nid_day, "snowf", "Snow fall", "kg/(m2*s)",
[776]	151	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	152	. "ave(X)", zstophy,zout)
[524]	153	c
[644]	154	c CALL histdef(nid_day, "snow_mass", "Snow Mass", "kg/m2",
[776]	155	c . iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	156	c . "ave(X)", zstophy,zout)
[524]	157	c
[644]	158	CALL histdef(nid_day, "evap", "Evaporation", "kg/(m2*s)",
[776]	159	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	160	. "ave(X)", zstophy,zout)
[524]	161	c
	162	CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
[776]	163	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	164	. "ave(X)", zstophy,zout)
[524]	165	c
	166	CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
[776]	167	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	168	. "ave(X)", zstophy,zout)
[524]	169	c
	170	CALL histdef(nid_day, "sols", "Net Solar rad. at surf.",
	171	. "W/m2",
[776]	172	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	173	. "ave(X)", zstophy,zout)
[524]	174	c
	175	CALL histdef(nid_day, "soll", "Net IR rad. at surface", "W/m2",
[776]	176	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	177	. "ave(X)", zstophy,zout)
[524]	178	c
[684]	179	CALL histdef(nid_day, "radsol", "Rayonnement au sol", "W/m2",
[776]	180	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[684]	181	. "ave(X)", zstophy,zout)
	182	c
[524]	183	CALL histdef(nid_day, "cldl", "Low-level cloudiness", "-",
[776]	184	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	185	. "ave(X)", zstophy,zout)
[524]	186	c
	187	CALL histdef(nid_day, "cldm", "Mid-level cloudiness", "-",
[776]	188	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	189	. "ave(X)", zstophy,zout)
[524]	190	c
	191	CALL histdef(nid_day, "cldh", "High-level cloudiness", "-",
[776]	192	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	193	. "ave(X)", zstophy,zout)
[524]	194	c
	195	CALL histdef(nid_day, "cldt", "Total cloudiness", "-",
[776]	196	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	197	. "ave(X)", zstophy,zout)
[524]	198	c
	199	CALL histdef(nid_day, "cldq", "Cloud liquid water path",
	200	. "kg/m2",
[776]	201	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	202	. "ave(X)", zstophy,zout)
[524]	203	c
	204	CALL histdef(nid_day, "prw", "Precipitable water", "kg/m2",
[776]	205	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	206	. "ave(X)", zstophy,zout)
[524]	207	c
	208	c Champs dynamiques sur niveaux de pression
	209	c
[644]	210	c DO k=1, nlevENS
	211	DO k=1, nlevSTD
	212	c
	213	c bb=clevSTD(k)
	214	c
	215	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
	216	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
	217	c aa=clevSTD(k)
	218	c bb2=clevSTD(k)
	219	c ENDIF
	220	c
	221	IF(bb2.EQ."850".OR.bb2.EQ."700"
	222	. .OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
	223	c
	224	CALL histdef(nid_day, "u"//bb2,
	225	. "Zonal wind "//bb2//"mb","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, "v"//bb2,
	230	. "Meridional wind "//bb2//"mb","m/s",
[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."850".OR.bb2.EQ."700".OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
[524]	235	c
[644]	236	IF(bb2.EQ."500") THEN
	237	c
	238	CALL histdef(nid_day, "w"//bb2, "Vertical wind "//bb2//"mb",
	239	. "m/s",
[776]	240	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	241	. "inst(X)", zout,zout)
[524]	242	c
[644]	243	CALL histdef(nid_day,"phi"//bb2, "Geopotential "//bb2//"mb",
	244	. "m",
[776]	245	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	246	. "inst(X)", zout,zout)
[524]	247	c
[644]	248	ENDIF !(bb2.EQ."500")
	249	c
	250	ENDDO !nlevSTD
	251	c
[524]	252	CALL histdef(nid_day, "slp", "Sea Level Pressure", "Pa",
[776]	253	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	254	. "ave(X)", zstophy,zout)
[524]	255	c
	256	CALL histdef(nid_day, "cape_max", "CAPE max.",
	257	. "J/kg",
[776]	258	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	259	. capemaxcels, zstophy,zout)
[524]	260	c
	261	CALL histdef(nid_day, "solldown", "Down. IR rad. at surface",
[776]	262	. "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	263	. "ave(X)", zstophy,zout)
[524]	264	c
	265	CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2",
[776]	266	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	267	. "ave(X)", zstophy,zout)
[524]	268	c
	269	CALL histdef(nid_day, "SWdnSFC", "SWdn at surface","W/m2",
[776]	270	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	271	. 32, "ave(X)", zstophy,zout)
[524]	272	c
[644]	273	cIM diagnostiques flux ocean-atm ou ocean-glace de mer
	274	cIM pour utilisation dans un modele de "slab" ocean
	275	c
	276	IF (OCEAN.EQ.'force ') THEN
	277	c
	278	CALL histdef(nid_day, "lmt_bils",
	279	. "Bilan au sol atmosphere forcee", "W/m2",
[776]	280	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	281	. "ave(X)", zstophy,zout)
	282	c
	283	ELSE IF (OCEAN.EQ.'slab ') THEN
	284	c
	285	CALL histdef(nid_day, "slab_bils",
	286	. "Bilan au sol slab", "W/m2",
[776]	287	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	288	. "ave(X)", zstophy,zout)
	289	c
	290	CALL histdef(nid_day, "tslab", "Slab SST ", "K",
[776]	291	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	292	. "ave(X)", zstophy,zout)
	293	c
	294	CALL histdef(nid_day, "seaice", "Slab seaice", "kg/m2",
[776]	295	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	296	. "ave(X)", zstophy,zout)
	297	c
	298	CALL histdef(nid_day, "siceh", "Slab seaice height", "m",
[776]	299	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	300	. "ave(X)", zstophy,zout)
	301	c
	302	ENDIF !OCEAN.EQ.force/slab
	303	c
	304	CALL histdef(nid_day, "fluxo",
	305	. "Flux turbulents ocean-atmosphere", "W/m2",
[776]	306	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	307	. "ave(X)", zstophy,zout)
	308	c
	309	CALL histdef(nid_day, "fluxg",
	310	. "Flux ocean-glace de mer", "W/m2",
[776]	311	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	312	. "ave(X)", zstophy,zout)
	313	c
[524]	314	ENDIF !lev_histday.GE.1
	315	c
	316	IF (lev_histday.GE.2) THEN
	317	c
	318	CALL histdef(nid_day, "bils", "Surf. total heat flux", "W/m2",
[776]	319	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	320	. "ave(X)", zstophy,zout)
[524]	321	c
	322	CALL histdef(nid_day, "fder", "Heat flux derivation", "W/m2/K",
[776]	323	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	324	. "ave(X)", zstophy,zout)
[524]	325	c
	326	ENDIF !lev_histday.GE.2
	327	c=================================================================
	328	IF(lev_histday.GE.3) THEN
	329	c=================================================================
	330	c INITIALISATION DES CHAMPS 3D
	331	c=================================================================
	332	c=================================================================
	333	c=================================================================
	334	c Champs 3D:
	335	c
	336	CALL histdef(nid_day, "temp", "Air temperature", "K",
[776]	337	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	338	. "ave(X)", zstophy,zout)
[524]	339	c
[879]	340	CALL histdef(nid_day, "rhum", "relative humidity", " ",
	341	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
	342	. "ave(X)", zstophy,zout)
	343	c
[524]	344	CALL histdef(nid_day, "ovap", "Specific humidity", "kg/kg",
[776]	345	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	346	. "ave(X)", zstophy,zout)
[524]	347	c
[644]	348	CALL histdef(nid_day, "geop", "Geopotential height", "m2/s2",
[776]	349	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	350	. "ave(X)", zstophy,zout)
[524]	351	c
	352	CALL histdef(nid_day, "vitu", "Zonal wind", "m/s",
[776]	353	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	354	. "ave(X)", zstophy,zout)
[524]	355	c
	356	CALL histdef(nid_day, "vitv", "Meridional wind", "m/s",
[776]	357	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	358	. "ave(X)", zstophy,zout)
[524]	359	c
[860]	360	CALL histdef(nid_day, "vitw", "Vertical wind", "Pa/s",
[776]	361	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	362	. "ave(X)", zstophy,zout)
[524]	363	c
	364	CALL histdef(nid_day, "pres", "Air pressure", "Pa",
[776]	365	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	366	. "ave(X)", zstophy,zout)
[524]	367	c
	368	c=================================================================
	369	c FIN INITIALISATION DES CHAMPS 3D
	370	c=================================================================
	371	ENDIF ! lev_histday.GE.3
	372	c
	373	IF (lev_histday.GE.4) THEN
	374	c=================================================================
	375	c
	376	c INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
	377	c
	378	c=================================================================
	379	c
	380	CALL histdef(nid_day, "SWupTOA", "SWup at TOA","W/m2",
[776]	381	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	382	. 32, "ave(X)", zstophy,zout)
[524]	383	c
	384	CALL histdef(nid_day, "SWupSFC", "SWup at surface","W/m2",
[776]	385	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	386	. 32, "ave(X)", zstophy,zout)
[524]	387	c
	388	CALL histdef(nid_day, "SWdnTOA", "SWdn at TOA","W/m2",
[776]	389	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	390	. 32, "ave(X)", zstophy,zout)
[524]	391	c
	392	CALL histdef(nid_day, "SWupTOAclr",
	393	. "SWup clear sky at TOA","W/m2",
[776]	394	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	395	. 32, "ave(X)", zstophy,zout)
[524]	396
	397	CALL histdef(nid_day, "SWupSFCclr",
	398	. "SWup clear sky at surface","W/m2",
[776]	399	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	400	. 32, "ave(X)", zstophy,zout)
[524]	401
	402	CALL histdef(nid_day, "SWdnTOAclr",
	403	. "SWdn clear sky at TOA","W/m2",
[776]	404	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	405	. 32, "ave(X)", zstophy,zout)
[524]	406
	407	CALL histdef(nid_day, "SWdnSFCclr",
	408	. "SWdn clear sky at surface","W/m2",
[776]	409	. iim,jj_nb,nhori, 1,1,1,-99,
[644]	410	. 32, "ave(X)", zstophy,zout)
[524]	411	c
[684]	412	cIM AMMA-MIP
	413	CALL histdef(nid_day, "LWdnSFC",
	414	. "LW down at surface","W/m2",
[776]	415	. iim,jj_nb,nhori, 1,1,1,-99,
[684]	416	. 32, "ave(X)", zstophy,zout)
	417	c
	418	CALL histdef(nid_day, "LWupSFC",
	419	. "LW down at surface","W/m2",
[776]	420	. iim,jj_nb,nhori, 1,1,1,-99,
[684]	421	. 32, "ave(X)", zstophy,zout)
	422	c
	423	cIM IF (iflag_con.GE.3) THEN
[890]	424	if (ok_ISCCP) then
	425	CALL histdef(nid_day, "ptop",
[684]	426	c . "Cld top pressure","mb",
	427	. "Cld top pressure (cf ISCCP simulator)","Pa",
[776]	428	. iim,jj_nb,nhori, 1,1,1,-99,
[684]	429	. 32, "ave(X)", zstophy,zout)
[890]	430	ENDIF
[684]	431	cIM ENDIF !(iflag_con.GE.3) THEN
	432	c
[524]	433	c=================================================================
	434	c=================================================================
	435	c=================================================================
	436	c INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
	437	c=================================================================
	438	c
[644]	439	CALL histdef(nid_day, "tter", "Surface Temperature ter", "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, "tlic", "Surface Temperature lic", "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, "toce", "Surface Temperature oce", "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, "tsic", "Surface Temperature sic", "K",
[776]	452	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	453	. "ave(X)", zstophy,zout)
[524]	454	c
[644]	455	CALL histdef(nid_day, "t2mter", "Temp.tere 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, "t2mlic", "Temp.lic 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, "t2moce", "Temp.oce 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, "t2msic", "Temp.sic 2m", "K",
[776]	468	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	469	. "ave(X)", zstophy,zout)
[524]	470	c
	471	CALL histdef(nid_day, "t2mter_min", "Temp.terre 2m min.",
	472	. "K",
[776]	473	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	474	. t2mincels, zstophy,zout)
[524]	475	c
	476	CALL histdef(nid_day, "t2mter_max", "Temp.terre 2m max.",
	477	. "K",
[776]	478	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	479	. t2maxcels, zstophy,zout)
[524]	480	c
	481	CALL histdef(nid_day, "u10mter", "Vent zonal ter 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, "u10mlic", "Vent zonal lic 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, "u10moce", "Vent zonal oce 10m", "m/s",
[776]	490	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	491	. "ave(X)", zstophy,zout)
[524]	492	c
	493	CALL histdef(nid_day, "u10msic", "Vent zonal sic 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, "v10mter", "Vent meridien ter 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, "v10mlic", "Vent meridien lic 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, "v10moce", "Vent meridien oce 10m",
[776]	506	. "m/s",iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	507	. "ave(X)", zstophy,zout)
[524]	508	c
	509	CALL histdef(nid_day, "v10msic", "Vent meridien sic 10m",
[776]	510	. "m/s",iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	511	. "ave(X)", zstophy,zout)
[524]	512	c
	513	DO nsrf = 1, nbsrf
	514	C
	515	call histdef(nid_day, "pourc_"//clnsurf(nsrf),
[644]	516	$ "% "//clnsurf(nsrf), "%",
[776]	517	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	518	$ "ave(X)", zstophy,zout)
[524]	519	C
[644]	520	call histdef(nid_day, "fract_"//clnsurf(nsrf),
	521	$ "Fraction "//clnsurf(nsrf), "1",
[776]	522	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	523	$ "ave(X)", zstophy,zout)
	524	C
[524]	525	call histdef(nid_day, "tsol_"//clnsurf(nsrf),
[644]	526	$ "Temperature "//clnsurf(nsrf), "K",
[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, "sens_"//clnsurf(nsrf),
	531	$ "Sensible heat flux "//clnsurf(nsrf), "W/m2",
[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, "lat_"//clnsurf(nsrf),
	536	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
[776]	537	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	538	$ "ave(X)", zstophy,zout)
[524]	539	C
	540	call histdef(nid_day, "taux_"//clnsurf(nsrf),
[644]	541	$ "Zonal wind stress "//clnsurf(nsrf),"Pa",
[776]	542	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	543	$ "ave(X)", zstophy,zout)
[524]	544
	545	call histdef(nid_day, "tauy_"//clnsurf(nsrf),
[644]	546	$ "Meridional wind stress "//clnsurf(nsrf), "Pa",
[776]	547	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	548	$ "ave(X)", zstophy,zout)
[524]	549	C
	550	call histdef(nid_day, "albe_"//clnsurf(nsrf),
[644]	551	$ "Albedo surf. "//clnsurf(nsrf), "-",
[776]	552	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	553	$ "ave(X)", zstophy,zout)
[524]	554	C
	555	call histdef(nid_day, "rugs_"//clnsurf(nsrf),
[644]	556	$ "Rugosite "//clnsurf(nsrf), "m",
[776]	557	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	558	$ "ave(X)", zstophy,zout)
[524]	559	C
[878]	560
	561
	562	! FH Sorties specifiques pour Mellor et Yamada
	563	if (iflag_pbl>1 .and. lev_histday.gt.10 ) then
	564	call histdef(nid_day, "tke_"//clnsurf(nsrf),
	565	$ "Max Turb. Kinetic Energy "//clnsurf(nsrf), "-",
	566	$ iim,jj_nb,nhori, klev,1,klev,nvert, 32,
	567	$ "ave(X)", zstophy,zout)
	568
	569	call histdef(nid_day, "tke_max_"//clnsurf(nsrf),
	570	$ "Max Turb. Kinetic Energy "//clnsurf(nsrf), "-",
	571	$ iim,jj_nb,nhori, klev,1,klev,nvert, 32,
	572	$ "t_max(X)", zstophy,zout)
	573	endif
	574
	575	C
[524]	576	END DO
	577	C
	578	c=================================================================
	579	c
	580	c FIN INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
	581	c
	582	c=================================================================
	583	ENDIF !lev_histday.GE.4
[644]	584	c
	585	IF(lev_histday.GE.5) THEN !lev_histday.GE.5
	586	c
	587	c rajout sorties F. Aires
	588	CALL histdef(nid_day, "tops0", "CS Solar rad. at TOA",
	589	. "W/m2",
[776]	590	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	591	. "ave(X)", zstophy,zout)
	592	c
	593	CALL histdef(nid_day, "topl0", "CS IR rad. at TOA",
	594	. "W/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, "sols0", "CS Net Solar rad. at surf.",
	599	. "W/m2",
[776]	600	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	601	. "ave(X)", zstophy,zout)
	602	c
	603	CALL histdef(nid_day, "soll0", "CS Net IR rad. at surface",
	604	. "W/m2",
[776]	605	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	606	. "ave(X)", zstophy,zout)
	607	c
	608	CALL histdef(nid_day, "rneb", "Cloud fraction", "-",
[776]	609	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	610	. "ave(X)", zstophy,zout)
	611	c
[879]	612	CALL histdef(nid_day, "dqpbl", "PBL dT", "K/s",
	613	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
	614	. "ave(X)", zstophy,zout)
	615	c
[644]	616	CALL histdef(nid_day, "rnebcon", "Convective Cloud Fraction"
	617	. , "-",
[776]	618	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[644]	619	. "ave(X)", zstophy,zout)
	620	c
[879]	621	CALL histdef(nid_day, "dtwak", "Wake dT", "K/s",
	622	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
	623	. "ave(X)", zstophy,zout)
	624	c
[644]	625	CALL histdef(nid_day,"lwcon","Cloud liquid water content",
	626	. "kg/kg",
[776]	627	. iim,jj_nb,nhori, klev,1,klev, nvert, 32,
[644]	628	. "ave(X)", zstophy,zout)
	629	c
	630	CALL histdef(nid_day,"iwcon","Cloud ice water content","kg/kg",
[776]	631	. iim,jj_nb,nhori, klev,1,klev, nvert, 32,
[644]	632	. "ave(X)", zstophy,zout)
	633	c
	634	CALL histdef(nid_day,"lwp","Cloud water path","kg/m2",
[776]	635	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	636	. "ave(X)", zstophy,zout)
	637	c
	638	CALL histdef(nid_day,"iwp","Cloud ice water path","kg/m2",
[776]	639	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	640	. "ave(X)", zstophy,zout)
	641	c
[879]	642
	643	CALL histdef(nid_day, "cape", "Conv avlbl pot ener", "J/kg",
	644	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	645	. "ave(X)", zstophy,zout)
	646
	647	CALL histdef(nid_day, "wh", "Wake height", "m",
	648	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	649	. "ave(X)", zstophy,zout)
	650	c
	651	CALL histdef(nid_day, "ws", "Wake surface", "m2",
	652	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	653	. "ave(X)", zstophy,zout)
	654	c
[644]	655	CALL histdef(nid_day,"meantaucld",
	656	. "ISCCP mean cloud optical thickness","1",
[776]	657	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	658	. "ave(X)", zstophy,zout)
	659	c
	660	CALL histdef(nid_day,"cldtau","Cloud optical thickness","1",
[776]	661	. iim,jj_nb,nhori, klev,1,klev, nvert, 32,
[644]	662	. "ave(X)", zstophy,zout)
	663	c
	664	CALL histdef(nid_day,"cldemi","Cloud optical emissivity","1",
[776]	665	. iim,jj_nb,nhori, klev,1,klev, nvert, 32,
[644]	666	. "ave(X)", zstophy,zout)
	667	c
[879]	668	c
	669	CALL histdef(nid_day, "cin", "Conv Inhibition", "J/kg",
	670	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	671	. "ave(X)", zstophy,zout)
	672	c
	673	CALL histdef(nid_day, "wale", "Wake Available Energy", "J/kg",
	674	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	675	. "ave(X)", zstophy,zout)
	676	c
	677	CALL histdef(nid_day, "walp",
	678	. "Available Lifting Energy due to wake", "W/m2",
	679	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	680	. "ave(X)", zstophy,zout)
	681	c
	682	CALL histdef(nid_day, "blale", "PBL Available Energy", "J/kg",
	683	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	684	. "ave(X)", zstophy,zout)
	685	c
	686	CALL histdef(nid_day, "blalp",
	687	. "Available Lifting Energy due to PBL", "W/m2",
	688	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	689	. "ave(X)", zstophy,zout)
	690	c
	691	CALL histdef(nid_day, "wdt1",
	692	. "Temp diff wake layer 1", "K",
	693	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	694	. "ave(X)", zstophy,zout)
	695	c
	696	CALL histdef(nid_day, "wdq1",
	697	. "Temp diff wake layer 1", "g/kg",
	698	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	699	. "ave(X)", zstophy,zout)
[684]	700	cIM rajout AMMA-MIP
	701	c
	702	CALL histdef(nid_day, "dtcon", "Convection dT", "K/s",
[776]	703	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[684]	704	. "ave(X)", zstophy,zout)
	705	c
	706	CALL histdef(nid_day, "dqcon", "Convection dQ", "(kg/kg)/s",
[776]	707	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[684]	708	. "ave(X)", zstophy,zout)
	709	c
	710	CALL histdef(nid_day, "rh2m", "Relative humidity at 2m", "%",
[776]	711	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[684]	712	. "ave(X)", zstophy,zout)
	713	c
	714	CALL histdef(nid_day, "qsat2m", "Saturant humidity at 2m", "%",
[776]	715	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[684]	716	. "ave(X)", zstophy,zout)
	717	c
	718	CALL histdef(nid_day, "tpot",
	719	. "Surface air potential temperature", "K",
[776]	720	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[684]	721	. "ave(X)", zstophy,zout)
	722	c
	723	CALL histdef(nid_day, "tpote",
	724	. "Surface air equivalent potential temperature", "K",
[776]	725	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[684]	726	. "ave(X)", zstophy,zout)
	727	c
[644]	728	ENDIF !lev_histday.GE.5
[524]	729	c=================================================================
	730	c
	731	CALL histend(nid_day)
	732	c
	733	ndex2d = 0
	734	ndex3d = 0
	735	c
	736	c=================================================================
[766]	737	c$OMP END MASTER
[524]	738	ENDIF ! fin de test sur ok_journe

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

Download in other formats: