Context Navigation

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

Last change on this file since 648 was 644, checked in by Laurent Fairhead, 20 years ago
Synchronisation avec tous les diagnostiques de Ionela IM Inclusion du slab ocean IM LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 10.9 KB

Line
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	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: