Context Navigation

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

Last change on this file since 891 was 890, checked in by Laurent Fairhead, 18 years ago
Mise sous flag ok_ISCCP d'une sortie d'une variable dependant du simulateur ISCCP Correction coquille sur des sorties sous-surfaces SD/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
	85	CALL histdef(nid_day, "plul", "Large-scale Precip.",
	86	. "kg/(s*m2)",
[776]	87	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[644]	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: