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

Last change on this file since 263 was 263, checked in by heinze, 16 years ago

Output of NetCDF messages with aid of message handling routine.

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