Home

Context Navigation

data_output_profiles.f90 @ 964

Last change on this file since 964 was 964, checked in by raasch, 12 years ago

old profil-parameters (cross_xtext, cross_normalized_x, etc. ) and respective code removed
(check_open, check_parameters, close_file, data_output_profiles, data_output_spectra, header, modules, parin)

reformatting (netcdf)

append feature removed from unit 14 (check_open)

Property svn:keywords set to Id

File size: 10.7 KB

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

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

Context Navigation

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

Download in other formats: