Context Navigation

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

Last change on this file since 525 was 524, checked in by lmdzadmin, 21 years ago
Initial revision
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 8.8 KB

Rev	Line
[524]	1	!
	2	! $Header$
	3	!
	4	IF (ok_instan) THEN
	5	c
	6	zsto = dtime * ecrit_ins
	7	zout = dtime * ecrit_ins
	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",
	26	. klev, presnivs, nvert)
	27	c call histvert(nid_ins, 'sig_s', 'Niveaux sigma','-',
	28	c . klev, znivsig, nvert)
	29	c
	30	c
	31	CALL histdef(nid_ins, "phis", "Surface geop. height", "-",
	32	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	33	. "once", zsto,zout)
	34	c
	35	CALL histdef(nid_ins, "aire", "Grid area", "-",
	36	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	37	. "once", zsto,zout)
	38	c
	39	c Champs 2D:
	40	c
	41	CALL histdef(nid_ins, "tsol", "Surface Temperature", "K",
	42	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	43	. "inst(X)", zsto,zout)
	44	c
	45	CALL histdef(nid_ins, "t2m", "Temperature 2m", "K",
	46	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	47	. "inst(X)", zsto,zout)
	48	c
	49	CALL histdef(nid_ins, "q2m", "Specific humidity 2m", "Kg/Kg",
	50	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	51	. "inst(X)", zsto,zout)
	52	c
	53	CALL histdef(nid_ins, "u10m", "Vent zonal 10m", "m/s",
	54	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	55	. "inst(X)", zsto,zout)
	56	c
	57	CALL histdef(nid_ins, "v10m", "Vent meridien 10m", "m/s",
	58	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	59	. "inst(X)", zsto,zout)
	60	c
	61	CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa",
	62	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	63	. "inst(X)", zsto,zout)
	64	c
	65	CALL histdef(nid_ins, "plul", "Large-scale Precip.", "mm/day",
	66	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	67	. "inst(X)", zsto,zout)
	68	c
	69	CALL histdef(nid_ins, "pluc", "Convective Precip.", "mm/day",
	70	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	71	. "inst(X)", zsto,zout)
	72	c
	73	CALL histdef(nid_ins, "cdrm", "Momentum drag coef.", "-",
	74	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	75	. "inst(X)", zsto,zout)
	76	c
	77	CALL histdef(nid_ins, "cdrh", "Heat drag coef.", "-",
	78	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	79	. "inst(X)", zsto,zout)
	80	c
	81	CALL histdef(nid_ins, "precip", "Precipitation Totale liq+sol",
	82	. "kg/(s*m2)",
	83	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	84	. "inst(X)", zsto,zout)
	85	c
	86	CALL histdef(nid_ins, "snow", "Snow fall", "kg/(s*m2)",
	87	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	88	. "inst(X)", zsto,zout)
	89	c
	90	CALL histdef(nid_ins, "snow_mass", "Snow Mass", "kg/m2",
	91	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	92	. "inst(X)", zsto,zout)
	93	c
	94	CALL histdef(nid_ins, "topl", "OLR", "W/m2",
	95	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	96	. "inst(X)", zsto,zout)
	97	c
	98	CALL histdef(nid_ins, "evap", "Evaporation", "kg/(s*m2)",
	99	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	100	. "inst(X)", zsto,zout)
	101	c
	102	CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2",
	103	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	104	. "inst(X)", zsto,zout)
	105	c
	106	CALL histdef(nid_ins, "soll", "IR rad. at surface", "W/m2",
	107	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	108	. "inst(X)", zsto,zout)
	109	c
	110	CALL histdef(nid_ins, "solldown", "Down. IR rad. at surface",
	111	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
	112	. "inst(X)", zsto,zout)
	113	c
	114	CALL histdef(nid_ins, "bils", "Surf. total heat flux", "W/m2",
	115	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	116	. "inst(X)", zsto,zout)
	117	c
	118	CALL histdef(nid_ins, "sens", "Sensible heat flux", "W/m2",
	119	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	120	. "inst(X)", zsto,zout)
	121	c
	122	CALL histdef(nid_ins, "fder", "Heat flux derivation", "W/m2",
	123	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	124	. "inst(X)", zsto,zout)
	125	c
	126	CALL histdef(nid_ins, "dtsvdfo", "Boundary-layer dTs(o)", "K/s",
	127	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	128	. "inst(X)", zsto,zout)
	129	c
	130	CALL histdef(nid_ins, "dtsvdft", "Boundary-layer dTs(t)", "K/s",
	131	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	132	. "inst(X)", zsto,zout)
	133	c
	134	CALL histdef(nid_ins, "dtsvdfg", "Boundary-layer dTs(g)", "K/s",
	135	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	136	. "inst(X)", zsto,zout)
	137	c
	138	CALL histdef(nid_ins, "dtsvdfi", "Boundary-layer dTs(g)", "K/s",
	139	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	140	. "inst(X)", zsto,zout)
	141
	142	DO nsrf = 1, nbsrf
	143	C
	144	call histdef(nid_ins, "pourc_"//clnsurf(nsrf),
	145	$ "Fraction"//clnsurf(nsrf), "W/m2",
	146	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	147	$ "inst(X)", zsto,zout)
	148
	149	call histdef(nid_ins, "sens_"//clnsurf(nsrf),
	150	$ "Sensible heat flux "//clnsurf(nsrf), "W/m2",
	151	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	152	$ "inst(X)", zsto,zout)
	153	c
	154	call histdef(nid_ins, "tsol_"//clnsurf(nsrf),
	155	$ "Surface Temperature"//clnsurf(nsrf), "W/m2",
	156	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	157	$ "inst(X)", zsto,zout)
	158	c
	159	call histdef(nid_ins, "lat_"//clnsurf(nsrf),
	160	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
	161	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	162	$ "inst(X)", zsto,zout)
	163	C
	164	call histdef(nid_ins, "taux_"//clnsurf(nsrf),
	165	$ "Zonal wind stress"//clnsurf(nsrf),"Pa",
	166	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	167	$ "inst(X)", zsto,zout)
	168
	169	call histdef(nid_ins, "tauy_"//clnsurf(nsrf),
	170	$ "Meridional xind stress "//clnsurf(nsrf), "Pa",
	171	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	172	$ "inst(X)", zsto,zout)
	173	c
	174	call histdef(nid_ins, "albe_"//clnsurf(nsrf),
	175	$ "Albedo "//clnsurf(nsrf), "-",
	176	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	177	$ "inst(X)", zsto,zout)
	178	c
	179	call histdef(nid_ins, "rugs_"//clnsurf(nsrf),
	180	$ "rugosite "//clnsurf(nsrf), "-",
	181	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	182	$ "inst(X)", zsto,zout)
	183	CXXX
	184	END DO
	185	CALL histdef(nid_ins, "rugs", "rugosity", "-",
	186	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	187	. "inst(X)", zsto,zout)
	188
	189	c
	190	CALL histdef(nid_ins, "albs", "Surface albedo", "-",
	191	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	192	. "inst(X)", zsto,zout)
	193	CALL histdef(nid_ins, "albslw", "Surface albedo LW", "-",
	194	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	195	. "inst(X)", zsto,zout)
	196	c
	197	c
	198	c Champs 3D:
	199	c
	200	CALL histdef(nid_ins, "temp", "Temperature", "K",
	201	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	202	. "inst(X)", zsto,zout)
	203	c
	204	CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s",
	205	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	206	. "inst(X)", zsto,zout)
	207	c
	208	CALL histdef(nid_ins, "vitv", "Merid wind", "m/s",
	209	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	210	. "inst(X)", zsto,zout)
	211	c
	212	CALL histdef(nid_ins, "geop", "Geopotential height", "m",
	213	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	214	. "inst(X)", zsto,zout)
	215	c
	216	CALL histdef(nid_ins, "pres", "Air pressure", "Pa",
	217	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	218	. "inst(X)", zsto,zout)
	219	c
	220	CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s",
	221	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	222	. "inst(X)", zsto,zout)
	223	c
	224	CALL histdef(nid_ins, "dqvdf", "Boundary-layer dQ", "Kg/Kg/s",
	225	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	226	. "inst(X)", zsto,zout)
	227	c
	228
	229	CALL histend(nid_ins)
	230	c
	231	ndex2d = 0
	232	ndex3d = 0
	233	c
	234	ENDIF

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

Download in other formats: