Context Navigation

source: LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histins.h @ 444

Last change on this file since 444 was 418, checked in by (none), 22 years ago
This commit was manufactured by cvs2svn to create branch 'rel-LF'.
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 8.9 KB

Rev	Line
[416]	1	IF (ok_instan) THEN
	2	c
	3	idayref = day_ref
	4	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
	5	c
	6	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
	7	DO i = 1, iim
	8	zx_lon(i,1) = rlon(i+1)
	9	zx_lon(i,jjmp1) = rlon(i+1)
	10	ENDDO
	11	DO ll=1,klev
	12	znivsig(ll)=float(ll)
	13	ENDDO
	14	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
	15	CALL histbeg("histins", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
	16	. 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
	17	. nhori, nid_ins)
	18	write(,)'Inst ', itau_phy, zjulian
	19	CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb",
	20	. klev, presnivs, nvert)
	21	c call histvert(nid_ins, 'sig_s', 'Niveaux sigma','-',
	22	c . klev, znivsig, nvert)
	23	c
	24	c
	25	zsto = dtime * ecrit_ins
	26	zout = dtime * ecrit_ins
	27	C
	28	CALL histdef(nid_ins, "phis", "Surface geop. height", "-",
	29	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	30	. "once", zsto,zout)
	31	c
	32	CALL histdef(nid_ins, "aire", "Grid area", "-",
	33	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	34	. "once", zsto,zout)
	35	c
	36	c Champs 2D:
	37	c
	38	CALL histdef(nid_ins, "tsol", "Surface Temperature", "K",
	39	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	40	. "inst(X)", zsto,zout)
	41	c
	42	CALL histdef(nid_ins, "t2m", "Temperature 2m", "K",
	43	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	44	. "inst(X)", zsto,zout)
	45	c
	46	CALL histdef(nid_ins, "q2m", "Specific humidity 2m", "Kg/Kg",
	47	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	48	. "inst(X)", zsto,zout)
	49	c
	50	CALL histdef(nid_ins, "u10m", "Vent zonal 10m", "m/s",
	51	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	52	. "inst(X)", zsto,zout)
	53	c
	54	CALL histdef(nid_ins, "v10m", "Vent meridien 10m", "m/s",
	55	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	56	. "inst(X)", zsto,zout)
	57	c
	58	CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa",
	59	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	60	. "inst(X)", zsto,zout)
	61	c
	62	CALL histdef(nid_ins, "plul", "Large-scale Precip.", "mm/day",
	63	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	64	. "inst(X)", zsto,zout)
	65	c
	66	CALL histdef(nid_ins, "pluc", "Convective Precip.", "mm/day",
	67	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	68	. "inst(X)", zsto,zout)
	69
	70	CALL histdef(nid_ins, "qsol", "Surface humidity", "mm",
	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",
	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",
	88	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	89	. "inst(X)", zsto,zout)
	90	c
	91	CALL histdef(nid_ins, "snow_mass", "Snow Mass", "kg/m2",
	92	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	93	. "inst(X)", zsto,zout)
	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",
	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),
	146	$ "Fraction"//clnsurf(nsrf), "W/m2",
	147	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	148	$ "inst(X)", zsto,zout)
	149
	150	call histdef(nid_ins, "sens_"//clnsurf(nsrf),
	151	$ "Sensible heat flux "//clnsurf(nsrf), "W/m2",
	152	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	153	$ "inst(X)", zsto,zout)
	154	c
	155	call histdef(nid_ins, "tsol_"//clnsurf(nsrf),
	156	$ "Surface Temperature"//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, "lat_"//clnsurf(nsrf),
	161	$ "Latent heat flux "//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, "taux_"//clnsurf(nsrf),
	166	$ "Zonal wind stress"//clnsurf(nsrf),"Pa",
	167	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	168	$ "inst(X)", zsto,zout)
	169
	170	call histdef(nid_ins, "tauy_"//clnsurf(nsrf),
	171	$ "Meridional xind stress "//clnsurf(nsrf), "Pa",
	172	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	173	$ "inst(X)", zsto,zout)
	174	c
	175	call histdef(nid_ins, "albe_"//clnsurf(nsrf),
	176	$ "Albedo "//clnsurf(nsrf), "-",
	177	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	178	$ "inst(X)", zsto,zout)
	179	c
	180	call histdef(nid_ins, "rugs_"//clnsurf(nsrf),
	181	$ "rugosite "//clnsurf(nsrf), "-",
	182	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
	183	$ "inst(X)", zsto,zout)
	184	CXXX
	185	END DO
	186	CALL histdef(nid_ins, "rugs", "rugosity", "-",
	187	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	188	. "inst(X)", zsto,zout)
	189
	190	c
	191	CALL histdef(nid_ins, "albs", "Surface albedo", "-",
	192	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	193	. "inst(X)", zsto,zout)
	194	CALL histdef(nid_ins, "albslw", "Surface albedo LW", "-",
	195	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	196	. "inst(X)", zsto,zout)
	197	c
	198	c
	199	c Champs 3D:
	200	c
	201	CALL histdef(nid_ins, "temp", "Temperature", "K",
	202	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	203	. "inst(X)", zsto,zout)
	204	c
	205	CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s",
	206	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	207	. "inst(X)", zsto,zout)
	208	c
	209	CALL histdef(nid_ins, "vitv", "Merid wind", "m/s",
	210	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	211	. "inst(X)", zsto,zout)
	212	c
	213	CALL histdef(nid_ins, "geop", "Geopotential height", "m",
	214	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	215	. "inst(X)", zsto,zout)
	216	c
	217	CALL histdef(nid_ins, "pres", "Air pressure", "Pa",
	218	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	219	. "inst(X)", zsto,zout)
	220	c
	221	CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s",
	222	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	223	. "inst(X)", zsto,zout)
	224	c
	225	CALL histdef(nid_ins, "dqvdf", "Boundary-layer dQ", "Kg/Kg/s",
	226	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	227	. "inst(X)", zsto,zout)
	228	c
	229
	230	CALL histend(nid_ins)
	231	c
	232	ndex2d = 0
	233	ndex3d = 0
	234	c
	235	ENDIF

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

Download in other formats: