Context Navigation

source: LMDZ4/trunk/libf/phylmd/ini_bilKP_ave.h @ 715

Last change on this file since 715 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:executable set to ``* Property svn:keywords set to `Author Date Id Revision`
File size: 8.6 KB

Rev	Line
[644]	1	c
	2	c $Header$
	3	c
	4	IF (ok_journe) THEN
	5	c
	6	zsto = dtime
[684]	7	zout = ecrit_day
[644]	8	typeval=tave
	9	c
	10	idayref = day_ref
	11	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
	12	c
	13	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
	14	DO i = 1, iim
	15	zx_lon(i,1) = rlon(i+1)
	16	zx_lon(i,jjmp1) = rlon(i+1)
	17	ENDDO
	18	DO ll=1,klev
	19	znivsig(ll)=float(ll)
	20	ENDDO
	21	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
	22	write(,)'zx_lon = ',zx_lon(:,1)
	23	write(,)'zx_lat = ',zx_lat(1,:)
	24	CALL histbeg("histbilKP_ave", iim,zx_lon(:,1), jjmp1,
	25	. zx_lat(1,:),
	26	. 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
	27	. nhori, nid_bilKPave)
	28	write(,)'Journee ', itau_phy, zjulian
	29	CALL histvert(nid_bilKPave, "presnivs",
	30	. "Vertical levels","mb",
	31	. klev, presnivs/100., nvert)
	32	c
	33	c
	34	c Champs 3D:
	35	c
	36	CALL histdef(nid_bilKPave,"ue",
	37	. "Zonal energy transport","-",
	38	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	39	. typeval, zsto,zout)
	40	c
	41	CALL histdef(nid_bilKPave,"ve",
	42	. "Merid energy transport","-",
	43	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	44	. typeval, zsto,zout)
	45	c
	46	CALL histdef(nid_bilKPave,"uq",
	47	. "Zonal humidity transport","-",
	48	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	49	. typeval, zsto,zout)
	50	c
	51	CALL histdef(nid_bilKPave,"vq",
	52	. "Merid humidity transport","-",
	53	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	54	. typeval, zsto,zout)
	55	c
	56	c Champs 3D:
	57	c
	58	CALL histdef(nid_bilKPave,"temp",
	59	. "Air temperature","K",
	60	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	61	. typeval, zsto,zout)
	62	c
	63	CALL histdef(nid_bilKPave,"ovap",
	64	. "Specific humidity","Kg/Kg",
	65	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	66	. typeval, zsto,zout)
	67	c
	68	CALL histdef(nid_bilKPave,"geop",
	69	. "Geopotential height","m",
	70	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	71	. typeval, zsto,zout)
	72	c
	73	CALL histdef(nid_bilKPave,"vitu",
	74	. "Zonal wind","m/s",
	75	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	76	. typeval, zsto,zout)
	77	c
	78	CALL histdef(nid_bilKPave,"vitv",
	79	. "Meridional wind","m/s",
	80	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	81	. typeval, zsto,zout)
	82	c
	83	CALL histdef(nid_bilKPave, "vitw",
	84	. "Vertical wind", "m/s",
	85	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	86	. typeval, zsto,zout)
	87	c
	88	CALL histdef(nid_bilKPave, "pres",
	89	. "Inter-Layer Air pressure",
	90	. "Pa",
	91	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	92	. typeval, zsto,zout)
	93	c
	94	CALL histdef(nid_bilKPave, "play",
	95	. "Mean-Layer Air pressure",
	96	. "Pa",
	97	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	98	. typeval, zsto,zout)
	99	c
	100	CALL histdef(nid_bilKPave, "oliq",
	101	. "Liquid water content",
	102	. "kg/kg",
	103	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	104	. typeval, zsto,zout)
	105	c
	106	CALL histdef(nid_bilKPave, "dtdyn",
	107	. "Dynamics dT", "K/s",
	108	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	109	. typeval, zsto,zout)
	110	c
	111	CALL histdef(nid_bilKPave, "dqdyn",
	112	. "Dynamics dQ", "Kg/Kg/s",
	113	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	114	. typeval, zsto,zout)
	115	c
	116	CALL histdef(nid_bilKPave, "dtcon",
	117	. "Convection dT", "K/s",
	118	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	119	. typeval, zsto,zout)
	120	c
	121	CALL histdef(nid_bilKPave, "ducon",
	122	. "Convection du", "m/s2",
	123	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	124	. typeval, zsto,zout)
	125	c
	126	CALL histdef(nid_bilKPave, "dvcon",
	127	. "Convection dv", "m/s2",
	128	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	129	. typeval, zsto,zout)
	130	c
	131	CALL histdef(nid_bilKPave,"dqcon",
	132	. "Convection dQ","Kg/Kg/s",
	133	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	134	. typeval, zsto,zout)
	135	c
	136	CALL histdef(nid_bilKPave, "dtlsc",
	137	. "Condensation dT", "K/s",
	138	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	139	. typeval, zsto,zout)
	140	c
	141	CALL histdef(nid_bilKPave,"dqlsc",
	142	. "Condensation dQ","Kg/Kg/s",
	143	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	144	. typeval, zsto,zout)
	145	c
	146	CALL histdef(nid_bilKPave,"dtvdf",
	147	. "Boundary-layer dT","K/s",
	148	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	149	. typeval, zsto,zout)
	150	c
	151	CALL histdef(nid_bilKPave, "dqvdf",
	152	. "Boundary-layer dQ",
	153	. "Kg/Kg/s",
	154	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	155	. typeval, zsto,zout)
	156	c
	157	CALL histdef(nid_bilKPave,"dtajs",
	158	. "Ajustement sec dT","K/s",
	159	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	160	. typeval, zsto,zout)
	161	c
	162	CALL histdef(nid_bilKPave, "dqajs",
	163	. "Ajustement sec dQ",
	164	. "Kg/Kg/s",
	165	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	166	. typeval, zsto,zout)
	167	c
	168	CALL histdef(nid_bilKPave, "dteva",
	169	. "Reevaporation dT", "K/s",
	170	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	171	. typeval, zsto,zout)
	172	c
	173	CALL histdef(nid_bilKPave,"dqeva",
	174	. "Reevaporation dQ",
	175	. "Kg/Kg/s",
	176	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	177	. typeval, zsto,zout)
	178
	179	c
	180	CALL histdef(nid_bilKPave, "dtswr",
	181	. "SW radiation dT", "K/s",
	182	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	183	. typeval, zsto,zout)
	184	c
	185	CALL histdef(nid_bilKPave, "dtsw0",
	186	. "SW radiation dT", "K/s",
	187	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	188	. typeval, zsto,zout)
	189	c
	190	CALL histdef(nid_bilKPave, "dtlwr",
	191	. "LW radiation dT", "K/s",
	192	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	193	. typeval, zsto,zout)
	194	c
	195	CALL histdef(nid_bilKPave, "dtlw0",
	196	. "LW radiation dT", "K/s",
	197	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	198	. typeval, zsto,zout)
	199	c
	200	CALL histdef(nid_bilKPave,"duvdf",
	201	. "Boundary-layer dU","m/s2",
	202	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	203	. typeval, zsto,zout)
	204	c
	205	CALL histdef(nid_bilKPave,"dvvdf",
	206	. "Boundary-layer dV","m/s2",
	207	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	208	. typeval, zsto,zout)
	209	c
	210	IF (ok_orodr) THEN
	211	IF (ok_orolf) THEN
	212	CALL histdef(nid_bilKPave, "duoli",
	213	. "Orography dU", "m/s2",
	214	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	215	. typeval, zsto,zout)
	216	c
	217	CALL histdef(nid_bilKPave, "dvoli",
	218	. "Orography dV", "m/s2",
	219	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	220	. typeval, zsto,zout)
	221	c
	222	ENDIF
	223	ENDIF
	224	C
	225	CALL histdef(nid_bilKPave, "duphy",
	226	. "Physiq dU","-",
	227	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	228	. typeval, zsto,zout)
	229	C
	230	CALL histdef(nid_bilKPave, "dvphy",
	231	. "Physiq dV","-",
	232	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	233	. typeval, zsto,zout)
	234	C
	235	CALL histdef(nid_bilKPave, "dtphy",
	236	. "Physiq dT","-",
	237	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	238	. typeval, zsto,zout)
	239	C
	240	CALL histdef(nid_bilKPave, "dqphy",
	241	. "Physiq dQ","-",
	242	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	243	. typeval, zsto,zout)
	244	C
	245	CALL histdef(nid_bilKPave, "dqlphy",
	246	. "Physiq dQl","-",
	247	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	248	. typeval, zsto,zout)
	249	C
	250	C
	251	CALL histend(nid_bilKPave)
	252	c
	253	ndex2d = 0
	254	ndex3d = 0
	255	c
	256	ENDIF ! fin de test sur ok_journe

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

Download in other formats: