source: palm/trunk/SOURCE/header.f90 @ 2299

Last change on this file since 2299 was 2299, checked in by maronga, 4 years ago

improvements for spinup mechanism

  • Property svn:keywords set to Id
File size: 90.7 KB
Line 
1!> @file header.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
6! terms of the GNU General Public License as published by the Free Software
7! Foundation, either version 3 of the License, or (at your option) any later
8! version.
9!
10! PALM is distributed in the hope that it will be useful, but WITHOUT ANY
11! WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
12! A PARTICULAR PURPOSE.  See the GNU General Public License for more details.
13!
14! You should have received a copy of the GNU General Public License along with
15! PALM. If not, see <http://www.gnu.org/licenses/>.
16!
17! Copyright 1997-2017 Leibniz Universitaet Hannover
18!------------------------------------------------------------------------------!
19!
20! Current revisions:
21! -----------------
22!
23!
24! Former revisions:
25! -----------------
26! $Id: header.f90 2299 2017-06-29 10:14:38Z maronga $
27! Modified output for spinups
28!
29! 2298 2017-06-29 09:28:18Z raasch
30! MPI2 related parts removed
31!
32! 2270 2017-06-09 12:18:47Z maronga
33! Renamed Prandtl layer to constant flux layer
34!
35! 2259 2017-06-08 09:09:11Z gronemeier
36! Implemented synthetic turbulence generator
37!
38! 2258 2017-06-08 07:55:13Z suehring
39! Bugfix, add pre-preprocessor directives to enable non-parrallel mode
40!
41! 2233 2017-05-30 18:08:54Z suehring
42!
43! 2232 2017-05-30 17:47:52Z suehring
44! Adjustments to new topography and surface concept
45! Generic tunnel setup added
46!
47! 2200 2017-04-11 11:37:51Z suehring
48! monotonic_adjustment removed
49!
50! 2118 2017-01-17 16:38:49Z raasch
51! OpenACC relatec code removed
52!
53! 2073 2016-11-30 14:34:05Z raasch
54! small bugfix concerning output of scalar profiles
55!
56! 2050 2016-11-08 15:00:55Z gronemeier
57! Implement turbulent outflow condition
58!
59! 2037 2016-10-26 11:15:40Z knoop
60! Anelastic approximation implemented
61!
62! 2000 2016-08-20 18:09:15Z knoop
63! Forced header and separation lines into 80 columns
64!
65! 1992 2016-08-12 15:14:59Z suehring
66! Adapted for top_scalarflux
67!
68! 1960 2016-07-12 16:34:24Z suehring
69! Treat humidity and passive scalar separately.
70! Modify misleading information concerning humidity.
71! Bugfix, change unit for humidity flux.
72!
73! 1957 2016-07-07 10:43:48Z suehring
74! flight module added
75!
76! 1931 2016-06-10 12:06:59Z suehring
77! Rename multigrid into multigrid_noopt
78!
79! 1902 2016-05-09 11:18:56Z suehring
80! Write information about masking_method only for multigrid solver
81!
82! 1849 2016-04-08 11:33:18Z hoffmann
83! Adapted for modularization of microphysics
84!
85! 1833 2016-04-07 14:23:03Z raasch
86! spectrum renamed spectra_mod, output of spectra related quantities moved to
87! spectra_mod
88!
89! 1831 2016-04-07 13:15:51Z hoffmann
90! turbulence renamed collision_turbulence,
91! drizzle renamed cloud_water_sedimentation
92!
93! 1826 2016-04-07 12:01:39Z maronga
94! Moved radiation model header output to the respective module.
95! Moved canopy model header output to the respective module.
96!
97! 1822 2016-04-07 07:49:42Z hoffmann
98! Tails removed. icloud_scheme replaced by microphysics_*
99!
100! 1817 2016-04-06 15:44:20Z maronga
101! Moved land_surface_model header output to the respective module.
102!
103! 1808 2016-04-05 19:44:00Z raasch
104! routine local_flush replaced by FORTRAN statement
105!
106! 1797 2016-03-21 16:50:28Z raasch
107! output of nesting datatransfer mode
108!
109! 1791 2016-03-11 10:41:25Z raasch
110! output of nesting informations of all domains
111!
112! 1788 2016-03-10 11:01:04Z maronga
113! Parameter dewfall removed
114!
115! 1786 2016-03-08 05:49:27Z raasch
116! cpp-direktives for spectra removed
117!
118! 1783 2016-03-06 18:36:17Z raasch
119! netcdf module and variable names changed, output of netcdf_deflate
120!
121! 1764 2016-02-28 12:45:19Z raasch
122! output of nesting informations
123!
124! 1697 2015-10-28 17:14:10Z raasch
125! small E- and F-FORMAT changes to avoid informative compiler messages about
126! insufficient field width
127!
128! 1691 2015-10-26 16:17:44Z maronga
129! Renamed prandtl_layer to constant_flux_layer, renames rif_min/rif_max to
130! zeta_min/zeta_max.
131!
132! 1682 2015-10-07 23:56:08Z knoop
133! Code annotations made doxygen readable
134!
135! 1675 2015-10-02 08:28:59Z gronemeier
136! Bugfix: Definition of topography grid levels
137!
138! 1660 2015-09-21 08:15:16Z gronemeier
139! Bugfix: Definition of building/street canyon height if vertical grid stretching
140!         starts below the maximum topography height.
141!
142! 1590 2015-05-08 13:56:27Z maronga
143! Bugfix: Added TRIM statements for character strings for LSM and radiation code
144!
145! 1585 2015-04-30 07:05:52Z maronga
146! Further output for radiation model(s).
147!
148! 1575 2015-03-27 09:56:27Z raasch
149! adjustments for psolver-queries, output of seed_follows_topography
150!
151! 1560 2015-03-06 10:48:54Z keck
152! output for recycling y shift
153!
154! 1557 2015-03-05 16:43:04Z suehring
155! output for monotonic limiter
156!
157! 1551 2015-03-03 14:18:16Z maronga
158! Added informal output for land surface model and radiation model. Removed typo.
159!
160! 1496 2014-12-02 17:25:50Z maronga
161! Renamed: "radiation -> "cloud_top_radiation"
162!
163! 1484 2014-10-21 10:53:05Z kanani
164! Changes due to new module structure of the plant canopy model:
165!   module plant_canopy_model_mod and output for new canopy model parameters
166!   (alpha_lad, beta_lad, lai_beta,...) added,
167!   drag_coefficient, leaf_surface_concentration and scalar_exchange_coefficient
168!   renamed to canopy_drag_coeff, leaf_surface_conc and leaf_scalar_exch_coeff,
169!   learde renamed leaf_area_density.
170! Bugfix: DO-WHILE-loop for lad header information additionally restricted
171! by maximum number of gradient levels (currently 10)
172!
173! 1482 2014-10-18 12:34:45Z raasch
174! information about calculated or predefined virtual processor topology adjusted
175!
176! 1468 2014-09-24 14:06:57Z maronga
177! Adapted for use on up to 6-digit processor cores
178!
179! 1429 2014-07-15 12:53:45Z knoop
180! header exended to provide ensemble_member_nr if specified
181!
182! 1376 2014-04-26 11:21:22Z boeske
183! Correction of typos
184!
185! 1365 2014-04-22 15:03:56Z boeske
186! New section 'Large scale forcing and nudging':
187! output of large scale forcing and nudging information,
188! new section for initial profiles created
189!
190! 1359 2014-04-11 17:15:14Z hoffmann
191! dt_sort_particles removed
192!
193! 1353 2014-04-08 15:21:23Z heinze
194! REAL constants provided with KIND-attribute
195!
196! 1327 2014-03-21 11:00:16Z raasch
197! parts concerning iso2d and avs output removed,
198! -netcdf output queries
199!
200! 1324 2014-03-21 09:13:16Z suehring
201! Bugfix: module spectrum added
202!
203! 1322 2014-03-20 16:38:49Z raasch
204! REAL functions provided with KIND-attribute,
205! some REAL constants defined as wp-kind
206!
207! 1320 2014-03-20 08:40:49Z raasch
208! ONLY-attribute added to USE-statements,
209! kind-parameters added to all INTEGER and REAL declaration statements,
210! kinds are defined in new module kinds,
211! revision history before 2012 removed,
212! comment fields (!:) to be used for variable explanations added to
213! all variable declaration statements
214!
215! 1308 2014-03-13 14:58:42Z fricke
216! output of the fixed number of output time levels
217! output_format adjusted for masked data if netcdf_data_format > 5
218!
219! 1299 2014-03-06 13:15:21Z heinze
220! output for using large_scale subsidence in combination
221! with large_scale_forcing
222! reformatting, more detailed explanations
223!
224! 1241 2013-10-30 11:36:58Z heinze
225! output for nudging + large scale forcing from external file
226!
227! 1216 2013-08-26 09:31:42Z raasch
228! output for transpose_compute_overlap
229!
230! 1212 2013-08-15 08:46:27Z raasch
231! output for poisfft_hybrid removed
232!
233! 1179 2013-06-14 05:57:58Z raasch
234! output of reference_state, use_reference renamed use_single_reference_value
235!
236! 1159 2013-05-21 11:58:22Z fricke
237! +use_cmax
238!
239! 1115 2013-03-26 18:16:16Z hoffmann
240! descriptions for Seifert-Beheng-cloud-physics-scheme added
241!
242! 1111 2013-03-08 23:54:10Z raasch
243! output of accelerator board information
244! ibc_p_b = 2 removed
245!
246! 1108 2013-03-05 07:03:32Z raasch
247! bugfix for r1106
248!
249! 1106 2013-03-04 05:31:38Z raasch
250! some format changes for coupled runs
251!
252! 1092 2013-02-02 11:24:22Z raasch
253! unused variables removed
254!
255! 1036 2012-10-22 13:43:42Z raasch
256! code put under GPL (PALM 3.9)
257!
258! 1031 2012-10-19 14:35:30Z raasch
259! output of netCDF data format modified
260!
261! 1015 2012-09-27 09:23:24Z raasch
262! output of Adjustment of mixing length to the Prandtl mixing length at first
263! grid point above ground removed
264!
265! 1003 2012-09-14 14:35:53Z raasch
266! output of information about equal/unequal subdomain size removed
267!
268! 1001 2012-09-13 14:08:46Z raasch
269! all actions concerning leapfrog- and upstream-spline-scheme removed
270!
271! 978 2012-08-09 08:28:32Z fricke
272! -km_damp_max, outflow_damping_width
273! +pt_damping_factor, pt_damping_width
274! +z0h
275!
276! 964 2012-07-26 09:14:24Z raasch
277! output of profil-related quantities removed
278!
279! 940 2012-07-09 14:31:00Z raasch
280! Output in case of simulations for pure neutral stratification (no pt-equation
281! solved)
282!
283! 927 2012-06-06 19:15:04Z raasch
284! output of masking_method for mg-solver
285!
286! 868 2012-03-28 12:21:07Z raasch
287! translation velocity in Galilean transformation changed to 0.6 * ug
288!
289! 833 2012-02-22 08:55:55Z maronga
290! Adjusted format for leaf area density
291!
292! 828 2012-02-21 12:00:36Z raasch
293! output of dissipation_classes + radius_classes
294!
295! 825 2012-02-19 03:03:44Z raasch
296! Output of cloud physics parameters/quantities complemented and restructured
297!
298! Revision 1.1  1997/08/11 06:17:20  raasch
299! Initial revision
300!
301!
302! Description:
303! ------------
304!> Writing a header with all important information about the current run.
305!> This subroutine is called three times, two times at the beginning
306!> (writing information on files RUN_CONTROL and HEADER) and one time at the
307!> end of the run, then writing additional information about CPU-usage on file
308!> header.
309!-----------------------------------------------------------------------------!
310 SUBROUTINE header
311 
312
313    USE arrays_3d,                                                             &
314        ONLY:  pt_init, q_init, s_init, sa_init, ug, vg, w_subs, zu, zw
315       
316    USE control_parameters
317       
318    USE cloud_parameters,                                                      &
319        ONLY:  cp, l_v, r_d
320
321    USE cpulog,                                                                &
322        ONLY:  log_point_s
323       
324    USE dvrp_variables,                                                        &
325        ONLY:  use_seperate_pe_for_dvrp_output
326       
327    USE flight_mod,                                                            &
328        ONLY:  flight_header
329       
330    USE grid_variables,                                                        &
331        ONLY:  dx, dy
332       
333    USE indices,                                                               &
334        ONLY:  mg_loc_ind, nnx, nny, nnz, nx, ny, nxl_mg, nxr_mg, nyn_mg,      &
335               nys_mg, nzt, nzt_mg
336       
337    USE kinds
338 
339    USE land_surface_model_mod,                                                &
340        ONLY: lsm_header
341
342    USE microphysics_mod,                                                      &
343        ONLY:  cloud_water_sedimentation, collision_turbulence,                &
344               c_sedimentation, limiter_sedimentation, nc_const,               &
345               ventilation_effect
346
347    USE model_1d,                                                              &
348        ONLY:  damp_level_ind_1d, dt_pr_1d, dt_run_control_1d, end_time_1d
349       
350    USE netcdf_interface,                                                      &
351        ONLY:  netcdf_data_format, netcdf_data_format_string, netcdf_deflate
352
353    USE particle_attributes,                                                   &
354        ONLY:  bc_par_b, bc_par_lr, bc_par_ns, bc_par_t, collision_kernel,     &
355               curvature_solution_effects,                                     &
356               density_ratio, dissipation_classes, dt_min_part, dt_prel,       &
357               dt_write_particle_data, end_time_prel,                          &
358               number_of_particle_groups, particle_advection,                  &
359               particle_advection_start,                                       &
360               particles_per_point, pdx, pdy, pdz,  psb, psl, psn, psr, pss,   &
361               pst, radius, radius_classes, random_start_position,             &
362               seed_follows_topography,                                        &
363               total_number_of_particles, use_sgs_for_particles,               &
364               vertical_particle_advection, write_particle_statistics
365       
366    USE pegrid
367
368    USE plant_canopy_model_mod,                                                &
369        ONLY:  pcm_header, plant_canopy
370
371    USE pmc_handle_communicator,                                               &
372        ONLY:  pmc_get_model_info
373
374    USE pmc_interface,                                                         &
375        ONLY:  nested_run, nesting_datatransfer_mode, nesting_mode
376
377    USE radiation_model_mod,                                                   &
378        ONLY:  radiation, radiation_header
379   
380    USE spectra_mod,                                                           &
381        ONLY:  calculate_spectra, spectra_header
382
383    USE surface_mod,                                                           &
384        ONLY:  surf_def_h
385
386    USE synthetic_turbulence_generator_mod,                                    &
387        ONLY:  stg_header
388
389    IMPLICIT NONE
390
391    CHARACTER (LEN=1)  ::  prec                !<
392   
393    CHARACTER (LEN=2)  ::  do2d_mode           !<
394   
395    CHARACTER (LEN=5)  ::  section_chr         !<
396   
397    CHARACTER (LEN=10) ::  coor_chr            !<
398    CHARACTER (LEN=10) ::  host_chr            !<
399   
400    CHARACTER (LEN=16) ::  begin_chr           !<
401   
402    CHARACTER (LEN=26) ::  ver_rev             !<
403
404    CHARACTER (LEN=32) ::  cpl_name            !<
405   
406    CHARACTER (LEN=40) ::  output_format       !<
407   
408    CHARACTER (LEN=70) ::  char1               !<
409    CHARACTER (LEN=70) ::  char2               !<
410    CHARACTER (LEN=70) ::  dopr_chr            !<
411    CHARACTER (LEN=70) ::  do2d_xy             !<
412    CHARACTER (LEN=70) ::  do2d_xz             !<
413    CHARACTER (LEN=70) ::  do2d_yz             !<
414    CHARACTER (LEN=70) ::  do3d_chr            !<
415    CHARACTER (LEN=70) ::  domask_chr          !<
416    CHARACTER (LEN=70) ::  run_classification  !<
417   
418    CHARACTER (LEN=85) ::  r_upper             !<
419    CHARACTER (LEN=85) ::  r_lower             !<
420   
421    CHARACTER (LEN=86) ::  coordinates         !<
422    CHARACTER (LEN=86) ::  gradients           !<
423    CHARACTER (LEN=86) ::  slices              !<
424    CHARACTER (LEN=86) ::  temperatures        !<
425    CHARACTER (LEN=86) ::  ugcomponent         !<
426    CHARACTER (LEN=86) ::  vgcomponent         !<
427
428    CHARACTER (LEN=1), DIMENSION(1:3) ::  dir = (/ 'x', 'y', 'z' /)  !<
429
430    INTEGER(iwp) ::  av             !<
431    INTEGER(iwp) ::  bh             !<
432    INTEGER(iwp) ::  blx            !<
433    INTEGER(iwp) ::  bly            !<
434    INTEGER(iwp) ::  bxl            !<
435    INTEGER(iwp) ::  bxr            !<
436    INTEGER(iwp) ::  byn            !<
437    INTEGER(iwp) ::  bys            !<
438    INTEGER(iwp) ::  ch             !<
439    INTEGER(iwp) ::  count          !<
440    INTEGER(iwp) ::  cpl_parent_id  !<
441    INTEGER(iwp) ::  cwx            !<
442    INTEGER(iwp) ::  cwy            !<
443    INTEGER(iwp) ::  cxl            !<
444    INTEGER(iwp) ::  cxr            !<
445    INTEGER(iwp) ::  cyn            !<
446    INTEGER(iwp) ::  cys            !<
447    INTEGER(iwp) ::  dim            !<
448    INTEGER(iwp) ::  i              !<
449    INTEGER(iwp) ::  io             !<
450    INTEGER(iwp) ::  j              !<
451    INTEGER(iwp) ::  k              !<
452    INTEGER(iwp) ::  l              !<
453    INTEGER(iwp) ::  ll             !<
454    INTEGER(iwp) ::  my_cpl_id      !<
455    INTEGER(iwp) ::  n              !<
456    INTEGER(iwp) ::  ncpl           !<
457    INTEGER(iwp) ::  npe_total      !<
458   
459
460    REAL(wp) ::  cpuseconds_per_simulated_second  !<
461    REAL(wp) ::  lower_left_coord_x               !< x-coordinate of nest domain
462    REAL(wp) ::  lower_left_coord_y               !< y-coordinate of nest domain
463
464!
465!-- Open the output file. At the end of the simulation, output is directed
466!-- to unit 19.
467    IF ( ( runnr == 0 .OR. force_print_header )  .AND. &
468         .NOT. simulated_time_at_begin /= simulated_time )  THEN
469       io = 15   !  header output on file RUN_CONTROL
470    ELSE
471       io = 19   !  header output on file HEADER
472    ENDIF
473    CALL check_open( io )
474
475!
476!-- At the end of the run, output file (HEADER) will be rewritten with
477!-- new information
478    IF ( io == 19 .AND. simulated_time_at_begin /= simulated_time ) REWIND( 19 )
479
480!
481!-- Determine kind of model run
482    IF ( TRIM( initializing_actions ) == 'read_restart_data' )  THEN
483       run_classification = 'restart run'
484    ELSEIF ( TRIM( initializing_actions ) == 'cyclic_fill' )  THEN
485       run_classification = 'run with cyclic fill of 3D - prerun data'
486    ELSEIF ( INDEX( initializing_actions, 'set_constant_profiles' ) /= 0 )  THEN
487       run_classification = 'run without 1D - prerun'
488    ELSEIF ( INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 )  THEN
489       run_classification = 'run with 1D - prerun'
490    ELSEIF ( INDEX( initializing_actions, 'by_user' ) /=0 )  THEN
491       run_classification = 'run initialized by user'
492    ELSE
493       message_string = ' unknown action(s): ' // TRIM( initializing_actions )
494       CALL message( 'header', 'PA0191', 0, 0, 0, 6, 0 )
495    ENDIF
496    IF ( nested_run )  run_classification = 'nested ' // run_classification
497    IF ( ocean )  THEN
498       run_classification = 'ocean - ' // run_classification
499    ELSE
500       run_classification = 'atmosphere - ' // run_classification
501    ENDIF
502
503!
504!-- Run-identification, date, time, host
505    host_chr = host(1:10)
506    ver_rev = TRIM( version ) // '  ' // TRIM( revision )
507    WRITE ( io, 100 )  ver_rev, TRIM( run_classification )
508    IF ( TRIM( coupling_mode ) /= 'uncoupled' )  THEN
509       WRITE ( io, 101 )  coupling_mode
510    ENDIF
511#if defined( __parallel )
512    IF ( coupling_start_time /= 0.0_wp  .AND. .NOT. spinup )  THEN
513       IF ( coupling_start_time > simulated_time_at_begin )  THEN
514          WRITE ( io, 109 )
515       ELSE
516          WRITE ( io, 114 )
517       ENDIF
518    ENDIF
519#endif
520    IF ( ensemble_member_nr /= 0 )  THEN
521       WRITE ( io, 512 )  run_date, run_identifier, run_time, runnr,           &
522                       ADJUSTR( host_chr ), ensemble_member_nr
523    ELSE
524       WRITE ( io, 102 )  run_date, run_identifier, run_time, runnr,           &
525                       ADJUSTR( host_chr )
526    ENDIF
527#if defined( __parallel )
528    IF ( npex == -1  .AND.  npey == -1 )  THEN
529       char1 = 'calculated'
530    ELSE
531       char1 = 'predefined'
532    ENDIF
533    IF ( threads_per_task == 1 )  THEN
534       WRITE ( io, 103 )  numprocs, pdims(1), pdims(2), TRIM( char1 )
535    ELSE
536       WRITE ( io, 104 )  numprocs*threads_per_task, numprocs, &
537                          threads_per_task, pdims(1), pdims(2), TRIM( char1 )
538    ENDIF
539    IF ( ( host(1:3) == 'ibm'  .OR.  host(1:3) == 'nec'  .OR.    &
540           host(1:2) == 'lc'   .OR.  host(1:3) == 'dec' )  .AND. &
541         npex == -1  .AND.  pdims(2) == 1 )                      &
542    THEN
543       WRITE ( io, 106 )
544    ELSEIF ( pdims(2) == 1 )  THEN
545       WRITE ( io, 107 )  'x'
546    ELSEIF ( pdims(1) == 1 )  THEN
547       WRITE ( io, 107 )  'y'
548    ENDIF
549    IF ( use_seperate_pe_for_dvrp_output )  WRITE ( io, 105 )
550    IF ( numprocs /= maximum_parallel_io_streams )  THEN
551       WRITE ( io, 108 )  maximum_parallel_io_streams
552    ENDIF
553#endif
554
555!
556!-- Nesting informations
557    IF ( nested_run )  THEN
558
559       WRITE ( io, 600 )  TRIM( nesting_mode ),                                &
560                          TRIM( nesting_datatransfer_mode )
561       CALL pmc_get_model_info( ncpl = ncpl, cpl_id = my_cpl_id )
562
563       DO  n = 1, ncpl
564          CALL pmc_get_model_info( request_for_cpl_id = n, cpl_name = cpl_name,&
565                                   cpl_parent_id = cpl_parent_id,              &
566                                   lower_left_x = lower_left_coord_x,          &
567                                   lower_left_y = lower_left_coord_y,          &
568                                   npe_total = npe_total )
569          IF ( n == my_cpl_id )  THEN
570             char1 = '*'
571          ELSE
572             char1 = ' '
573          ENDIF
574          WRITE ( io, 601 )  TRIM( char1 ), n, cpl_parent_id, npe_total,       &
575                             lower_left_coord_x, lower_left_coord_y,           &
576                             TRIM( cpl_name )
577       ENDDO
578    ENDIF
579    WRITE ( io, 99 )
580
581!
582!-- Numerical schemes
583    WRITE ( io, 110 )
584    WRITE ( io, 121 )  TRIM( approximation )
585    IF ( psolver(1:7) == 'poisfft' )  THEN
586       WRITE ( io, 111 )  TRIM( fft_method )
587       IF ( transpose_compute_overlap )  WRITE( io, 115 )
588    ELSEIF ( psolver == 'sor' )  THEN
589       WRITE ( io, 112 )  nsor_ini, nsor, omega_sor
590    ELSEIF ( psolver(1:9) == 'multigrid' )  THEN
591       WRITE ( io, 135 )  TRIM(psolver), cycle_mg, maximum_grid_level, ngsrb
592       IF ( mg_cycles == -1 )  THEN
593          WRITE ( io, 140 )  residual_limit
594       ELSE
595          WRITE ( io, 141 )  mg_cycles
596       ENDIF
597       IF ( mg_switch_to_pe0_level == 0 )  THEN
598          WRITE ( io, 136 )  nxr_mg(1)-nxl_mg(1)+1, nyn_mg(1)-nys_mg(1)+1, &
599                             nzt_mg(1)
600       ELSEIF (  mg_switch_to_pe0_level /= -1 )  THEN
601          WRITE ( io, 137 )  mg_switch_to_pe0_level,            &
602                             mg_loc_ind(2,0)-mg_loc_ind(1,0)+1, &
603                             mg_loc_ind(4,0)-mg_loc_ind(3,0)+1, &
604                             nzt_mg(mg_switch_to_pe0_level),    &
605                             nxr_mg(1)-nxl_mg(1)+1, nyn_mg(1)-nys_mg(1)+1, &
606                             nzt_mg(1)
607       ENDIF
608       IF ( psolver == 'multigrid_noopt' .AND. masking_method )  WRITE ( io, 144 )
609    ENDIF
610    IF ( call_psolver_at_all_substeps  .AND. timestep_scheme(1:5) == 'runge' ) &
611    THEN
612       WRITE ( io, 142 )
613    ENDIF
614
615    IF ( momentum_advec == 'pw-scheme' )  THEN
616       WRITE ( io, 113 )
617    ELSEIF (momentum_advec == 'ws-scheme' )  THEN
618       WRITE ( io, 503 )
619    ENDIF
620    IF ( scalar_advec == 'pw-scheme' )  THEN
621       WRITE ( io, 116 )
622    ELSEIF ( scalar_advec == 'ws-scheme' )  THEN
623       WRITE ( io, 504 )
624    ELSE
625       WRITE ( io, 118 )
626    ENDIF
627
628    WRITE ( io, 139 )  TRIM( loop_optimization )
629
630    IF ( galilei_transformation )  THEN
631       IF ( use_ug_for_galilei_tr )  THEN
632          char1 = '0.6 * geostrophic wind'
633       ELSE
634          char1 = 'mean wind in model domain'
635       ENDIF
636       IF ( simulated_time_at_begin == simulated_time )  THEN
637          char2 = 'at the start of the run'
638       ELSE
639          char2 = 'at the end of the run'
640       ENDIF
641       WRITE ( io, 119 )  TRIM( char1 ), TRIM( char2 ),                        &
642                          advected_distance_x/1000.0_wp,                       &
643                          advected_distance_y/1000.0_wp
644    ENDIF
645    WRITE ( io, 122 )  timestep_scheme
646    IF ( use_upstream_for_tke )  WRITE ( io, 143 )
647    IF ( rayleigh_damping_factor /= 0.0_wp )  THEN
648       IF ( .NOT. ocean )  THEN
649          WRITE ( io, 123 )  'above', rayleigh_damping_height, &
650               rayleigh_damping_factor
651       ELSE
652          WRITE ( io, 123 )  'below', rayleigh_damping_height, &
653               rayleigh_damping_factor
654       ENDIF
655    ENDIF
656    IF ( neutral )  WRITE ( io, 131 )  pt_surface
657    IF ( humidity )  THEN
658       IF ( .NOT. cloud_physics )  THEN
659          WRITE ( io, 129 )
660       ELSE
661          WRITE ( io, 130 )
662       ENDIF
663    ENDIF
664    IF ( passive_scalar )  WRITE ( io, 134 )
665    IF ( conserve_volume_flow )  THEN
666       WRITE ( io, 150 )  conserve_volume_flow_mode
667       IF ( TRIM( conserve_volume_flow_mode ) == 'bulk_velocity' )  THEN
668          WRITE ( io, 151 )  u_bulk, v_bulk
669       ENDIF
670    ELSEIF ( dp_external )  THEN
671       IF ( dp_smooth )  THEN
672          WRITE ( io, 152 )  dpdxy, dp_level_b, ', vertically smoothed.'
673       ELSE
674          WRITE ( io, 152 )  dpdxy, dp_level_b, '.'
675       ENDIF
676    ENDIF
677    WRITE ( io, 99 )
678
679!
680!-- Runtime and timestep information
681    WRITE ( io, 200 )
682    IF ( .NOT. dt_fixed )  THEN
683       WRITE ( io, 201 )  dt_max, cfl_factor
684    ELSE
685       WRITE ( io, 202 )  dt
686    ENDIF
687    WRITE ( io, 203 )  simulated_time_at_begin, end_time
688
689    IF ( time_restart /= 9999999.9_wp  .AND. &
690         simulated_time_at_begin == simulated_time )  THEN
691       IF ( dt_restart == 9999999.9_wp )  THEN
692          WRITE ( io, 204 )  ' Restart at:       ',time_restart
693       ELSE
694          WRITE ( io, 205 )  ' Restart at:       ',time_restart, dt_restart
695       ENDIF
696    ENDIF
697
698    IF ( simulated_time_at_begin /= simulated_time )  THEN
699       i = MAX ( log_point_s(10)%counts, 1 )
700       IF ( ( simulated_time - simulated_time_at_begin ) == 0.0_wp )  THEN
701          cpuseconds_per_simulated_second = 0.0_wp
702       ELSE
703          cpuseconds_per_simulated_second = log_point_s(10)%sum / &
704                                            ( simulated_time -    &
705                                              simulated_time_at_begin )
706       ENDIF
707       WRITE ( io, 206 )  simulated_time, log_point_s(10)%sum,      &
708                          log_point_s(10)%sum / REAL( i, KIND=wp ), &
709                          cpuseconds_per_simulated_second
710       IF ( time_restart /= 9999999.9_wp  .AND.  time_restart < end_time )  THEN
711          IF ( dt_restart == 9999999.9_wp )  THEN
712             WRITE ( io, 204 )  ' Next restart at:     ',time_restart
713          ELSE
714             WRITE ( io, 205 )  ' Next restart at:     ',time_restart, dt_restart
715          ENDIF
716       ENDIF
717    ENDIF
718
719
720!
721!-- Start time for coupled runs, if independent precursor runs for atmosphere
722!-- and ocean are used or have been used. In this case, coupling_start_time
723!-- defines the time when the coupling is switched on.
724    IF ( coupling_start_time /= 0.0_wp )  THEN
725       WRITE ( io, 207 )  coupling_start_time
726    ENDIF
727
728!
729!-- Computational grid
730    IF ( .NOT. ocean )  THEN
731       WRITE ( io, 250 )  dx, dy, dz, (nx+1)*dx, (ny+1)*dy, zu(nzt+1)
732       IF ( dz_stretch_level_index < nzt+1 )  THEN
733          WRITE ( io, 252 )  dz_stretch_level, dz_stretch_level_index, &
734                             dz_stretch_factor, dz_max
735       ENDIF
736    ELSE
737       WRITE ( io, 250 )  dx, dy, dz, (nx+1)*dx, (ny+1)*dy, zu(0)
738       IF ( dz_stretch_level_index > 0 )  THEN
739          WRITE ( io, 252 )  dz_stretch_level, dz_stretch_level_index, &
740                             dz_stretch_factor, dz_max
741       ENDIF
742    ENDIF
743    WRITE ( io, 254 )  nx, ny, nzt+1, MIN( nnx, nx+1 ), MIN( nny, ny+1 ), &
744                       MIN( nnz+2, nzt+2 )
745    IF ( sloping_surface )  WRITE ( io, 260 )  alpha_surface
746
747!
748!-- Large scale forcing and nudging
749    WRITE ( io, 160 )
750    IF ( large_scale_forcing )  THEN
751       WRITE ( io, 162 )
752       WRITE ( io, 163 )
753
754       IF ( large_scale_subsidence )  THEN
755          IF ( .NOT. use_subsidence_tendencies )  THEN
756             WRITE ( io, 164 )
757          ELSE
758             WRITE ( io, 165 )
759          ENDIF
760       ENDIF
761
762       IF ( bc_pt_b == 'dirichlet' )  THEN
763          WRITE ( io, 180 )
764       ELSEIF ( bc_pt_b == 'neumann' )  THEN
765          WRITE ( io, 181 )
766       ENDIF
767
768       IF ( bc_q_b == 'dirichlet' )  THEN
769          WRITE ( io, 182 )
770       ELSEIF ( bc_q_b == 'neumann' )  THEN
771          WRITE ( io, 183 )
772       ENDIF
773
774       WRITE ( io, 167 )
775       IF ( nudging )  THEN
776          WRITE ( io, 170 )
777       ENDIF
778    ELSE
779       WRITE ( io, 161 )
780       WRITE ( io, 171 )
781    ENDIF
782    IF ( large_scale_subsidence )  THEN
783       WRITE ( io, 168 )
784       WRITE ( io, 169 )
785    ENDIF
786
787!
788!-- Profile for the large scale vertial velocity
789!-- Building output strings, starting with surface value
790    IF ( large_scale_subsidence )  THEN
791       temperatures = '   0.0'
792       gradients = '------'
793       slices = '     0'
794       coordinates = '   0.0'
795       i = 1
796       DO  WHILE ( subs_vertical_gradient_level_i(i) /= -9999 )
797
798          WRITE (coor_chr,'(E10.2,7X)')  &
799                                w_subs(subs_vertical_gradient_level_i(i))
800          temperatures = TRIM( temperatures ) // ' ' // TRIM( coor_chr )
801
802          WRITE (coor_chr,'(E10.2,7X)')  subs_vertical_gradient(i)
803          gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
804
805          WRITE (coor_chr,'(I10,7X)')  subs_vertical_gradient_level_i(i)
806          slices = TRIM( slices ) // ' ' // TRIM( coor_chr )
807
808          WRITE (coor_chr,'(F10.2,7X)')  subs_vertical_gradient_level(i)
809          coordinates = TRIM( coordinates ) // ' '  // TRIM( coor_chr )
810
811          IF ( i == 10 )  THEN
812             EXIT
813          ELSE
814             i = i + 1
815          ENDIF
816
817       ENDDO
818
819 
820       IF ( .NOT. large_scale_forcing )  THEN
821          WRITE ( io, 426 )  TRIM( coordinates ), TRIM( temperatures ), &
822                             TRIM( gradients ), TRIM( slices )
823       ENDIF
824
825
826    ENDIF
827
828!-- Profile of the geostrophic wind (component ug)
829!-- Building output strings
830    WRITE ( ugcomponent, '(F6.2)' )  ug_surface
831    gradients = '------'
832    slices = '     0'
833    coordinates = '   0.0'
834    i = 1
835    DO  WHILE ( ug_vertical_gradient_level_ind(i) /= -9999 )
836     
837       WRITE (coor_chr,'(F6.2,1X)')  ug(ug_vertical_gradient_level_ind(i))
838       ugcomponent = TRIM( ugcomponent ) // '  ' // TRIM( coor_chr )
839
840       WRITE (coor_chr,'(F6.2,1X)')  ug_vertical_gradient(i)
841       gradients = TRIM( gradients ) // '  ' // TRIM( coor_chr )
842
843       WRITE (coor_chr,'(I6,1X)')  ug_vertical_gradient_level_ind(i)
844       slices = TRIM( slices ) // '  ' // TRIM( coor_chr )
845
846       WRITE (coor_chr,'(F6.1,1X)')  ug_vertical_gradient_level(i)
847       coordinates = TRIM( coordinates ) // '  ' // TRIM( coor_chr )
848
849       IF ( i == 10 )  THEN
850          EXIT
851       ELSE
852          i = i + 1
853       ENDIF
854
855    ENDDO
856
857    IF ( .NOT. large_scale_forcing )  THEN
858       WRITE ( io, 423 )  TRIM( coordinates ), TRIM( ugcomponent ), &
859                          TRIM( gradients ), TRIM( slices )
860    ENDIF
861
862!-- Profile of the geostrophic wind (component vg)
863!-- Building output strings
864    WRITE ( vgcomponent, '(F6.2)' )  vg_surface
865    gradients = '------'
866    slices = '     0'
867    coordinates = '   0.0'
868    i = 1
869    DO  WHILE ( vg_vertical_gradient_level_ind(i) /= -9999 )
870
871       WRITE (coor_chr,'(F6.2,1X)')  vg(vg_vertical_gradient_level_ind(i))
872       vgcomponent = TRIM( vgcomponent ) // '  ' // TRIM( coor_chr )
873
874       WRITE (coor_chr,'(F6.2,1X)')  vg_vertical_gradient(i)
875       gradients = TRIM( gradients ) // '  ' // TRIM( coor_chr )
876
877       WRITE (coor_chr,'(I6,1X)')  vg_vertical_gradient_level_ind(i)
878       slices = TRIM( slices ) // '  ' // TRIM( coor_chr )
879
880       WRITE (coor_chr,'(F6.1,1X)')  vg_vertical_gradient_level(i)
881       coordinates = TRIM( coordinates ) // '  ' // TRIM( coor_chr )
882
883       IF ( i == 10 )  THEN
884          EXIT
885       ELSE
886          i = i + 1
887       ENDIF
888 
889    ENDDO
890
891    IF ( .NOT. large_scale_forcing )  THEN
892       WRITE ( io, 424 )  TRIM( coordinates ), TRIM( vgcomponent ), &
893                          TRIM( gradients ), TRIM( slices )
894    ENDIF
895
896!
897!-- Topography
898    WRITE ( io, 270 )  topography
899    SELECT CASE ( TRIM( topography ) )
900
901       CASE ( 'flat' )
902          ! no actions necessary
903
904       CASE ( 'single_building' )
905          blx = INT( building_length_x / dx )
906          bly = INT( building_length_y / dy )
907          bh  = MINLOC( ABS( zw - building_height ), 1 ) - 1
908          IF ( ABS( zw(bh  ) - building_height ) == &
909               ABS( zw(bh+1) - building_height )    )  bh = bh + 1
910
911          IF ( building_wall_left == 9999999.9_wp )  THEN
912             building_wall_left = ( nx + 1 - blx ) / 2 * dx
913          ENDIF
914          bxl = INT ( building_wall_left / dx + 0.5_wp )
915          bxr = bxl + blx
916
917          IF ( building_wall_south == 9999999.9_wp )  THEN
918             building_wall_south = ( ny + 1 - bly ) / 2 * dy
919          ENDIF
920          bys = INT ( building_wall_south / dy + 0.5_wp )
921          byn = bys + bly
922
923          WRITE ( io, 271 )  building_length_x, building_length_y, &
924                             building_height, bxl, bxr, bys, byn
925
926       CASE ( 'single_street_canyon' )
927          ch  = MINLOC( ABS( zw - canyon_height ), 1 ) - 1
928          IF ( ABS( zw(ch  ) - canyon_height ) == &
929               ABS( zw(ch+1) - canyon_height )    )  ch = ch + 1
930          IF ( canyon_width_x /= 9999999.9_wp )  THEN
931!
932!--          Street canyon in y direction
933             cwx = NINT( canyon_width_x / dx )
934             IF ( canyon_wall_left == 9999999.9_wp )  THEN
935                canyon_wall_left = ( nx + 1 - cwx ) / 2 * dx
936             ENDIF
937             cxl = NINT( canyon_wall_left / dx )
938             cxr = cxl + cwx
939             WRITE ( io, 272 )  'y', canyon_height, ch, 'u', cxl, cxr
940
941          ELSEIF ( canyon_width_y /= 9999999.9_wp )  THEN
942!
943!--          Street canyon in x direction
944             cwy = NINT( canyon_width_y / dy )
945             IF ( canyon_wall_south == 9999999.9_wp )  THEN
946                canyon_wall_south = ( ny + 1 - cwy ) / 2 * dy
947             ENDIF
948             cys = NINT( canyon_wall_south / dy )
949             cyn = cys + cwy
950             WRITE ( io, 272 )  'x', canyon_height, ch, 'v', cys, cyn
951          ENDIF
952
953       CASE ( 'tunnel' )
954          IF ( tunnel_width_x /= 9999999.9_wp )  THEN
955!
956!--          Tunnel axis in y direction
957             IF ( tunnel_length == 9999999.9_wp  .OR.                          &
958                  tunnel_length >= ( nx + 1 ) * dx )  THEN
959                WRITE ( io, 273 )  'y', tunnel_height, tunnel_wall_depth,      &
960                                        tunnel_width_x
961             ELSE
962                WRITE ( io, 274 )  'y', tunnel_height, tunnel_wall_depth,      &
963                                        tunnel_width_x, tunnel_length
964             ENDIF
965
966          ELSEIF ( tunnel_width_y /= 9999999.9_wp )  THEN
967!
968!--          Tunnel axis in x direction
969             IF ( tunnel_length == 9999999.9_wp  .OR.                          &
970                  tunnel_length >= ( ny + 1 ) * dy )  THEN
971                WRITE ( io, 273 )  'x', tunnel_height, tunnel_wall_depth,      &
972                                        tunnel_width_y
973             ELSE
974                WRITE ( io, 274 )  'x', tunnel_height, tunnel_wall_depth,      &
975                                        tunnel_width_y, tunnel_length
976             ENDIF
977          ENDIF
978
979    END SELECT
980
981    IF ( TRIM( topography ) /= 'flat' )  THEN
982       IF ( TRIM( topography_grid_convention ) == ' ' )  THEN
983          IF ( TRIM( topography ) == 'single_building' .OR.  &
984               TRIM( topography ) == 'single_street_canyon' )  THEN
985             WRITE ( io, 278 )
986          ELSEIF ( TRIM( topography ) == 'read_from_file' )  THEN
987             WRITE ( io, 279 )
988          ENDIF
989       ELSEIF ( TRIM( topography_grid_convention ) == 'cell_edge' )  THEN
990          WRITE ( io, 278 )
991       ELSEIF ( TRIM( topography_grid_convention ) == 'cell_center' )  THEN
992          WRITE ( io, 279 )
993       ENDIF
994    ENDIF
995
996    IF ( synthetic_turbulence_generator )  CALL stg_header ( io )
997
998    IF ( plant_canopy )  CALL pcm_header ( io )
999
1000    IF ( land_surface )  CALL lsm_header ( io )
1001
1002    IF ( radiation )  CALL radiation_header ( io )
1003
1004!
1005!-- Boundary conditions
1006    IF ( ibc_p_b == 0 )  THEN
1007       r_lower = 'p(0)     = 0      |'
1008    ELSEIF ( ibc_p_b == 1 )  THEN
1009       r_lower = 'p(0)     = p(1)   |'
1010    ENDIF
1011    IF ( ibc_p_t == 0 )  THEN
1012       r_upper  = 'p(nzt+1) = 0      |'
1013    ELSE
1014       r_upper  = 'p(nzt+1) = p(nzt) |'
1015    ENDIF
1016
1017    IF ( ibc_uv_b == 0 )  THEN
1018       r_lower = TRIM( r_lower ) // ' uv(0)     = -uv(1)                |'
1019    ELSE
1020       r_lower = TRIM( r_lower ) // ' uv(0)     = uv(1)                 |'
1021    ENDIF
1022    IF ( TRIM( bc_uv_t ) == 'dirichlet_0' )  THEN
1023       r_upper  = TRIM( r_upper  ) // ' uv(nzt+1) = 0                     |'
1024    ELSEIF ( ibc_uv_t == 0 )  THEN
1025       r_upper  = TRIM( r_upper  ) // ' uv(nzt+1) = ug(nzt+1), vg(nzt+1)  |'
1026    ELSE
1027       r_upper  = TRIM( r_upper  ) // ' uv(nzt+1) = uv(nzt)               |'
1028    ENDIF
1029
1030    IF ( ibc_pt_b == 0 )  THEN
1031       IF ( land_surface )  THEN
1032          r_lower = TRIM( r_lower ) // ' pt(0)     = from soil model'
1033       ELSE
1034          r_lower = TRIM( r_lower ) // ' pt(0)     = pt_surface'
1035       ENDIF
1036    ELSEIF ( ibc_pt_b == 1 )  THEN
1037       r_lower = TRIM( r_lower ) // ' pt(0)     = pt(1)'
1038    ELSEIF ( ibc_pt_b == 2 )  THEN
1039       r_lower = TRIM( r_lower ) // ' pt(0)     = from coupled model'
1040    ENDIF
1041    IF ( ibc_pt_t == 0 )  THEN
1042       r_upper  = TRIM( r_upper  ) // ' pt(nzt+1) = pt_top'
1043    ELSEIF( ibc_pt_t == 1 )  THEN
1044       r_upper  = TRIM( r_upper  ) // ' pt(nzt+1) = pt(nzt)'
1045    ELSEIF( ibc_pt_t == 2 )  THEN
1046       r_upper  = TRIM( r_upper  ) // ' pt(nzt+1) = pt(nzt) + dpt/dz_ini'
1047
1048    ENDIF
1049
1050    WRITE ( io, 300 )  r_lower, r_upper
1051
1052    IF ( .NOT. constant_diffusion )  THEN
1053       IF ( ibc_e_b == 1 )  THEN
1054          r_lower = 'e(0)     = e(1)'
1055       ELSE
1056          r_lower = 'e(0)     = e(1) = (u*/0.1)**2'
1057       ENDIF
1058       r_upper = 'e(nzt+1) = e(nzt) = e(nzt-1)'
1059
1060       WRITE ( io, 301 )  'e', r_lower, r_upper       
1061
1062    ENDIF
1063
1064    IF ( ocean )  THEN
1065       r_lower = 'sa(0)    = sa(1)'
1066       IF ( ibc_sa_t == 0 )  THEN
1067          r_upper =  'sa(nzt+1) = sa_surface'
1068       ELSE
1069          r_upper =  'sa(nzt+1) = sa(nzt)'
1070       ENDIF
1071       WRITE ( io, 301 ) 'sa', r_lower, r_upper
1072    ENDIF
1073
1074    IF ( humidity )  THEN
1075       IF ( ibc_q_b == 0 )  THEN
1076          IF ( land_surface )  THEN
1077             r_lower = 'q(0)     = from soil model'
1078          ELSE
1079             r_lower = 'q(0)     = q_surface'
1080          ENDIF
1081
1082       ELSE
1083          r_lower = 'q(0)      = q(1)'
1084       ENDIF
1085       IF ( ibc_q_t == 0 )  THEN
1086          r_upper =  'q(nzt+1) = q_top'
1087       ELSE
1088          r_upper =  'q(nzt+1) = q(nzt) + dq/dz'
1089       ENDIF
1090       WRITE ( io, 301 ) 'q', r_lower, r_upper
1091    ENDIF
1092
1093    IF ( passive_scalar )  THEN
1094       IF ( ibc_s_b == 0 )  THEN
1095          r_lower = 's(0)      = s_surface'
1096       ELSE
1097          r_lower = 's(0)      = s(1)'
1098       ENDIF
1099       IF ( ibc_s_t == 0 )  THEN
1100          r_upper =  's(nzt+1) = s_top'
1101       ELSEIF ( ibc_s_t == 1 )  THEN
1102          r_upper =  's(nzt+1) = s(nzt)'
1103       ELSEIF ( ibc_s_t == 2 )  THEN
1104          r_upper =  's(nzt+1) = s(nzt) + ds/dz'
1105       ENDIF
1106       WRITE ( io, 301 ) 's', r_lower, r_upper
1107    ENDIF
1108
1109    IF ( use_surface_fluxes )  THEN
1110       WRITE ( io, 303 )
1111       IF ( constant_heatflux )  THEN
1112          IF ( large_scale_forcing .AND. lsf_surf )  THEN
1113             IF ( surf_def_h(0)%ns >= 1 )  WRITE ( io, 306 )  surf_def_h(0)%shf(1)
1114          ELSE
1115             WRITE ( io, 306 )  surface_heatflux
1116          ENDIF
1117          IF ( random_heatflux )  WRITE ( io, 307 )
1118       ENDIF
1119       IF ( humidity  .AND.  constant_waterflux )  THEN
1120          IF ( large_scale_forcing .AND. lsf_surf )  THEN
1121             WRITE ( io, 311 ) surf_def_h(0)%qsws(1)
1122          ELSE
1123             WRITE ( io, 311 ) surface_waterflux
1124          ENDIF
1125       ENDIF
1126       IF ( passive_scalar  .AND.  constant_scalarflux )  THEN
1127          WRITE ( io, 313 ) surface_scalarflux
1128       ENDIF
1129    ENDIF
1130
1131    IF ( use_top_fluxes )  THEN
1132       WRITE ( io, 304 )
1133       IF ( coupling_mode == 'uncoupled' )  THEN
1134          WRITE ( io, 320 )  top_momentumflux_u, top_momentumflux_v
1135          IF ( constant_top_heatflux )  THEN
1136             WRITE ( io, 306 )  top_heatflux
1137          ENDIF
1138       ELSEIF ( coupling_mode == 'ocean_to_atmosphere' )  THEN
1139          WRITE ( io, 316 )
1140       ENDIF
1141       IF ( ocean  .AND.  constant_top_salinityflux )                          &
1142          WRITE ( io, 309 )  top_salinityflux
1143       IF ( humidity       )  WRITE ( io, 315 )
1144       IF ( passive_scalar .AND.  constant_top_scalarflux )                    &
1145          WRITE ( io, 302 ) top_scalarflux
1146    ENDIF
1147
1148    IF ( constant_flux_layer )  THEN
1149       WRITE ( io, 305 )  (zu(1)-zu(0)), roughness_length,                     &
1150                          z0h_factor*roughness_length, kappa,                  &
1151                          zeta_min, zeta_max
1152       IF ( .NOT. constant_heatflux )  WRITE ( io, 308 )
1153       IF ( humidity  .AND.  .NOT. constant_waterflux )  THEN
1154          WRITE ( io, 312 )
1155       ENDIF
1156       IF ( passive_scalar  .AND.  .NOT. constant_scalarflux )  THEN
1157          WRITE ( io, 314 )
1158       ENDIF
1159    ELSE
1160       IF ( INDEX(initializing_actions, 'set_1d-model_profiles') /= 0 )  THEN
1161          WRITE ( io, 310 )  zeta_min, zeta_max
1162       ENDIF
1163    ENDIF
1164
1165    WRITE ( io, 317 )  bc_lr, bc_ns
1166    IF ( .NOT. bc_lr_cyc  .OR.  .NOT. bc_ns_cyc )  THEN
1167       WRITE ( io, 318 )  use_cmax, pt_damping_width, pt_damping_factor       
1168       IF ( turbulent_inflow )  THEN
1169          IF ( .NOT. recycling_yshift ) THEN
1170             WRITE ( io, 319 )  recycling_width, recycling_plane, &
1171                                inflow_damping_height, inflow_damping_width
1172          ELSE
1173             WRITE ( io, 322 )  recycling_width, recycling_plane, &
1174                                inflow_damping_height, inflow_damping_width
1175          END IF
1176       ENDIF
1177       IF ( turbulent_outflow )  THEN
1178          WRITE ( io, 323 )  outflow_source_plane, INT(outflow_source_plane/dx)
1179       ENDIF
1180    ENDIF
1181
1182!
1183!-- Initial Profiles
1184    WRITE ( io, 321 )
1185!
1186!-- Initial wind profiles
1187    IF ( u_profile(1) /= 9999999.9_wp )  WRITE ( io, 427 )
1188
1189!
1190!-- Initial temperature profile
1191!-- Building output strings, starting with surface temperature
1192    WRITE ( temperatures, '(F6.2)' )  pt_surface
1193    gradients = '------'
1194    slices = '     0'
1195    coordinates = '   0.0'
1196    i = 1
1197    DO  WHILE ( pt_vertical_gradient_level_ind(i) /= -9999 )
1198
1199       WRITE (coor_chr,'(F7.2)')  pt_init(pt_vertical_gradient_level_ind(i))
1200       temperatures = TRIM( temperatures ) // ' ' // TRIM( coor_chr )
1201
1202       WRITE (coor_chr,'(F7.2)')  pt_vertical_gradient(i)
1203       gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
1204
1205       WRITE (coor_chr,'(I7)')  pt_vertical_gradient_level_ind(i)
1206       slices = TRIM( slices ) // ' ' // TRIM( coor_chr )
1207
1208       WRITE (coor_chr,'(F7.1)')  pt_vertical_gradient_level(i)
1209       coordinates = TRIM( coordinates ) // ' '  // TRIM( coor_chr )
1210
1211       IF ( i == 10 )  THEN
1212          EXIT
1213       ELSE
1214          i = i + 1
1215       ENDIF
1216
1217    ENDDO
1218
1219    IF ( .NOT. nudging )  THEN
1220       WRITE ( io, 420 )  TRIM( coordinates ), TRIM( temperatures ), &
1221                          TRIM( gradients ), TRIM( slices )
1222    ELSE
1223       WRITE ( io, 428 ) 
1224    ENDIF
1225
1226!
1227!-- Initial humidity profile
1228!-- Building output strings, starting with surface humidity
1229    IF ( humidity )  THEN
1230       WRITE ( temperatures, '(E8.1)' )  q_surface
1231       gradients = '--------'
1232       slices = '       0'
1233       coordinates = '     0.0'
1234       i = 1
1235       DO  WHILE ( q_vertical_gradient_level_ind(i) /= -9999 )
1236         
1237          WRITE (coor_chr,'(E8.1,4X)')  q_init(q_vertical_gradient_level_ind(i))
1238          temperatures = TRIM( temperatures ) // '  ' // TRIM( coor_chr )
1239
1240          WRITE (coor_chr,'(E8.1,4X)')  q_vertical_gradient(i)
1241          gradients = TRIM( gradients ) // '  ' // TRIM( coor_chr )
1242         
1243          WRITE (coor_chr,'(I8,4X)')  q_vertical_gradient_level_ind(i)
1244          slices = TRIM( slices ) // '  ' // TRIM( coor_chr )
1245         
1246          WRITE (coor_chr,'(F8.1,4X)')  q_vertical_gradient_level(i)
1247          coordinates = TRIM( coordinates ) // '  '  // TRIM( coor_chr )
1248
1249          IF ( i == 10 )  THEN
1250             EXIT
1251          ELSE
1252             i = i + 1
1253          ENDIF
1254
1255       ENDDO
1256
1257       IF ( .NOT. nudging )  THEN
1258          WRITE ( io, 421 )  TRIM( coordinates ), TRIM( temperatures ),        &
1259                             TRIM( gradients ), TRIM( slices )
1260       ENDIF
1261    ENDIF
1262!
1263!-- Initial scalar profile
1264!-- Building output strings, starting with surface humidity
1265    IF ( passive_scalar )  THEN
1266       WRITE ( temperatures, '(E8.1)' )  s_surface
1267       gradients = '--------'
1268       slices = '       0'
1269       coordinates = '     0.0'
1270       i = 1
1271       DO  WHILE ( s_vertical_gradient_level_ind(i) /= -9999 )
1272         
1273          WRITE (coor_chr,'(E8.1,4X)')  s_init(s_vertical_gradient_level_ind(i))
1274          temperatures = TRIM( temperatures ) // '  ' // TRIM( coor_chr )
1275
1276          WRITE (coor_chr,'(E8.1,4X)')  s_vertical_gradient(i)
1277          gradients = TRIM( gradients ) // '  ' // TRIM( coor_chr )
1278         
1279          WRITE (coor_chr,'(I8,4X)')  s_vertical_gradient_level_ind(i)
1280          slices = TRIM( slices ) // '  ' // TRIM( coor_chr )
1281         
1282          WRITE (coor_chr,'(F8.1,4X)')  s_vertical_gradient_level(i)
1283          coordinates = TRIM( coordinates ) // '  '  // TRIM( coor_chr )
1284
1285          IF ( i == 10 )  THEN
1286             EXIT
1287          ELSE
1288             i = i + 1
1289          ENDIF
1290
1291       ENDDO
1292
1293       WRITE ( io, 422 )  TRIM( coordinates ), TRIM( temperatures ),           &
1294                          TRIM( gradients ), TRIM( slices )
1295    ENDIF   
1296
1297!
1298!-- Initial salinity profile
1299!-- Building output strings, starting with surface salinity
1300    IF ( ocean )  THEN
1301       WRITE ( temperatures, '(F6.2)' )  sa_surface
1302       gradients = '------'
1303       slices = '     0'
1304       coordinates = '   0.0'
1305       i = 1
1306       DO  WHILE ( sa_vertical_gradient_level_ind(i) /= -9999 )
1307
1308          WRITE (coor_chr,'(F7.2)')  sa_init(sa_vertical_gradient_level_ind(i))
1309          temperatures = TRIM( temperatures ) // ' ' // TRIM( coor_chr )
1310
1311          WRITE (coor_chr,'(F7.2)')  sa_vertical_gradient(i)
1312          gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
1313
1314          WRITE (coor_chr,'(I7)')  sa_vertical_gradient_level_ind(i)
1315          slices = TRIM( slices ) // ' ' // TRIM( coor_chr )
1316
1317          WRITE (coor_chr,'(F7.1)')  sa_vertical_gradient_level(i)
1318          coordinates = TRIM( coordinates ) // ' '  // TRIM( coor_chr )
1319
1320          IF ( i == 10 )  THEN
1321             EXIT
1322          ELSE
1323             i = i + 1
1324          ENDIF
1325
1326       ENDDO
1327
1328       WRITE ( io, 425 )  TRIM( coordinates ), TRIM( temperatures ), &
1329                          TRIM( gradients ), TRIM( slices )
1330    ENDIF
1331
1332
1333!
1334!-- Listing of 1D-profiles
1335    WRITE ( io, 325 )  dt_dopr_listing
1336    IF ( averaging_interval_pr /= 0.0_wp )  THEN
1337       WRITE ( io, 326 )  averaging_interval_pr, dt_averaging_input_pr
1338    ENDIF
1339
1340!
1341!-- DATA output
1342    WRITE ( io, 330 )
1343    IF ( averaging_interval_pr /= 0.0_wp )  THEN
1344       WRITE ( io, 326 )  averaging_interval_pr, dt_averaging_input_pr
1345    ENDIF
1346
1347!
1348!-- 1D-profiles
1349    dopr_chr = 'Profile:'
1350    IF ( dopr_n /= 0 )  THEN
1351       WRITE ( io, 331 )
1352
1353       output_format = ''
1354       output_format = netcdf_data_format_string
1355       IF ( netcdf_deflate == 0 )  THEN
1356          WRITE ( io, 344 )  output_format
1357       ELSE
1358          WRITE ( io, 354 )  TRIM( output_format ), netcdf_deflate
1359       ENDIF
1360
1361       DO  i = 1, dopr_n
1362          dopr_chr = TRIM( dopr_chr ) // ' ' // TRIM( data_output_pr(i) ) // ','
1363          IF ( LEN_TRIM( dopr_chr ) >= 60 )  THEN
1364             WRITE ( io, 332 )  dopr_chr
1365             dopr_chr = '       :'
1366          ENDIF
1367       ENDDO
1368
1369       IF ( dopr_chr /= '' )  THEN
1370          WRITE ( io, 332 )  dopr_chr
1371       ENDIF
1372       WRITE ( io, 333 )  dt_dopr, averaging_interval_pr, dt_averaging_input_pr
1373       IF ( skip_time_dopr /= 0.0_wp )  WRITE ( io, 339 )  skip_time_dopr
1374    ENDIF
1375
1376!
1377!-- 2D-arrays
1378    DO  av = 0, 1
1379
1380       i = 1
1381       do2d_xy = ''
1382       do2d_xz = ''
1383       do2d_yz = ''
1384       DO  WHILE ( do2d(av,i) /= ' ' )
1385
1386          l = MAX( 2, LEN_TRIM( do2d(av,i) ) )
1387          do2d_mode = do2d(av,i)(l-1:l)
1388
1389          SELECT CASE ( do2d_mode )
1390             CASE ( 'xy' )
1391                ll = LEN_TRIM( do2d_xy )
1392                do2d_xy = do2d_xy(1:ll) // ' ' // do2d(av,i)(1:l-3) // ','
1393             CASE ( 'xz' )
1394                ll = LEN_TRIM( do2d_xz )
1395                do2d_xz = do2d_xz(1:ll) // ' ' // do2d(av,i)(1:l-3) // ','
1396             CASE ( 'yz' )
1397                ll = LEN_TRIM( do2d_yz )
1398                do2d_yz = do2d_yz(1:ll) // ' ' // do2d(av,i)(1:l-3) // ','
1399          END SELECT
1400
1401          i = i + 1
1402
1403       ENDDO
1404
1405       IF ( ( ( do2d_xy /= ''  .AND.  section(1,1) /= -9999 )  .OR.    &
1406              ( do2d_xz /= ''  .AND.  section(1,2) /= -9999 )  .OR.    &
1407              ( do2d_yz /= ''  .AND.  section(1,3) /= -9999 ) ) )  THEN
1408
1409          IF (  av == 0 )  THEN
1410             WRITE ( io, 334 )  ''
1411          ELSE
1412             WRITE ( io, 334 )  '(time-averaged)'
1413          ENDIF
1414
1415          IF ( do2d_at_begin )  THEN
1416             begin_chr = 'and at the start'
1417          ELSE
1418             begin_chr = ''
1419          ENDIF
1420
1421          output_format = ''
1422          output_format = netcdf_data_format_string
1423          IF ( netcdf_deflate == 0 )  THEN
1424             WRITE ( io, 344 )  output_format
1425          ELSE
1426             WRITE ( io, 354 )  TRIM( output_format ), netcdf_deflate
1427          ENDIF
1428
1429          IF ( do2d_xy /= ''  .AND.  section(1,1) /= -9999 )  THEN
1430             i = 1
1431             slices = '/'
1432             coordinates = '/'
1433!
1434!--          Building strings with index and coordinate information of the
1435!--          slices
1436             DO  WHILE ( section(i,1) /= -9999 )
1437
1438                WRITE (section_chr,'(I5)')  section(i,1)
1439                section_chr = ADJUSTL( section_chr )
1440                slices = TRIM( slices ) // TRIM( section_chr ) // '/'
1441
1442                IF ( section(i,1) == -1 )  THEN
1443                   WRITE (coor_chr,'(F10.1)')  -1.0_wp
1444                ELSE
1445                   WRITE (coor_chr,'(F10.1)')  zu(section(i,1))
1446                ENDIF
1447                coor_chr = ADJUSTL( coor_chr )
1448                coordinates = TRIM( coordinates ) // TRIM( coor_chr ) // '/'
1449
1450                i = i + 1
1451             ENDDO
1452             IF ( av == 0 )  THEN
1453                WRITE ( io, 335 )  'XY', do2d_xy, dt_do2d_xy, &
1454                                   TRIM( begin_chr ), 'k', TRIM( slices ), &
1455                                   TRIM( coordinates )
1456                IF ( skip_time_do2d_xy /= 0.0_wp )  THEN
1457                   WRITE ( io, 339 )  skip_time_do2d_xy
1458                ENDIF
1459             ELSE
1460                WRITE ( io, 342 )  'XY', do2d_xy, dt_data_output_av, &
1461                                   TRIM( begin_chr ), averaging_interval, &
1462                                   dt_averaging_input, 'k', TRIM( slices ), &
1463                                   TRIM( coordinates )
1464                IF ( skip_time_data_output_av /= 0.0_wp )  THEN
1465                   WRITE ( io, 339 )  skip_time_data_output_av
1466                ENDIF
1467             ENDIF
1468             IF ( netcdf_data_format > 4 )  THEN
1469                WRITE ( io, 352 )  ntdim_2d_xy(av)
1470             ELSE
1471                WRITE ( io, 353 )
1472             ENDIF
1473          ENDIF
1474
1475          IF ( do2d_xz /= ''  .AND.  section(1,2) /= -9999 )  THEN
1476             i = 1
1477             slices = '/'
1478             coordinates = '/'
1479!
1480!--          Building strings with index and coordinate information of the
1481!--          slices
1482             DO  WHILE ( section(i,2) /= -9999 )
1483
1484                WRITE (section_chr,'(I5)')  section(i,2)
1485                section_chr = ADJUSTL( section_chr )
1486                slices = TRIM( slices ) // TRIM( section_chr ) // '/'
1487
1488                WRITE (coor_chr,'(F10.1)')  section(i,2) * dy
1489                coor_chr = ADJUSTL( coor_chr )
1490                coordinates = TRIM( coordinates ) // TRIM( coor_chr ) // '/'
1491
1492                i = i + 1
1493             ENDDO
1494             IF ( av == 0 )  THEN
1495                WRITE ( io, 335 )  'XZ', do2d_xz, dt_do2d_xz, &
1496                                   TRIM( begin_chr ), 'j', TRIM( slices ), &
1497                                   TRIM( coordinates )
1498                IF ( skip_time_do2d_xz /= 0.0_wp )  THEN
1499                   WRITE ( io, 339 )  skip_time_do2d_xz
1500                ENDIF
1501             ELSE
1502                WRITE ( io, 342 )  'XZ', do2d_xz, dt_data_output_av, &
1503                                   TRIM( begin_chr ), averaging_interval, &
1504                                   dt_averaging_input, 'j', TRIM( slices ), &
1505                                   TRIM( coordinates )
1506                IF ( skip_time_data_output_av /= 0.0_wp )  THEN
1507                   WRITE ( io, 339 )  skip_time_data_output_av
1508                ENDIF
1509             ENDIF
1510             IF ( netcdf_data_format > 4 )  THEN
1511                WRITE ( io, 352 )  ntdim_2d_xz(av)
1512             ELSE
1513                WRITE ( io, 353 )
1514             ENDIF
1515          ENDIF
1516
1517          IF ( do2d_yz /= ''  .AND.  section(1,3) /= -9999 )  THEN
1518             i = 1
1519             slices = '/'
1520             coordinates = '/'
1521!
1522!--          Building strings with index and coordinate information of the
1523!--          slices
1524             DO  WHILE ( section(i,3) /= -9999 )
1525
1526                WRITE (section_chr,'(I5)')  section(i,3)
1527                section_chr = ADJUSTL( section_chr )
1528                slices = TRIM( slices ) // TRIM( section_chr ) // '/'
1529
1530                WRITE (coor_chr,'(F10.1)')  section(i,3) * dx
1531                coor_chr = ADJUSTL( coor_chr )
1532                coordinates = TRIM( coordinates ) // TRIM( coor_chr ) // '/'
1533
1534                i = i + 1
1535             ENDDO
1536             IF ( av == 0 )  THEN
1537                WRITE ( io, 335 )  'YZ', do2d_yz, dt_do2d_yz, &
1538                                   TRIM( begin_chr ), 'i', TRIM( slices ), &
1539                                   TRIM( coordinates )
1540                IF ( skip_time_do2d_yz /= 0.0_wp )  THEN
1541                   WRITE ( io, 339 )  skip_time_do2d_yz
1542                ENDIF
1543             ELSE
1544                WRITE ( io, 342 )  'YZ', do2d_yz, dt_data_output_av, &
1545                                   TRIM( begin_chr ), averaging_interval, &
1546                                   dt_averaging_input, 'i', TRIM( slices ), &
1547                                   TRIM( coordinates )
1548                IF ( skip_time_data_output_av /= 0.0_wp )  THEN
1549                   WRITE ( io, 339 )  skip_time_data_output_av
1550                ENDIF
1551             ENDIF
1552             IF ( netcdf_data_format > 4 )  THEN
1553                WRITE ( io, 352 )  ntdim_2d_yz(av)
1554             ELSE
1555                WRITE ( io, 353 )
1556             ENDIF
1557          ENDIF
1558
1559       ENDIF
1560
1561    ENDDO
1562
1563!
1564!-- 3d-arrays
1565    DO  av = 0, 1
1566
1567       i = 1
1568       do3d_chr = ''
1569       DO  WHILE ( do3d(av,i) /= ' ' )
1570
1571          do3d_chr = TRIM( do3d_chr ) // ' ' // TRIM( do3d(av,i) ) // ','
1572          i = i + 1
1573
1574       ENDDO
1575
1576       IF ( do3d_chr /= '' )  THEN
1577          IF ( av == 0 )  THEN
1578             WRITE ( io, 336 )  ''
1579          ELSE
1580             WRITE ( io, 336 )  '(time-averaged)'
1581          ENDIF
1582
1583          output_format = netcdf_data_format_string
1584          IF ( netcdf_deflate == 0 )  THEN
1585             WRITE ( io, 344 )  output_format
1586          ELSE
1587             WRITE ( io, 354 )  TRIM( output_format ), netcdf_deflate
1588          ENDIF
1589
1590          IF ( do3d_at_begin )  THEN
1591             begin_chr = 'and at the start'
1592          ELSE
1593             begin_chr = ''
1594          ENDIF
1595          IF ( av == 0 )  THEN
1596             WRITE ( io, 337 )  do3d_chr, dt_do3d, TRIM( begin_chr ), &
1597                                zu(nz_do3d), nz_do3d
1598          ELSE
1599             WRITE ( io, 343 )  do3d_chr, dt_data_output_av,           &
1600                                TRIM( begin_chr ), averaging_interval, &
1601                                dt_averaging_input, zu(nz_do3d), nz_do3d
1602          ENDIF
1603
1604          IF ( netcdf_data_format > 4 )  THEN
1605             WRITE ( io, 352 )  ntdim_3d(av)
1606          ELSE
1607             WRITE ( io, 353 )
1608          ENDIF
1609
1610          IF ( av == 0 )  THEN
1611             IF ( skip_time_do3d /= 0.0_wp )  THEN
1612                WRITE ( io, 339 )  skip_time_do3d
1613             ENDIF
1614          ELSE
1615             IF ( skip_time_data_output_av /= 0.0_wp )  THEN
1616                WRITE ( io, 339 )  skip_time_data_output_av
1617             ENDIF
1618          ENDIF
1619
1620       ENDIF
1621
1622    ENDDO
1623
1624!
1625!-- masked arrays
1626    IF ( masks > 0 )  WRITE ( io, 345 )  &
1627         mask_scale_x, mask_scale_y, mask_scale_z
1628    DO  mid = 1, masks
1629       DO  av = 0, 1
1630
1631          i = 1
1632          domask_chr = ''
1633          DO  WHILE ( domask(mid,av,i) /= ' ' )
1634             domask_chr = TRIM( domask_chr ) // ' ' //  &
1635                          TRIM( domask(mid,av,i) ) // ','
1636             i = i + 1
1637          ENDDO
1638
1639          IF ( domask_chr /= '' )  THEN
1640             IF ( av == 0 )  THEN
1641                WRITE ( io, 346 )  '', mid
1642             ELSE
1643                WRITE ( io, 346 )  ' (time-averaged)', mid
1644             ENDIF
1645
1646             output_format = netcdf_data_format_string
1647!--          Parallel output not implemented for mask data, hence
1648!--          output_format must be adjusted.
1649             IF ( netcdf_data_format == 5 ) output_format = 'netCDF4/HDF5'
1650             IF ( netcdf_data_format == 6 ) output_format = 'netCDF4/HDF5 classic'
1651             IF ( netcdf_deflate == 0 )  THEN
1652                WRITE ( io, 344 )  output_format
1653             ELSE
1654                WRITE ( io, 354 )  TRIM( output_format ), netcdf_deflate
1655             ENDIF
1656
1657             IF ( av == 0 )  THEN
1658                WRITE ( io, 347 )  domask_chr, dt_domask(mid)
1659             ELSE
1660                WRITE ( io, 348 )  domask_chr, dt_data_output_av, &
1661                                   averaging_interval, dt_averaging_input
1662             ENDIF
1663
1664             IF ( av == 0 )  THEN
1665                IF ( skip_time_domask(mid) /= 0.0_wp )  THEN
1666                   WRITE ( io, 339 )  skip_time_domask(mid)
1667                ENDIF
1668             ELSE
1669                IF ( skip_time_data_output_av /= 0.0_wp )  THEN
1670                   WRITE ( io, 339 )  skip_time_data_output_av
1671                ENDIF
1672             ENDIF
1673!
1674!--          output locations
1675             DO  dim = 1, 3
1676                IF ( mask(mid,dim,1) >= 0.0_wp )  THEN
1677                   count = 0
1678                   DO  WHILE ( mask(mid,dim,count+1) >= 0.0_wp )
1679                      count = count + 1
1680                   ENDDO
1681                   WRITE ( io, 349 )  dir(dim), dir(dim), mid, dir(dim), &
1682                                      mask(mid,dim,:count)
1683                ELSEIF ( mask_loop(mid,dim,1) < 0.0_wp .AND.  &
1684                         mask_loop(mid,dim,2) < 0.0_wp .AND.  &
1685                         mask_loop(mid,dim,3) == 0.0_wp )  THEN
1686                   WRITE ( io, 350 )  dir(dim), dir(dim)
1687                ELSEIF ( mask_loop(mid,dim,3) == 0.0_wp )  THEN
1688                   WRITE ( io, 351 )  dir(dim), dir(dim), mid, dir(dim), &
1689                                      mask_loop(mid,dim,1:2)
1690                ELSE
1691                   WRITE ( io, 351 )  dir(dim), dir(dim), mid, dir(dim), &
1692                                      mask_loop(mid,dim,1:3)
1693                ENDIF
1694             ENDDO
1695          ENDIF
1696
1697       ENDDO
1698    ENDDO
1699
1700!
1701!-- Timeseries
1702    IF ( dt_dots /= 9999999.9_wp )  THEN
1703       WRITE ( io, 340 )
1704
1705       output_format = netcdf_data_format_string
1706       IF ( netcdf_deflate == 0 )  THEN
1707          WRITE ( io, 344 )  output_format
1708       ELSE
1709          WRITE ( io, 354 )  TRIM( output_format ), netcdf_deflate
1710       ENDIF
1711       WRITE ( io, 341 )  dt_dots
1712    ENDIF
1713
1714#if defined( __dvrp_graphics )
1715!
1716!-- Dvrp-output
1717    IF ( dt_dvrp /= 9999999.9_wp )  THEN
1718       WRITE ( io, 360 )  dt_dvrp, TRIM( dvrp_output ), TRIM( dvrp_host ), &
1719                          TRIM( dvrp_username ), TRIM( dvrp_directory )
1720       i = 1
1721       l = 0
1722       m = 0
1723       DO WHILE ( mode_dvrp(i) /= ' ' )
1724          IF ( mode_dvrp(i)(1:10) == 'isosurface' )  THEN
1725             READ ( mode_dvrp(i), '(10X,I2)' )  j
1726             l = l + 1
1727             IF ( do3d(0,j) /= ' ' )  THEN
1728                WRITE ( io, 361 )  TRIM( do3d(0,j) ), threshold(l), &
1729                                   isosurface_color(:,l)
1730             ENDIF
1731          ELSEIF ( mode_dvrp(i)(1:6) == 'slicer' )  THEN
1732             READ ( mode_dvrp(i), '(6X,I2)' )  j
1733             m = m + 1
1734             IF ( do2d(0,j) /= ' ' )  THEN
1735                WRITE ( io, 362 )  TRIM( do2d(0,j) ), &
1736                                   slicer_range_limits_dvrp(:,m)
1737             ENDIF
1738          ENDIF
1739          i = i + 1
1740       ENDDO
1741
1742       WRITE ( io, 365 )  groundplate_color, superelevation_x, &
1743                          superelevation_y, superelevation, clip_dvrp_l, &
1744                          clip_dvrp_r, clip_dvrp_s, clip_dvrp_n
1745
1746       IF ( TRIM( topography ) /= 'flat' )  THEN
1747          WRITE ( io, 366 )  topography_color
1748          IF ( cluster_size > 1 )  THEN
1749             WRITE ( io, 367 )  cluster_size
1750          ENDIF
1751       ENDIF
1752
1753    ENDIF
1754#endif
1755!
1756!-- Output of virtual flight information
1757    IF ( virtual_flight )  CALL flight_header( io )
1758
1759!
1760!-- Output of spectra related quantities
1761    IF ( calculate_spectra )  CALL spectra_header( io )
1762
1763    WRITE ( io, 99 )
1764
1765!
1766!-- Physical quantities
1767    WRITE ( io, 400 )
1768
1769!
1770!-- Geostrophic parameters
1771    WRITE ( io, 410 )  phi, omega, f, fs
1772
1773!
1774!-- Other quantities
1775    WRITE ( io, 411 )  g
1776
1777    WRITE ( io, 412 )  TRIM( reference_state )
1778    IF ( use_single_reference_value )  THEN
1779       IF ( ocean )  THEN
1780          WRITE ( io, 413 )  prho_reference
1781       ELSE
1782          WRITE ( io, 414 )  pt_reference
1783       ENDIF
1784    ENDIF
1785
1786!
1787!-- Cloud physics parameters
1788    IF ( cloud_physics )  THEN
1789       WRITE ( io, 415 )
1790       WRITE ( io, 416 ) surface_pressure, r_d, rho_surface, cp, l_v
1791       IF ( microphysics_seifert )  THEN
1792          WRITE ( io, 510 ) 1.0E-6_wp * nc_const
1793          WRITE ( io, 511 ) c_sedimentation
1794       ENDIF
1795    ENDIF
1796
1797!
1798!-- Cloud physcis parameters / quantities / numerical methods
1799    WRITE ( io, 430 )
1800    IF ( humidity .AND. .NOT. cloud_physics .AND. .NOT. cloud_droplets)  THEN
1801       WRITE ( io, 431 )
1802    ELSEIF ( humidity  .AND.  cloud_physics )  THEN
1803       WRITE ( io, 432 )
1804       IF ( cloud_top_radiation )  WRITE ( io, 132 )
1805       IF ( microphysics_kessler )  THEN
1806          WRITE ( io, 133 )
1807       ELSEIF ( microphysics_seifert )  THEN
1808          IF ( cloud_water_sedimentation )  WRITE ( io, 506 )
1809          WRITE ( io, 505 )
1810          IF ( collision_turbulence )  WRITE ( io, 507 )
1811          IF ( ventilation_effect )  WRITE ( io, 508 )
1812          IF ( limiter_sedimentation )  WRITE ( io, 509 )
1813       ENDIF
1814    ELSEIF ( humidity  .AND.  cloud_droplets )  THEN
1815       WRITE ( io, 433 )
1816       IF ( curvature_solution_effects )  WRITE ( io, 434 )
1817       IF ( collision_kernel /= 'none' )  THEN
1818          WRITE ( io, 435 )  TRIM( collision_kernel )
1819          IF ( collision_kernel(6:9) == 'fast' )  THEN
1820             WRITE ( io, 436 )  radius_classes, dissipation_classes
1821          ENDIF
1822       ELSE
1823          WRITE ( io, 437 )
1824       ENDIF
1825    ENDIF
1826
1827!
1828!-- LES / turbulence parameters
1829    WRITE ( io, 450 )
1830
1831!--
1832! ... LES-constants used must still be added here
1833!--
1834    IF ( constant_diffusion )  THEN
1835       WRITE ( io, 451 )  km_constant, km_constant/prandtl_number, &
1836                          prandtl_number
1837    ENDIF
1838    IF ( .NOT. constant_diffusion)  THEN
1839       IF ( e_init > 0.0_wp )  WRITE ( io, 455 )  e_init
1840       IF ( e_min > 0.0_wp )  WRITE ( io, 454 )  e_min
1841       IF ( wall_adjustment )  WRITE ( io, 453 )  wall_adjustment_factor
1842    ENDIF
1843
1844!
1845!-- Special actions during the run
1846    WRITE ( io, 470 )
1847    IF ( create_disturbances )  THEN
1848       WRITE ( io, 471 )  dt_disturb, disturbance_amplitude,                   &
1849                          zu(disturbance_level_ind_b), disturbance_level_ind_b,&
1850                          zu(disturbance_level_ind_t), disturbance_level_ind_t
1851       IF ( .NOT. bc_lr_cyc  .OR.  .NOT. bc_ns_cyc )  THEN
1852          WRITE ( io, 472 )  inflow_disturbance_begin, inflow_disturbance_end
1853       ELSE
1854          WRITE ( io, 473 )  disturbance_energy_limit
1855       ENDIF
1856       WRITE ( io, 474 )  TRIM( random_generator )
1857    ENDIF
1858    IF ( pt_surface_initial_change /= 0.0_wp )  THEN
1859       WRITE ( io, 475 )  pt_surface_initial_change
1860    ENDIF
1861    IF ( humidity  .AND.  q_surface_initial_change /= 0.0_wp )  THEN
1862       WRITE ( io, 476 )  q_surface_initial_change       
1863    ENDIF
1864    IF ( passive_scalar  .AND.  q_surface_initial_change /= 0.0_wp )  THEN
1865       WRITE ( io, 477 )  q_surface_initial_change       
1866    ENDIF
1867
1868    IF ( particle_advection )  THEN
1869!
1870!--    Particle attributes
1871       WRITE ( io, 480 )  particle_advection_start, dt_prel, bc_par_lr, &
1872                          bc_par_ns, bc_par_b, bc_par_t, particle_maximum_age, &
1873                          end_time_prel
1874       IF ( use_sgs_for_particles )  WRITE ( io, 488 )  dt_min_part
1875       IF ( random_start_position )  WRITE ( io, 481 )
1876       IF ( seed_follows_topography )  WRITE ( io, 496 )
1877       IF ( particles_per_point > 1 )  WRITE ( io, 489 )  particles_per_point
1878       WRITE ( io, 495 )  total_number_of_particles
1879       IF ( dt_write_particle_data /= 9999999.9_wp )  THEN
1880          WRITE ( io, 485 )  dt_write_particle_data
1881          IF ( netcdf_data_format > 1 )  THEN
1882             output_format = 'netcdf (64 bit offset) and binary'
1883          ELSE
1884             output_format = 'netcdf and binary'
1885          ENDIF
1886          IF ( netcdf_deflate == 0 )  THEN
1887             WRITE ( io, 344 )  output_format
1888          ELSE
1889             WRITE ( io, 354 )  TRIM( output_format ), netcdf_deflate
1890          ENDIF
1891       ENDIF
1892       IF ( dt_dopts /= 9999999.9_wp )  WRITE ( io, 494 )  dt_dopts
1893       IF ( write_particle_statistics )  WRITE ( io, 486 )
1894
1895       WRITE ( io, 487 )  number_of_particle_groups
1896
1897       DO  i = 1, number_of_particle_groups
1898          IF ( i == 1  .AND.  density_ratio(i) == 9999999.9_wp )  THEN
1899             WRITE ( io, 490 )  i, 0.0_wp
1900             WRITE ( io, 492 )
1901          ELSE
1902             WRITE ( io, 490 )  i, radius(i)
1903             IF ( density_ratio(i) /= 0.0_wp )  THEN
1904                WRITE ( io, 491 )  density_ratio(i)
1905             ELSE
1906                WRITE ( io, 492 )
1907             ENDIF
1908          ENDIF
1909          WRITE ( io, 493 )  psl(i), psr(i), pss(i), psn(i), psb(i), pst(i), &
1910                             pdx(i), pdy(i), pdz(i)
1911          IF ( .NOT. vertical_particle_advection(i) )  WRITE ( io, 482 )
1912       ENDDO
1913
1914    ENDIF
1915
1916
1917!
1918!-- Parameters of 1D-model
1919    IF ( INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 )  THEN
1920       WRITE ( io, 500 )  end_time_1d, dt_run_control_1d, dt_pr_1d, &
1921                          mixing_length_1d, dissipation_1d
1922       IF ( damp_level_ind_1d /= nzt+1 )  THEN
1923          WRITE ( io, 502 )  zu(damp_level_ind_1d), damp_level_ind_1d
1924       ENDIF
1925    ENDIF
1926
1927!
1928!-- User-defined information
1929    CALL user_header( io )
1930
1931    WRITE ( io, 99 )
1932
1933!
1934!-- Write buffer contents to disc immediately
1935    FLUSH( io )
1936
1937!
1938!-- Here the FORMATs start
1939
1940 99 FORMAT (1X,78('-'))
1941100 FORMAT (/1X,'******************************',4X,44('-')/        &
1942            1X,'* ',A,' *',4X,A/                               &
1943            1X,'******************************',4X,44('-'))
1944101 FORMAT (35X,'coupled run: ',A/ &
1945            35X,42('-'))
1946102 FORMAT (/' Date:                 ',A8,4X,'Run:       ',A20/      &
1947            ' Time:                 ',A8,4X,'Run-No.:   ',I2.2/     &
1948            ' Run on host:        ',A10)
1949#if defined( __parallel )
1950103 FORMAT (' Number of PEs:',10X,I6,4X,'Processor grid (x,y): (',I4,',',I4, &
1951              ')',1X,A)
1952104 FORMAT (' Number of PEs:',10X,I6,4X,'Tasks:',I4,'   threads per task:',I4/ &
1953              35X,'Processor grid (x,y): (',I4,',',I4,')',1X,A)
1954105 FORMAT (35X,'One additional PE is used to handle'/37X,'the dvrp output!')
1955106 FORMAT (35X,'A 1d-decomposition along x is forced'/ &
1956            35X,'because the job is running on an SMP-cluster')
1957107 FORMAT (35X,'A 1d-decomposition along ',A,' is used')
1958108 FORMAT (35X,'Max. # of parallel I/O streams is ',I5)
1959109 FORMAT (35X,'Precursor run for coupled atmos-ocean run'/ &
1960            35X,42('-'))
1961114 FORMAT (35X,'Coupled atmosphere-ocean run following'/ &
1962            35X,'independent precursor runs'/             &
1963            35X,42('-'))
1964#endif
1965110 FORMAT (/' Numerical Schemes:'/ &
1966             ' -----------------'/)
1967121 FORMAT (' --> Use the ',A,' approximation for the model equations.')
1968111 FORMAT (' --> Solve perturbation pressure via FFT using ',A,' routines')
1969112 FORMAT (' --> Solve perturbation pressure via SOR-Red/Black-Schema'/ &
1970            '     Iterations (initial/other): ',I3,'/',I3,'  omega =',F6.3)
1971113 FORMAT (' --> Momentum advection via Piascek-Williams-Scheme (Form C3)', &
1972                  ' or Upstream')
1973115 FORMAT ('     FFT and transpositions are overlapping')
1974116 FORMAT (' --> Scalar advection via Piascek-Williams-Scheme (Form C3)', &
1975                  ' or Upstream')
1976118 FORMAT (' --> Scalar advection via Bott-Chlond-Scheme')
1977119 FORMAT (' --> Galilei-Transform applied to horizontal advection:'/ &
1978            '     translation velocity = ',A/ &
1979            '     distance advected ',A,':  ',F8.3,' km(x)  ',F8.3,' km(y)')
1980122 FORMAT (' --> Time differencing scheme: ',A)
1981123 FORMAT (' --> Rayleigh-Damping active, starts ',A,' z = ',F8.2,' m'/ &
1982            '     maximum damping coefficient:',F6.3, ' 1/s')
1983129 FORMAT (' --> Additional prognostic equation for the specific humidity')
1984130 FORMAT (' --> Additional prognostic equation for the total water content')
1985131 FORMAT (' --> No pt-equation solved. Neutral stratification with pt = ', &
1986                  F6.2, ' K assumed')
1987132 FORMAT ('     Parameterization of long-wave radiation processes via'/ &
1988            '     effective emissivity scheme')
1989133 FORMAT ('     Precipitation parameterization via Kessler-Scheme')
1990134 FORMAT (' --> Additional prognostic equation for a passive scalar')
1991135 FORMAT (' --> Solve perturbation pressure via ',A,' method (', &
1992                  A,'-cycle)'/ &
1993            '     number of grid levels:                   ',I2/ &
1994            '     Gauss-Seidel red/black iterations:       ',I2)
1995136 FORMAT ('     gridpoints of coarsest subdomain (x,y,z): (',I3,',',I3,',', &
1996                  I3,')')
1997137 FORMAT ('     level data gathered on PE0 at level:     ',I2/ &
1998            '     gridpoints of coarsest subdomain (x,y,z): (',I3,',',I3,',', &
1999                  I3,')'/ &
2000            '     gridpoints of coarsest domain (x,y,z):    (',I3,',',I3,',', &
2001                  I3,')')
2002139 FORMAT (' --> Loop optimization method: ',A)
2003140 FORMAT ('     maximum residual allowed:                ',E10.3)
2004141 FORMAT ('     fixed number of multigrid cycles:        ',I4)
2005142 FORMAT ('     perturbation pressure is calculated at every Runge-Kutta ', &
2006                  'step')
2007143 FORMAT ('     Euler/upstream scheme is used for the SGS turbulent ', &
2008                  'kinetic energy')
2009144 FORMAT ('     masking method is used')
2010150 FORMAT (' --> Volume flow at the right and north boundary will be ', &
2011                  'conserved'/ &
2012            '     using the ',A,' mode')
2013151 FORMAT ('     with u_bulk = ',F7.3,' m/s and v_bulk = ',F7.3,' m/s')
2014152 FORMAT (' --> External pressure gradient directly prescribed by the user:',&
2015           /'     ',2(1X,E12.5),'Pa/m in x/y direction', &
2016           /'     starting from dp_level_b =', F8.3, 'm', A /)
2017160 FORMAT (//' Large scale forcing and nudging:'/ &
2018              ' -------------------------------'/)
2019161 FORMAT (' --> No large scale forcing from external is used (default) ')
2020162 FORMAT (' --> Large scale forcing from external file LSF_DATA is used: ')
2021163 FORMAT ('     - large scale advection tendencies ')
2022164 FORMAT ('     - large scale subsidence velocity w_subs ')
2023165 FORMAT ('     - large scale subsidence tendencies ')
2024167 FORMAT ('     - and geostrophic wind components ug and vg')
2025168 FORMAT (' --> Large-scale vertical motion is used in the ', &
2026                  'prognostic equation(s) for')
2027169 FORMAT ('     the scalar(s) only')
2028170 FORMAT (' --> Nudging is used')
2029171 FORMAT (' --> No nudging is used (default) ')
2030180 FORMAT ('     - prescribed surface values for temperature')
2031181 FORMAT ('     - prescribed surface fluxes for temperature')
2032182 FORMAT ('     - prescribed surface values for humidity')
2033183 FORMAT ('     - prescribed surface fluxes for humidity')
2034200 FORMAT (//' Run time and time step information:'/ &
2035             ' ----------------------------------'/)
2036201 FORMAT ( ' Timestep:             variable     maximum value: ',F6.3,' s', &
2037             '    CFL-factor:',F5.2)
2038202 FORMAT ( ' Timestep:          dt = ',F6.3,' s'/)
2039203 FORMAT ( ' Start time:          ',F9.3,' s'/ &
2040             ' End time:            ',F9.3,' s')
2041204 FORMAT ( A,F9.3,' s')
2042205 FORMAT ( A,F9.3,' s',5X,'restart every',17X,F9.3,' s')
2043206 FORMAT (/' Time reached:        ',F9.3,' s'/ &
2044             ' CPU-time used:       ',F9.3,' s     per timestep:               ', &
2045               '  ',F9.3,' s'/                                                    &
2046             '                                      per second of simulated tim', &
2047               'e: ',F9.3,' s')
2048207 FORMAT ( ' Coupling start time: ',F9.3,' s')
2049250 FORMAT (//' Computational grid and domain size:'/ &
2050              ' ----------------------------------'// &
2051              ' Grid length:      dx =    ',F7.3,' m    dy =    ',F7.3, &
2052              ' m    dz =    ',F7.3,' m'/ &
2053              ' Domain size:       x = ',F10.3,' m     y = ',F10.3, &
2054              ' m  z(u) = ',F10.3,' m'/)
2055252 FORMAT (' dz constant up to ',F10.3,' m (k=',I4,'), then stretched by', &
2056              ' factor:',F6.3/ &
2057            ' maximum dz not to be exceeded is dz_max = ',F10.3,' m'/)
2058254 FORMAT (' Number of gridpoints (x,y,z):  (0:',I4,', 0:',I4,', 0:',I4,')'/ &
2059            ' Subdomain size (x,y,z):        (  ',I4,',   ',I4,',   ',I4,')'/)
2060260 FORMAT (/' The model has a slope in x-direction. Inclination angle: ',F6.2,&
2061             ' degrees')
2062270 FORMAT (//' Topography information:'/ &
2063              ' ----------------------'// &
2064              1X,'Topography: ',A)
2065271 FORMAT (  ' Building size (x/y/z) in m: ',F5.1,' / ',F5.1,' / ',F5.1/ &
2066              ' Horizontal index bounds (l/r/s/n): ',I4,' / ',I4,' / ',I4, &
2067                ' / ',I4)
2068272 FORMAT (  ' Single quasi-2D street canyon of infinite length in ',A, &
2069              ' direction' / &
2070              ' Canyon height: ', F6.2, 'm, ch = ', I4, '.'      / &
2071              ' Canyon position (',A,'-walls): cxl = ', I4,', cxr = ', I4, '.')
2072273 FORMAT (  ' Tunnel of infinite length in ',A, &
2073              ' direction' / &
2074              ' Tunnel height: ', F6.2, / &
2075              ' Tunnel-wall depth: ', F6.2      / &
2076              ' Tunnel width: ', F6.2 )
2077274 FORMAT (  ' Tunnel in ', A, ' direction.' / &
2078              ' Tunnel height: ', F6.2, / &   
2079              ' Tunnel-wall depth: ', F6.2      / &
2080              ' Tunnel width: ', F6.2, / &
2081              ' Tunnel length: ', F6.2 )
2082278 FORMAT (' Topography grid definition convention:'/ &
2083            ' cell edge (staggered grid points'/  &
2084            ' (u in x-direction, v in y-direction))' /)
2085279 FORMAT (' Topography grid definition convention:'/ &
2086            ' cell center (scalar grid points)' /)
2087300 FORMAT (//' Boundary conditions:'/ &
2088             ' -------------------'// &
2089             '                     p                    uv             ', &
2090             '                     pt'// &
2091             ' B. bound.: ',A/ &
2092             ' T. bound.: ',A)
2093301 FORMAT (/'                     ',A// &
2094             ' B. bound.: ',A/ &
2095             ' T. bound.: ',A)
2096303 FORMAT (/' Bottom surface fluxes are used in diffusion terms at k=1')
2097304 FORMAT (/' Top surface fluxes are used in diffusion terms at k=nzt')
2098305 FORMAT (//'    Constant flux layer between bottom surface and first ',     &
2099              'computational u,v-level:'// &
2100             '       z_mo = ',F6.2,' m   z0 =',F7.4,' m   z0h =',F8.5,&
2101             ' m   kappa =',F5.2/ &
2102             '       Rif value range:   ',F8.2,' <= rif <=',F6.2)
2103306 FORMAT ('       Predefined constant heatflux:   ',F9.6,' K m/s')
2104307 FORMAT ('       Heatflux has a random normal distribution')
2105308 FORMAT ('       Predefined surface temperature')
2106309 FORMAT ('       Predefined constant salinityflux:   ',F9.6,' psu m/s')
2107310 FORMAT (//'    1D-Model:'// &
2108             '       Rif value range:   ',F6.2,' <= rif <=',F6.2)
2109311 FORMAT ('       Predefined constant humidity flux: ',E10.3,' kg/kg m/s')
2110312 FORMAT ('       Predefined surface humidity')
2111313 FORMAT ('       Predefined constant scalar flux: ',E10.3,' kg/(m**2 s)')
2112314 FORMAT ('       Predefined scalar value at the surface')
2113302 FORMAT ('       Predefined constant scalarflux:   ',F9.6,' kg/(m**2 s)')
2114315 FORMAT ('       Humidity flux at top surface is 0.0')
2115316 FORMAT ('       Sensible heatflux and momentum flux from coupled ', &
2116                    'atmosphere model')
2117317 FORMAT (//' Lateral boundaries:'/ &
2118            '       left/right:  ',A/    &
2119            '       north/south: ',A)
2120318 FORMAT (/'       use_cmax: ',L1 / &
2121            '       pt damping layer width = ',F8.2,' m, pt ', &
2122                    'damping factor =',F7.4)
2123319 FORMAT ('       turbulence recycling at inflow switched on'/ &
2124            '       width of recycling domain: ',F7.1,' m   grid index: ',I4/ &
2125            '       inflow damping height: ',F6.1,' m   width: ',F6.1,' m')
2126320 FORMAT ('       Predefined constant momentumflux:  u: ',F9.6,' m**2/s**2'/ &
2127            '                                          v: ',F9.6,' m**2/s**2')
2128321 FORMAT (//' Initial profiles:'/ &
2129              ' ----------------')
2130322 FORMAT ('       turbulence recycling at inflow switched on'/ &
2131            '       y shift of the recycled inflow turbulence switched on'/ &
2132            '       width of recycling domain: ',F7.1,' m   grid index: ',I4/ &
2133            '       inflow damping height: ',F6.1,' m   width: ',F6.1,' m'/)
2134323 FORMAT ('       turbulent outflow conditon switched on'/ &
2135            '       position of outflow source plane: ',F7.1,' m   ', &
2136                    'grid index: ', I4)
2137325 FORMAT (//' List output:'/ &
2138             ' -----------'//  &
2139            '    1D-Profiles:'/    &
2140            '       Output every             ',F8.2,' s')
2141326 FORMAT ('       Time averaged over       ',F8.2,' s'/ &
2142            '       Averaging input every    ',F8.2,' s')
2143330 FORMAT (//' Data output:'/ &
2144             ' -----------'/)
2145331 FORMAT (/'    1D-Profiles:')
2146332 FORMAT (/'       ',A)
2147333 FORMAT ('       Output every             ',F8.2,' s',/ &
2148            '       Time averaged over       ',F8.2,' s'/ &
2149            '       Averaging input every    ',F8.2,' s')
2150334 FORMAT (/'    2D-Arrays',A,':')
2151335 FORMAT (/'       ',A2,'-cross-section  Arrays: ',A/ &
2152            '       Output every             ',F8.2,' s  ',A/ &
2153            '       Cross sections at ',A1,' = ',A/ &
2154            '       scalar-coordinates:   ',A,' m'/)
2155336 FORMAT (/'    3D-Arrays',A,':')
2156337 FORMAT (/'       Arrays: ',A/ &
2157            '       Output every             ',F8.2,' s  ',A/ &
2158            '       Upper output limit at    ',F8.2,' m  (GP ',I4,')'/)
2159339 FORMAT ('       No output during initial ',F8.2,' s')
2160340 FORMAT (/'    Time series:')
2161341 FORMAT ('       Output every             ',F8.2,' s'/)
2162342 FORMAT (/'       ',A2,'-cross-section  Arrays: ',A/ &
2163            '       Output every             ',F8.2,' s  ',A/ &
2164            '       Time averaged over       ',F8.2,' s'/ &
2165            '       Averaging input every    ',F8.2,' s'/ &
2166            '       Cross sections at ',A1,' = ',A/ &
2167            '       scalar-coordinates:   ',A,' m'/)
2168343 FORMAT (/'       Arrays: ',A/ &
2169            '       Output every             ',F8.2,' s  ',A/ &
2170            '       Time averaged over       ',F8.2,' s'/ &
2171            '       Averaging input every    ',F8.2,' s'/ &
2172            '       Upper output limit at    ',F8.2,' m  (GP ',I4,')'/)
2173344 FORMAT ('       Output format: ',A/)
2174345 FORMAT (/'    Scaling lengths for output locations of all subsequent mask IDs:',/ &
2175            '       mask_scale_x (in x-direction): ',F9.3, ' m',/ &
2176            '       mask_scale_y (in y-direction): ',F9.3, ' m',/ &
2177            '       mask_scale_z (in z-direction): ',F9.3, ' m' )
2178346 FORMAT (/'    Masked data output',A,' for mask ID ',I2, ':')
2179347 FORMAT ('       Variables: ',A/ &
2180            '       Output every             ',F8.2,' s')
2181348 FORMAT ('       Variables: ',A/ &
2182            '       Output every             ',F8.2,' s'/ &
2183            '       Time averaged over       ',F8.2,' s'/ &
2184            '       Averaging input every    ',F8.2,' s')
2185349 FORMAT (/'       Output locations in ',A,'-direction in multiples of ', &
2186            'mask_scale_',A,' predefined by array mask_',I2.2,'_',A,':'/ &
2187            13('       ',8(F8.2,',')/) )
2188350 FORMAT (/'       Output locations in ',A,'-direction: ', &
2189            'all gridpoints along ',A,'-direction (default).' )
2190351 FORMAT (/'       Output locations in ',A,'-direction in multiples of ', &
2191            'mask_scale_',A,' constructed from array mask_',I2.2,'_',A,'_loop:'/ &
2192            '          loop begin:',F8.2,', end:',F8.2,', stride:',F8.2 )
2193352 FORMAT  (/'       Number of output time levels allowed: ',I3 /)
2194353 FORMAT  (/'       Number of output time levels allowed: unlimited' /)
2195354 FORMAT ('       Output format: ',A, '   compressed with level: ',I1/)
2196#if defined( __dvrp_graphics )
2197360 FORMAT ('    Plot-Sequence with dvrp-software:'/ &
2198            '       Output every      ',F7.1,' s'/ &
2199            '       Output mode:      ',A/ &
2200            '       Host / User:      ',A,' / ',A/ &
2201            '       Directory:        ',A// &
2202            '       The sequence contains:')
2203361 FORMAT (/'       Isosurface of "',A,'"    Threshold value: ', E12.3/ &
2204            '          Isosurface color: (',F4.2,',',F4.2,',',F4.2,') (R,G,B)')
2205362 FORMAT (/'       Slicer plane ',A/ &
2206            '       Slicer limits: [',F6.2,',',F6.2,']')
2207365 FORMAT (/'       Groundplate color: (',F4.2,',',F4.2,',',F4.2,') (R,G,B)'/ &
2208            '       Superelevation along (x,y,z): (',F4.1,',',F4.1,',',F4.1, &
2209                     ')'/ &
2210            '       Clipping limits: from x = ',F9.1,' m to x = ',F9.1,' m'/ &
2211            '                        from y = ',F9.1,' m to y = ',F9.1,' m')
2212366 FORMAT (/'       Topography color: (',F4.2,',',F4.2,',',F4.2,') (R,G,B)')
2213367 FORMAT ('       Polygon reduction for topography: cluster_size = ', I1)
2214#endif
2215400 FORMAT (//' Physical quantities:'/ &
2216              ' -------------------'/)
2217410 FORMAT ('    Geograph. latitude  :   phi    = ',F4.1,' degr'/   &
2218            '    Angular velocity    :   omega  =',E10.3,' rad/s'/  &
2219            '    Coriolis parameter  :   f      = ',F9.6,' 1/s'/    &
2220            '                            f*     = ',F9.6,' 1/s')
2221411 FORMAT (/'    Gravity             :   g      = ',F4.1,' m/s**2')
2222412 FORMAT (/'    Reference state used in buoyancy terms: ',A)
2223413 FORMAT ('       Reference density in buoyancy terms: ',F8.3,' kg/m**3')
2224414 FORMAT ('       Reference temperature in buoyancy terms: ',F8.4,' K')
2225415 FORMAT (/' Cloud physics parameters:'/ &
2226             ' ------------------------'/)
2227416 FORMAT ('    Surface pressure   :   p_0   = ',F7.2,' hPa'/      &
2228            '    Gas constant       :   R     = ',F5.1,' J/(kg K)'/ &
2229            '    Density of air     :   rho_0 =',F6.3,' kg/m**3'/  &
2230            '    Specific heat cap. :   c_p   = ',F6.1,' J/(kg K)'/ &
2231            '    Vapourization heat :   L_v   =',E9.2,' J/kg')
2232417 FORMAT ('    Geograph. longitude :   lambda = ',F4.1,' degr')
2233418 FORMAT (/'    Day of the year at model start :   day_init      =     ',I3 &
2234            /'    UTC time at model start        :   time_utc_init = ',F7.1' s')
2235420 FORMAT (/'    Characteristic levels of the initial temperature profile:'// &
2236            '       Height:        ',A,'  m'/ &
2237            '       Temperature:   ',A,'  K'/ &
2238            '       Gradient:      ',A,'  K/100m'/ &
2239            '       Gridpoint:     ',A)
2240421 FORMAT (/'    Characteristic levels of the initial humidity profile:'// &
2241            '       Height:      ',A,'  m'/ &
2242            '       Humidity:    ',A,'  kg/kg'/ &
2243            '       Gradient:    ',A,'  (kg/kg)/100m'/ &
2244            '       Gridpoint:   ',A)
2245422 FORMAT (/'    Characteristic levels of the initial scalar profile:'// &
2246            '       Height:                  ',A,'  m'/ &
2247            '       Scalar concentration:    ',A,'  kg/m**3'/ &
2248            '       Gradient:                ',A,'  (kg/m**3)/100m'/ &
2249            '       Gridpoint:               ',A)
2250423 FORMAT (/'    Characteristic levels of the geo. wind component ug:'// &
2251            '       Height:      ',A,'  m'/ &
2252            '       ug:          ',A,'  m/s'/ &
2253            '       Gradient:    ',A,'  1/100s'/ &
2254            '       Gridpoint:   ',A)
2255424 FORMAT (/'    Characteristic levels of the geo. wind component vg:'// &
2256            '       Height:      ',A,'  m'/ &
2257            '       vg:          ',A,'  m/s'/ &
2258            '       Gradient:    ',A,'  1/100s'/ &
2259            '       Gridpoint:   ',A)
2260425 FORMAT (/'    Characteristic levels of the initial salinity profile:'// &
2261            '       Height:     ',A,'  m'/ &
2262            '       Salinity:   ',A,'  psu'/ &
2263            '       Gradient:   ',A,'  psu/100m'/ &
2264            '       Gridpoint:  ',A)
2265426 FORMAT (/'    Characteristic levels of the subsidence/ascent profile:'// &
2266            '       Height:      ',A,'  m'/ &
2267            '       w_subs:      ',A,'  m/s'/ &
2268            '       Gradient:    ',A,'  (m/s)/100m'/ &
2269            '       Gridpoint:   ',A)
2270427 FORMAT (/'    Initial wind profiles (u,v) are interpolated from given'// &
2271                  ' profiles')
2272428 FORMAT (/'    Initial profiles (u, v, pt, q) are taken from file '/ &
2273             '    NUDGING_DATA')
2274430 FORMAT (//' Cloud physics quantities / methods:'/ &
2275              ' ----------------------------------'/)
2276431 FORMAT ('    Humidity is considered, bu no condensation')
2277432 FORMAT ('    Bulk scheme with liquid water potential temperature and'/ &
2278            '    total water content is used.'/ &
2279            '    Condensation is parameterized via 0% - or 100% scheme.')
2280433 FORMAT ('    Cloud droplets treated explicitly using the Lagrangian part', &
2281                 'icle model')
2282434 FORMAT ('    Curvature and solution effecs are considered for growth of', &
2283                 ' droplets < 1.0E-6 m')
2284435 FORMAT ('    Droplet collision is handled by ',A,'-kernel')
2285436 FORMAT ('       Fast kernel with fixed radius- and dissipation classes ', &
2286                    'are used'/ &
2287            '          number of radius classes:       ',I3,'    interval ', &
2288                       '[1.0E-6,2.0E-4] m'/ &
2289            '          number of dissipation classes:   ',I2,'    interval ', &
2290                       '[0,1000] cm**2/s**3')
2291437 FORMAT ('    Droplet collision is switched off')
2292450 FORMAT (//' LES / Turbulence quantities:'/ &
2293              ' ---------------------------'/)
2294451 FORMAT ('    Diffusion coefficients are constant:'/ &
2295            '    Km = ',F6.2,' m**2/s   Kh = ',F6.2,' m**2/s   Pr = ',F5.2)
2296453 FORMAT ('    Mixing length is limited to',F5.2,' * z')
2297454 FORMAT ('    TKE is not allowed to fall below ',E9.2,' (m/s)**2')
2298455 FORMAT ('    initial TKE is prescribed as ',E9.2,' (m/s)**2')
2299470 FORMAT (//' Actions during the simulation:'/ &
2300              ' -----------------------------'/)
2301471 FORMAT ('    Disturbance impulse (u,v) every :   ',F6.2,' s'/            &
2302            '    Disturbance amplitude           :    ',F5.2, ' m/s'/       &
2303            '    Lower disturbance level         : ',F8.2,' m (GP ',I4,')'/  &
2304            '    Upper disturbance level         : ',F8.2,' m (GP ',I4,')')
2305472 FORMAT ('    Disturbances continued during the run from i/j =',I4, &
2306                 ' to i/j =',I4)
2307473 FORMAT ('    Disturbances cease as soon as the disturbance energy exceeds',&
2308                 F6.3, ' m**2/s**2')
2309474 FORMAT ('    Random number generator used    : ',A/)
2310475 FORMAT ('    The surface temperature is increased (or decreased, ', &
2311                 'respectively, if'/ &
2312            '    the value is negative) by ',F5.2,' K at the beginning of the',&
2313                 ' 3D-simulation'/)
2314476 FORMAT ('    The surface humidity is increased (or decreased, ',&
2315                 'respectively, if the'/ &
2316            '    value is negative) by ',E8.1,' kg/kg at the beginning of', &
2317                 ' the 3D-simulation'/)
2318477 FORMAT ('    The scalar value is increased at the surface (or decreased, ',&
2319                 'respectively, if the'/ &
2320            '    value is negative) by ',E8.1,' kg/m**3 at the beginning of', &
2321                 ' the 3D-simulation'/)
2322480 FORMAT ('    Particles:'/ &
2323            '    ---------'// &
2324            '       Particle advection is active (switched on at t = ', F7.1, &
2325                    ' s)'/ &
2326            '       Start of new particle generations every  ',F6.1,' s'/ &
2327            '       Boundary conditions: left/right: ', A, ' north/south: ', A/&
2328            '                            bottom:     ', A, ' top:         ', A/&
2329            '       Maximum particle age:                 ',F9.1,' s'/ &
2330            '       Advection stopped at t = ',F9.1,' s'/)
2331481 FORMAT ('       Particles have random start positions'/)
2332482 FORMAT ('          Particles are advected only horizontally'/)
2333485 FORMAT ('       Particle data are written on file every ', F9.1, ' s')
2334486 FORMAT ('       Particle statistics are written on file'/)
2335487 FORMAT ('       Number of particle groups: ',I2/)
2336488 FORMAT ('       SGS velocity components are used for particle advection'/ &
2337            '          minimum timestep for advection:', F8.5/)
2338489 FORMAT ('       Number of particles simultaneously released at each ', &
2339                    'point: ', I5/)
2340490 FORMAT ('       Particle group ',I2,':'/ &
2341            '          Particle radius: ',E10.3, 'm')
2342491 FORMAT ('          Particle inertia is activated'/ &
2343            '             density_ratio (rho_fluid/rho_particle) =',F6.3/)
2344492 FORMAT ('          Particles are advected only passively (no inertia)'/)
2345493 FORMAT ('          Boundaries of particle source: x:',F8.1,' - ',F8.1,' m'/&
2346            '                                         y:',F8.1,' - ',F8.1,' m'/&
2347            '                                         z:',F8.1,' - ',F8.1,' m'/&
2348            '          Particle distances:  dx = ',F8.1,' m  dy = ',F8.1, &
2349                       ' m  dz = ',F8.1,' m'/)
2350494 FORMAT ('       Output of particle time series in NetCDF format every ', &
2351                    F8.2,' s'/)
2352495 FORMAT ('       Number of particles in total domain: ',I10/)
2353496 FORMAT ('       Initial vertical particle positions are interpreted ', &
2354                    'as relative to the given topography')
2355500 FORMAT (//' 1D-Model parameters:'/                           &
2356              ' -------------------'//                           &
2357            '    Simulation time:                   ',F8.1,' s'/ &
2358            '    Run-controll output every:         ',F8.1,' s'/ &
2359            '    Vertical profile output every:     ',F8.1,' s'/ &
2360            '    Mixing length calculation:         ',A/         &
2361            '    Dissipation calculation:           ',A/)
2362502 FORMAT ('    Damping layer starts from ',F7.1,' m (GP ',I4,')'/)
2363503 FORMAT (' --> Momentum advection via Wicker-Skamarock-Scheme 5th order')
2364504 FORMAT (' --> Scalar advection via Wicker-Skamarock-Scheme 5th order')
2365505 FORMAT ('    Precipitation parameterization via Seifert-Beheng-Scheme')
2366506 FORMAT ('    Cloud water sedimentation parameterization via Stokes law')
2367507 FORMAT ('    Turbulence effects on precipitation process')
2368508 FORMAT ('    Ventilation effects on evaporation of rain drops')
2369509 FORMAT ('    Slope limiter used for sedimentation process')
2370510 FORMAT ('    Droplet density    :   N_c   = ',F6.1,' 1/cm**3')
2371511 FORMAT ('    Sedimentation Courant number:                  '/&
2372            '                               C_s   =',F4.1,'        ')
2373512 FORMAT (/' Date:                 ',A8,6X,'Run:       ',A20/      &
2374            ' Time:                 ',A8,6X,'Run-No.:   ',I2.2/     &
2375            ' Run on host:        ',A10,6X,'En-No.:    ',I2.2)
2376600 FORMAT (/' Nesting informations:'/ &
2377            ' --------------------'/ &
2378            ' Nesting mode:                     ',A/ &
2379            ' Nesting-datatransfer mode:        ',A// &
2380            ' Nest id  parent  number   lower left coordinates   name'/ &
2381            ' (*=me)     id    of PEs      x (m)     y (m)' )
2382601 FORMAT (2X,A1,1X,I2.2,6X,I2.2,5X,I5,5X,F8.2,2X,F8.2,5X,A)
2383
2384 END SUBROUTINE header
Note: See TracBrowser for help on using the repository browser.