Home

Context Navigation

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

Last change on this file since 197 was 197, checked in by raasch, 16 years ago
further adjustments for SGI and other small changes
Property svn:keywords set to `Id`
File size: 24.5 KB

Line
1	SUBROUTINE data_output_profiles
2
3	!------------------------------------------------------------------------------!
4	! Actual revisions:
5	! -----------------
6	! Time coordinate t=0 stored on netcdf-file only if an output is required for
7	! this time for at least one of the profiles
8	!
9	! Former revisions:
10	! -----------------
11	! $Id: data_output_profiles.f90 197 2008-09-16 15:29:03Z raasch $
12	! RCS Log replace by Id keyword, revision history cleaned up
13	!
14	! 87 2007-05-22 15:46:47Z raasch
15	! var_hom renamed pr_palm
16	!
17	! Revision 1.18 2006/08/16 14:27:04 raasch
18	! PRINT* statements for testing removed
19	!
20	! Revision 1.1 1997/09/12 06:28:48 raasch
21	! Initial revision
22	!
23	!
24	! Description:
25	! ------------
26	! Plot output of 1D-profiles for PROFIL
27	!------------------------------------------------------------------------------!
28
29	USE control_parameters
30	USE cpulog
31	USE indices
32	USE interfaces
33	USE netcdf_control
34	USE pegrid
35	USE profil_parameter
36	USE statistics
37
38	IMPLICIT NONE
39
40
41	INTEGER :: i, id, ilc, ils, j, k, sr
42	LOGICAL :: output_for_t0
43	REAL :: uxma, uxmi
44
45
46	!
47	!-- If required, compute statistics
48	IF ( .NOT. flow_statistics_called ) CALL flow_statistics
49
50	!
51	!-- Flow_statistics has its own CPU time measurement
52	CALL cpu_log( log_point(15), 'data_output_profiles', 'start' )
53
54	!
55	!-- If required, compute temporal average
56	IF ( averaging_interval_pr == 0.0 ) THEN
57	hom_sum(:,:,:) = hom(:,1,:,:)
58	ELSE
59	IF ( average_count_pr > 0 ) THEN
60	hom_sum = hom_sum / REAL( average_count_pr )
61	ELSE
62	!
63	!-- This case may happen if dt_dopr is changed in the d3par-list of
64	!-- a restart run
65	RETURN
66	ENDIF
67	ENDIF
68
69
70	IF ( myid == 0 ) THEN
71
72	!
73	!-- Plot-output for each (sub-)region
74
75	!
76	!-- Open file for profile output in NetCDF format
77	IF ( netcdf_output ) THEN
78	CALL check_open( 104 )
79	ENDIF
80
81	!
82	!-- Open PROFIL-output files for each (sub-)region
83	IF ( profil_output ) THEN
84	DO sr = 0, statistic_regions
85	CALL check_open( 40 + sr )
86	ENDDO
87	ENDIF
88
89	!
90	!-- Increment the counter for number of output times
91	dopr_time_count = dopr_time_count + 1
92
93	!
94	!-- Re-set to zero the counter for the number of profiles already written
95	!-- at the current output time into the respective crosses
96	cross_pnc_local = 0
97
98	!
99	!-- Output of initial profiles
100	IF ( dopr_time_count == 1 ) THEN
101
102	IF ( netcdf_output ) THEN
103	#if defined( __netcdf )
104	!
105	!-- Store initial time (t=0) to time axis, but only if an output
106	!-- is required for at least one of the profiles
107	output_for_t0 = .FALSE.
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	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 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	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 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	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 331 )
130
131	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(3), & ! tsw2
132	(/ hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), &
133	start = (/ 1 /), count = (/ 1 /) )
134	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 332 )
135
136	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(4), & ! ws3
137	(/ hom_sum(nzb+8,pr_palm,normalizing_region)**3 /), &
138	start = (/ 1 /), count = (/ 1 /) )
139	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 333 )
140
141	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(5), & ! ws2tsw
142	(/ hom_sum(nzb+8,pr_palm,normalizing_region)*3 &
143	hom_sum(nzb+3,pr_palm,normalizing_region) /), &
144	start = (/ 1 /), count = (/ 1 /) )
145	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 334 )
146
147	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(6), & ! wstsw2
148	(/ hom_sum(nzb+8,pr_palm,normalizing_region) * &
149	hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), &
150	start = (/ 1 /), count = (/ 1 /) )
151	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 335 )
152
153	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(7), & ! z_i
154	(/ hom_sum(nzb+6,pr_palm,normalizing_region) /), &
155	start = (/ 1 /), count = (/ 1 /) )
156	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 336 )
157	#endif
158	ENDIF
159	!
160	!-- Loop over all 1D variables
161	DO i = 1, dopr_n
162
163	IF ( dopr_initial_index(i) /= 0 ) THEN
164
165	!
166	!-- Output for the individual (sub-)regions
167	DO sr = 0, statistic_regions
168
169	IF ( profil_output ) THEN
170	id = 40 + sr
171	!
172	!-- Write Label-Header
173	WRITE ( id, 100 ) TRIM( data_output_pr(i) ), '(t=0)'
174	!
175	!-- Write total profile
176	DO k = nzb, nzt+1
177	WRITE ( id, 101 ) hom(k,2,dopr_initial_index(i),sr), &
178	hom(k,1,dopr_initial_index(i),sr)
179	ENDDO
180	!
181	!-- Write separation label
182	WRITE ( id, 102 )
183	ENDIF
184
185	IF ( netcdf_output ) THEN
186	#if defined( __netcdf )
187	!
188	!-- Write data to netcdf file
189	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_dopr(i,sr), &
190	hom(nzb:nzt+1,1,dopr_initial_index(i),sr), &
191	start = (/ 1, 1 /), &
192	count = (/ nzt-nzb+2, 1 /) )
193	IF ( nc_stat /= NF90_NOERR ) &
194	CALL handle_netcdf_error( 337 )
195	#endif
196	ENDIF
197
198	ENDDO
199
200	IF ( profil_output ) THEN
201	!
202	!-- Determine indices for later NAMELIST-output (s. below)
203	profile_number = profile_number + 1
204	j = dopr_crossindex(i)
205	IF ( j /= 0 ) THEN
206	cross_profile_number_count(j) = &
207	cross_profile_number_count(j) + 1
208	k = cross_profile_number_count(j)
209	cross_profile_numbers(k,j) = profile_number
210	!
211	!-- Initial profiles are always drawn as solid lines in
212	!-- anti-background colour.
213	cross_linecolors(k,j) = 1
214	cross_linestyles(k,j) = 0
215	!
216	!-- If required, extend x-value range of the respective
217	!-- cross, provided it has not been specified in &
218	!-- check_parameters. Determination over all (sub-)regions.
219	IF ( cross_uxmin(j) == 0.0 .AND. &
220	cross_uxmax(j) == 0.0 ) THEN
221
222	DO sr = 0, statistic_regions
223
224	uxmi = &
225	MINVAL( hom(:nz_do1d,1,dopr_initial_index(i),sr) )
226
227	uxma = &
228	MAXVAL( hom(:nz_do1d,1,dopr_initial_index(i),sr) )
229	!
230	!-- When the value range of the first line in the
231	!-- corresponding cross is determined, its value range
232	!-- is simply adopted.
233	IF ( cross_uxmin_computed(j) > &
234	cross_uxmax_computed(j) ) THEN
235	cross_uxmin_computed(j) = uxmi
236	cross_uxmax_computed(j) = uxma
237	ELSE
238	cross_uxmin_computed(j) = &
239	MIN( cross_uxmin_computed(j), uxmi )
240	cross_uxmax_computed(j) = &
241	MAX( cross_uxmax_computed(j), uxma )
242	ENDIF
243
244	ENDDO
245
246	ENDIF
247	!
248	!-- If required, determine and note normalizing factors
249	SELECT CASE ( cross_normalized_x(j) )
250
251	CASE ( 'ts2' )
252	cross_normx_factor(k,j) = &
253	( hom_sum(nzb+3,pr_palm,normalizing_region) )**2
254	CASE ( 'wpt0' )
255	cross_normx_factor(k,j) = &
256	hom_sum(nzb,18,normalizing_region)
257	CASE ( 'wsts2' )
258	cross_normx_factor(k,j) = &
259	hom_sum(nzb+8,pr_palm,normalizing_region) &
260	* ( hom_sum(nzb+3,pr_palm,normalizing_region) )**2
261	CASE ( 'ws2' )
262	cross_normx_factor(k,j) = &
263	( hom_sum(nzb+8,pr_palm,normalizing_region) )**2
264	CASE ( 'ws2ts' )
265	cross_normx_factor(k,j) = &
266	( hom_sum(nzb+8,pr_palm,normalizing_region) )**2 &
267	* hom_sum(nzb+3,pr_palm,normalizing_region)
268	CASE ( 'ws3' )
269	cross_normx_factor(k,j) = &
270	( hom_sum(nzb+8,pr_palm,normalizing_region) )**3
271
272	END SELECT
273
274	SELECT CASE ( cross_normalized_y(j) )
275
276	CASE ( 'z_i' )
277	cross_normy_factor(k,j) = &
278	hom_sum(nzb+6,pr_palm,normalizing_region)
279
280	END SELECT
281
282	!
283	!-- Check the normalizing factors for zeros and deactivate
284	!-- the normalization, if required.
285	IF ( cross_normx_factor(k,j) == 0.0 .OR. &
286	cross_normy_factor(k,j) == 0.0 ) THEN
287	PRINT*,'+++ WARNING data_output_profiles: normalizi', &
288	'ng cross ',j, ' is not possible since one o', &
289	'f the'
290	PRINT*,' normalizing factors is ',&
291	'zero!'
292	PRINT*,' cross_normx_factor(',k,',',j,') = ', &
293	cross_normx_factor(k,j)
294	PRINT*,' cross_normy_factor(',k,',',j,') = ', &
295	cross_normy_factor(k,j)
296	cross_normx_factor(k,j) = 1.0
297	cross_normy_factor(k,j) = 1.0
298	cross_normalized_x(j) = ' '
299	cross_normalized_y(j) = ' '
300	ENDIF
301
302	!
303	!-- If required, extend normalized x-value range of the
304	!-- respective cross, provided it has not been specified in
305	!-- check_parameters. Determination over all (sub-)regions.
306	IF ( cross_uxmin_normalized(j) == 0.0 .AND. &
307	cross_uxmax_normalized(j) == 0.0 ) THEN
308
309	DO sr = 0, statistic_regions
310
311	uxmi = MINVAL( hom(:nz_do1d,1, &
312	dopr_initial_index(i),sr) ) / &
313	cross_normx_factor(k,j)
314	uxma = MAXVAL( hom(:nz_do1d,1, &
315	dopr_initial_index(i),sr) ) / &
316	cross_normx_factor(k,j)
317	!
318	!-- When the value range of the first line in the
319	!-- corresponding cross is determined, its value range
320	!-- is simply adopted.
321	IF ( cross_uxmin_normalized_computed(j) > &
322	cross_uxmax_normalized_computed(j) ) THEN
323	cross_uxmin_normalized_computed(j) = uxmi
324	cross_uxmax_normalized_computed(j) = uxma
325	ELSE
326	cross_uxmin_normalized_computed(j) = &
327	MIN( cross_uxmin_normalized_computed(j), uxmi )
328	cross_uxmax_normalized_computed(j) = &
329	MAX( cross_uxmax_normalized_computed(j), uxma )
330	ENDIF
331
332	ENDDO
333
334	ENDIF
335
336	ENDIF ! Index determination
337
338	ENDIF ! profil output
339
340	ENDIF ! Initial profile available
341
342	ENDDO ! Loop over dopr_n for initial profiles
343
344	IF ( netcdf_output .AND. output_for_t0 ) THEN
345	dopr_time_count = dopr_time_count + 1
346	ENDIF
347
348	ENDIF ! Initial profiles
349
350	IF ( netcdf_output ) THEN
351	#if defined( __netcdf )
352	!
353	!-- Store time to time axis
354	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_time_pr, &
355	(/ simulated_time /), &
356	start = (/ dopr_time_count /), &
357	count = (/ 1 /) )
358	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 338 )
359
360	!
361	!-- Store normalization factors
362	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(1), & ! wpt0
363	(/ hom_sum(nzb,18,normalizing_region) /), &
364	start = (/ dopr_time_count /), &
365	count = (/ 1 /) )
366	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 339 )
367
368	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(2), & ! ws2
369	(/ hom_sum(nzb+8,pr_palm,normalizing_region)**2 /), &
370	start = (/ dopr_time_count /), &
371	count = (/ 1 /) )
372	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 340 )
373
374	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(3), & ! tsw2
375	(/ hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), &
376	start = (/ dopr_time_count /), &
377	count = (/ 1 /) )
378	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 341 )
379
380	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(4), & ! ws3
381	(/ hom_sum(nzb+8,pr_palm,normalizing_region)**3 /), &
382	start = (/ dopr_time_count /), &
383	count = (/ 1 /) )
384	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 342 )
385
386	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(5), & ! ws2tsw
387	(/ hom_sum(nzb+8,pr_palm,normalizing_region)*3 &
388	hom_sum(nzb+3,pr_palm,normalizing_region) /), &
389	start = (/ dopr_time_count /), &
390	count = (/ 1 /) )
391	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 343 )
392
393	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(6), & ! wstsw2
394	(/ hom_sum(nzb+8,pr_palm,normalizing_region) * &
395	hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), &
396	start = (/ dopr_time_count /), &
397	count = (/ 1 /) )
398	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 344 )
399
400	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(7), & ! z_i
401	(/ hom_sum(nzb+6,pr_palm,normalizing_region) /), &
402	start = (/ dopr_time_count /), &
403	count = (/ 1 /) )
404	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 345 )
405	#endif
406	ENDIF
407
408	!
409	!-- Output of the individual (non-initial) profiles
410	DO i = 1, dopr_n
411
412	!
413	!-- Output for the individual (sub-)domains
414	DO sr = 0, statistic_regions
415
416	IF ( profil_output ) THEN
417	id = 40 + sr
418	!
419	!-- Write Label-Header
420	WRITE ( id, 100 ) TRIM( dopr_label(i) ), simulated_time_chr
421	!
422	!-- Output of total profile
423	DO k = nzb, nzt+1
424	WRITE ( id, 101 ) hom(k,2,dopr_index(i),sr), &
425	hom_sum(k,dopr_index(i),sr)
426	ENDDO
427	!
428	!-- Write separation label
429	WRITE ( id, 102 )
430	ENDIF
431
432	IF ( netcdf_output ) THEN
433	#if defined( __netcdf )
434	!
435	!-- Write data to netcdf file
436	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_dopr(i,sr), &
437	hom_sum(nzb:nzt+1,dopr_index(i),sr),&
438	start = (/ 1, dopr_time_count /), &
439	count = (/ nzt-nzb+2, 1 /) )
440	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 346 )
441	#endif
442	ENDIF
443
444	ENDDO
445
446	IF ( profil_output ) THEN
447	!
448	!-- Determine profile number on file and note the data for later
449	!-- NAMELIST output, if the respective profile is to be drawn by
450	!-- PROFIL (if it shall not be drawn, the variable dopr_crossindex has
451	!-- the value 0, otherwise the number of the coordinate cross)
452	profile_number = profile_number + 1
453	j = dopr_crossindex(i)
454
455	IF ( j /= 0 ) THEN
456	cross_profile_number_count(j) = cross_profile_number_count(j) +1
457	k = cross_profile_number_count(j)
458	cross_pnc_local(j) = cross_pnc_local(j) +1
459	cross_profile_numbers(k,j) = profile_number
460	ilc = MOD( dopr_time_count, 10 )
461	IF ( ilc == 0 ) ilc = 10
462	cross_linecolors(k,j) = linecolors(ilc)
463	ils = MOD( cross_pnc_local(j), 11 )
464	IF ( ils == 0 ) ils = 11
465	cross_linestyles(k,j) = linestyles(ils)
466	!
467	!-- If required, extend x-value range of the respective coordinate
468	!-- cross, provided it has not been specified in check_parameters.
469	!-- Determination over all (sub-)regions.
470	IF ( cross_uxmin(j) == 0.0 .AND. cross_uxmax(j) == 0.0 ) THEN
471
472	DO sr = 0, statistic_regions
473
474	uxmi = MINVAL( hom_sum(:nz_do1d,dopr_index(i),sr) )
475	uxma = MAXVAL( hom_sum(:nz_do1d,dopr_index(i),sr) )
476	!
477	!-- When the value range of the first line in the
478	!-- corresponding cross is determined, its value range is
479	!-- simply adopted.
480	IF ( cross_uxmin_computed(j) > cross_uxmax_computed(j) ) &
481	THEN
482	cross_uxmin_computed(j) = uxmi
483	cross_uxmax_computed(j) = uxma
484	ELSE
485	cross_uxmin_computed(j) = &
486	MIN( cross_uxmin_computed(j), uxmi )
487	cross_uxmax_computed(j) = &
488	MAX( cross_uxmax_computed(j), uxma )
489	ENDIF
490
491	ENDDO
492
493	ENDIF
494	!
495	!-- If required, store the normalizing factors
496	SELECT CASE ( cross_normalized_x(j) )
497
498	CASE ( 'tsw2' )
499	cross_normx_factor(k,j) = &
500	( hom_sum(nzb+11,pr_palm,normalizing_region) )**2
501	CASE ( 'wpt0' )
502	cross_normx_factor(k,j) = &
503	hom_sum(nzb,18,normalizing_region)
504	CASE ( 'wstsw2' )
505	cross_normx_factor(k,j) = &
506	hom_sum(nzb+8,pr_palm,normalizing_region) &
507	* ( hom_sum(nzb+11,pr_palm,normalizing_region) )**2
508	CASE ( 'ws2' )
509	cross_normx_factor(k,j) = &
510	( hom_sum(nzb+8,pr_palm,normalizing_region) )**2
511	CASE ( 'ws2tsw' )
512	cross_normx_factor(k,j) = &
513	( hom_sum(nzb+8,pr_palm,normalizing_region) )**2&
514	* hom_sum(nzb+11,pr_palm,normalizing_region)
515	CASE ( 'ws3' )
516	cross_normx_factor(k,j) = &
517	( hom_sum(nzb+8,pr_palm,normalizing_region) )**3
518
519	END SELECT
520	SELECT CASE ( cross_normalized_y(j) )
521
522	CASE ( 'z_i' )
523	cross_normy_factor(k,j) = &
524	hom_sum(nzb+6,pr_palm,normalizing_region)
525
526	END SELECT
527
528	!
529	!-- Check the normalizing factors for zeros and deactivate the
530	!-- normalization, if required.
531	IF ( cross_normx_factor(k,j) == 0.0 .OR. &
532	cross_normy_factor(k,j) == 0.0 ) THEN
533	PRINT*,'+++ WARNING data_output_profiles: normalizing ',j, &
534	' cross is not possible since one of the'
535	PRINT*,' normalizing factors is zero!'
536	PRINT*,' cross_normx_factor(',k,',',j,') = ', &
537	cross_normx_factor(k,j)
538	PRINT*,' cross_normy_factor(',k,',',j,') = ', &
539	cross_normy_factor(k,j)
540	cross_normx_factor(k,j) = 1.0
541	cross_normy_factor(k,j) = 1.0
542	cross_normalized_x(j) = ' '
543	cross_normalized_y(j) = ' '
544	ENDIF
545
546	!
547	!-- If required, extend normalized x-value range of the respective
548	!-- cross, provided it has not been specified in check_parameters.
549	!-- Determination over all (sub-)regions.
550	IF ( cross_uxmin_normalized(j) == 0.0 .AND. &
551	cross_uxmax_normalized(j) == 0.0 ) THEN
552
553	DO sr = 0, statistic_regions
554
555	uxmi = MINVAL( hom(:nz_do1d,1,dopr_index(i),sr) ) / &
556	cross_normx_factor(k,j)
557	uxma = MAXVAL( hom(:nz_do1d,1,dopr_index(i),sr) ) / &
558	cross_normx_factor(k,j)
559	!
560	!-- When the value range of the first line in the
561	!-- corresponding cross is determined, its value range is
562	!-- simply adopted.
563	IF ( cross_uxmin_normalized_computed(j) > &
564	cross_uxmax_normalized_computed(j) ) THEN
565	cross_uxmin_normalized_computed(j) = uxmi
566	cross_uxmax_normalized_computed(j) = uxma
567	ELSE
568	cross_uxmin_normalized_computed(j) = &
569	MIN( cross_uxmin_normalized_computed(j), uxmi )
570	cross_uxmax_normalized_computed(j) = &
571	MAX( cross_uxmax_normalized_computed(j), uxma )
572	ENDIF
573
574	ENDDO
575
576	ENDIF
577
578	ENDIF ! Index determination
579
580	ENDIF ! profil output
581
582	ENDDO ! Loop over dopr_n
583
584	ENDIF ! Output on PE0
585
586	!
587	!-- If averaging has been done above, the summation counter must be re-set.
588	IF ( averaging_interval_pr /= 0.0 ) THEN
589	average_count_pr = 0
590	ENDIF
591
592	CALL cpu_log( log_point(15), 'data_output_profiles','stop', 'nobarrier' )
593
594	!
595	!-- Formats
596	100 FORMAT ('#1 ',A,1X,A)
597	101 FORMAT (E15.7,1X,E15.7)
598	102 FORMAT ('NEXT')
599
600	END SUBROUTINE data_output_profiles

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |