source: palm/tags/release-3.4/SOURCE/data_output_profiles.f90 @ 2718

Last change on this file since 2718 was 90, checked in by raasch, 17 years ago

New:
---
Calculation and output of user-defined profiles. New &userpar parameters data_output_pr_user and max_pr_user.

check_parameters, flow_statistics, modules, parin, read_var_list, user_interface, write_var_list

Changed:


Division through dt_3d replaced by multiplication of the inverse. For performance optimisation, this is done in the loop calculating the divergence instead of using a seperate loop. (pres.f90) var_hom and var_sum renamed pr_palm.

data_output_profiles, flow_statistics, init_3d_model, modules, parin, pres, read_var_list, run_control, time_integration

Errors:


Bugfix: work_fft*_vec removed from some PRIVATE-declarations (poisfft).

Bugfix: field_chr renamed field_char (user_interface).

Bugfix: output of use_upstream_for_tke (header).

header, poisfft, user_interface

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