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

Last change on this file since 328 was 291, checked in by raasch, 16 years ago

changes for coupling with independent precursor runs; z_i calculation with Sullivan criterion

  • Property svn:keywords set to Id
File size: 25.0 KB
RevLine 
[1]1 SUBROUTINE data_output_profiles
2
3!------------------------------------------------------------------------------!
[254]4! Current revisions:
[1]5! -----------------
[291]6! simulated_time in NetCDF output replaced by time_since_reference_point.
[263]7! Output of NetCDF messages with aid of message handling routine.
[254]8! Output of messages replaced by message handling routine.
[1]9!
[254]10!
[1]11! Former revisions:
12! -----------------
[3]13! $Id: data_output_profiles.f90 291 2009-04-16 12:07:26Z letzel $
[198]14!
15! 197 2008-09-16 15:29:03Z raasch
16! Time coordinate t=0 stored on netcdf-file only if an output is required for
17! this time for at least one of the profiles
18!
19! February 2007
[3]20! RCS Log replace by Id keyword, revision history cleaned up
21!
[90]22! 87 2007-05-22 15:46:47Z raasch
23! var_hom renamed pr_palm
24!
[1]25! Revision 1.18  2006/08/16 14:27:04  raasch
26! PRINT* statements for testing removed
27!
28! Revision 1.1  1997/09/12 06:28:48  raasch
29! Initial revision
30!
31!
32! Description:
33! ------------
34! Plot output of 1D-profiles for PROFIL
35!------------------------------------------------------------------------------!
36
37    USE control_parameters
38    USE cpulog
39    USE indices
40    USE interfaces
41    USE netcdf_control
42    USE pegrid
43    USE profil_parameter
44    USE statistics
45
46    IMPLICIT NONE
47
48
49    INTEGER ::  i, id, ilc, ils, j, k, sr
[197]50    LOGICAL ::  output_for_t0
[1]51    REAL    ::  uxma, uxmi
52
53
54!
55!-- If required, compute statistics
56    IF ( .NOT. flow_statistics_called )  CALL flow_statistics
57
58!
59!-- Flow_statistics has its own CPU time measurement
60    CALL cpu_log( log_point(15), 'data_output_profiles', 'start' )
61
62!
63!-- If required, compute temporal average
64    IF ( averaging_interval_pr == 0.0 )  THEN
65       hom_sum(:,:,:) = hom(:,1,:,:)
66    ELSE
67       IF ( average_count_pr > 0 )  THEN
68          hom_sum = hom_sum / REAL( average_count_pr )
69       ELSE
70!
71!--       This case may happen if dt_dopr is changed in the d3par-list of
72!--       a restart run
73          RETURN
74       ENDIF
75    ENDIF
76
77   
78    IF ( myid == 0 )  THEN
79
80!
81!--    Plot-output for each (sub-)region
82
83!
84!--    Open file for profile output in NetCDF format
85       IF ( netcdf_output )  THEN
86          CALL check_open( 104 )
87       ENDIF
88
89!
90!--    Open PROFIL-output files for each (sub-)region
91       IF ( profil_output )  THEN
92          DO  sr = 0, statistic_regions
93             CALL check_open( 40 + sr )
94          ENDDO
95       ENDIF
96
97!
98!--    Increment the counter for number of output times
99       dopr_time_count = dopr_time_count + 1
100
101!
102!--    Re-set to zero the counter for the number of profiles already written
103!--    at the current output time into the respective crosses
104       cross_pnc_local = 0
105
106!
107!--    Output of initial profiles
108       IF ( dopr_time_count == 1 )  THEN
109
110          IF ( netcdf_output )  THEN
111#if defined( __netcdf )
112!
[197]113!--          Store initial time (t=0) to time axis, but only if an output
114!--          is required for at least one of the profiles
115             output_for_t0 = .FALSE.
116             DO  i = 1, dopr_n
117                IF ( dopr_initial_index(i) /= 0 )  THEN
118                   nc_stat = NF90_PUT_VAR( id_set_pr, id_var_time_pr,  &
119                                           (/ 0.0 /), start = (/ 1 /), &
120                                           count = (/ 1 /) )
[263]121                   CALL handle_netcdf_error( 'data_output_profiles', 329 )
[197]122                   output_for_t0 = .TRUE.
123                   EXIT
124                ENDIF
125             ENDDO
[1]126
127!
128!--          Store normalization factors
129             nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(1), &  ! wpt0
130                                  (/ hom_sum(nzb,18,normalizing_region) /), &
131                                     start = (/ 1 /), count = (/ 1 /) )
[263]132             CALL handle_netcdf_error( 'data_output_profiles', 330 )
[1]133
134             nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(2), &  ! ws2
[87]135                        (/ hom_sum(nzb+8,pr_palm,normalizing_region)**2 /), &
[1]136                                     start = (/ 1 /), count = (/ 1 /) )
[263]137             CALL handle_netcdf_error( 'data_output_profiles', 331 )
[1]138             nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(3), &  ! tsw2
[87]139                        (/ hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), &
[1]140                                     start = (/ 1 /), count = (/ 1 /) )
[263]141             CALL handle_netcdf_error( 'data_output_profiles', 332 )
[1]142             nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(4), &  ! ws3
[87]143                        (/ hom_sum(nzb+8,pr_palm,normalizing_region)**3 /), &
[1]144                                     start = (/ 1 /), count = (/ 1 /) )
[263]145             CALL handle_netcdf_error( 'data_output_profiles', 333 )
[1]146
147             nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(5), &  ! ws2tsw
[87]148                        (/ hom_sum(nzb+8,pr_palm,normalizing_region)**3 *   &
149                           hom_sum(nzb+3,pr_palm,normalizing_region)    /), &
[1]150                                     start = (/ 1 /), count = (/ 1 /) )
[263]151             CALL handle_netcdf_error( 'data_output_profiles', 334 )
[1]152
153             nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(6), &  ! wstsw2
[87]154                        (/ hom_sum(nzb+8,pr_palm,normalizing_region) *      &
155                           hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), &
[1]156                                     start = (/ 1 /), count = (/ 1 /) )
[263]157             CALL handle_netcdf_error( 'data_output_profiles', 335 )
[1]158
159             nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(7), &  ! z_i
[87]160                           (/ hom_sum(nzb+6,pr_palm,normalizing_region) /), &
[1]161                                     start = (/ 1 /), count = (/ 1 /) )
[263]162             CALL handle_netcdf_error( 'data_output_profiles', 336 )
[1]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 /) )
[263]199                      CALL handle_netcdf_error( 'data_output_profiles', 337 )
[1]200#endif
201                   ENDIF
202
203                ENDDO
204
205                IF ( profil_output )  THEN
206!
207!--                Determine indices for later NAMELIST-output (s. below)
208                   profile_number = profile_number + 1
209                   j = dopr_crossindex(i)
210                   IF ( j /= 0 )  THEN
211                      cross_profile_number_count(j) = &
212                                               cross_profile_number_count(j) + 1
213                      k = cross_profile_number_count(j)
214                      cross_profile_numbers(k,j) = profile_number
215!
216!--                   Initial profiles are always drawn as solid lines in
217!--                   anti-background colour.
218                      cross_linecolors(k,j) = 1
219                      cross_linestyles(k,j) = 0
220!
221!--                   If required, extend x-value range of the respective
222!--                   cross, provided it has not been specified in &
223!--                   check_parameters. Determination over all (sub-)regions.
224                      IF ( cross_uxmin(j) == 0.0  .AND. &
225                           cross_uxmax(j) == 0.0 )  THEN
226
227                         DO  sr = 0, statistic_regions
228
229                            uxmi = &
230                            MINVAL( hom(:nz_do1d,1,dopr_initial_index(i),sr) )
231
232                            uxma = &
233                            MAXVAL( hom(:nz_do1d,1,dopr_initial_index(i),sr) )
234!
235!--                         When the value range of the first line in the
236!--                         corresponding cross is determined, its value range
237!--                         is simply adopted.
238                            IF ( cross_uxmin_computed(j) > &
239                                 cross_uxmax_computed(j) )  THEN
240                               cross_uxmin_computed(j) = uxmi
241                               cross_uxmax_computed(j) = uxma
242                            ELSE
243                               cross_uxmin_computed(j) = &
244                                            MIN( cross_uxmin_computed(j), uxmi )
245                               cross_uxmax_computed(j) = &
246                                            MAX( cross_uxmax_computed(j), uxma )
247                            ENDIF
248
249                         ENDDO
250
251                      ENDIF
252!
253!--                   If required, determine and note normalizing factors
254                      SELECT CASE ( cross_normalized_x(j) )
255
256                         CASE ( 'ts2' )
257                            cross_normx_factor(k,j) = &
[87]258                             ( hom_sum(nzb+3,pr_palm,normalizing_region) )**2
[1]259                         CASE ( 'wpt0' )
260                            cross_normx_factor(k,j) = &
261                             hom_sum(nzb,18,normalizing_region)
262                         CASE ( 'wsts2' )
263                            cross_normx_factor(k,j) = &
[87]264                             hom_sum(nzb+8,pr_palm,normalizing_region)  &
265                           * ( hom_sum(nzb+3,pr_palm,normalizing_region) )**2
[1]266                         CASE ( 'ws2' )
267                            cross_normx_factor(k,j) = &
[87]268                             ( hom_sum(nzb+8,pr_palm,normalizing_region) )**2
[1]269                         CASE ( 'ws2ts' )
270                            cross_normx_factor(k,j) = &
[87]271                           ( hom_sum(nzb+8,pr_palm,normalizing_region) )**2 &
272                           * hom_sum(nzb+3,pr_palm,normalizing_region)
[1]273                         CASE ( 'ws3' )
274                            cross_normx_factor(k,j) = &
[87]275                             ( hom_sum(nzb+8,pr_palm,normalizing_region) )**3
[1]276
277                      END SELECT
278
279                      SELECT CASE ( cross_normalized_y(j) )
280
281                         CASE ( 'z_i' )
282                            cross_normy_factor(k,j) = &
[87]283                                    hom_sum(nzb+6,pr_palm,normalizing_region)
[1]284
285                      END SELECT
286
287!
288!--                   Check the normalizing factors for zeros and deactivate
289!--                   the normalization, if required.
290                      IF ( cross_normx_factor(k,j) == 0.0  .OR. &
291                           cross_normy_factor(k,j) == 0.0 )  THEN
[274]292                         WRITE( message_string, * ) 'data_output_profiles: ',  &
293                              'normalizing cross ',j, ' is not possible ',     &
294                              'since one of the & normalizing factors ',       &
295                              'is zero! & cross_normx_factor(',k,',',j,') = ', &
296                                                      cross_normx_factor(k,j), &
297                              ' & cross_normy_factor(',k,',',j,') = ',         &
298                                                      cross_normy_factor(k,j)
299                         CALL message( 'data_output_profiles', 'PA0185',&
300                                                                0, 1, 0, 6, 0 )
[1]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
[197]349          IF ( netcdf_output  .AND.  output_for_t0 )  THEN
350             dopr_time_count = dopr_time_count + 1
351          ENDIF
[1]352
353       ENDIF   ! Initial profiles
354
355       IF ( netcdf_output )  THEN
356#if defined( __netcdf )
357!
358!--       Store time to time axis         
[291]359          nc_stat = NF90_PUT_VAR( id_set_pr, id_var_time_pr,        &
360                                  (/ time_since_reference_point /), &
361                                  start = (/ dopr_time_count /),    &
[1]362                                  count = (/ 1 /) )
[263]363          CALL handle_netcdf_error( 'data_output_profiles', 338 )
[1]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 /) )
[263]371          CALL handle_netcdf_error( 'data_output_profiles', 339 )
[1]372
373          nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(2), &  ! ws2
[87]374                        (/ hom_sum(nzb+8,pr_palm,normalizing_region)**2 /), &
[1]375                                  start = (/ dopr_time_count /),               &
376                                  count = (/ 1 /) )
[263]377          CALL handle_netcdf_error( 'data_output_profiles', 340 )
[1]378
379          nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(3), &  ! tsw2
[87]380                        (/ hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), &
[1]381                                  start = (/ dopr_time_count /),               &
382                                  count = (/ 1 /) )
[263]383          CALL handle_netcdf_error( 'data_output_profiles', 341 )
[1]384
385          nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(4), &  ! ws3
[87]386                        (/ hom_sum(nzb+8,pr_palm,normalizing_region)**3 /), &
[1]387                                  start = (/ dopr_time_count /),               &
388                                  count = (/ 1 /) )
[263]389          CALL handle_netcdf_error( 'data_output_profiles', 342 )
[1]390
391          nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(5), &  ! ws2tsw
[87]392                        (/ hom_sum(nzb+8,pr_palm,normalizing_region)**3 *   &
393                           hom_sum(nzb+3,pr_palm,normalizing_region)    /), &
[1]394                                  start = (/ dopr_time_count /),               &
395                                  count = (/ 1 /) )
[263]396          CALL handle_netcdf_error( 'data_output_profiles', 343 )
397         
[1]398          nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(6), &  ! wstsw2
[87]399                        (/ hom_sum(nzb+8,pr_palm,normalizing_region) *      &
400                           hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), &
[1]401                                  start = (/ dopr_time_count /),               &
402                                  count = (/ 1 /) )
[263]403          CALL handle_netcdf_error( 'data_output_profiles', 344 )
[1]404
405          nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(7), &  ! z_i
[87]406                           (/ hom_sum(nzb+6,pr_palm,normalizing_region) /), &
[1]407                                  start = (/ dopr_time_count /),               &
408                                  count = (/ 1 /) )
[263]409          CALL handle_netcdf_error( 'data_output_profiles', 345 )
[1]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 /) )
[263]445                CALL handle_netcdf_error( 'data_output_profiles', 346 )
[1]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) = &
[87]505                            ( hom_sum(nzb+11,pr_palm,normalizing_region) )**2
[1]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) = &
[87]511                              hom_sum(nzb+8,pr_palm,normalizing_region)  &
512                          * ( hom_sum(nzb+11,pr_palm,normalizing_region) )**2
[1]513                   CASE ( 'ws2' )
514                      cross_normx_factor(k,j) = &
[87]515                            ( hom_sum(nzb+8,pr_palm,normalizing_region) )**2
[1]516                   CASE ( 'ws2tsw' )
517                      cross_normx_factor(k,j) = &
[87]518                            ( hom_sum(nzb+8,pr_palm,normalizing_region) )**2&
519                            * hom_sum(nzb+11,pr_palm,normalizing_region)
[1]520                   CASE ( 'ws3' )
521                      cross_normx_factor(k,j) = &
[87]522                            ( hom_sum(nzb+8,pr_palm,normalizing_region) )**3
[1]523
524                END SELECT
525                SELECT CASE ( cross_normalized_y(j) )
526
527                   CASE ( 'z_i' )
528                      cross_normy_factor(k,j) = &
[87]529                                   hom_sum(nzb+6,pr_palm,normalizing_region)
[1]530
531                END SELECT
532
533!
[254]534!--             Check the normalizing factors for zeros and deactivate
535!--             the normalization, if required.
[1]536                IF ( cross_normx_factor(k,j) == 0.0  .OR. &
537                     cross_normy_factor(k,j) == 0.0 )  THEN
[274]538                   WRITE( message_string, * ) 'data_output_profiles: ',  &
539                              'normalizing cross ',j, ' is not possible ',     &
540                              'since one of the & normalizing factors ',       &
541                              'is zero! & 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',&
546                                                                0, 1, 0, 6, 0 )
[254]547                    cross_normx_factor(k,j) = 1.0
548                    cross_normy_factor(k,j) = 1.0
549                    cross_normalized_x(j) = ' '
550                    cross_normalized_y(j) = ' '
[1]551                ENDIF
552
553!
554!--             If required, extend normalized x-value range of the respective 
555!--             cross, provided it has not been specified in check_parameters.
556!--             Determination over all (sub-)regions.
557                IF ( cross_uxmin_normalized(j) == 0.0  .AND. &
558                     cross_uxmax_normalized(j) == 0.0 )  THEN
559
560                   DO  sr = 0, statistic_regions
561
562                      uxmi = MINVAL( hom(:nz_do1d,1,dopr_index(i),sr) ) / &
563                             cross_normx_factor(k,j)
564                      uxma = MAXVAL( hom(:nz_do1d,1,dopr_index(i),sr) ) / &
565                             cross_normx_factor(k,j)
566!
567!--                   When the value range of the first line in the
568!--                   corresponding cross is determined, its value range is
569!--                   simply adopted.
570                      IF ( cross_uxmin_normalized_computed(j) > &
571                           cross_uxmax_normalized_computed(j) )  THEN
572                         cross_uxmin_normalized_computed(j) = uxmi
573                         cross_uxmax_normalized_computed(j) = uxma
574                      ELSE
575                         cross_uxmin_normalized_computed(j) = &
576                                MIN( cross_uxmin_normalized_computed(j), uxmi )
577                         cross_uxmax_normalized_computed(j) = &
578                                MAX( cross_uxmax_normalized_computed(j), uxma )
579                      ENDIF
580
581                   ENDDO
582
583                ENDIF
584
585             ENDIF   ! Index determination
586
587          ENDIF   ! profil output
588
589       ENDDO   ! Loop over dopr_n
590
591    ENDIF  ! Output on PE0
592
593!
594!-- If averaging has been done above, the summation counter must be re-set.
595    IF ( averaging_interval_pr /= 0.0 )  THEN
596       average_count_pr = 0
597    ENDIF
598
599    CALL cpu_log( log_point(15), 'data_output_profiles','stop', 'nobarrier' )
600
601!
602!-- Formats
603100 FORMAT ('#1 ',A,1X,A)
604101 FORMAT (E15.7,1X,E15.7)
605102 FORMAT ('NEXT')
606
607 END SUBROUTINE data_output_profiles
Note: See TracBrowser for help on using the repository browser.