Context Navigation

source: LMDZ4/branches/V3_test/libf/phylmd/ini_histhf.h @ 704

Last change on this file since 704 was 704, checked in by Laurent Fairhead, 19 years ago
Inclusion des modifs de Y. Meurdesoif pour la version V3 LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 12.9 KB

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

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

Download in other formats: