Context Navigation

source: LMDZ4/trunk/libf/phylmd/ini_histins.h @ 758

Last change on this file since 758 was 684, checked in by lmdzadmin, 19 years ago
Ajout nouveaux flags (ip_ebil_phy, ok_slab_sicOBS) ds physiq.def + Lecture freq. ecriture en nombre de jours IM
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 10.8 KB

Rev	Line
[524]	1	!
	2	! $Header$
	3	!
	4	IF (ok_instan) THEN
	5	c
[684]	6	zsto = ecrit_ins
	7	zout = ecrit_ins
[524]	8	c
	9	idayref = day_ref
	10	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
	11	c
	12	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
	13	DO i = 1, iim
	14	zx_lon(i,1) = rlon(i+1)
	15	zx_lon(i,jjmp1) = rlon(i+1)
	16	ENDDO
	17	DO ll=1,klev
	18	znivsig(ll)=float(ll)
	19	ENDDO
	20	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
	21	CALL histbeg("histins", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
	22	. 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
	23	. nhori, nid_ins)
	24	write(,)'Inst ', itau_phy, zjulian
	25	CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb",
[644]	26	cIM . klev, presnivs, nvert)
	27	. klev, presnivs/100., nvert)
[524]	28	c call histvert(nid_ins, 'sig_s', 'Niveaux sigma','-',
	29	c . klev, znivsig, nvert)
	30	c
	31	c
	32	CALL histdef(nid_ins, "phis", "Surface geop. height", "-",
	33	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	34	. "once", zsto,zout)
	35	c
	36	CALL histdef(nid_ins, "aire", "Grid area", "-",
	37	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	38	. "once", zsto,zout)
	39	c
	40	c Champs 2D:
	41	c
	42	CALL histdef(nid_ins, "tsol", "Surface Temperature", "K",
	43	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	44	. "inst(X)", zsto,zout)
	45	c
	46	CALL histdef(nid_ins, "t2m", "Temperature 2m", "K",
	47	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	48	. "inst(X)", zsto,zout)
	49	c
	50	CALL histdef(nid_ins, "q2m", "Specific humidity 2m", "Kg/Kg",
	51	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	52	. "inst(X)", zsto,zout)
	53	c
	54	CALL histdef(nid_ins, "u10m", "Vent zonal 10m", "m/s",
	55	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	56	. "inst(X)", zsto,zout)
	57	c
	58	CALL histdef(nid_ins, "v10m", "Vent meridien 10m", "m/s",
	59	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	60	. "inst(X)", zsto,zout)
	61	c
	62	CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa",
	63	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	64	. "inst(X)", zsto,zout)
	65	c
	66	CALL histdef(nid_ins, "plul", "Large-scale Precip.", "mm/day",
	67	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	68	. "inst(X)", zsto,zout)
	69	c
	70	CALL histdef(nid_ins, "pluc", "Convective Precip.", "mm/day",
	71	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	72	. "inst(X)", zsto,zout)
	73	c
	74	CALL histdef(nid_ins, "cdrm", "Momentum drag coef.", "-",
	75	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	76	. "inst(X)", zsto,zout)
	77	c
	78	CALL histdef(nid_ins, "cdrh", "Heat drag coef.", "-",
	79	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	80	. "inst(X)", zsto,zout)
	81	c
	82	CALL histdef(nid_ins, "precip", "Precipitation Totale liq+sol",
	83	. "kg/(s*m2)",
	84	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	85	. "inst(X)", zsto,zout)
	86	c
	87	CALL histdef(nid_ins, "snow", "Snow fall", "kg/(s*m2)",
	88	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	89	. "inst(X)", zsto,zout)
	90	c
[644]	91	c CALL histdef(nid_ins, "snow_mass", "Snow Mass", "kg/m2",
	92	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	93	c . "inst(X)", zsto,zout)
[524]	94	c
	95	CALL histdef(nid_ins, "topl", "OLR", "W/m2",
	96	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	97	. "inst(X)", zsto,zout)
	98	c
	99	CALL histdef(nid_ins, "evap", "Evaporation", "kg/(s*m2)",
	100	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	101	. "inst(X)", zsto,zout)
	102	c
	103	CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2",
	104	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	105	. "inst(X)", zsto,zout)
	106	c
	107	CALL histdef(nid_ins, "soll", "IR rad. at surface", "W/m2",
	108	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	109	. "inst(X)", zsto,zout)
	110	c
	111	CALL histdef(nid_ins, "solldown", "Down. IR rad. at surface",
	112	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
	113	. "inst(X)", zsto,zout)
	114	c
	115	CALL histdef(nid_ins, "bils", "Surf. total heat flux", "W/m2",
	116	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	117	. "inst(X)", zsto,zout)
	118	c
	119	CALL histdef(nid_ins, "sens", "Sensible heat flux", "W/m2",
	120	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	121	. "inst(X)", zsto,zout)
	122	c
	123	CALL histdef(nid_ins, "fder", "Heat flux derivation", "W/m2",
	124	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	125	. "inst(X)", zsto,zout)
	126	c
	127	CALL histdef(nid_ins, "dtsvdfo", "Boundary-layer dTs(o)", "K/s",
	128	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	129	. "inst(X)", zsto,zout)
	130	c
	131	CALL histdef(nid_ins, "dtsvdft", "Boundary-layer dTs(t)", "K/s",
	132	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	133	. "inst(X)", zsto,zout)
	134	c
	135	CALL histdef(nid_ins, "dtsvdfg", "Boundary-layer dTs(g)", "K/s",
	136	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	137	. "inst(X)", zsto,zout)
	138	c
	139	CALL histdef(nid_ins, "dtsvdfi", "Boundary-layer dTs(g)", "K/s",
	140	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	141	. "inst(X)", zsto,zout)
	142
	143	DO nsrf = 1, nbsrf
	144	C
	145	call histdef(nid_ins, "pourc_"//clnsurf(nsrf),
[644]	146	$ "% "//clnsurf(nsrf), "%",
[524]	147	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	148	$ "inst(X)", zsto,zout)
	149
[644]	150	call histdef(nid_ins, "fract_"//clnsurf(nsrf),
	151	$ "Fraction "//clnsurf(nsrf), "1",
	152	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	153	$ "inst(X)", zsto,zout)
	154
[524]	155	call histdef(nid_ins, "sens_"//clnsurf(nsrf),
	156	$ "Sensible heat flux "//clnsurf(nsrf), "W/m2",
	157	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	158	$ "inst(X)", zsto,zout)
	159	c
	160	call histdef(nid_ins, "tsol_"//clnsurf(nsrf),
	161	$ "Surface Temperature"//clnsurf(nsrf), "W/m2",
	162	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	163	$ "inst(X)", zsto,zout)
	164	c
	165	call histdef(nid_ins, "lat_"//clnsurf(nsrf),
	166	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
	167	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	168	$ "inst(X)", zsto,zout)
	169	C
	170	call histdef(nid_ins, "taux_"//clnsurf(nsrf),
	171	$ "Zonal wind stress"//clnsurf(nsrf),"Pa",
	172	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	173	$ "inst(X)", zsto,zout)
	174
	175	call histdef(nid_ins, "tauy_"//clnsurf(nsrf),
	176	$ "Meridional xind stress "//clnsurf(nsrf), "Pa",
	177	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	178	$ "inst(X)", zsto,zout)
	179	c
	180	call histdef(nid_ins, "albe_"//clnsurf(nsrf),
	181	$ "Albedo "//clnsurf(nsrf), "-",
	182	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	183	$ "inst(X)", zsto,zout)
	184	c
	185	call histdef(nid_ins, "rugs_"//clnsurf(nsrf),
	186	$ "rugosite "//clnsurf(nsrf), "-",
	187	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	188	$ "inst(X)", zsto,zout)
	189	CXXX
	190	END DO
	191	CALL histdef(nid_ins, "rugs", "rugosity", "-",
	192	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	193	. "inst(X)", zsto,zout)
	194
	195	c
	196	CALL histdef(nid_ins, "albs", "Surface albedo", "-",
	197	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	198	. "inst(X)", zsto,zout)
	199	CALL histdef(nid_ins, "albslw", "Surface albedo LW", "-",
	200	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	201	. "inst(X)", zsto,zout)
	202	c
	203	c
[644]	204	cIM cf. AM 081204 BEG
	205	c HBTM2
	206	CALL histdef(nid_ins, "s_pblh", "Boundary Layer Height", "m",
	207	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	208	. "inst(X)", zsto,zout)
	209	c
	210	CALL histdef(nid_ins, "s_pblt", "T at Boundary Layer Height",
	211	. "K",
	212	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	213	. "inst(X)", zsto,zout)
	214	c
	215	CALL histdef(nid_ins, "s_lcl", "Condensation level", "m",
	216	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	217	. "inst(X)", zsto,zout)
	218	c
	219	CALL histdef(nid_ins, "s_capCL", "Conv avlbl pot ener for
	220	. ABL", "J/m2",
	221	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	222	. "inst(X)", zsto,zout)
	223	c
	224	CALL histdef(nid_ins, "s_oliqCL", "Liq Water in BL", "kg/m2",
	225	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	226	. "inst(X)", zsto,zout)
	227	c
	228	CALL histdef(nid_ins, "s_cteiCL", "Instability criteria
	229	. (ABL)", "K",
	230	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	231	. "inst(X)", zsto,zout)
	232	c
	233	CALL histdef(nid_ins, "s_therm", "Exces du thermique", "K",
	234	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	235	. "inst(X)", zsto,zout)
	236	c
	237	CALL histdef(nid_ins, "s_trmb1", "deep_cape(HBTM2)", "J/m2",
	238	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	239	. "inst(X)", zsto,zout)
	240	c
	241	CALL histdef(nid_ins, "s_trmb2", "inhibition (HBTM2)", "J/m2",
	242	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	243	. "inst(X)", zsto,zout)
	244	c
	245	CALL histdef(nid_ins, "s_trmb3", "Point Omega (HBTM2)", "m",
	246	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	247	. "inst(X)", zsto,zout)
	248	c
	249	cIM cf. AM 081204 END
	250	c
[524]	251	c Champs 3D:
	252	c
	253	CALL histdef(nid_ins, "temp", "Temperature", "K",
	254	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	255	. "inst(X)", zsto,zout)
	256	c
	257	CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s",
	258	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	259	. "inst(X)", zsto,zout)
	260	c
	261	CALL histdef(nid_ins, "vitv", "Merid wind", "m/s",
	262	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	263	. "inst(X)", zsto,zout)
	264	c
	265	CALL histdef(nid_ins, "geop", "Geopotential height", "m",
	266	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	267	. "inst(X)", zsto,zout)
	268	c
	269	CALL histdef(nid_ins, "pres", "Air pressure", "Pa",
	270	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	271	. "inst(X)", zsto,zout)
	272	c
	273	CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s",
	274	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	275	. "inst(X)", zsto,zout)
	276	c
	277	CALL histdef(nid_ins, "dqvdf", "Boundary-layer dQ", "Kg/Kg/s",
	278	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	279	. "inst(X)", zsto,zout)
	280	c
	281
	282	CALL histend(nid_ins)
	283	c
	284	ndex2d = 0
	285	ndex3d = 0
	286	c
	287	ENDIF

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

Download in other formats: