Home

Context Navigation

source: palm/trunk/SOURCE/data_output_profiles.f90 @ 1026

Last change on this file since 1026 was 965, checked in by raasch, 12 years ago
last commit documented
Property svn:keywords set to `Id`
File size: 10.8 KB

Line
1	SUBROUTINE data_output_profiles
2
3	!------------------------------------------------------------------------------!
4	! Current revisions:
5	! -----------------
6	!
7	!
8	! Former revisions:
9	! -----------------
10	! $Id: data_output_profiles.f90 965 2012-07-26 09:29:52Z letzel $
11	!
12	! 964 2012-07-26 09:14:24Z raasch
13	! code for profil-output removed
14	!
15	! 345 2009-07-01 14:37:56Z heinze
16	! In case of restart runs without extension, initial profiles are not written
17	! to NetCDF-file anymore.
18	! simulated_time in NetCDF output replaced by time_since_reference_point.
19	! Output of NetCDF messages with aid of message handling routine.
20	! Output of messages replaced by message handling routine.
21	!
22	! 197 2008-09-16 15:29:03Z raasch
23	! Time coordinate t=0 stored on netcdf-file only if an output is required for
24	! this time for at least one of the profiles
25	!
26	! February 2007
27	! RCS Log replace by Id keyword, revision history cleaned up
28	!
29	! 87 2007-05-22 15:46:47Z raasch
30	! var_hom renamed pr_palm
31	!
32	! Revision 1.18 2006/08/16 14:27:04 raasch
33	! PRINT* statements for testing removed
34	!
35	! Revision 1.1 1997/09/12 06:28:48 raasch
36	! Initial revision
37	!
38	!
39	! Description:
40	! ------------
41	! Plot output of 1D-profiles for PROFIL
42	!------------------------------------------------------------------------------!
43
44	USE control_parameters
45	USE cpulog
46	USE indices
47	USE interfaces
48	USE netcdf_control
49	USE pegrid
50	USE profil_parameter
51	USE statistics
52
53	IMPLICIT NONE
54
55
56	INTEGER :: i, id, ilc, ils, j, k, sr
57
58	!
59	!-- If required, compute statistics
60	IF ( .NOT. flow_statistics_called ) CALL flow_statistics
61
62	!
63	!-- Flow_statistics has its own CPU time measurement
64	CALL cpu_log( log_point(15), 'data_output_profiles', 'start' )
65
66	!
67	!-- If required, compute temporal average
68	IF ( averaging_interval_pr == 0.0 ) THEN
69	hom_sum(:,:,:) = hom(:,1,:,:)
70	ELSE
71	IF ( average_count_pr > 0 ) THEN
72	hom_sum = hom_sum / REAL( average_count_pr )
73	ELSE
74	!
75	!-- This case may happen if dt_dopr is changed in the d3par-list of
76	!-- a restart run
77	RETURN
78	ENDIF
79	ENDIF
80
81
82	IF ( myid == 0 ) THEN
83
84	!
85	!-- Plot-output for each (sub-)region
86
87	!
88	!-- Open file for profile output in NetCDF format
89	IF ( netcdf_output ) THEN
90	CALL check_open( 104 )
91	ENDIF
92
93	!
94	!-- Increment the counter for number of output times
95	dopr_time_count = dopr_time_count + 1
96
97	!
98	!-- Output of initial profiles
99	IF ( dopr_time_count == 1 ) THEN
100
101	IF ( .NOT. output_for_t0 ) THEN
102
103	IF ( netcdf_output ) THEN
104	#if defined( __netcdf )
105	!
106	!-- Store initial time (t=0) to time axis, but only if an output
107	!-- is required for at least one of the profiles
108	DO i = 1, dopr_n
109	IF ( dopr_initial_index(i) /= 0 ) THEN
110	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_time_pr, &
111	(/ 0.0 /), start = (/ 1 /), &
112	count = (/ 1 /) )
113	CALL handle_netcdf_error( 'data_output_profiles', 329 )
114	output_for_t0 = .TRUE.
115	EXIT
116	ENDIF
117	ENDDO
118
119	!
120	!-- Store normalization factors
121	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(1), & ! wpt0
122	(/ hom_sum(nzb,18,normalizing_region) /), &
123	start = (/ 1 /), count = (/ 1 /) )
124	CALL handle_netcdf_error( 'data_output_profiles', 330 )
125
126	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(2), & ! ws2
127	(/ hom_sum(nzb+8,pr_palm,normalizing_region)**2 /), &
128	start = (/ 1 /), count = (/ 1 /) )
129	CALL handle_netcdf_error( 'data_output_profiles', 331 )
130	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(3), & ! tsw2
131	(/ hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), &
132	start = (/ 1 /), count = (/ 1 /) )
133	CALL handle_netcdf_error( 'data_output_profiles', 332 )
134	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(4), & ! ws3
135	(/ hom_sum(nzb+8,pr_palm,normalizing_region)**3 /), &
136	start = (/ 1 /), count = (/ 1 /) )
137	CALL handle_netcdf_error( 'data_output_profiles', 333 )
138
139	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(5), &!ws2tsw
140	(/ hom_sum(nzb+8,pr_palm,normalizing_region)*3 &
141	hom_sum(nzb+3,pr_palm,normalizing_region) /), &
142	start = (/ 1 /), count = (/ 1 /) )
143	CALL handle_netcdf_error( 'data_output_profiles', 334 )
144
145	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(6), &!wstsw2
146	(/ hom_sum(nzb+8,pr_palm,normalizing_region) * &
147	hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), &
148	start = (/ 1 /), count = (/ 1 /) )
149	CALL handle_netcdf_error( 'data_output_profiles', 335 )
150
151	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(7), & ! z_i
152	(/ hom_sum(nzb+6,pr_palm,normalizing_region) /), &
153	start = (/ 1 /), count = (/ 1 /) )
154	CALL handle_netcdf_error( 'data_output_profiles', 336 )
155
156	#endif
157	ENDIF
158	!
159	!-- Loop over all 1D variables
160	DO i = 1, dopr_n
161
162	IF ( dopr_initial_index(i) /= 0 ) THEN
163
164	!
165	!-- Output for the individual (sub-)regions
166	DO sr = 0, statistic_regions
167
168	IF ( netcdf_output ) THEN
169	#if defined( __netcdf )
170	!
171	!-- Write data to netcdf file
172	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_dopr(i,sr), &
173	hom(nzb:nzt+1,1,dopr_initial_index(i),sr), &
174	start = (/ 1, 1 /), &
175	count = (/ nzt-nzb+2, 1 /) )
176	CALL handle_netcdf_error( 'data_output_profiles', 337 )
177	#endif
178	ENDIF
179
180	ENDDO
181
182	ENDIF ! Initial profile available
183
184	ENDDO ! Loop over dopr_n for initial profiles
185
186	IF ( netcdf_output .AND. output_for_t0 ) THEN
187	dopr_time_count = dopr_time_count + 1
188	ENDIF
189
190	END IF
191	ENDIF ! Initial profiles
192
193	IF ( netcdf_output ) THEN
194	#if defined( __netcdf )
195
196	!
197	!-- Store time to time axis
198	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_time_pr, &
199	(/ time_since_reference_point /), &
200	start = (/ dopr_time_count /), &
201	count = (/ 1 /) )
202	CALL handle_netcdf_error( 'data_output_profiles', 338 )
203
204	!
205	!-- Store normalization factors
206	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(1), & ! wpt0
207	(/ hom_sum(nzb,18,normalizing_region) /), &
208	start = (/ dopr_time_count /), &
209	count = (/ 1 /) )
210	CALL handle_netcdf_error( 'data_output_profiles', 339 )
211
212	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(2), & ! ws2
213	(/ hom_sum(nzb+8,pr_palm,normalizing_region)**2 /), &
214	start = (/ dopr_time_count /), &
215	count = (/ 1 /) )
216	CALL handle_netcdf_error( 'data_output_profiles', 340 )
217
218	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(3), & ! tsw2
219	(/ hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), &
220	start = (/ dopr_time_count /), &
221	count = (/ 1 /) )
222	CALL handle_netcdf_error( 'data_output_profiles', 341 )
223
224	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(4), & ! ws3
225	(/ hom_sum(nzb+8,pr_palm,normalizing_region)**3 /), &
226	start = (/ dopr_time_count /), &
227	count = (/ 1 /) )
228	CALL handle_netcdf_error( 'data_output_profiles', 342 )
229
230	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(5), & ! ws2tsw
231	(/ hom_sum(nzb+8,pr_palm,normalizing_region)*3 &
232	hom_sum(nzb+3,pr_palm,normalizing_region) /), &
233	start = (/ dopr_time_count /), &
234	count = (/ 1 /) )
235	CALL handle_netcdf_error( 'data_output_profiles', 343 )
236
237	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(6), & ! wstsw2
238	(/ hom_sum(nzb+8,pr_palm,normalizing_region) * &
239	hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), &
240	start = (/ dopr_time_count /), &
241	count = (/ 1 /) )
242	CALL handle_netcdf_error( 'data_output_profiles', 344 )
243
244	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(7), & ! z_i
245	(/ hom_sum(nzb+6,pr_palm,normalizing_region) /), &
246	start = (/ dopr_time_count /), &
247	count = (/ 1 /) )
248	CALL handle_netcdf_error( 'data_output_profiles', 345 )
249	#endif
250	ENDIF
251
252	!
253	!-- Output of the individual (non-initial) profiles
254	DO i = 1, dopr_n
255
256	!
257	!-- Output for the individual (sub-)domains
258	DO sr = 0, statistic_regions
259
260	IF ( netcdf_output ) THEN
261	#if defined( __netcdf )
262	!
263	!-- Write data to netcdf file
264	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_dopr(i,sr), &
265	hom_sum(nzb:nzt+1,dopr_index(i),sr),&
266	start = (/ 1, dopr_time_count /), &
267	count = (/ nzt-nzb+2, 1 /) )
268	CALL handle_netcdf_error( 'data_output_profiles', 346 )
269	#endif
270	ENDIF
271
272	ENDDO
273
274	ENDDO
275
276	ENDIF ! Output on PE0
277
278	!
279	!-- If averaging has been done above, the summation counter must be re-set.
280	IF ( averaging_interval_pr /= 0.0 ) THEN
281	average_count_pr = 0
282	ENDIF
283
284	CALL cpu_log( log_point(15), 'data_output_profiles','stop', 'nobarrier' )
285
286	!
287	!-- Formats
288	100 FORMAT ('#1 ',A,1X,A)
289	101 FORMAT (E15.7,1X,E15.7)
290	102 FORMAT ('NEXT')
291
292	END SUBROUTINE data_output_profiles

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |