Context Navigation

source: LMDZ4/branches/pre_V3/libf/phylmd/ini_histins.h @ 4096

Last change on this file since 4096 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

Line
1	!
2	! $Header$
3	!
4	IF (ok_instan) THEN
5	c
6	zsto = ecrit_ins
7	zout = 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	cIM . klev, presnivs, nvert)
27	. klev, presnivs/100., nvert)
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
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)
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),
146	$ "% "//clnsurf(nsrf), "%",
147	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
148	$ "inst(X)", zsto,zout)
149
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
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
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
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: