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

Last change on this file since 254 was 254, checked in by heinze, 12 years ago

Output of messages replaced by message handling routine.

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