source: palm/trunk/SOURCE/data_output_spectra.f90 @ 1729

Last change on this file since 1729 was 1683, checked in by knoop, 9 years ago

last commit documented

  • Property svn:keywords set to Id
File size: 21.8 KB
Line 
1!> @file data_output_spectra.f90
2!--------------------------------------------------------------------------------!
3! This file is part of PALM.
4!
5! PALM is free software: you can redistribute it and/or modify it under the terms
6! of the GNU General Public License as published by the Free Software Foundation,
7! either version 3 of the License, or (at your option) any later version.
8!
9! PALM is distributed in the hope that it will be useful, but WITHOUT ANY
10! WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
11! A PARTICULAR PURPOSE.  See the GNU General Public License for more details.
12!
13! You should have received a copy of the GNU General Public License along with
14! PALM. If not, see <http://www.gnu.org/licenses/>.
15!
16! Copyright 1997-2014 Leibniz Universitaet Hannover
17!--------------------------------------------------------------------------------!
18!
19! Current revisions:
20! ------------------
21!
22!
23! Former revisions:
24! -----------------
25! $Id: data_output_spectra.f90 1683 2015-10-07 23:57:51Z gronemeier $
26!
27! 1682 2015-10-07 23:56:08Z knoop
28! Code annotations made doxygen readable
29!
30! 1353 2014-04-08 15:21:23Z heinze
31! REAL constants provided with KIND-attribute
32!
33! 1327 2014-03-21 11:00:16Z raasch
34! -netcdf output queries
35!
36! 1324 2014-03-21 09:13:16Z suehring
37! Bugfix: module statistics and module spectrum added, missing variables in ONLY
38! arguments added
39!
40! 1322 2014-03-20 16:38:49Z raasch
41! REAL functions provided with KIND-attribute
42!
43! 1320 2014-03-20 08:40:49Z raasch
44! ONLY-attribute added to USE-statements,
45! kind-parameters added to all INTEGER and REAL declaration statements,
46! kinds are defined in new module kinds,
47! revision history before 2012 removed,
48! comment fields (!:) to be used for variable explanations added to
49! all variable declaration statements
50!
51! 1318 2014-03-17 13:35:16Z raasch
52! module interfaces removed
53!
54! 1036 2012-10-22 13:43:42Z raasch
55! code put under GPL (PALM 3.9)
56!
57! 964 2012-07-26 09:14:24Z raasch
58! code for profil-output removed
59!
60! Revision 1.1  2001/01/05 15:14:20  raasch
61! Initial revision
62!
63!
64! Description:
65! ------------
66!> Writing spectra data on file, using a special format which allows
67!> plotting of these data with PROFIL-graphic-software
68!------------------------------------------------------------------------------!
69 SUBROUTINE data_output_spectra
70 
71#if defined( __spectra )
72
73    USE control_parameters,                                                    &
74        ONLY:  average_count_sp, averaging_interval_sp, dosp_time_count,       &
75               message_string, run_description_header,                         &
76               time_since_reference_point
77
78    USE cpulog,                                                                &
79        ONLY:  cpu_log, log_point
80
81    USE kinds
82
83    USE netcdf_control
84
85    USE pegrid
86
87    USE spectrum,                                                              &
88        ONLY:  data_output_sp, spectra_direction
89
90    USE statistics,                                                            &
91        ONLY:  spectrum_x, spectrum_y
92
93
94    IMPLICIT NONE
95
96    INTEGER(iwp) ::  cranz_x !<
97    INTEGER(iwp) ::  cranz_y !<
98    INTEGER(iwp) ::  m       !<
99    INTEGER(iwp) ::  pr      !<
100   
101    LOGICAL      ::  frame_x !<
102    LOGICAL      ::  frame_y !<
103
104    CALL cpu_log( log_point(31), 'data_output_spectra', 'start' )
105
106!
107!-- Output is only performed on PE0
108    IF ( myid == 0 )  THEN
109
110!
111!--    Open file for spectra output in NetCDF format
112       CALL check_open( 107 )
113
114!
115!--    Increment the counter for number of output times
116       dosp_time_count = dosp_time_count + 1
117
118#if defined( __netcdf )
119!
120!--    Update the spectra time axis
121       nc_stat = NF90_PUT_VAR( id_set_sp, id_var_time_sp,        &
122                               (/ time_since_reference_point /), &
123                               start = (/ dosp_time_count /), count = (/ 1 /) )
124       CALL handle_netcdf_error( 'data_output_spectra', 47 )
125#endif
126
127!
128!--    If necessary, calculate time average and reset average counter
129       IF ( average_count_sp == 0 )  THEN
130           message_string = 'no spectra data available'
131           CALL message( 'data_output_spectra', 'PA0186', 0, 0, 0, 6, 0 )
132       ENDIF
133       IF ( average_count_sp /= 1 )  THEN
134          spectrum_x = spectrum_x / REAL( average_count_sp, KIND=wp )
135          spectrum_y = spectrum_y / REAL( average_count_sp, KIND=wp )
136          average_count_sp = 0
137       ENDIF
138
139!
140!--    Loop over all spectra defined by the user
141       m = 1
142       DO WHILE ( data_output_sp(m) /= ' '  .AND.  m <= 10 )
143
144          SELECT CASE ( TRIM( data_output_sp(m) ) )
145
146             CASE ( 'u' )
147                pr = 1
148
149             CASE ( 'v' )
150                pr = 2
151
152             CASE ( 'w' )
153                pr = 3
154
155             CASE ( 'pt' )
156                pr = 4
157
158             CASE ( 'q' )
159                pr = 5
160
161             CASE DEFAULT
162!
163!--             The DEFAULT case is reached either if the parameter
164!--             data_output_sp(m) contains a wrong character string or if the
165!--             user has coded a special case in the user interface. There, the
166!--             subroutine user_spectra checks which of these two conditions
167!--             applies.
168                CALL user_spectra( 'data_output', m, pr )
169
170          END SELECT
171
172!
173!--       Output of spectra in NetCDF format
174!--       Output of x-spectra
175          IF ( INDEX( spectra_direction(m), 'x' ) /= 0 ) THEN
176             CALL output_spectra_netcdf( m, 'x' )
177          ENDIF
178!
179!--       Output of y-spectra
180          IF ( INDEX( spectra_direction(m), 'y' ) /= 0 ) THEN
181             CALL output_spectra_netcdf( m, 'y' )
182          ENDIF
183
184!
185!--       Increase counter for next spectrum
186          m = m + 1
187
188       ENDDO
189
190!
191!--    Reset spectra values
192       spectrum_x = 0.0_wp; spectrum_y = 0.0_wp
193
194    ENDIF
195
196    CALL cpu_log( log_point(31), 'data_output_spectra', 'stop' )
197
198#if defined( __parallel )
199!    CALL MPI_BARRIER( comm2d, ierr )  ! really necessary
200#endif
201
202#endif
203 END SUBROUTINE data_output_spectra
204
205
206!------------------------------------------------------------------------------!
207! Description:
208! ------------
209!> @todo Missing subroutine description.
210!------------------------------------------------------------------------------!
211 SUBROUTINE output_spectra_netcdf( nsp, direction )
212#if defined( __netcdf )
213
214    USE constants,                                                             &
215        ONLY:  pi
216
217    USE control_parameters,                                                    &
218        ONLY:  dosp_time_count
219
220    USE grid_variables,                                                        &
221        ONLY:  dx, dy
222
223    USE indices,                                                               &
224        ONLY:  nx, ny
225
226    USE kinds
227
228    USE netcdf_control
229
230    USE spectrum,                                                              &
231        ONLY:  n_sp_x, n_sp_y
232
233    USE statistics,                                                            &
234        ONLY:  spectrum_x, spectrum_y
235
236    IMPLICIT NONE
237
238    CHARACTER (LEN=1), INTENT(IN) ::  direction     !<
239
240    INTEGER(iwp), INTENT(IN)      ::  nsp           !<
241
242    INTEGER(iwp)                  ::  i             !<
243    INTEGER(iwp)                  ::  k             !<
244
245    REAL(wp)                      ::  frequency     !<
246
247    REAL(wp), DIMENSION(nx/2)     ::  netcdf_data_x !<
248    REAL(wp), DIMENSION(ny/2)     ::  netcdf_data_y !<
249
250
251    IF ( direction == 'x' )  THEN
252
253       DO  k = 1, n_sp_x
254
255          DO  i = 1, nx/2
256             frequency = 2.0_wp * pi * i / ( dx * ( nx + 1 ) )
257             netcdf_data_x(i) = frequency * spectrum_x(i,k,nsp)
258          ENDDO
259
260          nc_stat = NF90_PUT_VAR( id_set_sp, id_var_dospx(nsp), netcdf_data_x, &
261                                  start = (/ 1, k, dosp_time_count /), &
262                                  count = (/ nx/2, 1, 1 /) )
263          CALL handle_netcdf_error( 'data_output_spectra', 348 )
264
265       ENDDO
266
267    ENDIF
268
269    IF ( direction == 'y' )  THEN
270
271       DO  k = 1, n_sp_y
272
273          DO  i = 1, ny/2
274             frequency = 2.0_wp * pi * i / ( dy * ( ny + 1 ) )
275             netcdf_data_y(i) = frequency * spectrum_y(i,k,nsp)
276          ENDDO
277
278          nc_stat = NF90_PUT_VAR( id_set_sp, id_var_dospy(nsp), netcdf_data_y, &
279                                  start = (/ 1, k, dosp_time_count /), &
280                                  count = (/ ny/2, 1, 1 /) )
281          CALL handle_netcdf_error( 'data_output_spectra', 349 )
282
283       ENDDO
284
285    ENDIF
286
287#endif
288 END SUBROUTINE output_spectra_netcdf
289
290
291#if defined( __spectra )
292!------------------------------------------------------------------------------!
293! Description:
294! ------------
295!> @todo Missing subroutine description.
296!------------------------------------------------------------------------------!
297 SUBROUTINE data_output_spectra_x( m, cranz, pr, frame_written )
298
299    USE arrays_3d,                                                             &
300        ONLY:  zu, zw
301    USE constants,                                                             &
302        ONLY:  pi
303
304    USE control_parameters,                                                    &
305        ONLY:  averaging_interval_sp, run_description_header, simulated_time_chr
306
307    USE grid_variables,                                                        &
308        ONLY:  dx
309
310    USE indices,                                                               &
311        ONLY:  nx
312
313    USE kinds
314
315    USE pegrid
316
317    USE statistics,                                                            &
318        ONLY:  spectrum_x
319
320    USE spectrum,                                                              &
321        ONLY:  comp_spectra_level, header_char, lstyles, klist_x, n_sp_x,      &
322               plot_spectra_level, utext_char, ytext_char
323
324    IMPLICIT NONE
325
326    CHARACTER (LEN=30) ::  atext !<
327   
328    INTEGER(iwp)       ::  i     !<
329    INTEGER(iwp)       ::  j     !<
330    INTEGER(iwp)       ::  k     !<
331    INTEGER(iwp)       ::  m     !<
332    INTEGER(iwp)       ::  pr    !<
333   
334    LOGICAL            ::  frame_written   !<
335   
336    REAL(wp)           ::  frequency = 0.0_wp !<
337!
338!-- Variables needed for PROFIL-namelist
339    CHARACTER (LEN=80) ::  rtext                !<
340    CHARACTER (LEN=80) ::  utext                !<
341    CHARACTER (LEN=80) ::  xtext = 'k in m>->1' !<
342    CHARACTER (LEN=80) ::  ytext                !<
343
344    INTEGER(iwp)       ::  cranz       !<
345    INTEGER(iwp)       ::  labforx = 3 !<
346    INTEGER(iwp)       ::  labfory = 3 !<
347    INTEGER(iwp)       ::  legpos  = 3 !<
348    INTEGER(iwp)       ::  timodex = 1 !<
349   
350    INTEGER(iwp), DIMENSION(1:100) ::  cucol  = 1      !<
351    INTEGER(iwp), DIMENSION(1:100) ::  klist  = 999999 !<
352    INTEGER(iwp), DIMENSION(1:100) ::  lstyle = 0      !<
353   
354    LOGICAL ::  datleg = .TRUE. !<
355    LOGICAL ::  grid = .TRUE.   !<
356    LOGICAL ::  lclose = .TRUE. !<
357    LOGICAL ::  rand = .TRUE.   !<
358    LOGICAL ::  swap = .TRUE.   !<
359    LOGICAL ::  twoxa = .TRUE.  !<
360    LOGICAL ::  xlog = .TRUE.   !<
361    LOGICAL ::  ylog = .TRUE.   !<
362   
363    REAL(wp) ::  gwid = 0.1_wp    !<
364    REAL(wp) ::  rlegfak = 0.7_wp !<
365    REAL(wp) ::  uxmin            !<
366    REAL(wp) ::  uxmax            !<
367    REAL(wp) ::  uymin            !<
368    REAL(wp) ::  uymax            !<
369     
370    REAL(wp), DIMENSION(1:100) ::  lwid = 0.6_wp !<
371    REAL(wp), DIMENSION(100)   ::  uyma          !<
372    REAL(wp), DIMENSION(100)   ::  uymi          !<
373
374    NAMELIST /RAHMEN/  cranz, datleg, rtext, swap
375    NAMELIST /CROSS/   rand, cucol, grid, gwid, klist, labforx, labfory,      &
376                       legpos, lclose, lstyle, lwid, rlegfak, timodex, utext, &
377                       uxmin, uxmax, uymin, uymax, twoxa, xlog, xtext, ylog,  &
378                       ytext
379
380
381    rtext = '\0.5 ' // run_description_header
382
383!
384!-- Open parameter- and data-file
385    CALL check_open( 81 )
386    CALL check_open( 82 )
387
388!
389!-- Write file header,
390!-- write RAHMEN-parameters (pr=3: w-array is on zw, other arrys on zu,
391!-- pr serves as an index for output of strings (axis-labels) of the
392!-- different quantities u, v, w, pt and q)
393    DO  k = 1, n_sp_x
394       IF ( k < 100 )  THEN
395          IF ( pr == 3 )  THEN
396             WRITE ( 82, 100 )  '#', k, header_char( pr ),        &
397                                INT( zw(comp_spectra_level(k)) ), &
398                                simulated_time_chr
399          ELSE
400             WRITE ( 82, 100 )  '#', k, header_char( pr ),        &
401                                INT( zu(comp_spectra_level(k)) ), &
402                                simulated_time_chr
403          ENDIF
404       ELSE
405          IF ( pr == 3 )  THEN
406             WRITE ( 82, 101 )  '#', k, header_char( pr ),        &
407                                INT( zw(comp_spectra_level(k)) ), &
408                                simulated_time_chr
409          ELSE
410             WRITE ( 82, 101 )  '#', k, header_char( pr ),        &
411                                INT( zu(comp_spectra_level(k)) ), &
412                                simulated_time_chr
413          ENDIF
414       ENDIF
415    ENDDO
416
417    IF ( .NOT. frame_written )  THEN
418       WRITE ( 81, RAHMEN )
419       frame_written = .TRUE.
420    ENDIF
421
422!
423!-- Write all data and calculate uymi and uyma. They serve to calculate
424!-- the CROSS-parameters uymin and uymax
425    uymi = 999.999_wp; uyma = -999.999_wp
426    DO  i = 1, nx/2
427       frequency = 2.0_wp * pi * i / ( dx * ( nx + 1 ) )
428       WRITE ( 82, 102 )  frequency, ( frequency * spectrum_x(i,k,m), k = 1, &
429                          n_sp_x )
430       DO  k = 1, n_sp_x
431          uymi(k) = MIN( uymi(k), frequency * spectrum_x(i,k,m) )
432          uyma(k) = MAX( uyma(k), frequency * spectrum_x(i,k,m) )
433       ENDDO
434    ENDDO
435
436!
437!-- Determine CROSS-parameters
438    cucol(1:n_sp_x)  = (/ ( k, k = 1, n_sp_x ) /)
439    lstyle(1:n_sp_x) = (/ ( lstyles(k), k = 1, n_sp_x ) /)
440
441!
442!-- Calculate klist-values from the available comp_spectra_level values
443    i = 1; k = 1
444    DO WHILE ( i <= 100  .AND.  plot_spectra_level(i) /= 999999 )
445       DO WHILE ( k <= n_sp_x  .AND. &
446                  plot_spectra_level(i) >= comp_spectra_level(k) )
447          IF ( plot_spectra_level(i) == comp_spectra_level(k) )  THEN
448             klist(i) = k + klist_x
449          ELSE
450             uymi(k) =  999.999_wp
451             uyma(k) = -999.999_wp
452          ENDIF
453          k = k + 1
454       ENDDO
455       i = i + 1
456    ENDDO
457    uymi(k:n_sp_x) =  999.999_wp
458    uyma(k:n_sp_x) = -999.999_wp
459    utext = 'x'//utext_char( pr )
460    IF ( averaging_interval_sp /= 0.0_wp ) THEN
461       WRITE ( atext, 104 )  averaging_interval_sp
462       utext = TRIM(utext) // ',  ' // TRIM( atext )
463    ENDIF
464    uxmin = 0.8_wp * 2.0_wp * pi        / ( dx * ( nx + 1 ) )
465    uxmax = 1.2_wp * 2.0_wp * pi * nx/2 / ( dx * ( nx + 1 ) )
466    uymin = 0.8_wp * MIN (  999.999_wp, MINVAL ( uymi ) )
467    uymax = 1.2_wp * MAX ( -999.999_wp, MAXVAL ( uyma ) )
468    ytext = ytext_char( pr )
469
470!
471!-- Output of CROSS-parameters
472    WRITE ( 81, CROSS )
473
474!
475!-- Increase counter by the number of profiles written in the actual block
476    klist_x = klist_x + n_sp_x
477
478!
479!-- Write end-mark
480    WRITE ( 82, 103 )
481
482!
483!-- Close parameter- and data-file
484    CALL close_file( 81 )
485    CALL close_file( 82 )
486
487!
488!-- Formats
489100 FORMAT (A,I1,1X,A,1X,I4,'m ',A)
490101 FORMAT (A,I2,1X,A,1X,I4,'m ',A)
491102 FORMAT (E15.7,100(1X,E15.7))
492103 FORMAT ('NEXT')
493104 FORMAT ('time averaged over',F7.1,' s')
494
495 END SUBROUTINE data_output_spectra_x
496
497
498!------------------------------------------------------------------------------!
499! Description:
500! ------------
501!> @todo Missing subroutine description.
502!------------------------------------------------------------------------------!
503 SUBROUTINE data_output_spectra_y( m, cranz, pr, frame_written )
504
505    USE arrays_3d,                                                             &
506        ONLY:  zu, zw
507 
508    USE constants,                                                             &
509        ONLY:  pi
510
511    USE control_parameters,                                                    &
512        ONLY:  averaging_interval_sp, run_description_header, simulated_time_chr
513
514    USE grid_variables,                                                        &
515        ONLY:  dy
516
517    USE indices,                                                               &
518        ONLY:  ny
519
520    USE kinds
521
522    USE pegrid
523   
524    USE statistics,                                                             &
525        ONLY:  spectrum_y
526
527    USE spectrum,                                                              &
528        ONLY:  comp_spectra_level, header_char, klist_y, lstyles, n_sp_y,      &
529               plot_spectra_level, utext_char, ytext_char
530
531    IMPLICIT NONE
532
533   
534    CHARACTER (LEN=30) ::  atext !<
535   
536    INTEGER(iwp)       ::  i     !<
537    INTEGER(iwp)       ::  j     !<
538    INTEGER(iwp)       ::  k     !<
539    INTEGER(iwp)       ::  m     !<
540    INTEGER(iwp)       ::  pr    !<
541   
542    LOGICAL            :: frame_written   !<
543   
544    REAL(wp)           :: frequency = 0.0_wp !<
545
546!
547!-- Variables needed for PROFIL-namelist
548    CHARACTER (LEN=80) ::  rtext                !<
549    CHARACTER (LEN=80) ::  utext                !<
550    CHARACTER (LEN=80) ::  xtext = 'k in m>->1' !<
551    CHARACTER (LEN=80) ::  ytext                !<
552
553    INTEGER(iwp) ::  cranz       !<
554    INTEGER(iwp) ::  labforx = 3 !<
555    INTEGER(iwp) ::  labfory = 3 !<
556    INTEGER(iwp) ::  legpos  = 3 !<
557    INTEGER(iwp) ::  timodex = 1 !<
558   
559    INTEGER(iwp), DIMENSION(1:100) ::  cucol  = 1      !<
560    INTEGER(iwp), DIMENSION(1:100) ::  klist  = 999999 !<
561    INTEGER(iwp), DIMENSION(1:100) ::  lstyle = 0      !<
562   
563    LOGICAL ::  datleg = .TRUE. !<
564    LOGICAL ::  grid = .TRUE.   !<
565    LOGICAL ::  lclose = .TRUE. !<
566    LOGICAL ::  rand = .TRUE.   !<
567    LOGICAL ::  swap = .TRUE.   !<
568    LOGICAL ::  twoxa = .TRUE.  !<
569    LOGICAL ::  xlog = .TRUE.   !<
570    LOGICAL ::  ylog = .TRUE.   !<
571   
572    REAL(wp) ::  gwid = 0.1_wp     !<
573    REAL(wp) ::  rlegfak = 0.7_wp  !<
574    REAL(wp) ::  uxmin             !<
575    REAL(wp) ::  uxmax             !<
576    REAL(wp) ::  uymin             !<
577    REAL(wp) ::  uymax             !<
578   
579    REAL(wp), DIMENSION(1:100) ::  lwid = 0.6_wp !<
580   
581    REAL(wp), DIMENSION(100)   ::  uyma          !<
582    REAL(wp), DIMENSION(100)   ::  uymi          !<
583
584    NAMELIST /RAHMEN/  cranz, datleg, rtext, swap
585    NAMELIST /CROSS/   rand, cucol, grid, gwid, klist, labforx, labfory,      &
586                       legpos, lclose, lstyle, lwid, rlegfak, timodex, utext, &
587                       uxmin, uxmax, uymin, uymax, twoxa, xlog, xtext, ylog,  &
588                       ytext
589
590
591    rtext = '\0.5 ' // run_description_header
592
593!
594!-- Open parameter- and data-file
595    CALL check_open( 83 )
596    CALL check_open( 84 )
597
598!
599!-- Write file header,
600!-- write RAHMEN-parameters (pr=3: w-array is on zw, other arrys on zu,
601!-- pr serves as an index for output of strings (axis-labels) of the
602!-- different quantities u, v, w, pt and q)
603    DO  k = 1, n_sp_y
604       IF ( k < 100 )  THEN
605          IF ( pr == 3 ) THEN
606             WRITE ( 84, 100 )  '#', k, header_char( pr ),        &
607                                INT( zw(comp_spectra_level(k)) ), &
608                                simulated_time_chr
609          ELSE
610             WRITE ( 84, 100 )  '#', k, header_char( pr ),        &
611                                INT( zu(comp_spectra_level(k)) ), &
612                                simulated_time_chr
613          ENDIF
614       ELSE
615          IF ( pr == 3 )  THEN
616             WRITE ( 84, 101 )  '#', k, header_char( pr ),        &
617                                INT( zw(comp_spectra_level(k)) ), &
618                                simulated_time_chr
619          ELSE
620             WRITE ( 84, 101 )  '#', k, header_char( pr ),        &
621                                INT( zu(comp_spectra_level(k)) ), &
622                                simulated_time_chr
623          ENDIF
624       ENDIF
625    ENDDO
626
627    IF ( .NOT. frame_written )  THEN
628       WRITE ( 83, RAHMEN )
629       frame_written = .TRUE.
630    ENDIF
631
632!
633!-- Write all data and calculate uymi and uyma. They serve to calculate
634!-- the CROSS-parameters uymin and uymax
635    uymi = 999.999_wp; uyma = -999.999_wp
636    DO  j = 1, ny/2
637       frequency = 2.0_wp * pi * j / ( dy * ( ny + 1 ) )
638       WRITE ( 84, 102 ) frequency, ( frequency * spectrum_y(j,k,m), &
639                                      k = 1, n_sp_y ) 
640       DO k = 1, n_sp_y
641          uymi(k) = MIN( uymi(k), frequency * spectrum_y(j,k,m) )
642          uyma(k) = MAX( uyma(k), frequency * spectrum_y(j,k,m) )
643       ENDDO
644    ENDDO
645
646!
647!-- Determine CROSS-parameters
648    cucol(1:n_sp_y)  = (/ ( k, k = 1, n_sp_y ) /)
649    lstyle(1:n_sp_y) = (/ ( lstyles(k), k = 1, n_sp_y ) /)
650
651!
652!-- Calculate klist-values from the available comp_spectra_level values
653    j = 1; k = 1
654    DO WHILE ( j <= 100  .AND.  plot_spectra_level(j) /= 999999 )
655       DO WHILE ( k <= n_sp_y  .AND. &
656                  plot_spectra_level(j) >= comp_spectra_level(k) )
657          IF ( plot_spectra_level(j) == comp_spectra_level(k) )  THEN
658             klist(j) = k + klist_y
659          ELSE
660             uymi(k) =  999.999_wp
661             uyma(k) = -999.999_wp
662          ENDIF
663          k = k + 1
664       ENDDO
665       j = j + 1
666    ENDDO
667    uymi(k:n_sp_y) =  999.999_wp
668    uyma(k:n_sp_y) = -999.999_wp
669    utext = 'y'//utext_char( pr )
670    IF ( averaging_interval_sp /= 0.0_wp )  THEN
671       WRITE ( atext, 104 )  averaging_interval_sp
672       utext = TRIM(utext) // ',  ' // TRIM( atext )
673    ENDIF
674    uxmin = 0.8_wp * 2.0_wp * pi        / ( dy * ( ny + 1 ) )
675    uxmax = 1.2_wp * 2.0_wp * pi * ny/2 / ( dy * ( ny + 1 ) )
676    uymin = 0.8_wp * MIN (  999.999_wp, MINVAL ( uymi ) )
677    uymax = 1.2_wp * MAX ( -999.999_wp, MAXVAL ( uyma ) )
678    ytext = ytext_char( pr )
679
680!
681!-- Output CROSS-parameters
682    WRITE ( 83, CROSS )
683
684!
685!-- Increase counter by the number of profiles written in the actual block
686    klist_y = klist_y + n_sp_y
687
688!
689!-- Write end-mark
690    WRITE ( 84, 103 ) 
691
692!
693!-- Close parameter- and data-file
694    CALL close_file( 83 )
695    CALL close_file( 84 )
696
697!
698!-- Formats
699100 FORMAT (A,I1,1X,A,1X,I4,'m ',A)
700101 FORMAT (A,I2,1X,A,1X,I4,'m ',A)
701102 FORMAT (E15.7,100(1X,E15.7))
702103 FORMAT ('NEXT')
703104 FORMAT ('time averaged over',F7.1,' s')
704
705 END SUBROUTINE data_output_spectra_y
706#endif
Note: See TracBrowser for help on using the repository browser.