Context Navigation

module_diagnostics.f90 @ 220

Last change on this file since 220 was 207, checked in by aslmd, 14 years ago
MESOSCALE: A GENERAL CLEAN-UP FOLLOWING UPDATING THE USER MANUAL. EVERYTHING ESSENTIAL IS IN MESOSCALE (much lighter than before). EVERYTHING FOR DEVELOPPERS OR EXPERTS IS IN MESOSCALE_DEV.
File size: 20.6 KB

Line
1	!! Diagnostics
2	!! See moddule_diagnostic for which fields are available
3
4	MODULE module_diagnostics
5
6	USE output_module
7	USE gridinfo_module
8	USE module_interp
9	USE module_arrays
10	USE constants_module
11	USE module_pressure
12
13
14	USE module_calc_pressure
15	USE module_calc_height
16	USE module_calc_theta
17	USE module_calc_tk
18	USE module_calc_tc
19	USE module_calc_td
20	USE module_calc_rh
21	USE module_calc_uvmet
22	USE module_calc_wdir
23	USE module_calc_wspd
24	USE module_calc_slp
25	USE module_calc_dbz
26	USE module_calc_cape
27	USE module_calc_clfr
28
29	CONTAINS
30
31	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
32	! Name: process_diagnostics
33	! Purpose: All new calls to diagnostics routines go here
34	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
35	SUBROUTINE process_diagnostics ()
36
37	IMPLICIT NONE
38
39	! Local variables
40	character (len=128) :: cname, cunits, cdesc, c_nm
41	real, allocatable, dimension(:,:) :: SCR2
42	real, pointer, dimension(:,:,:) :: SCR3
43	real, allocatable, dimension(:,:,:,:) :: SCR4
44	real, pointer, dimension(:,:,:) :: SCR
45	real, pointer, dimension(:,:,:) :: data_out
46	integer :: nxout, nyout, nzout
47	integer :: ii, jj, kk
48
49
50
51	!!! If we don't have PRES (Pa) already - calculate and keep
52	IF ( .not. have_PRES ) THEN
53	IF ( iprogram == 8 .AND. have_P .AND. have_PB ) THEN
54	ALLOCATE(SCR3(west_east_dim,south_north_dim,bottom_top_dim))
55	SCR3 = P + PB
56	c_nm = 'PRES'
57	CALL keep_arrays(c_nm, SCR3)
58	END IF
59	IF ( iprogram == 6 .AND. have_MU .AND. have_MUB .AND. have_ZNU .AND. have_ZNW .AND. have_PTOP) THEN
60	CALL pressure(SCR3)
61	c_nm = 'PRES'
62	CALL keep_arrays(c_nm, SCR3)
63	END IF
64	DEALLOCATE(SCR3)
65	END IF
66
67
68	!!! If we don't have TK (K) already - calculate and keep
69	IF ( .not. have_TK ) THEN
70	IF ( have_T .AND. have_PRES ) THEN
71	ALLOCATE(SCR3(west_east_dim,south_north_dim,bottom_top_dim))
72	c_nm = 'TK'
73	!! SCR3 = (T+300.) * ( PRES / p0 )**RCP
74	SCR3 = (T+220.) * ( PRES / p0 )**RCP
75	CALL keep_arrays(c_nm, SCR3)
76	DEALLOCATE(SCR3)
77	END IF
78	END IF
79
80
81
82	!!
83	!! DIAGNOSTICS GO HERE
84	!!
85
86
87
88	!!! Calculate pressure in hPA
89	IF ( INDEX(plot_these_fields,",pressure,") /= 0 ) THEN
90	IF ( have_PRES ) THEN
91	IF ( ASSOCIATED(SCR) ) DEALLOCATE(SCR)
92	ALLOCATE(SCR(west_east_dim,south_north_dim,bottom_top_dim))
93	CALL calc_pressure(SCR, cname, cdesc, cunits)
94	CALL interp (SCR, west_east_dim, south_north_dim, bottom_top_dim, &
95	data_out, nxout, nyout, nzout, &
96	vert_array, interp_levels, number_of_zlevs)
97	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
98	DEALLOCATE(SCR)
99	END IF
100	END IF
101
102
103	!!! Calculate height from PH and PHB
104	IF ( INDEX(plot_these_fields,",height,") /= 0 ) THEN
105	IF ( have_PH .AND. have_PHB ) THEN
106	IF ( ASSOCIATED(SCR) ) DEALLOCATE(SCR)
107	ALLOCATE(SCR(west_east_dim,south_north_dim,bottom_top_dim))
108	CALL calc_height(SCR, cname, cdesc, cunits)
109	CALL interp (SCR, west_east_dim, south_north_dim, bottom_top_dim, &
110	data_out, nxout, nyout, nzout, &
111	vert_array, interp_levels, number_of_zlevs)
112	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
113	DEALLOCATE(SCR)
114	END IF
115	END IF
116
117
118	!!! Calculate temperature
119	IF ( INDEX(plot_these_fields,",tk,") /= 0 ) THEN
120	IF ( have_TK ) THEN
121	IF ( ASSOCIATED(SCR) ) DEALLOCATE(SCR)
122	ALLOCATE(SCR(west_east_dim,south_north_dim,bottom_top_dim))
123	CALL calc_tk(SCR, cname, cdesc, cunits)
124	CALL interp (SCR, west_east_dim, south_north_dim, bottom_top_dim, &
125	data_out, nxout, nyout, nzout, &
126	vert_array, interp_levels, number_of_zlevs)
127	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
128	DEALLOCATE(SCR)
129	END IF
130	END IF
131	IF ( INDEX(plot_these_fields,",tc,") /= 0 ) THEN
132	IF ( have_TK ) THEN
133	IF ( ASSOCIATED(SCR) ) DEALLOCATE(SCR)
134	ALLOCATE(SCR(west_east_dim,south_north_dim,bottom_top_dim))
135	CALL calc_tc(SCR, cname, cdesc, cunits)
136	CALL interp (SCR, west_east_dim, south_north_dim, bottom_top_dim, &
137	data_out, nxout, nyout, nzout, &
138	vert_array, interp_levels, number_of_zlevs)
139	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
140	DEALLOCATE(SCR)
141	END IF
142	END IF
143
144
145	!!! Calculate theta
146	IF ( INDEX(plot_these_fields,",theta,") /= 0 ) THEN
147	IF ( have_T ) THEN
148	IF ( ASSOCIATED(SCR) ) DEALLOCATE(SCR)
149	ALLOCATE(SCR(west_east_dim,south_north_dim,bottom_top_dim))
150	CALL calc_theta(SCR, cname, cdesc, cunits)
151	CALL interp (SCR, west_east_dim, south_north_dim, bottom_top_dim, &
152	data_out, nxout, nyout, nzout, &
153	vert_array, interp_levels, number_of_zlevs)
154	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
155	DEALLOCATE(SCR)
156	END IF
157	END IF
158
159
160	!!! Dewpoint temperature
161	IF ( INDEX(plot_these_fields,",td,") /= 0 ) THEN
162	IF ( have_QV .AND. have_PRES ) THEN
163	IF ( ASSOCIATED(SCR) ) DEALLOCATE(SCR)
164	ALLOCATE(SCR(west_east_dim,south_north_dim,bottom_top_dim))
165	CALL calc_td(SCR, cname, cdesc, cunits)
166	CALL interp (SCR, west_east_dim, south_north_dim, bottom_top_dim, &
167	data_out, nxout, nyout, nzout, &
168	vert_array, interp_levels, number_of_zlevs)
169	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
170	DEALLOCATE(SCR)
171	END IF
172	END IF
173
174
175	!!! Relative Humidity
176	IF ( INDEX(plot_these_fields,",rh,") /= 0 .OR. INDEX(plot_these_fields,",clfr,") /= 0 ) THEN
177	IF ( have_TK .AND. have_QV .AND. have_PRES ) THEN
178	IF ( ASSOCIATED(SCR) ) DEALLOCATE(SCR)
179	ALLOCATE(SCR(west_east_dim,south_north_dim,bottom_top_dim))
180	CALL calc_rh(SCR, cname, cdesc, cunits)
181	CALL interp (SCR, west_east_dim, south_north_dim, bottom_top_dim, &
182	data_out, nxout, nyout, nzout, &
183	vert_array, interp_levels, number_of_zlevs)
184	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
185
186	!!! Cloud fractions
187	IF ( INDEX(plot_these_fields,",clfr,") /= 0 ) THEN
188	IF ( ALLOCATED(SCR4) ) DEALLOCATE(SCR4)
189	ALLOCATE(SCR4(west_east_dim,south_north_dim,bottom_top_dim,4))
190	SCR4 = 0.0
191	CALL calc_clfr(SCR4, cname, cdesc, cunits, SCR)
192	DEALLOCATE(SCR)
193	ALLOCATE(SCR(west_east_dim,south_north_dim,1))
194	cname = "clflo"
195	cdesc = "Low Cloud Fraction"
196	SCR(:,:,1) = SCR4(:,:,1,1)
197	CALL interp (SCR, west_east_dim, south_north_dim, 1, &
198	data_out, nxout, nyout, nzout, &
199	vert_array, interp_levels, number_of_zlevs)
200	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
201	cname = "clfmi"
202	cdesc = "Mid Cloud Fraction"
203	SCR(:,:,1) = SCR4(:,:,1,2)
204	CALL interp (SCR, west_east_dim, south_north_dim, 1, &
205	data_out, nxout, nyout, nzout, &
206	vert_array, interp_levels, number_of_zlevs)
207	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
208	cname = "clfhi"
209	cdesc = "High Cloud Fraction"
210	SCR(:,:,1) = SCR4(:,:,1,3)
211	CALL interp (SCR, west_east_dim, south_north_dim, 1, &
212	data_out, nxout, nyout, nzout, &
213	vert_array, interp_levels, number_of_zlevs)
214	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
215	END IF
216
217	DEALLOCATE(SCR)
218	END IF
219	END IF
220
221
222	!!! Wind speed
223	IF ( INDEX(plot_these_fields,",wspd,") /= 0 ) THEN
224	IF ( have_UUU .AND. have_VVV ) THEN
225	IF ( ASSOCIATED(SCR) ) DEALLOCATE(SCR)
226	ALLOCATE(SCR(west_east_dim,south_north_dim,bottom_top_dim))
227	CALL calc_wspd(SCR, cname, cdesc, cunits, 1)
228	CALL interp (SCR, west_east_dim, south_north_dim, bottom_top_dim, &
229	data_out, nxout, nyout, nzout, &
230	vert_array, interp_levels, number_of_zlevs)
231	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
232	DEALLOCATE(SCR)
233	END IF
234	END IF
235
236
237	!!! Wind direction
238	IF ( INDEX(plot_these_fields,",wdir,") /= 0 ) THEN
239	IF ( have_UUU .AND. have_VVV ) THEN
240	IF ( ASSOCIATED(SCR) ) DEALLOCATE(SCR)
241	ALLOCATE(SCR(west_east_dim,south_north_dim,bottom_top_dim))
242	CALL calc_wdir(SCR, cname, cdesc, cunits, 1)
243	CALL interp (SCR, west_east_dim, south_north_dim, bottom_top_dim, &
244	data_out, nxout, nyout, nzout, &
245	vert_array, interp_levels, number_of_zlevs)
246	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
247	DEALLOCATE(SCR)
248	END IF
249	END IF
250
251
252	!!! Wind speed - 10m
253	IF ( INDEX(plot_these_fields,",ws10,") /= 0 ) THEN
254	IF ( have_U10 .AND. have_V10 ) THEN
255	IF ( ASSOCIATED(SCR) ) DEALLOCATE(SCR)
256	ALLOCATE(SCR(west_east_dim,south_north_dim,1))
257	CALL calc_wspd(SCR, cname, cdesc, cunits, 0)
258	CALL interp (SCR, west_east_dim, south_north_dim, 1, &
259	data_out, nxout, nyout, nzout, &
260	vert_array, interp_levels, number_of_zlevs)
261	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
262	DEALLOCATE(SCR)
263	END IF
264	END IF
265
266
267	!!! Wind direction - 10m
268	IF ( INDEX(plot_these_fields,",wd10,") /= 0 ) THEN
269	IF ( have_U10 .AND. have_V10 ) THEN
270	IF ( ASSOCIATED(SCR) ) DEALLOCATE(SCR)
271	ALLOCATE(SCR(west_east_dim,south_north_dim,1))
272	CALL calc_wdir(SCR, cname, cdesc, cunits, 0)
273	CALL interp (SCR, west_east_dim, south_north_dim, 1, &
274	data_out, nxout, nyout, nzout, &
275	vert_array, interp_levels, number_of_zlevs)
276	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
277	DEALLOCATE(SCR)
278	END IF
279	END IF
280
281
282	!!! Rotated U and V
283	IF ( INDEX(plot_these_fields,",umet,") /= 0 .OR. &
284	INDEX(plot_these_fields,",vmet,") /= 0) THEN
285	IF ( have_UUU .AND. have_VVV .AND. have_XLAT .AND. have_XLONG ) THEN
286
287	IF ( ALLOCATED(SCR4) ) DEALLOCATE(SCR4)
288	ALLOCATE(SCR4(west_east_dim,south_north_dim,bottom_top_dim,4))
289	SCR4 = 0.0
290	CALL calc_uvmet(SCR4, cname, cdesc, cunits)
291
292	ALLOCATE(SCR(west_east_dim,south_north_dim,bottom_top_dim))
293	IF ( INDEX(plot_these_fields,",umet,") /= 0 ) THEN
294	cname = "umet"
295	SCR(:,:,:) = SCR4(:,:,:,1)
296	CALL interp (SCR, west_east_dim, south_north_dim, bottom_top_dim, &
297	data_out, nxout, nyout, nzout, &
298	vert_array, interp_levels, number_of_zlevs)
299	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
300	END IF
301	IF ( INDEX(plot_these_fields,",vmet,") /= 0) THEN
302	cname = "vmet"
303	SCR(:,:,:) = SCR4(:,:,:,2)
304	CALL interp (SCR, west_east_dim, south_north_dim, bottom_top_dim, &
305	data_out, nxout, nyout, nzout, &
306	vert_array, interp_levels, number_of_zlevs)
307	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
308	END IF
309	DEALLOCATE(SCR)
310	DEALLOCATE(SCR4)
311
312	END IF
313	END IF
314
315
316	!!! Rotated U10 and V10
317	IF ( INDEX(plot_these_fields,",u10m,") /= 0 .OR. &
318	INDEX(plot_these_fields,",v10m,") /= 0) THEN
319	IF ( have_U10 .AND. have_V10 .AND. have_XLAT .AND. have_XLONG ) THEN
320
321	IF ( ALLOCATED(SCR4) ) DEALLOCATE(SCR4)
322	ALLOCATE(SCR4(west_east_dim,south_north_dim,bottom_top_dim,4))
323	SCR4 = 0.0
324	CALL calc_uvmet(SCR4, cname, cdesc, cunits)
325
326	ALLOCATE(SCR(west_east_dim,south_north_dim,1))
327	IF ( INDEX(plot_these_fields,",u10m,") /= 0 ) THEN
328	cname = "u10m"
329	SCR(:,:,1) = SCR4(:,:,1,1)
330	CALL interp (SCR, west_east_dim, south_north_dim, 1, &
331	data_out, nxout, nyout, nzout, &
332	vert_array, interp_levels, number_of_zlevs)
333	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
334	END IF
335	IF ( INDEX(plot_these_fields,",v10m,") /= 0) THEN
336	cname = "v10m"
337	SCR(:,:,1) = SCR4(:,:,1,2)
338	CALL interp( SCR, west_east_dim, south_north_dim, 1, &
339	data_out, nxout, nyout, nzout, &
340	vert_array, interp_levels, number_of_zlevs)
341	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
342	END IF
343	DEALLOCATE(SCR)
344	DEALLOCATE(SCR4)
345
346	END IF
347	END IF
348
349
350	!!! Sea Level Pressure
351	IF ( INDEX(plot_these_fields,",slp,") /= 0 ) THEN
352	IF ( have_TK .AND. have_QV .AND. have_PRES .AND. have_PH .AND. have_PHB) THEN
353	IF ( ASSOCIATED(SCR) ) DEALLOCATE(SCR)
354	ALLOCATE(SCR(west_east_dim,south_north_dim,bottom_top_dim))
355	CALL calc_slp(SCR, cname, cdesc, cunits)
356	CALL interp (SCR, west_east_dim, south_north_dim, 1, &
357	data_out, nxout, nyout, nzout, &
358	vert_array, interp_levels, number_of_zlevs)
359	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
360	DEALLOCATE(SCR)
361	END IF
362	END IF
363
364
365	!!! Reflectivity
366	IF ( INDEX(plot_these_fields,",dbz,") /= 0 .OR. &
367	INDEX(plot_these_fields,",max_dbz,") /= 0) THEN
368	IF ( have_QV .AND. have_QR .AND. have_TK .AND. have_PRES ) THEN
369	CALL calc_dbz(SCR, cname, cdesc, cunits)
370
371	IF ( INDEX(plot_these_fields,",dbz,") /= 0 ) THEN
372	CALL interp (SCR, west_east_dim, south_north_dim, bottom_top_dim, &
373	data_out, nxout, nyout, nzout, &
374	vert_array, interp_levels, number_of_zlevs)
375	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
376	END IF
377
378	IF ( INDEX(plot_these_fields,",max_dbz,") /= 0) THEN
379	ALLOCATE(SCR2(west_east_dim,south_north_dim))
380	SCR2 = 0.0
381	DO jj = 1,south_north_dim
382	DO ii = 1,west_east_dim
383	DO kk = 1,bottom_top_dim
384	SCR2(ii,jj) = MAX( SCR2(ii,jj) , SCR(ii,jj,kk) )
385	END DO
386	END DO
387	END DO
388	IF (ASSOCIATED(SCR)) DEALLOCATE(SCR)
389	ALLOCATE(SCR(west_east_dim,south_north_dim,1))
390	SCR(:,:,1) = SCR2(:,:)
391	cname = "max_dbz"
392	cdesc = "Max Reflectivity"
393	CALL interp (SCR, west_east_dim, south_north_dim, 1, &
394	data_out, nxout, nyout, nzout, &
395	vert_array, interp_levels, number_of_zlevs)
396	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
397	DEALLOCATE(SCR2)
398	DEALLOCATE(SCR)
399	END IF
400
401	END IF
402	END IF
403
404
405
406	!!! 3D CAPE/CIN
407	IF ( INDEX(plot_these_fields,",cape,") /= 0 .OR. &
408	INDEX(plot_these_fields,",cin,") /= 0) THEN
409	IF ( have_HGT .AND. have_QV .AND. have_PRES .AND. have_TK .AND. &
410	have_PH .AND. have_PHB) THEN
411
412	IF ( ALLOCATED(SCR4) ) DEALLOCATE(SCR4)
413	ALLOCATE(SCR4(west_east_dim,south_north_dim,bottom_top_dim,4))
414	SCR4 = 0.0
415	CALL calc_cape(SCR4, cname, cdesc, cunits,1)
416
417	ALLOCATE(SCR(west_east_dim,south_north_dim,bottom_top_dim))
418	IF ( INDEX(plot_these_fields,",cape,") /= 0 ) THEN
419	cname = "cape"
420	cdesc = "CAPE"
421	cunits = "J/kg"
422	SCR(:,:,:) = SCR4(:,:,:,1)
423	CALL interp (SCR, west_east_dim, south_north_dim, bottom_top_dim, &
424	data_out, nxout, nyout, nzout, &
425	vert_array, interp_levels, number_of_zlevs)
426	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
427	END IF
428	IF ( INDEX(plot_these_fields,",cin,") /= 0) THEN
429	cname = "cin"
430	cdesc = "CIN"
431	cunits = "J/kg"
432	SCR(:,:,:) = SCR4(:,:,:,2)
433	CALL interp (SCR, west_east_dim, south_north_dim, bottom_top_dim, &
434	data_out, nxout, nyout, nzout, &
435	vert_array, interp_levels, number_of_zlevs)
436	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
437	END IF
438	DEALLOCATE(SCR)
439	DEALLOCATE(SCR4)
440
441	END IF
442	END IF
443
444
445
446	!!! 2D CAPE/CIN & LCL/LFC
447	IF ( INDEX(plot_these_fields,",mcape,") /= 0 .OR. &
448	INDEX(plot_these_fields,",mcin,") /= 0 .OR. &
449	INDEX(plot_these_fields,",lcl,") /= 0 .OR. &
450	INDEX(plot_these_fields,",lfc,") /= 0) THEN
451	IF ( have_HGT .AND. have_QV .AND. have_PRES .AND. have_TK .AND. &
452	have_PH .AND. have_PHB) THEN
453
454	IF ( ALLOCATED(SCR4) ) DEALLOCATE(SCR4)
455	ALLOCATE(SCR4(west_east_dim,south_north_dim,bottom_top_dim,4))
456	SCR4 = 0.0
457	CALL calc_cape(SCR4, cname, cdesc, cunits,0)
458
459	ALLOCATE(SCR(west_east_dim,south_north_dim,1))
460	IF ( INDEX(plot_these_fields,",mcape,") /= 0 ) THEN
461	cname = "mcape"
462	cdesc = "MCAPE"
463	cunits = "J/kg"
464	SCR(:,:,1) = SCR4(:,:,1,1)
465	CALL interp (SCR, west_east_dim, south_north_dim, 1, &
466	data_out, nxout, nyout, nzout, &
467	vert_array, interp_levels, number_of_zlevs)
468	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
469	END IF
470	IF ( INDEX(plot_these_fields,",mcin,") /= 0) THEN
471	cname = "mcin"
472	cdesc = "MCIN"
473	cunits = "J/kg"
474	SCR(:,:,1) = SCR4(:,:,1,2)
475	CALL interp (SCR, west_east_dim, south_north_dim, 1, &
476	data_out, nxout, nyout, nzout, &
477	vert_array, interp_levels, number_of_zlevs)
478	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
479	END IF
480	IF ( INDEX(plot_these_fields,",lcl,") /= 0) THEN
481	cname = "lcl"
482	cdesc = "LCL"
483	cunits = "meters AGL"
484	SCR(:,:,1) = SCR4(:,:,1,3)
485	CALL interp (SCR, west_east_dim, south_north_dim, 1, &
486	data_out, nxout, nyout, nzout, &
487	vert_array, interp_levels, number_of_zlevs)
488	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
489	END IF
490	IF ( INDEX(plot_these_fields,",lfc,") /= 0) THEN
491	cname = "lfc"
492	cdesc = "LFC"
493	cunits = "meters AGL"
494	SCR(:,:,1) = SCR4(:,:,1,4)
495	CALL interp (SCR, west_east_dim, south_north_dim, 1, &
496	data_out, nxout, nyout, nzout, &
497	vert_array, interp_levels, number_of_zlevs)
498	CALL write_dat (data_out, nxout, nyout, nzout, cname, cdesc, cunits)
499	END IF
500	DEALLOCATE(SCR)
501	DEALLOCATE(SCR4)
502
503	END IF
504	END IF
505
506
507	IF ( ALLOCATED(SCR4) ) DEALLOCATE(SCR4)
508	IF ( ASSOCIATED(SCR) ) DEALLOCATE(SCR)
509
510	END SUBROUTINE process_diagnostics
511
512	END MODULE module_diagnostics

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

Download in other formats:

Original Format