Context Navigation

source: LMDZ4/trunk/libf/phylmd/ini_histhf.h @ 665

Last change on this file since 665 was 644, checked in by Laurent Fairhead, 20 years ago
Synchronisation avec tous les diagnostiques de Ionela IM Inclusion du slab ocean IM LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 12.6 KB

Rev	Line
[524]	1	!
	2	! $Header$
	3	!
	4
	5	IF (ok_hf) THEN
	6	c
[644]	7	zstohf = dtime * FLOAT(ecrit_hf)
	8	zstophy = dtime
	9	zstoday = dtime * FLOAT(ecrit_day)
	10	zout = dtime * FLOAT(ecrit_hf)
[524]	11	c
	12	idayref = day_ref
	13	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
	14
	15	c
	16	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
	17	DO i = 1, iim
	18	zx_lon(i,1) = rlon(i+1)
	19	zx_lon(i,jjmp1) = rlon(i+1)
	20	ENDDO
	21
	22	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
	23	CALL histbeg("histhf", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
	24	. 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
	25	. nhori, nid_hf)
	26
	27	CALL histvert(nid_hf, "presnivs", "Vertical levels", "mb",
[644]	28	. klev, presnivs/100., nvert)
[524]	29	c
	30
[644]	31	IF(type_run.EQ."CLIM".OR.type_run.EQ."ENSP") THEN
	32	c
[524]	33	IF(lev_histhf.GE.1) THEN
	34	c
	35	CALL histdef(nid_hf, "aireTER","Grid area CONT","-",
	36	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	37	. "once", zstohf,zout)
[524]	38	c
	39	CALL histdef(nid_hf, "contfracATM","% sfce ter+lic ","-",
	40	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	41	. "once", zstohf,zout)
[524]	42	c
	43	CALL histdef(nid_hf, "contfracOR","% sfce terre OR", "-",
	44	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	45	. "once", zstohf,zout)
[524]	46	c
	47	c Champs 2D:
	48	c
	49	CALL histdef(nid_hf, "t2m", "Temperature 2m", "K",
	50	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	51	. "inst(X)", zstohf,zout)
[524]	52	c
	53	CALL histdef(nid_hf, "q2m", "Specific humidity", "kg/kg",
	54	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	55	. "inst(X)", zstohf,zout)
[524]	56	c
	57	CALL histdef(nid_hf, "psol", "Surface Pressure", "Pa",
	58	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	59	. "inst(X)", zstohf,zout)
[524]	60	c
[644]	61	CALL histdef(nid_hf, "precip", "Precipitation", "kg/(s*m2)",
[524]	62	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	63	. "ave(X)", zstophy,zout)
[524]	64	c
	65	cIM ENSEMBLES BEG
	66	CALL histdef(nid_hf, "tsol", "Surface Temperature", "K",
	67	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	68	. "inst(X)", zstohf,zout)
[524]	69	c
	70	CALL histdef(nid_hf, "slp", "Sea Level Pressure", "Pa",
	71	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	72	. "inst(X)", zstohf,zout)
[524]	73	c
	74	CALL histdef(nid_hf, "u10m", "Vent zonal 10m", "m/s",
	75	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	76	. "inst(X)", zstohf,zout)
[524]	77	c
	78	CALL histdef(nid_hf, "v10m", "Vent meridien 10m", "m/s",
	79	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	80	. "inst(X)", zstohf,zout)
[524]	81	c
	82	CALL histdef(nid_hf, "wind10m","10-m wind speed","m/s",
	83	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	84	. "inst(X)", zstohf,zout)
[524]	85	c
[644]	86	DO k=1, nlevSTD
	87	c
	88	c bb=clevSTD(k)
	89	c
	90	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
	91	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
	92	c
	93	c IF(k.GE.2) THEN
	94	c aa=clevSTD(k)
	95	c bb=aa(1:lnblnk1(aa))
	96	c ENDIF
	97	c
	98	IF(bb2.EQ."500") THEN
	99	CALL histdef(nid_hf, "phi"//bb2,
	100	. "Geopotential "//bb2//"mb", "m",
[524]	101	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	102	. "inst(X)", zstohf,zout)
	103	ENDIF !bb2.EQ."500"
	104	ENDDO !k=1, nlevSTD
[524]	105	c
	106	ENDIF !lev_histhf.GE.1
	107	c
	108	IF(lev_histhf.GE.2) THEN
	109	c
[644]	110	CALL histdef(nid_hf, "cldt", "Total cloudiness", "%",
[524]	111	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	112	. "ave(X)", zstophy,zout)
[524]	113	c
[644]	114	CALL histdef(nid_hf, "SWnetOR",
	115	. "Sfce net SW radiation OR", "W/m2",
[524]	116	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	117	. "ave(X)", zstophy,zout)
[524]	118	c
[644]	119	CALL histdef(nid_hf, "SWdownOR",
	120	. "Sfce incident SW radiation OR", "W/m2",
[590]	121	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	122	. "ave(X)", zstophy,zout)
[590]	123	c
[524]	124	CALL histdef(nid_hf, "LWdownOR",
[644]	125	. "Sfce incident LW radiation OR", "W/m2",
[524]	126	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	127	. "ave(X)", zstophy,zout)
[524]	128	ENDIF !lev_histhf.GE.2
	129	c
	130	cIM ENSEMBLES END
	131	c
	132	IF(lev_histhf.GE.3) THEN
	133	c
[644]	134	DO k=1, nlevSTD
[524]	135	c
[644]	136	c bb=clevSTD(k)
	137	c
	138	c IF(k.GE.2) THEN
	139	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
	140	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
	141	c aa=clevSTD(k)
	142	c bb=aa(1:lnblnk1(aa))
	143	c ENDIF
	144	c
	145	IF(bb2.EQ."850".or.bb2.EQ."700".OR.
	146	$ bb2.EQ."500".or.bb2.EQ."200") THEN
	147	c
	148	CALL histdef(nid_hf, "t"//bb2,
	149	. "Temperature "//bb2//"mb","K",
[524]	150	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	151	. "inst(X)", zstohf,zout)
[524]	152	c
[644]	153	IF(bb2.NE."500") THEN
	154	CALL histdef(nid_hf, "phi"//bb2,
	155	. "Geopotential "//bb2//"mb", "m",
[524]	156	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	157	. "inst(X)", zstohf,zout)
[524]	158	ENDIF
	159	c
[644]	160	CALL histdef(nid_hf, "q"//bb2,
	161	. "Specific humidity "//bb2//"mb","kg/kg",
[524]	162	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	163	. "inst(X)", zstohf,zout)
[524]	164	c
[644]	165	CALL histdef(nid_hf, "u"//bb2,
	166	. "Zonal wind "//bb2//"mb","m/s",
[524]	167	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	168	. "inst(X)", zstohf,zout)
[524]	169	c
[644]	170	CALL histdef(nid_hf, "v"//bb2,
	171	. "Meridional wind "//bb2//"mb","m/s",
[524]	172	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[644]	173	. "inst(X)", zstohf,zout)
[524]	174	c
[644]	175	ENDIF !IF(bb2.EQ."850".or.bb2.EQ."700".OR.
[524]	176	c
[644]	177	ENDDO !nlevSTD
[524]	178	c
	179	ENDIF !lev_histhf.GE.3
	180	c
	181	IF(lev_histhf.GE.4) THEN
[644]	182	c
	183	cIM sorties hf 3D
	184	c
	185	#undef histhf3d
[524]	186	#define histhf3d
	187	#ifdef histhf3d
	188	#include "ini_histhf3d.h"
	189	#endif
	190	ENDIF !lev_histhf.GE.4
	191	c
[644]	192	ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
[524]	193	c
[644]	194	IF(lev_histhf.GE.1) THEN
[524]	195	c
[644]	196	DO k=1, nlevSTD
	197	c
	198	c
	199	c bb=clevSTD(k)
	200	c
	201	c IF(k.GE.2) THEN
	202	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
	203	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
	204	c aa=clevSTD(k)
	205	c bb=aa(1:lnblnk1(aa))
	206	c ENDIF
	207	c
	208	IF(bb2.EQ."850".or.bb2.EQ."200") THEN
	209	c
	210	CALL histdef(nid_hf, "u"//bb2,
	211	. "Zonal wind "//bb2//"mb","m/s",
	212	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	213	. "inst(X)", zstohf,zstohf)
	214	c
	215	CALL histdef(nid_hf, "v"//bb2,
	216	. "Meridional wind "//bb2//"mb","m/s",
	217	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	218	. "inst(X)", zstohf,zstohf)
	219	c
	220	ENDIF !bb2.EQ."850".or.bb2.EQ."200"
	221	ENDDO !k=1, nlevSTD
	222	c
	223	CALL histdef(nid_hf, "topl",
	224	. "Outgoing longwave radiation", "W/m2",
	225	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	226	. "ave(X)", zstophy,zout)
	227	c
	228	CALL histdef(nid_hf, "precip","Total precipitation rate",
	229	. "kg/(s*m2)",iim,jjmp1,nhori, 1,1,1, -99, 32,
	230	. "ave(X)", zstophy,zout)
	231	c
	232	CALL histdef(nid_hf, "slp",
	233	. "Mean sea-level pressure", "Pa",
	234	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	235	. "inst(X)", zstohf,zout)
	236	c
	237	cIM 280504 BEG
	238	ENDIF !(lev_histhf.GE.1) THEN
	239	c
	240	IF(lev_histhf.GE.2) THEN
	241	c
	242	DO k=1, nlevSTD
	243	c
	244	c bb=clevSTD(k)
	245	c
	246	c IF(k.GE.2) THEN
	247	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
	248	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
	249	c aa=clevSTD(k)
	250	c bb=aa(1:lnblnk1(aa))
	251	c ENDIF
	252	c
	253	IF(bb2.EQ."850".OR.bb2.EQ."500") THEN
	254	c
	255	CALL histdef(nid_hf, "t"//bb2,
	256	. "Temperature "//bb2//"mb","K",
	257	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	258	. "inst(X)", zstohf,zstohf)
	259	c
	260	CALL histdef(nid_hf, "q"//bb2,
	261	. "Specific humidity "//bb2//"mb","kg/kg",
	262	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	263	. "inst(X)", zstohf,zstohf)
	264	c
	265	ENDIF !bb2
	266	c
	267	IF(bb2.EQ."500") THEN
	268	c
	269	CALL histdef(nid_hf, "phi"//bb2,
	270	. "Geopotential "//bb2//"mb", "m",
	271	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	272	. "inst(X)", zstohf,zstohf)
	273	c
	274	CALL histdef(nid_hf, "w"//bb2,
	275	. "Vertical motion "//bb2//"mb","Pa/s",
	276	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	277	. "inst(X)", zstohf,zstohf)
	278	ENDIF !bb2.EQ."500"
	279	c
	280	IF(bb3.EQ."50") THEN
	281	c
	282	CALL histdef(nid_hf, "t"//bb3,
	283	. "Temperature "//bb3//"mb","K",
	284	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	285	. "inst(X)", zstohf,zstohf)
	286	c
	287	CALL histdef(nid_hf, "u"//bb3,
	288	. "Zonal wind "//bb3//"mb","m/s",
	289	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	290	. "inst(X)", zstohf,zstohf)
	291	c
	292	CALL histdef(nid_hf, "v"//bb3,
	293	. "Meridional wind "//bb3//"mb","m/s",
	294	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	295	. "inst(X)", zstohf,zstohf)
	296	c
	297	ENDIF !bb3
	298	c
	299	ENDDO !k=1, nlevSTD
	300	c
	301	CALL histdef(nid_hf, "u10m", "Vent zonal 10m", "m/s",
	302	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	303	. "inst(X)", zstohf,zout)
	304	c
	305	CALL histdef(nid_hf, "v10m", "Vent meridien 10m", "m/s",
	306	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	307	. "inst(X)", zstohf,zout)
	308	c
	309	CALL histdef(nid_hf, "t2m", "Temperature 2m", "K",
	310	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	311	. "inst(X)", zstohf,zout)
	312	c
	313	CALL histdef(nid_hf, "q2m", "Specific humidity", "kg/kg",
	314	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	315	. "inst(X)", zstohf,zout)
	316	c
	317	CALL histdef(nid_hf, "psol", "Surface Pressure", "Pa",
	318	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	319	. "inst(X)", zstohf,zout)
	320	c
	321	CALL histdef(nid_hf, "cldt", "Total cloudiness", "%",
	322	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	323	. "ave(X)", zstohf,zout)
	324	c
	325	CALL histdef(nid_hf, "taux",
	326	$ "Zonal wind stress","Pa",
	327	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	328	$ "ave(X)", zstohf,zout)
	329	c
	330	CALL histdef(nid_hf, "tauy",
	331	$ "Meridional wind stress ", "Pa",
	332	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	333	$ "ave(X)", zstohf,zout)
	334	c
	335	CALL histdef(nid_hf, "prw", "Precipitable water", "kg/m2",
	336	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	337	cIM 090904 . "ave(X)", zstohf,zout)
	338	. "ave(X)", zstophy,zout)
	339	c
	340	CALL histdef(nid_hf, "sens", "Sensible heat flux", "W/m2",
	341	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	342	cIM 090904 . "ave(X)", zstohf,zout)
	343	. "ave(X)", zstophy,zout)
	344	c
	345	CALL histdef(nid_hf, "flat", "Latent heat flux", "W/m2",
	346	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	347	. "ave(X)", zstohf,zout)
	348	c
	349	c CALL histdef(nid_hf, "snow_mass", "Snow Mass", "kg/m2",
	350	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	351	c . "ave(X)", zstohf,zout)
	352	c
	353	CALL histdef(nid_hf, "SWdnSFC", "SWdn at surface","W/m2",
	354	. iim,jjmp1,nhori, 1,1,1,-99,
	355	. 32, "ave(X)", zstohf,zout)
	356	c
	357	CALL histdef(nid_hf, "SWupSFC", "SWup at surface","W/m2",
	358	. iim,jjmp1,nhori, 1,1,1,-99,
	359	. 32, "ave(X)", zstohf,zout)
	360	c
	361	CALL histdef(nid_hf,"LWdnSFC","Down. IR rad. at surface",
	362	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
	363	. "ave(X)", zstohf,zout)
	364	c
	365	CALL histdef(nid_hf,"LWupSFC","Upwd. IR rad. at surface",
	366	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
	367	. "ave(X)", zstohf,zout)
	368	c
	369	CALL histdef(nid_hf, "SWdnTOA", "SWdn at TOA","W/m2",
	370	. iim,jjmp1,nhori, 1,1,1,-99,
	371	. 32, "ave(X)", zstohf,zout)
	372	c
	373	CALL histdef(nid_hf, "SWupTOA", "SWup at TOA","W/m2",
	374	. iim,jjmp1,nhori, 1,1,1,-99,
	375	. 32, "ave(X)", zstohf,zout)
	376	c
	377	DO k=1, nbteta
	378	CALL histdef(nid_hf, "PV"//ctetaSTD(k),
	379	$ "Vorticite potentielle "//ctetaSTD(k)//"K",
	380	$ "1/(Pa*s)",iim,jjmp1,nhori, 1,1,1, -99,
	381	$ 32, "inst(X)", zstohf,zout)
	382	ENDDO !k=1, nbteta
	383	c
	384	c planetary boundary layer height
	385	c
	386	c ENDIF !1.EQ.0
	387	c
	388	ENDIF !(lev_histhf.GE.2) THEN
	389	cIM 280504 END
	390	c
	391	ENDIF !type_run
	392	c
	393	CALL histend(nid_hf)
	394	c
	395	ENDIF ! ok_hf

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

Download in other formats: