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

Last change on this file since 1783 was 1783, checked in by raasch, 6 years ago

NetCDF routines modularized; new parameter netcdf_deflate; further changes in the pmc

  • Property svn:keywords set to Id
File size: 22.2 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! name change of netcdf routines and module + related changes
22!
23! Former revisions:
24! -----------------
25! $Id: data_output_spectra.f90 1783 2016-03-06 18:36:17Z raasch $
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( __netcdf )
72#if defined( __spectra )
73
74    USE control_parameters,                                                    &
75        ONLY:  average_count_sp, averaging_interval_sp, dosp_time_count,       &
76               message_string, run_description_header,                         &
77               time_since_reference_point
78
79    USE cpulog,                                                                &
80        ONLY:  cpu_log, log_point
81
82    USE kinds
83
84    USE NETCDF
85
86    USE netcdf_interface,                                                      &
87        ONLY:  id_set_sp, id_var_time_sp, nc_stat, netcdf_handle_error
88
89    USE pegrid
90
91    USE spectrum,                                                              &
92        ONLY:  data_output_sp, spectra_direction
93
94    USE statistics,                                                            &
95        ONLY:  spectrum_x, spectrum_y
96
97
98    IMPLICIT NONE
99
100    INTEGER(iwp) ::  cranz_x !<
101    INTEGER(iwp) ::  cranz_y !<
102    INTEGER(iwp) ::  m       !<
103    INTEGER(iwp) ::  pr      !<
104   
105    LOGICAL      ::  frame_x !<
106    LOGICAL      ::  frame_y !<
107
108    CALL cpu_log( log_point(31), 'data_output_spectra', 'start' )
109
110!
111!-- Output is only performed on PE0
112    IF ( myid == 0 )  THEN
113
114!
115!--    Open file for spectra output in NetCDF format
116       CALL check_open( 107 )
117
118!
119!--    Increment the counter for number of output times
120       dosp_time_count = dosp_time_count + 1
121
122!
123!--    Update the spectra time axis
124       nc_stat = NF90_PUT_VAR( id_set_sp, id_var_time_sp,        &
125                               (/ time_since_reference_point /), &
126                               start = (/ dosp_time_count /), count = (/ 1 /) )
127       CALL netcdf_handle_error( 'data_output_spectra', 47 )
128
129!
130!--    If necessary, calculate time average and reset average counter
131       IF ( average_count_sp == 0 )  THEN
132           message_string = 'no spectra data available'
133           CALL message( 'data_output_spectra', 'PA0186', 0, 0, 0, 6, 0 )
134       ENDIF
135       IF ( average_count_sp /= 1 )  THEN
136          spectrum_x = spectrum_x / REAL( average_count_sp, KIND=wp )
137          spectrum_y = spectrum_y / REAL( average_count_sp, KIND=wp )
138          average_count_sp = 0
139       ENDIF
140
141!
142!--    Loop over all spectra defined by the user
143       m = 1
144       DO WHILE ( data_output_sp(m) /= ' '  .AND.  m <= 10 )
145
146          SELECT CASE ( TRIM( data_output_sp(m) ) )
147
148             CASE ( 'u' )
149                pr = 1
150
151             CASE ( 'v' )
152                pr = 2
153
154             CASE ( 'w' )
155                pr = 3
156
157             CASE ( 'pt' )
158                pr = 4
159
160             CASE ( 'q' )
161                pr = 5
162
163             CASE DEFAULT
164!
165!--             The DEFAULT case is reached either if the parameter
166!--             data_output_sp(m) contains a wrong character string or if the
167!--             user has coded a special case in the user interface. There, the
168!--             subroutine user_spectra checks which of these two conditions
169!--             applies.
170                CALL user_spectra( 'data_output', m, pr )
171
172          END SELECT
173
174!
175!--       Output of spectra in NetCDF format
176!--       Output of x-spectra
177          IF ( INDEX( spectra_direction(m), 'x' ) /= 0 ) THEN
178             CALL output_spectra_netcdf( m, 'x' )
179          ENDIF
180!
181!--       Output of y-spectra
182          IF ( INDEX( spectra_direction(m), 'y' ) /= 0 ) THEN
183             CALL output_spectra_netcdf( m, 'y' )
184          ENDIF
185
186!
187!--       Increase counter for next spectrum
188          m = m + 1
189
190       ENDDO
191
192!
193!--    Reset spectra values
194       spectrum_x = 0.0_wp; spectrum_y = 0.0_wp
195
196    ENDIF
197
198    CALL cpu_log( log_point(31), 'data_output_spectra', 'stop' )
199
200#if defined( __parallel )
201!    CALL MPI_BARRIER( comm2d, ierr )  ! really necessary
202#endif
203
204#endif
205#endif
206 END SUBROUTINE data_output_spectra
207
208
209!------------------------------------------------------------------------------!
210! Description:
211! ------------
212!> @todo Missing subroutine description.
213!------------------------------------------------------------------------------!
214 SUBROUTINE output_spectra_netcdf( nsp, direction )
215#if defined( __netcdf )
216
217    USE constants,                                                             &
218        ONLY:  pi
219
220    USE control_parameters,                                                    &
221        ONLY:  dosp_time_count
222
223    USE grid_variables,                                                        &
224        ONLY:  dx, dy
225
226    USE indices,                                                               &
227        ONLY:  nx, ny
228
229    USE kinds
230
231    USE NETCDF
232
233    USE netcdf_interface,                                                      &
234        ONLY:  id_set_sp, id_var_dospx, id_var_dospy, nc_stat,                 &
235               netcdf_handle_error
236
237    USE spectrum,                                                              &
238        ONLY:  n_sp_x, n_sp_y
239
240    USE statistics,                                                            &
241        ONLY:  spectrum_x, spectrum_y
242
243    IMPLICIT NONE
244
245    CHARACTER (LEN=1), INTENT(IN) ::  direction     !<
246
247    INTEGER(iwp), INTENT(IN)      ::  nsp           !<
248
249    INTEGER(iwp)                  ::  i             !<
250    INTEGER(iwp)                  ::  k             !<
251
252    REAL(wp)                      ::  frequency     !<
253
254    REAL(wp), DIMENSION(nx/2)     ::  netcdf_data_x !<
255    REAL(wp), DIMENSION(ny/2)     ::  netcdf_data_y !<
256
257
258    IF ( direction == 'x' )  THEN
259
260       DO  k = 1, n_sp_x
261
262          DO  i = 1, nx/2
263             frequency = 2.0_wp * pi * i / ( dx * ( nx + 1 ) )
264             netcdf_data_x(i) = frequency * spectrum_x(i,k,nsp)
265          ENDDO
266
267          nc_stat = NF90_PUT_VAR( id_set_sp, id_var_dospx(nsp), netcdf_data_x, &
268                                  start = (/ 1, k, dosp_time_count /), &
269                                  count = (/ nx/2, 1, 1 /) )
270          CALL netcdf_handle_error( 'data_output_spectra', 348 )
271
272       ENDDO
273
274    ENDIF
275
276    IF ( direction == 'y' )  THEN
277
278       DO  k = 1, n_sp_y
279
280          DO  i = 1, ny/2
281             frequency = 2.0_wp * pi * i / ( dy * ( ny + 1 ) )
282             netcdf_data_y(i) = frequency * spectrum_y(i,k,nsp)
283          ENDDO
284
285          nc_stat = NF90_PUT_VAR( id_set_sp, id_var_dospy(nsp), netcdf_data_y, &
286                                  start = (/ 1, k, dosp_time_count /), &
287                                  count = (/ ny/2, 1, 1 /) )
288          CALL netcdf_handle_error( 'data_output_spectra', 349 )
289
290       ENDDO
291
292    ENDIF
293
294#endif
295 END SUBROUTINE output_spectra_netcdf
296
297
298#if defined( __spectra )
299!------------------------------------------------------------------------------!
300! Description:
301! ------------
302!> @todo Missing subroutine description.
303!------------------------------------------------------------------------------!
304 SUBROUTINE data_output_spectra_x( m, cranz, pr, frame_written )
305
306    USE arrays_3d,                                                             &
307        ONLY:  zu, zw
308    USE constants,                                                             &
309        ONLY:  pi
310
311    USE control_parameters,                                                    &
312        ONLY:  averaging_interval_sp, run_description_header, simulated_time_chr
313
314    USE grid_variables,                                                        &
315        ONLY:  dx
316
317    USE indices,                                                               &
318        ONLY:  nx
319
320    USE kinds
321
322    USE pegrid
323
324    USE statistics,                                                            &
325        ONLY:  spectrum_x
326
327    USE spectrum,                                                              &
328        ONLY:  comp_spectra_level, header_char, lstyles, klist_x, n_sp_x,      &
329               plot_spectra_level, utext_char, ytext_char
330
331    IMPLICIT NONE
332
333    CHARACTER (LEN=30) ::  atext !<
334   
335    INTEGER(iwp)       ::  i     !<
336    INTEGER(iwp)       ::  j     !<
337    INTEGER(iwp)       ::  k     !<
338    INTEGER(iwp)       ::  m     !<
339    INTEGER(iwp)       ::  pr    !<
340   
341    LOGICAL            ::  frame_written   !<
342   
343    REAL(wp)           ::  frequency = 0.0_wp !<
344!
345!-- Variables needed for PROFIL-namelist
346    CHARACTER (LEN=80) ::  rtext                !<
347    CHARACTER (LEN=80) ::  utext                !<
348    CHARACTER (LEN=80) ::  xtext = 'k in m>->1' !<
349    CHARACTER (LEN=80) ::  ytext                !<
350
351    INTEGER(iwp)       ::  cranz       !<
352    INTEGER(iwp)       ::  labforx = 3 !<
353    INTEGER(iwp)       ::  labfory = 3 !<
354    INTEGER(iwp)       ::  legpos  = 3 !<
355    INTEGER(iwp)       ::  timodex = 1 !<
356   
357    INTEGER(iwp), DIMENSION(1:100) ::  cucol  = 1      !<
358    INTEGER(iwp), DIMENSION(1:100) ::  klist  = 999999 !<
359    INTEGER(iwp), DIMENSION(1:100) ::  lstyle = 0      !<
360   
361    LOGICAL ::  datleg = .TRUE. !<
362    LOGICAL ::  grid = .TRUE.   !<
363    LOGICAL ::  lclose = .TRUE. !<
364    LOGICAL ::  rand = .TRUE.   !<
365    LOGICAL ::  swap = .TRUE.   !<
366    LOGICAL ::  twoxa = .TRUE.  !<
367    LOGICAL ::  xlog = .TRUE.   !<
368    LOGICAL ::  ylog = .TRUE.   !<
369   
370    REAL(wp) ::  gwid = 0.1_wp    !<
371    REAL(wp) ::  rlegfak = 0.7_wp !<
372    REAL(wp) ::  uxmin            !<
373    REAL(wp) ::  uxmax            !<
374    REAL(wp) ::  uymin            !<
375    REAL(wp) ::  uymax            !<
376     
377    REAL(wp), DIMENSION(1:100) ::  lwid = 0.6_wp !<
378    REAL(wp), DIMENSION(100)   ::  uyma          !<
379    REAL(wp), DIMENSION(100)   ::  uymi          !<
380
381    NAMELIST /RAHMEN/  cranz, datleg, rtext, swap
382    NAMELIST /CROSS/   rand, cucol, grid, gwid, klist, labforx, labfory,      &
383                       legpos, lclose, lstyle, lwid, rlegfak, timodex, utext, &
384                       uxmin, uxmax, uymin, uymax, twoxa, xlog, xtext, ylog,  &
385                       ytext
386
387
388    rtext = '\0.5 ' // run_description_header
389
390!
391!-- Open parameter- and data-file
392    CALL check_open( 81 )
393    CALL check_open( 82 )
394
395!
396!-- Write file header,
397!-- write RAHMEN-parameters (pr=3: w-array is on zw, other arrys on zu,
398!-- pr serves as an index for output of strings (axis-labels) of the
399!-- different quantities u, v, w, pt and q)
400    DO  k = 1, n_sp_x
401       IF ( k < 100 )  THEN
402          IF ( pr == 3 )  THEN
403             WRITE ( 82, 100 )  '#', k, header_char( pr ),        &
404                                INT( zw(comp_spectra_level(k)) ), &
405                                simulated_time_chr
406          ELSE
407             WRITE ( 82, 100 )  '#', k, header_char( pr ),        &
408                                INT( zu(comp_spectra_level(k)) ), &
409                                simulated_time_chr
410          ENDIF
411       ELSE
412          IF ( pr == 3 )  THEN
413             WRITE ( 82, 101 )  '#', k, header_char( pr ),        &
414                                INT( zw(comp_spectra_level(k)) ), &
415                                simulated_time_chr
416          ELSE
417             WRITE ( 82, 101 )  '#', k, header_char( pr ),        &
418                                INT( zu(comp_spectra_level(k)) ), &
419                                simulated_time_chr
420          ENDIF
421       ENDIF
422    ENDDO
423
424    IF ( .NOT. frame_written )  THEN
425       WRITE ( 81, RAHMEN )
426       frame_written = .TRUE.
427    ENDIF
428
429!
430!-- Write all data and calculate uymi and uyma. They serve to calculate
431!-- the CROSS-parameters uymin and uymax
432    uymi = 999.999_wp; uyma = -999.999_wp
433    DO  i = 1, nx/2
434       frequency = 2.0_wp * pi * i / ( dx * ( nx + 1 ) )
435       WRITE ( 82, 102 )  frequency, ( frequency * spectrum_x(i,k,m), k = 1, &
436                          n_sp_x )
437       DO  k = 1, n_sp_x
438          uymi(k) = MIN( uymi(k), frequency * spectrum_x(i,k,m) )
439          uyma(k) = MAX( uyma(k), frequency * spectrum_x(i,k,m) )
440       ENDDO
441    ENDDO
442
443!
444!-- Determine CROSS-parameters
445    cucol(1:n_sp_x)  = (/ ( k, k = 1, n_sp_x ) /)
446    lstyle(1:n_sp_x) = (/ ( lstyles(k), k = 1, n_sp_x ) /)
447
448!
449!-- Calculate klist-values from the available comp_spectra_level values
450    i = 1; k = 1
451    DO WHILE ( i <= 100  .AND.  plot_spectra_level(i) /= 999999 )
452       DO WHILE ( k <= n_sp_x  .AND. &
453                  plot_spectra_level(i) >= comp_spectra_level(k) )
454          IF ( plot_spectra_level(i) == comp_spectra_level(k) )  THEN
455             klist(i) = k + klist_x
456          ELSE
457             uymi(k) =  999.999_wp
458             uyma(k) = -999.999_wp
459          ENDIF
460          k = k + 1
461       ENDDO
462       i = i + 1
463    ENDDO
464    uymi(k:n_sp_x) =  999.999_wp
465    uyma(k:n_sp_x) = -999.999_wp
466    utext = 'x'//utext_char( pr )
467    IF ( averaging_interval_sp /= 0.0_wp ) THEN
468       WRITE ( atext, 104 )  averaging_interval_sp
469       utext = TRIM(utext) // ',  ' // TRIM( atext )
470    ENDIF
471    uxmin = 0.8_wp * 2.0_wp * pi        / ( dx * ( nx + 1 ) )
472    uxmax = 1.2_wp * 2.0_wp * pi * nx/2 / ( dx * ( nx + 1 ) )
473    uymin = 0.8_wp * MIN (  999.999_wp, MINVAL ( uymi ) )
474    uymax = 1.2_wp * MAX ( -999.999_wp, MAXVAL ( uyma ) )
475    ytext = ytext_char( pr )
476
477!
478!-- Output of CROSS-parameters
479    WRITE ( 81, CROSS )
480
481!
482!-- Increase counter by the number of profiles written in the actual block
483    klist_x = klist_x + n_sp_x
484
485!
486!-- Write end-mark
487    WRITE ( 82, 103 )
488
489!
490!-- Close parameter- and data-file
491    CALL close_file( 81 )
492    CALL close_file( 82 )
493
494!
495!-- Formats
496100 FORMAT (A,I1,1X,A,1X,I4,'m ',A)
497101 FORMAT (A,I2,1X,A,1X,I4,'m ',A)
498102 FORMAT (E15.7,100(1X,E15.7))
499103 FORMAT ('NEXT')
500104 FORMAT ('time averaged over',F7.1,' s')
501
502 END SUBROUTINE data_output_spectra_x
503
504
505!------------------------------------------------------------------------------!
506! Description:
507! ------------
508!> @todo Missing subroutine description.
509!------------------------------------------------------------------------------!
510 SUBROUTINE data_output_spectra_y( m, cranz, pr, frame_written )
511
512    USE arrays_3d,                                                             &
513        ONLY:  zu, zw
514 
515    USE constants,                                                             &
516        ONLY:  pi
517
518    USE control_parameters,                                                    &
519        ONLY:  averaging_interval_sp, run_description_header, simulated_time_chr
520
521    USE grid_variables,                                                        &
522        ONLY:  dy
523
524    USE indices,                                                               &
525        ONLY:  ny
526
527    USE kinds
528
529    USE pegrid
530   
531    USE statistics,                                                             &
532        ONLY:  spectrum_y
533
534    USE spectrum,                                                              &
535        ONLY:  comp_spectra_level, header_char, klist_y, lstyles, n_sp_y,      &
536               plot_spectra_level, utext_char, ytext_char
537
538    IMPLICIT NONE
539
540   
541    CHARACTER (LEN=30) ::  atext !<
542   
543    INTEGER(iwp)       ::  i     !<
544    INTEGER(iwp)       ::  j     !<
545    INTEGER(iwp)       ::  k     !<
546    INTEGER(iwp)       ::  m     !<
547    INTEGER(iwp)       ::  pr    !<
548   
549    LOGICAL            :: frame_written   !<
550   
551    REAL(wp)           :: frequency = 0.0_wp !<
552
553!
554!-- Variables needed for PROFIL-namelist
555    CHARACTER (LEN=80) ::  rtext                !<
556    CHARACTER (LEN=80) ::  utext                !<
557    CHARACTER (LEN=80) ::  xtext = 'k in m>->1' !<
558    CHARACTER (LEN=80) ::  ytext                !<
559
560    INTEGER(iwp) ::  cranz       !<
561    INTEGER(iwp) ::  labforx = 3 !<
562    INTEGER(iwp) ::  labfory = 3 !<
563    INTEGER(iwp) ::  legpos  = 3 !<
564    INTEGER(iwp) ::  timodex = 1 !<
565   
566    INTEGER(iwp), DIMENSION(1:100) ::  cucol  = 1      !<
567    INTEGER(iwp), DIMENSION(1:100) ::  klist  = 999999 !<
568    INTEGER(iwp), DIMENSION(1:100) ::  lstyle = 0      !<
569   
570    LOGICAL ::  datleg = .TRUE. !<
571    LOGICAL ::  grid = .TRUE.   !<
572    LOGICAL ::  lclose = .TRUE. !<
573    LOGICAL ::  rand = .TRUE.   !<
574    LOGICAL ::  swap = .TRUE.   !<
575    LOGICAL ::  twoxa = .TRUE.  !<
576    LOGICAL ::  xlog = .TRUE.   !<
577    LOGICAL ::  ylog = .TRUE.   !<
578   
579    REAL(wp) ::  gwid = 0.1_wp     !<
580    REAL(wp) ::  rlegfak = 0.7_wp  !<
581    REAL(wp) ::  uxmin             !<
582    REAL(wp) ::  uxmax             !<
583    REAL(wp) ::  uymin             !<
584    REAL(wp) ::  uymax             !<
585   
586    REAL(wp), DIMENSION(1:100) ::  lwid = 0.6_wp !<
587   
588    REAL(wp), DIMENSION(100)   ::  uyma          !<
589    REAL(wp), DIMENSION(100)   ::  uymi          !<
590
591    NAMELIST /RAHMEN/  cranz, datleg, rtext, swap
592    NAMELIST /CROSS/   rand, cucol, grid, gwid, klist, labforx, labfory,      &
593                       legpos, lclose, lstyle, lwid, rlegfak, timodex, utext, &
594                       uxmin, uxmax, uymin, uymax, twoxa, xlog, xtext, ylog,  &
595                       ytext
596
597
598    rtext = '\0.5 ' // run_description_header
599
600!
601!-- Open parameter- and data-file
602    CALL check_open( 83 )
603    CALL check_open( 84 )
604
605!
606!-- Write file header,
607!-- write RAHMEN-parameters (pr=3: w-array is on zw, other arrys on zu,
608!-- pr serves as an index for output of strings (axis-labels) of the
609!-- different quantities u, v, w, pt and q)
610    DO  k = 1, n_sp_y
611       IF ( k < 100 )  THEN
612          IF ( pr == 3 ) THEN
613             WRITE ( 84, 100 )  '#', k, header_char( pr ),        &
614                                INT( zw(comp_spectra_level(k)) ), &
615                                simulated_time_chr
616          ELSE
617             WRITE ( 84, 100 )  '#', k, header_char( pr ),        &
618                                INT( zu(comp_spectra_level(k)) ), &
619                                simulated_time_chr
620          ENDIF
621       ELSE
622          IF ( pr == 3 )  THEN
623             WRITE ( 84, 101 )  '#', k, header_char( pr ),        &
624                                INT( zw(comp_spectra_level(k)) ), &
625                                simulated_time_chr
626          ELSE
627             WRITE ( 84, 101 )  '#', k, header_char( pr ),        &
628                                INT( zu(comp_spectra_level(k)) ), &
629                                simulated_time_chr
630          ENDIF
631       ENDIF
632    ENDDO
633
634    IF ( .NOT. frame_written )  THEN
635       WRITE ( 83, RAHMEN )
636       frame_written = .TRUE.
637    ENDIF
638
639!
640!-- Write all data and calculate uymi and uyma. They serve to calculate
641!-- the CROSS-parameters uymin and uymax
642    uymi = 999.999_wp; uyma = -999.999_wp
643    DO  j = 1, ny/2
644       frequency = 2.0_wp * pi * j / ( dy * ( ny + 1 ) )
645       WRITE ( 84, 102 ) frequency, ( frequency * spectrum_y(j,k,m), &
646                                      k = 1, n_sp_y ) 
647       DO k = 1, n_sp_y
648          uymi(k) = MIN( uymi(k), frequency * spectrum_y(j,k,m) )
649          uyma(k) = MAX( uyma(k), frequency * spectrum_y(j,k,m) )
650       ENDDO
651    ENDDO
652
653!
654!-- Determine CROSS-parameters
655    cucol(1:n_sp_y)  = (/ ( k, k = 1, n_sp_y ) /)
656    lstyle(1:n_sp_y) = (/ ( lstyles(k), k = 1, n_sp_y ) /)
657
658!
659!-- Calculate klist-values from the available comp_spectra_level values
660    j = 1; k = 1
661    DO WHILE ( j <= 100  .AND.  plot_spectra_level(j) /= 999999 )
662       DO WHILE ( k <= n_sp_y  .AND. &
663                  plot_spectra_level(j) >= comp_spectra_level(k) )
664          IF ( plot_spectra_level(j) == comp_spectra_level(k) )  THEN
665             klist(j) = k + klist_y
666          ELSE
667             uymi(k) =  999.999_wp
668             uyma(k) = -999.999_wp
669          ENDIF
670          k = k + 1
671       ENDDO
672       j = j + 1
673    ENDDO
674    uymi(k:n_sp_y) =  999.999_wp
675    uyma(k:n_sp_y) = -999.999_wp
676    utext = 'y'//utext_char( pr )
677    IF ( averaging_interval_sp /= 0.0_wp )  THEN
678       WRITE ( atext, 104 )  averaging_interval_sp
679       utext = TRIM(utext) // ',  ' // TRIM( atext )
680    ENDIF
681    uxmin = 0.8_wp * 2.0_wp * pi        / ( dy * ( ny + 1 ) )
682    uxmax = 1.2_wp * 2.0_wp * pi * ny/2 / ( dy * ( ny + 1 ) )
683    uymin = 0.8_wp * MIN (  999.999_wp, MINVAL ( uymi ) )
684    uymax = 1.2_wp * MAX ( -999.999_wp, MAXVAL ( uyma ) )
685    ytext = ytext_char( pr )
686
687!
688!-- Output CROSS-parameters
689    WRITE ( 83, CROSS )
690
691!
692!-- Increase counter by the number of profiles written in the actual block
693    klist_y = klist_y + n_sp_y
694
695!
696!-- Write end-mark
697    WRITE ( 84, 103 ) 
698
699!
700!-- Close parameter- and data-file
701    CALL close_file( 83 )
702    CALL close_file( 84 )
703
704!
705!-- Formats
706100 FORMAT (A,I1,1X,A,1X,I4,'m ',A)
707101 FORMAT (A,I2,1X,A,1X,I4,'m ',A)
708102 FORMAT (E15.7,100(1X,E15.7))
709103 FORMAT ('NEXT')
710104 FORMAT ('time averaged over',F7.1,' s')
711
712 END SUBROUTINE data_output_spectra_y
713#endif
Note: See TracBrowser for help on using the repository browser.