Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Normal
Revision Log

ini_histhf.h @ 4999

Last change on this file since 4999 was 912, checked in by (none), 17 years ago
This commit was manufactured by cvs2svn to create tag 'LMDZ4_V3_2'.
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 13.2 KB

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

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

Context Navigation

source: LMDZ4/tags/LMDZ4_V3_2/libf/phylmd/ini_histhf.h @ 4999

Download in other formats: