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

Last change on this file since 1682 was 1682, checked in by knoop, 6 years ago

Code annotations made doxygen readable

  • Property svn:keywords set to Id
File size: 92.2 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 terms
6! of the GNU General Public License as published by the Free Software Foundation,
7! either version 3 of the License, or (at your option) any later version.
8!
9! PALM is distributed in the hope that it will be useful, but WITHOUT ANY
10! WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
11! A PARTICULAR PURPOSE.  See the GNU General Public License for more details.
12!
13! You should have received a copy of the GNU General Public License along with
14! PALM. If not, see <http://www.gnu.org/licenses/>.
15!
16! Copyright 1997-2014 Leibniz Universitaet Hannover
17!--------------------------------------------------------------------------------!
18!
19! Current revisions:
20! -----------------
21! Code annotations made doxygen readable
22!
23! Former revisions:
24! -----------------
25! $Id: header.f90 1682 2015-10-07 23:56:08Z knoop $
26!
27! 1675 2015-10-02 08:28:59Z gronemeier
28! Bugfix: Definition of topography grid levels
29!
30! 1660 2015-09-21 08:15:16Z gronemeier
31! Bugfix: Definition of building/street canyon height if vertical grid stretching
32!         starts below the maximum topography height.
33!
34! 1590 2015-05-08 13:56:27Z maronga
35! Bugfix: Added TRIM statements for character strings for LSM and radiation code
36!
37! 1585 2015-04-30 07:05:52Z maronga
38! Further output for radiation model(s).
39!
40! 1575 2015-03-27 09:56:27Z raasch
41! adjustments for psolver-queries, output of seed_follows_topography
42!
43! 1560 2015-03-06 10:48:54Z keck
44! output for recycling y shift
45!
46! 1557 2015-03-05 16:43:04Z suehring
47! output for monotonic limiter
48!
49! 1551 2015-03-03 14:18:16Z maronga
50! Added informal output for land surface model and radiation model. Removed typo.
51!
52! 1496 2014-12-02 17:25:50Z maronga
53! Renamed: "radiation -> "cloud_top_radiation"
54!
55! 1484 2014-10-21 10:53:05Z kanani
56! Changes due to new module structure of the plant canopy model:
57!   module plant_canopy_model_mod and output for new canopy model parameters
58!   (alpha_lad, beta_lad, lai_beta,...) added,
59!   drag_coefficient, leaf_surface_concentration and scalar_exchange_coefficient
60!   renamed to canopy_drag_coeff, leaf_surface_conc and leaf_scalar_exch_coeff,
61!   learde renamed leaf_area_density.
62! Bugfix: DO-WHILE-loop for lad header information additionally restricted
63! by maximum number of gradient levels (currently 10)
64!
65! 1482 2014-10-18 12:34:45Z raasch
66! information about calculated or predefined virtual processor topology adjusted
67!
68! 1468 2014-09-24 14:06:57Z maronga
69! Adapted for use on up to 6-digit processor cores
70!
71! 1429 2014-07-15 12:53:45Z knoop
72! header exended to provide ensemble_member_nr if specified
73!
74! 1376 2014-04-26 11:21:22Z boeske
75! Correction of typos
76!
77! 1365 2014-04-22 15:03:56Z boeske
78! New section 'Large scale forcing and nudging':
79! output of large scale forcing and nudging information,
80! new section for initial profiles created
81!
82! 1359 2014-04-11 17:15:14Z hoffmann
83! dt_sort_particles removed
84!
85! 1353 2014-04-08 15:21:23Z heinze
86! REAL constants provided with KIND-attribute
87!
88! 1327 2014-03-21 11:00:16Z raasch
89! parts concerning iso2d and avs output removed,
90! -netcdf output queries
91!
92! 1324 2014-03-21 09:13:16Z suehring
93! Bugfix: module spectrum added
94!
95! 1322 2014-03-20 16:38:49Z raasch
96! REAL functions provided with KIND-attribute,
97! some REAL constants defined as wp-kind
98!
99! 1320 2014-03-20 08:40:49Z raasch
100! ONLY-attribute added to USE-statements,
101! kind-parameters added to all INTEGER and REAL declaration statements,
102! kinds are defined in new module kinds,
103! revision history before 2012 removed,
104! comment fields (!:) to be used for variable explanations added to
105! all variable declaration statements
106!
107! 1308 2014-03-13 14:58:42Z fricke
108! output of the fixed number of output time levels
109! output_format adjusted for masked data if netcdf_data_format > 5
110!
111! 1299 2014-03-06 13:15:21Z heinze
112! output for using large_scale subsidence in combination
113! with large_scale_forcing
114! reformatting, more detailed explanations
115!
116! 1241 2013-10-30 11:36:58Z heinze
117! output for nudging + large scale forcing from external file
118!
119! 1216 2013-08-26 09:31:42Z raasch
120! output for transpose_compute_overlap
121!
122! 1212 2013-08-15 08:46:27Z raasch
123! output for poisfft_hybrid removed
124!
125! 1179 2013-06-14 05:57:58Z raasch
126! output of reference_state, use_reference renamed use_single_reference_value
127!
128! 1159 2013-05-21 11:58:22Z fricke
129! +use_cmax
130!
131! 1115 2013-03-26 18:16:16Z hoffmann
132! descriptions for Seifert-Beheng-cloud-physics-scheme added
133!
134! 1111 2013-03-08 23:54:10Z raasch
135! output of accelerator board information
136! ibc_p_b = 2 removed
137!
138! 1108 2013-03-05 07:03:32Z raasch
139! bugfix for r1106
140!
141! 1106 2013-03-04 05:31:38Z raasch
142! some format changes for coupled runs
143!
144! 1092 2013-02-02 11:24:22Z raasch
145! unused variables removed
146!
147! 1036 2012-10-22 13:43:42Z raasch
148! code put under GPL (PALM 3.9)
149!
150! 1031 2012-10-19 14:35:30Z raasch
151! output of netCDF data format modified
152!
153! 1015 2012-09-27 09:23:24Z raasch
154! output of Adjustment of mixing length to the Prandtl mixing length at first
155! grid point above ground removed
156!
157! 1003 2012-09-14 14:35:53Z raasch
158! output of information about equal/unequal subdomain size removed
159!
160! 1001 2012-09-13 14:08:46Z raasch
161! all actions concerning leapfrog- and upstream-spline-scheme removed
162!
163! 978 2012-08-09 08:28:32Z fricke
164! -km_damp_max, outflow_damping_width
165! +pt_damping_factor, pt_damping_width
166! +z0h
167!
168! 964 2012-07-26 09:14:24Z raasch
169! output of profil-related quantities removed
170!
171! 940 2012-07-09 14:31:00Z raasch
172! Output in case of simulations for pure neutral stratification (no pt-equation
173! solved)
174!
175! 927 2012-06-06 19:15:04Z raasch
176! output of masking_method for mg-solver
177!
178! 868 2012-03-28 12:21:07Z raasch
179! translation velocity in Galilean transformation changed to 0.6 * ug
180!
181! 833 2012-02-22 08:55:55Z maronga
182! Adjusted format for leaf area density
183!
184! 828 2012-02-21 12:00:36Z raasch
185! output of dissipation_classes + radius_classes
186!
187! 825 2012-02-19 03:03:44Z raasch
188! Output of cloud physics parameters/quantities complemented and restructured
189!
190! Revision 1.1  1997/08/11 06:17:20  raasch
191! Initial revision
192!
193!
194! Description:
195! ------------
196!> Writing a header with all important information about the actual run.
197!> This subroutine is called three times, two times at the beginning
198!> (writing information on files RUN_CONTROL and HEADER) and one time at the
199!> end of the run, then writing additional information about CPU-usage on file
200!> header.
201!-----------------------------------------------------------------------------!
202 SUBROUTINE header
203 
204
205    USE arrays_3d,                                                             &
206        ONLY:  pt_init, qsws, q_init, sa_init, shf, ug, vg, w_subs, zu, zw
207       
208    USE control_parameters
209       
210    USE cloud_parameters,                                                      &
211        ONLY:  cp, curvature_solution_effects, c_sedimentation,                &
212               limiter_sedimentation, l_v, nc_const, r_d, ventilation_effect
213       
214    USE cpulog,                                                                &
215        ONLY:  log_point_s
216       
217    USE dvrp_variables,                                                        &
218        ONLY:  use_seperate_pe_for_dvrp_output
219       
220    USE grid_variables,                                                        &
221        ONLY:  dx, dy
222       
223    USE indices,                                                               &
224        ONLY:  mg_loc_ind, nnx, nny, nnz, nx, ny, nxl_mg, nxr_mg, nyn_mg,      &
225               nys_mg, nzt, nzt_mg
226       
227    USE kinds
228   
229    USE land_surface_model_mod,                                                &
230        ONLY:  conserve_water_content, dewfall, land_surface, nzb_soil,        &
231               nzt_soil, root_fraction, soil_moisture, soil_temperature,       &
232               soil_type, soil_type_name, veg_type, veg_type_name, zs
233 
234    USE model_1d,                                                              &
235        ONLY:  damp_level_ind_1d, dt_pr_1d, dt_run_control_1d, end_time_1d
236       
237    USE particle_attributes,                                                   &
238        ONLY:  bc_par_b, bc_par_lr, bc_par_ns, bc_par_t, collision_kernel,     &
239               density_ratio, dissipation_classes, dt_min_part, dt_prel,       &
240               dt_write_particle_data, end_time_prel,                          &
241               maximum_number_of_tailpoints, maximum_tailpoint_age,            &
242               minimum_tailpoint_distance, number_of_particle_groups,          &
243               particle_advection, particle_advection_start,                   &
244               particles_per_point, pdx, pdy, pdz,  psb, psl, psn, psr, pss,   &
245               pst, radius, radius_classes, random_start_position,             &
246               seed_follows_topography,                                        &
247               total_number_of_particles, use_particle_tails,                  &
248               use_sgs_for_particles, total_number_of_tails,                   &
249               vertical_particle_advection, write_particle_statistics
250       
251    USE pegrid
252
253    USE plant_canopy_model_mod,                                                &
254        ONLY:  alpha_lad, beta_lad, calc_beta_lad_profile, canopy_drag_coeff,  &
255               canopy_mode, cthf, lad, lad_surface, lad_vertical_gradient,     &
256               lad_vertical_gradient_level, lad_vertical_gradient_level_ind,   &
257               lai_beta, leaf_scalar_exch_coeff, leaf_surface_conc, pch_index, &
258               plant_canopy
259
260    USE radiation_model_mod,                                                   &
261        ONLY:  albedo, albedo_type, albedo_type_name, constant_albedo,         &
262               day_init, dt_radiation, lambda, lw_radiation, net_radiation,    &
263               radiation, radiation_scheme, sw_radiation, time_utc_init
264   
265    USE spectrum,                                                              &
266        ONLY:  comp_spectra_level, data_output_sp, plot_spectra_level,         &
267               spectra_direction
268
269    IMPLICIT NONE
270
271    CHARACTER (LEN=1)  ::  prec                !<
272   
273    CHARACTER (LEN=2)  ::  do2d_mode           !<
274   
275    CHARACTER (LEN=5)  ::  section_chr         !<
276   
277    CHARACTER (LEN=10) ::  coor_chr            !<
278    CHARACTER (LEN=10) ::  host_chr            !<
279   
280    CHARACTER (LEN=16) ::  begin_chr           !<
281   
282    CHARACTER (LEN=26) ::  ver_rev             !<
283   
284    CHARACTER (LEN=40) ::  output_format       !<
285   
286    CHARACTER (LEN=70) ::  char1               !<
287    CHARACTER (LEN=70) ::  char2               !<
288    CHARACTER (LEN=70) ::  dopr_chr            !<
289    CHARACTER (LEN=70) ::  do2d_xy             !<
290    CHARACTER (LEN=70) ::  do2d_xz             !<
291    CHARACTER (LEN=70) ::  do2d_yz             !<
292    CHARACTER (LEN=70) ::  do3d_chr            !<
293    CHARACTER (LEN=70) ::  domask_chr          !<
294    CHARACTER (LEN=70) ::  run_classification  !<
295   
296    CHARACTER (LEN=85) ::  roben               !<
297    CHARACTER (LEN=85) ::  runten              !<
298   
299    CHARACTER (LEN=86) ::  coordinates         !<
300    CHARACTER (LEN=86) ::  gradients           !<
301    CHARACTER (LEN=86) ::  leaf_area_density   !<
302    CHARACTER (LEN=86) ::  roots               !<
303    CHARACTER (LEN=86) ::  slices              !<
304    CHARACTER (LEN=86) ::  temperatures        !<
305    CHARACTER (LEN=86) ::  ugcomponent         !<
306    CHARACTER (LEN=86) ::  vgcomponent         !<
307
308    CHARACTER (LEN=1), DIMENSION(1:3) ::  dir = (/ 'x', 'y', 'z' /)  !<
309
310    INTEGER(iwp) ::  av        !<
311    INTEGER(iwp) ::  bh        !<
312    INTEGER(iwp) ::  blx       !<
313    INTEGER(iwp) ::  bly       !<
314    INTEGER(iwp) ::  bxl       !<
315    INTEGER(iwp) ::  bxr       !<
316    INTEGER(iwp) ::  byn       !<
317    INTEGER(iwp) ::  bys       !<
318    INTEGER(iwp) ::  ch        !<
319    INTEGER(iwp) ::  count     !<
320    INTEGER(iwp) ::  cwx       !<
321    INTEGER(iwp) ::  cwy       !<
322    INTEGER(iwp) ::  cxl       !<
323    INTEGER(iwp) ::  cxr       !<
324    INTEGER(iwp) ::  cyn       !<
325    INTEGER(iwp) ::  cys       !<
326    INTEGER(iwp) ::  dim       !<
327    INTEGER(iwp) ::  i         !<
328    INTEGER(iwp) ::  io        !<
329    INTEGER(iwp) ::  j         !<
330    INTEGER(iwp) ::  k         !<
331    INTEGER(iwp) ::  l         !<
332    INTEGER(iwp) ::  ll        !<
333    INTEGER(iwp) ::  mpi_type  !<
334   
335    REAL(wp) ::  canopy_height                    !< canopy height (in m)
336    REAL(wp) ::  cpuseconds_per_simulated_second  !<
337
338!
339!-- Open the output file. At the end of the simulation, output is directed
340!-- to unit 19.
341    IF ( ( runnr == 0 .OR. force_print_header )  .AND. &
342         .NOT. simulated_time_at_begin /= simulated_time )  THEN
343       io = 15   !  header output on file RUN_CONTROL
344    ELSE
345       io = 19   !  header output on file HEADER
346    ENDIF
347    CALL check_open( io )
348
349!
350!-- At the end of the run, output file (HEADER) will be rewritten with
351!-- new information
352    IF ( io == 19 .AND. simulated_time_at_begin /= simulated_time ) REWIND( 19 )
353
354!
355!-- Determine kind of model run
356    IF ( TRIM( initializing_actions ) == 'read_restart_data' )  THEN
357       run_classification = '3D - restart run'
358    ELSEIF ( TRIM( initializing_actions ) == 'cyclic_fill' )  THEN
359       run_classification = '3D - run with cyclic fill of 3D - prerun data'
360    ELSEIF ( INDEX( initializing_actions, 'set_constant_profiles' ) /= 0 )  THEN
361       run_classification = '3D - run without 1D - prerun'
362    ELSEIF ( INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 )  THEN
363       run_classification = '3D - run with 1D - prerun'
364    ELSEIF ( INDEX( initializing_actions, 'by_user' ) /=0 )  THEN
365       run_classification = '3D - run initialized by user'
366    ELSE
367       message_string = ' unknown action(s): ' // TRIM( initializing_actions )
368       CALL message( 'header', 'PA0191', 0, 0, 0, 6, 0 )
369    ENDIF
370    IF ( ocean )  THEN
371       run_classification = 'ocean - ' // run_classification
372    ELSE
373       run_classification = 'atmosphere - ' // run_classification
374    ENDIF
375
376!
377!-- Run-identification, date, time, host
378    host_chr = host(1:10)
379    ver_rev = TRIM( version ) // '  ' // TRIM( revision )
380    WRITE ( io, 100 )  ver_rev, TRIM( run_classification )
381    IF ( TRIM( coupling_mode ) /= 'uncoupled' )  THEN
382#if defined( __mpi2 )
383       mpi_type = 2
384#else
385       mpi_type = 1
386#endif
387       WRITE ( io, 101 )  mpi_type, coupling_mode
388    ENDIF
389#if defined( __parallel )
390    IF ( coupling_start_time /= 0.0_wp )  THEN
391       IF ( coupling_start_time > simulated_time_at_begin )  THEN
392          WRITE ( io, 109 )
393       ELSE
394          WRITE ( io, 114 )
395       ENDIF
396    ENDIF
397#endif
398    IF ( ensemble_member_nr /= 0 )  THEN
399       WRITE ( io, 512 )  run_date, run_identifier, run_time, runnr,           &
400                       ADJUSTR( host_chr ), ensemble_member_nr
401    ELSE
402       WRITE ( io, 102 )  run_date, run_identifier, run_time, runnr,           &
403                       ADJUSTR( host_chr )
404    ENDIF
405#if defined( __parallel )
406    IF ( npex == -1  .AND.  npey == -1 )  THEN
407       char1 = 'calculated'
408    ELSE
409       char1 = 'predefined'
410    ENDIF
411    IF ( threads_per_task == 1 )  THEN
412       WRITE ( io, 103 )  numprocs, pdims(1), pdims(2), TRIM( char1 )
413    ELSE
414       WRITE ( io, 104 )  numprocs*threads_per_task, numprocs, &
415                          threads_per_task, pdims(1), pdims(2), TRIM( char1 )
416    ENDIF
417    IF ( num_acc_per_node /= 0 )  WRITE ( io, 117 )  num_acc_per_node   
418    IF ( ( host(1:3) == 'ibm'  .OR.  host(1:3) == 'nec'  .OR.    &
419           host(1:2) == 'lc'   .OR.  host(1:3) == 'dec' )  .AND. &
420         npex == -1  .AND.  pdims(2) == 1 )                      &
421    THEN
422       WRITE ( io, 106 )
423    ELSEIF ( pdims(2) == 1 )  THEN
424       WRITE ( io, 107 )  'x'
425    ELSEIF ( pdims(1) == 1 )  THEN
426       WRITE ( io, 107 )  'y'
427    ENDIF
428    IF ( use_seperate_pe_for_dvrp_output )  WRITE ( io, 105 )
429    IF ( numprocs /= maximum_parallel_io_streams )  THEN
430       WRITE ( io, 108 )  maximum_parallel_io_streams
431    ENDIF
432#else
433    IF ( num_acc_per_node /= 0 )  WRITE ( io, 120 )  num_acc_per_node
434#endif
435    WRITE ( io, 99 )
436
437!
438!-- Numerical schemes
439    WRITE ( io, 110 )
440    IF ( psolver(1:7) == 'poisfft' )  THEN
441       WRITE ( io, 111 )  TRIM( fft_method )
442       IF ( transpose_compute_overlap )  WRITE( io, 115 )
443    ELSEIF ( psolver == 'sor' )  THEN
444       WRITE ( io, 112 )  nsor_ini, nsor, omega_sor
445    ELSEIF ( psolver(1:9) == 'multigrid' )  THEN
446       WRITE ( io, 135 )  TRIM(psolver), cycle_mg, maximum_grid_level, ngsrb
447       IF ( mg_cycles == -1 )  THEN
448          WRITE ( io, 140 )  residual_limit
449       ELSE
450          WRITE ( io, 141 )  mg_cycles
451       ENDIF
452       IF ( mg_switch_to_pe0_level == 0 )  THEN
453          WRITE ( io, 136 )  nxr_mg(1)-nxl_mg(1)+1, nyn_mg(1)-nys_mg(1)+1, &
454                             nzt_mg(1)
455       ELSEIF (  mg_switch_to_pe0_level /= -1 )  THEN
456          WRITE ( io, 137 )  mg_switch_to_pe0_level,            &
457                             mg_loc_ind(2,0)-mg_loc_ind(1,0)+1, &
458                             mg_loc_ind(4,0)-mg_loc_ind(3,0)+1, &
459                             nzt_mg(mg_switch_to_pe0_level),    &
460                             nxr_mg(1)-nxl_mg(1)+1, nyn_mg(1)-nys_mg(1)+1, &
461                             nzt_mg(1)
462       ENDIF
463       IF ( masking_method )  WRITE ( io, 144 )
464    ENDIF
465    IF ( call_psolver_at_all_substeps  .AND. timestep_scheme(1:5) == 'runge' ) &
466    THEN
467       WRITE ( io, 142 )
468    ENDIF
469
470    IF ( momentum_advec == 'pw-scheme' )  THEN
471       WRITE ( io, 113 )
472    ELSEIF (momentum_advec == 'ws-scheme' )  THEN
473       WRITE ( io, 503 )
474    ENDIF
475    IF ( scalar_advec == 'pw-scheme' )  THEN
476       WRITE ( io, 116 )
477    ELSEIF ( scalar_advec == 'ws-scheme' )  THEN
478       WRITE ( io, 504 )
479    ELSEIF ( scalar_advec == 'ws-scheme-mono' )  THEN
480       WRITE ( io, 513 )
481    ELSE
482       WRITE ( io, 118 )
483    ENDIF
484
485    WRITE ( io, 139 )  TRIM( loop_optimization )
486
487    IF ( galilei_transformation )  THEN
488       IF ( use_ug_for_galilei_tr )  THEN
489          char1 = '0.6 * geostrophic wind'
490       ELSE
491          char1 = 'mean wind in model domain'
492       ENDIF
493       IF ( simulated_time_at_begin == simulated_time )  THEN
494          char2 = 'at the start of the run'
495       ELSE
496          char2 = 'at the end of the run'
497       ENDIF
498       WRITE ( io, 119 )  TRIM( char1 ), TRIM( char2 ),                        &
499                          advected_distance_x/1000.0_wp,                       &
500                          advected_distance_y/1000.0_wp
501    ENDIF
502    WRITE ( io, 122 )  timestep_scheme
503    IF ( use_upstream_for_tke )  WRITE ( io, 143 )
504    IF ( rayleigh_damping_factor /= 0.0_wp )  THEN
505       IF ( .NOT. ocean )  THEN
506          WRITE ( io, 123 )  'above', rayleigh_damping_height, &
507               rayleigh_damping_factor
508       ELSE
509          WRITE ( io, 123 )  'below', rayleigh_damping_height, &
510               rayleigh_damping_factor
511       ENDIF
512    ENDIF
513    IF ( neutral )  WRITE ( io, 131 )  pt_surface
514    IF ( humidity )  THEN
515       IF ( .NOT. cloud_physics )  THEN
516          WRITE ( io, 129 )
517       ELSE
518          WRITE ( io, 130 )
519       ENDIF
520    ENDIF
521    IF ( passive_scalar )  WRITE ( io, 134 )
522    IF ( conserve_volume_flow )  THEN
523       WRITE ( io, 150 )  conserve_volume_flow_mode
524       IF ( TRIM( conserve_volume_flow_mode ) == 'bulk_velocity' )  THEN
525          WRITE ( io, 151 )  u_bulk, v_bulk
526       ENDIF
527    ELSEIF ( dp_external )  THEN
528       IF ( dp_smooth )  THEN
529          WRITE ( io, 152 )  dpdxy, dp_level_b, ', vertically smoothed.'
530       ELSE
531          WRITE ( io, 152 )  dpdxy, dp_level_b, '.'
532       ENDIF
533    ENDIF
534    WRITE ( io, 99 )
535
536!
537!-- Runtime and timestep information
538    WRITE ( io, 200 )
539    IF ( .NOT. dt_fixed )  THEN
540       WRITE ( io, 201 )  dt_max, cfl_factor
541    ELSE
542       WRITE ( io, 202 )  dt
543    ENDIF
544    WRITE ( io, 203 )  simulated_time_at_begin, end_time
545
546    IF ( time_restart /= 9999999.9_wp  .AND. &
547         simulated_time_at_begin == simulated_time )  THEN
548       IF ( dt_restart == 9999999.9_wp )  THEN
549          WRITE ( io, 204 )  ' Restart at:       ',time_restart
550       ELSE
551          WRITE ( io, 205 )  ' Restart at:       ',time_restart, dt_restart
552       ENDIF
553    ENDIF
554
555    IF ( simulated_time_at_begin /= simulated_time )  THEN
556       i = MAX ( log_point_s(10)%counts, 1 )
557       IF ( ( simulated_time - simulated_time_at_begin ) == 0.0_wp )  THEN
558          cpuseconds_per_simulated_second = 0.0_wp
559       ELSE
560          cpuseconds_per_simulated_second = log_point_s(10)%sum / &
561                                            ( simulated_time -    &
562                                              simulated_time_at_begin )
563       ENDIF
564       WRITE ( io, 206 )  simulated_time, log_point_s(10)%sum,      &
565                          log_point_s(10)%sum / REAL( i, KIND=wp ), &
566                          cpuseconds_per_simulated_second
567       IF ( time_restart /= 9999999.9_wp  .AND.  time_restart < end_time )  THEN
568          IF ( dt_restart == 9999999.9_wp )  THEN
569             WRITE ( io, 204 )  ' Next restart at:     ',time_restart
570          ELSE
571             WRITE ( io, 205 )  ' Next restart at:     ',time_restart, dt_restart
572          ENDIF
573       ENDIF
574    ENDIF
575
576
577!
578!-- Start time for coupled runs, if independent precursor runs for atmosphere
579!-- and ocean are used or have been used. In this case, coupling_start_time
580!-- defines the time when the coupling is switched on.
581    IF ( coupling_start_time /= 0.0_wp )  THEN
582       WRITE ( io, 207 )  coupling_start_time
583    ENDIF
584
585!
586!-- Computational grid
587    IF ( .NOT. ocean )  THEN
588       WRITE ( io, 250 )  dx, dy, dz, (nx+1)*dx, (ny+1)*dy, zu(nzt+1)
589       IF ( dz_stretch_level_index < nzt+1 )  THEN
590          WRITE ( io, 252 )  dz_stretch_level, dz_stretch_level_index, &
591                             dz_stretch_factor, dz_max
592       ENDIF
593    ELSE
594       WRITE ( io, 250 )  dx, dy, dz, (nx+1)*dx, (ny+1)*dy, zu(0)
595       IF ( dz_stretch_level_index > 0 )  THEN
596          WRITE ( io, 252 )  dz_stretch_level, dz_stretch_level_index, &
597                             dz_stretch_factor, dz_max
598       ENDIF
599    ENDIF
600    WRITE ( io, 254 )  nx, ny, nzt+1, MIN( nnx, nx+1 ), MIN( nny, ny+1 ), &
601                       MIN( nnz+2, nzt+2 )
602    IF ( sloping_surface )  WRITE ( io, 260 )  alpha_surface
603
604!
605!-- Large scale forcing and nudging
606    WRITE ( io, 160 )
607    IF ( large_scale_forcing )  THEN
608       WRITE ( io, 162 )
609       WRITE ( io, 163 )
610
611       IF ( large_scale_subsidence )  THEN
612          IF ( .NOT. use_subsidence_tendencies )  THEN
613             WRITE ( io, 164 )
614          ELSE
615             WRITE ( io, 165 )
616          ENDIF
617       ENDIF
618
619       IF ( bc_pt_b == 'dirichlet' )  THEN
620          WRITE ( io, 180 )
621       ELSEIF ( bc_pt_b == 'neumann' )  THEN
622          WRITE ( io, 181 )
623       ENDIF
624
625       IF ( bc_q_b == 'dirichlet' )  THEN
626          WRITE ( io, 182 )
627       ELSEIF ( bc_q_b == 'neumann' )  THEN
628          WRITE ( io, 183 )
629       ENDIF
630
631       WRITE ( io, 167 )
632       IF ( nudging )  THEN
633          WRITE ( io, 170 )
634       ENDIF
635    ELSE
636       WRITE ( io, 161 )
637       WRITE ( io, 171 )
638    ENDIF
639    IF ( large_scale_subsidence )  THEN
640       WRITE ( io, 168 )
641       WRITE ( io, 169 )
642    ENDIF
643
644!
645!-- Profile for the large scale vertial velocity
646!-- Building output strings, starting with surface value
647    IF ( large_scale_subsidence )  THEN
648       temperatures = '   0.0'
649       gradients = '------'
650       slices = '     0'
651       coordinates = '   0.0'
652       i = 1
653       DO  WHILE ( subs_vertical_gradient_level_i(i) /= -9999 )
654
655          WRITE (coor_chr,'(E10.2,7X)')  &
656                                w_subs(subs_vertical_gradient_level_i(i))
657          temperatures = TRIM( temperatures ) // ' ' // TRIM( coor_chr )
658
659          WRITE (coor_chr,'(E10.2,7X)')  subs_vertical_gradient(i)
660          gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
661
662          WRITE (coor_chr,'(I10,7X)')  subs_vertical_gradient_level_i(i)
663          slices = TRIM( slices ) // ' ' // TRIM( coor_chr )
664
665          WRITE (coor_chr,'(F10.2,7X)')  subs_vertical_gradient_level(i)
666          coordinates = TRIM( coordinates ) // ' '  // TRIM( coor_chr )
667
668          IF ( i == 10 )  THEN
669             EXIT
670          ELSE
671             i = i + 1
672          ENDIF
673
674       ENDDO
675
676 
677       IF ( .NOT. large_scale_forcing )  THEN
678          WRITE ( io, 426 )  TRIM( coordinates ), TRIM( temperatures ), &
679                             TRIM( gradients ), TRIM( slices )
680       ENDIF
681
682
683    ENDIF
684
685!-- Profile of the geostrophic wind (component ug)
686!-- Building output strings
687    WRITE ( ugcomponent, '(F6.2)' )  ug_surface
688    gradients = '------'
689    slices = '     0'
690    coordinates = '   0.0'
691    i = 1
692    DO  WHILE ( ug_vertical_gradient_level_ind(i) /= -9999 )
693     
694       WRITE (coor_chr,'(F6.2,1X)')  ug(ug_vertical_gradient_level_ind(i))
695       ugcomponent = TRIM( ugcomponent ) // '  ' // TRIM( coor_chr )
696
697       WRITE (coor_chr,'(F6.2,1X)')  ug_vertical_gradient(i)
698       gradients = TRIM( gradients ) // '  ' // TRIM( coor_chr )
699
700       WRITE (coor_chr,'(I6,1X)')  ug_vertical_gradient_level_ind(i)
701       slices = TRIM( slices ) // '  ' // TRIM( coor_chr )
702
703       WRITE (coor_chr,'(F6.1,1X)')  ug_vertical_gradient_level(i)
704       coordinates = TRIM( coordinates ) // '  ' // TRIM( coor_chr )
705
706       IF ( i == 10 )  THEN
707          EXIT
708       ELSE
709          i = i + 1
710       ENDIF
711
712    ENDDO
713
714    IF ( .NOT. large_scale_forcing )  THEN
715       WRITE ( io, 423 )  TRIM( coordinates ), TRIM( ugcomponent ), &
716                          TRIM( gradients ), TRIM( slices )
717    ENDIF
718
719!-- Profile of the geostrophic wind (component vg)
720!-- Building output strings
721    WRITE ( vgcomponent, '(F6.2)' )  vg_surface
722    gradients = '------'
723    slices = '     0'
724    coordinates = '   0.0'
725    i = 1
726    DO  WHILE ( vg_vertical_gradient_level_ind(i) /= -9999 )
727
728       WRITE (coor_chr,'(F6.2,1X)')  vg(vg_vertical_gradient_level_ind(i))
729       vgcomponent = TRIM( vgcomponent ) // '  ' // TRIM( coor_chr )
730
731       WRITE (coor_chr,'(F6.2,1X)')  vg_vertical_gradient(i)
732       gradients = TRIM( gradients ) // '  ' // TRIM( coor_chr )
733
734       WRITE (coor_chr,'(I6,1X)')  vg_vertical_gradient_level_ind(i)
735       slices = TRIM( slices ) // '  ' // TRIM( coor_chr )
736
737       WRITE (coor_chr,'(F6.1,1X)')  vg_vertical_gradient_level(i)
738       coordinates = TRIM( coordinates ) // '  ' // TRIM( coor_chr )
739
740       IF ( i == 10 )  THEN
741          EXIT
742       ELSE
743          i = i + 1
744       ENDIF
745 
746    ENDDO
747
748    IF ( .NOT. large_scale_forcing )  THEN
749       WRITE ( io, 424 )  TRIM( coordinates ), TRIM( vgcomponent ), &
750                          TRIM( gradients ), TRIM( slices )
751    ENDIF
752
753!
754!-- Topography
755    WRITE ( io, 270 )  topography
756    SELECT CASE ( TRIM( topography ) )
757
758       CASE ( 'flat' )
759          ! no actions necessary
760
761       CASE ( 'single_building' )
762          blx = INT( building_length_x / dx )
763          bly = INT( building_length_y / dy )
764          bh  = MINLOC( ABS( zw - building_height ), 1 ) - 1
765          IF ( ABS( zw(bh  ) - building_height ) == &
766               ABS( zw(bh+1) - building_height )    )  bh = bh + 1
767
768          IF ( building_wall_left == 9999999.9_wp )  THEN
769             building_wall_left = ( nx + 1 - blx ) / 2 * dx
770          ENDIF
771          bxl = INT ( building_wall_left / dx + 0.5_wp )
772          bxr = bxl + blx
773
774          IF ( building_wall_south == 9999999.9_wp )  THEN
775             building_wall_south = ( ny + 1 - bly ) / 2 * dy
776          ENDIF
777          bys = INT ( building_wall_south / dy + 0.5_wp )
778          byn = bys + bly
779
780          WRITE ( io, 271 )  building_length_x, building_length_y, &
781                             building_height, bxl, bxr, bys, byn
782
783       CASE ( 'single_street_canyon' )
784          ch  = MINLOC( ABS( zw - canyon_height ), 1 ) - 1
785          IF ( ABS( zw(ch  ) - canyon_height ) == &
786               ABS( zw(ch+1) - canyon_height )    )  ch = ch + 1
787          IF ( canyon_width_x /= 9999999.9_wp )  THEN
788!
789!--          Street canyon in y direction
790             cwx = NINT( canyon_width_x / dx )
791             IF ( canyon_wall_left == 9999999.9_wp )  THEN
792                canyon_wall_left = ( nx + 1 - cwx ) / 2 * dx
793             ENDIF
794             cxl = NINT( canyon_wall_left / dx )
795             cxr = cxl + cwx
796             WRITE ( io, 272 )  'y', canyon_height, ch, 'u', cxl, cxr
797
798          ELSEIF ( canyon_width_y /= 9999999.9_wp )  THEN
799!
800!--          Street canyon in x direction
801             cwy = NINT( canyon_width_y / dy )
802             IF ( canyon_wall_south == 9999999.9_wp )  THEN
803                canyon_wall_south = ( ny + 1 - cwy ) / 2 * dy
804             ENDIF
805             cys = NINT( canyon_wall_south / dy )
806             cyn = cys + cwy
807             WRITE ( io, 272 )  'x', canyon_height, ch, 'v', cys, cyn
808          ENDIF
809
810    END SELECT
811
812    IF ( TRIM( topography ) /= 'flat' )  THEN
813       IF ( TRIM( topography_grid_convention ) == ' ' )  THEN
814          IF ( TRIM( topography ) == 'single_building' .OR.  &
815               TRIM( topography ) == 'single_street_canyon' )  THEN
816             WRITE ( io, 278 )
817          ELSEIF ( TRIM( topography ) == 'read_from_file' )  THEN
818             WRITE ( io, 279 )
819          ENDIF
820       ELSEIF ( TRIM( topography_grid_convention ) == 'cell_edge' )  THEN
821          WRITE ( io, 278 )
822       ELSEIF ( TRIM( topography_grid_convention ) == 'cell_center' )  THEN
823          WRITE ( io, 279 )
824       ENDIF
825    ENDIF
826
827    IF ( plant_canopy )  THEN
828   
829       canopy_height = pch_index * dz
830
831       WRITE ( io, 280 )  canopy_mode, canopy_height, pch_index,               &
832                          canopy_drag_coeff
833       IF ( passive_scalar )  THEN
834          WRITE ( io, 281 )  leaf_scalar_exch_coeff,                           &
835                             leaf_surface_conc
836       ENDIF
837
838!
839!--    Heat flux at the top of vegetation
840       WRITE ( io, 282 )  cthf
841
842!
843!--    Leaf area density profile, calculated either from given vertical
844!--    gradients or from beta probability density function.
845       IF (  .NOT.  calc_beta_lad_profile )  THEN
846
847!--       Building output strings, starting with surface value
848          WRITE ( leaf_area_density, '(F6.4)' )  lad_surface
849          gradients = '------'
850          slices = '     0'
851          coordinates = '   0.0'
852          i = 1
853          DO  WHILE ( i < 11  .AND.  lad_vertical_gradient_level_ind(i) /= -9999 )
854
855             WRITE (coor_chr,'(F7.2)')  lad(lad_vertical_gradient_level_ind(i))
856             leaf_area_density = TRIM( leaf_area_density ) // ' ' // TRIM( coor_chr )
857 
858             WRITE (coor_chr,'(F7.2)')  lad_vertical_gradient(i)
859             gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
860
861             WRITE (coor_chr,'(I7)')  lad_vertical_gradient_level_ind(i)
862             slices = TRIM( slices ) // ' ' // TRIM( coor_chr )
863
864             WRITE (coor_chr,'(F7.1)')  lad_vertical_gradient_level(i)
865             coordinates = TRIM( coordinates ) // ' '  // TRIM( coor_chr )
866
867             i = i + 1
868          ENDDO
869
870          WRITE ( io, 283 )  TRIM( coordinates ), TRIM( leaf_area_density ),              &
871                             TRIM( gradients ), TRIM( slices )
872
873       ELSE
874       
875          WRITE ( leaf_area_density, '(F6.4)' )  lad_surface
876          coordinates = '   0.0'
877         
878          DO  k = 1, pch_index
879
880             WRITE (coor_chr,'(F7.2)')  lad(k)
881             leaf_area_density = TRIM( leaf_area_density ) // ' ' // TRIM( coor_chr )
882 
883             WRITE (coor_chr,'(F7.1)')  zu(k)
884             coordinates = TRIM( coordinates ) // ' '  // TRIM( coor_chr )
885
886          ENDDO       
887
888          WRITE ( io, 284 ) TRIM( coordinates ), TRIM( leaf_area_density ), alpha_lad,    &
889                            beta_lad, lai_beta
890
891       ENDIF 
892
893    ENDIF
894
895
896    IF ( land_surface )  THEN
897
898       temperatures = ''
899       gradients    = '' ! use for humidity here
900       coordinates  = '' ! use for height
901       roots        = '' ! use for root fraction
902       slices       = '' ! use for index
903
904       i = 1
905       DO i = nzb_soil, nzt_soil
906          WRITE (coor_chr,'(F10.2,7X)') soil_temperature(i)
907          temperatures = TRIM( temperatures ) // ' ' // TRIM( coor_chr )
908
909          WRITE (coor_chr,'(F10.2,7X)') soil_moisture(i)
910          gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
911
912          WRITE (coor_chr,'(F10.2,7X)')  - zs(i)
913          coordinates = TRIM( coordinates ) // ' '  // TRIM( coor_chr )
914
915          WRITE (coor_chr,'(F10.2,7X)')  root_fraction(i)
916          roots = TRIM( roots ) // ' '  // TRIM( coor_chr )
917
918          WRITE (coor_chr,'(I10,7X)')  i
919          slices = TRIM( slices ) // ' '  // TRIM( coor_chr )
920
921
922       ENDDO
923
924!
925!--    Write land surface model header
926       WRITE( io, 419 )
927       IF ( conserve_water_content )  THEN
928          WRITE( io, 440 )
929       ELSE
930          WRITE( io, 441 )
931       ENDIF
932
933       IF ( dewfall )  THEN
934          WRITE( io, 442 )
935       ELSE
936          WRITE( io, 443 )
937       ENDIF
938
939       WRITE( io, 438 ) TRIM( veg_type_name(veg_type) ),                       &
940                        TRIM (soil_type_name(soil_type) )
941       WRITE( io, 439 ) TRIM( coordinates ), TRIM( temperatures ),             &
942                        TRIM( gradients ), TRIM( roots ), TRIM( slices )
943
944
945    ENDIF
946
947    IF ( radiation )  THEN
948!
949!--    Write radiation model header
950       WRITE( io, 444 )
951
952       IF ( radiation_scheme == "constant" )  THEN
953          WRITE( io, 445 ) net_radiation
954       ELSEIF ( radiation_scheme == "clear-sky" )  THEN
955          WRITE( io, 446 )
956       ELSEIF ( radiation_scheme == "rrtmg" )  THEN
957          WRITE( io, 447 )
958          IF ( .NOT. lw_radiation )  WRITE( io, 458 )
959          IF ( .NOT. sw_radiation )  WRITE( io, 459 )
960       ENDIF
961
962       IF ( albedo_type == 0 )  THEN
963          WRITE( io, 448 ) albedo
964       ELSE
965          WRITE( io, 456 ) TRIM( albedo_type_name(albedo_type) )
966       ENDIF
967       IF ( constant_albedo )  THEN
968          WRITE( io, 457 )
969       ENDIF
970       WRITE( io, 449 ) dt_radiation
971    ENDIF
972
973
974!
975!-- Boundary conditions
976    IF ( ibc_p_b == 0 )  THEN
977       runten = 'p(0)     = 0      |'
978    ELSEIF ( ibc_p_b == 1 )  THEN
979       runten = 'p(0)     = p(1)   |'
980    ENDIF
981    IF ( ibc_p_t == 0 )  THEN
982       roben  = 'p(nzt+1) = 0      |'
983    ELSE
984       roben  = 'p(nzt+1) = p(nzt) |'
985    ENDIF
986
987    IF ( ibc_uv_b == 0 )  THEN
988       runten = TRIM( runten ) // ' uv(0)     = -uv(1)                |'
989    ELSE
990       runten = TRIM( runten ) // ' uv(0)     = uv(1)                 |'
991    ENDIF
992    IF ( TRIM( bc_uv_t ) == 'dirichlet_0' )  THEN
993       roben  = TRIM( roben  ) // ' uv(nzt+1) = 0                     |'
994    ELSEIF ( ibc_uv_t == 0 )  THEN
995       roben  = TRIM( roben  ) // ' uv(nzt+1) = ug(nzt+1), vg(nzt+1)  |'
996    ELSE
997       roben  = TRIM( roben  ) // ' uv(nzt+1) = uv(nzt)               |'
998    ENDIF
999
1000    IF ( ibc_pt_b == 0 )  THEN
1001       IF ( land_surface )  THEN
1002          runten = TRIM( runten ) // ' pt(0)     = from soil model'
1003       ELSE
1004          runten = TRIM( runten ) // ' pt(0)     = pt_surface'
1005       ENDIF
1006    ELSEIF ( ibc_pt_b == 1 )  THEN
1007       runten = TRIM( runten ) // ' pt(0)     = pt(1)'
1008    ELSEIF ( ibc_pt_b == 2 )  THEN
1009       runten = TRIM( runten ) // ' pt(0)     = from coupled model'
1010    ENDIF
1011    IF ( ibc_pt_t == 0 )  THEN
1012       roben  = TRIM( roben  ) // ' pt(nzt+1) = pt_top'
1013    ELSEIF( ibc_pt_t == 1 )  THEN
1014       roben  = TRIM( roben  ) // ' pt(nzt+1) = pt(nzt)'
1015    ELSEIF( ibc_pt_t == 2 )  THEN
1016       roben  = TRIM( roben  ) // ' pt(nzt+1) = pt(nzt) + dpt/dz_ini'
1017
1018    ENDIF
1019
1020    WRITE ( io, 300 )  runten, roben
1021
1022    IF ( .NOT. constant_diffusion )  THEN
1023       IF ( ibc_e_b == 1 )  THEN
1024          runten = 'e(0)     = e(1)'
1025       ELSE
1026          runten = 'e(0)     = e(1) = (u*/0.1)**2'
1027       ENDIF
1028       roben = 'e(nzt+1) = e(nzt) = e(nzt-1)'
1029
1030       WRITE ( io, 301 )  'e', runten, roben       
1031
1032    ENDIF
1033
1034    IF ( ocean )  THEN
1035       runten = 'sa(0)    = sa(1)'
1036       IF ( ibc_sa_t == 0 )  THEN
1037          roben =  'sa(nzt+1) = sa_surface'
1038       ELSE
1039          roben =  'sa(nzt+1) = sa(nzt)'
1040       ENDIF
1041       WRITE ( io, 301 ) 'sa', runten, roben
1042    ENDIF
1043
1044    IF ( humidity )  THEN
1045       IF ( ibc_q_b == 0 )  THEN
1046          IF ( land_surface )  THEN
1047             runten = 'q(0)     = from soil model'
1048          ELSE
1049             runten = 'q(0)     = q_surface'
1050          ENDIF
1051
1052       ELSE
1053          runten = 'q(0)     = q(1)'
1054       ENDIF
1055       IF ( ibc_q_t == 0 )  THEN
1056          roben =  'q(nzt)   = q_top'
1057       ELSE
1058          roben =  'q(nzt)   = q(nzt-1) + dq/dz'
1059       ENDIF
1060       WRITE ( io, 301 ) 'q', runten, roben
1061    ENDIF
1062
1063    IF ( passive_scalar )  THEN
1064       IF ( ibc_q_b == 0 )  THEN
1065          runten = 's(0)     = s_surface'
1066       ELSE
1067          runten = 's(0)     = s(1)'
1068       ENDIF
1069       IF ( ibc_q_t == 0 )  THEN
1070          roben =  's(nzt)   = s_top'
1071       ELSE
1072          roben =  's(nzt)   = s(nzt-1) + ds/dz'
1073       ENDIF
1074       WRITE ( io, 301 ) 's', runten, roben
1075    ENDIF
1076
1077    IF ( use_surface_fluxes )  THEN
1078       WRITE ( io, 303 )
1079       IF ( constant_heatflux )  THEN
1080          IF ( large_scale_forcing .AND. lsf_surf )  THEN
1081             WRITE ( io, 306 )  shf(0,0)
1082          ELSE
1083             WRITE ( io, 306 )  surface_heatflux
1084          ENDIF
1085          IF ( random_heatflux )  WRITE ( io, 307 )
1086       ENDIF
1087       IF ( humidity  .AND.  constant_waterflux )  THEN
1088          IF ( large_scale_forcing .AND. lsf_surf )  THEN
1089             WRITE ( io, 311 ) qsws(0,0)
1090          ELSE
1091             WRITE ( io, 311 ) surface_waterflux
1092          ENDIF
1093       ENDIF
1094       IF ( passive_scalar  .AND.  constant_waterflux )  THEN
1095          WRITE ( io, 313 ) surface_waterflux
1096       ENDIF
1097    ENDIF
1098
1099    IF ( use_top_fluxes )  THEN
1100       WRITE ( io, 304 )
1101       IF ( coupling_mode == 'uncoupled' )  THEN
1102          WRITE ( io, 320 )  top_momentumflux_u, top_momentumflux_v
1103          IF ( constant_top_heatflux )  THEN
1104             WRITE ( io, 306 )  top_heatflux
1105          ENDIF
1106       ELSEIF ( coupling_mode == 'ocean_to_atmosphere' )  THEN
1107          WRITE ( io, 316 )
1108       ENDIF
1109       IF ( ocean  .AND.  constant_top_salinityflux )  THEN
1110          WRITE ( io, 309 )  top_salinityflux
1111       ENDIF
1112       IF ( humidity  .OR.  passive_scalar )  THEN
1113          WRITE ( io, 315 )
1114       ENDIF
1115    ENDIF
1116
1117    IF ( prandtl_layer )  THEN
1118       WRITE ( io, 305 )  (zu(1)-zu(0)), roughness_length, &
1119                          z0h_factor*roughness_length, kappa, &
1120                          rif_min, rif_max
1121       IF ( .NOT. constant_heatflux )  WRITE ( io, 308 )
1122       IF ( humidity  .AND.  .NOT. constant_waterflux )  THEN
1123          WRITE ( io, 312 )
1124       ENDIF
1125       IF ( passive_scalar  .AND.  .NOT. constant_waterflux )  THEN
1126          WRITE ( io, 314 )
1127       ENDIF
1128    ELSE
1129       IF ( INDEX(initializing_actions, 'set_1d-model_profiles') /= 0 )  THEN
1130          WRITE ( io, 310 )  rif_min, rif_max
1131       ENDIF
1132    ENDIF
1133
1134    WRITE ( io, 317 )  bc_lr, bc_ns
1135    IF ( .NOT. bc_lr_cyc  .OR.  .NOT. bc_ns_cyc )  THEN
1136       WRITE ( io, 318 )  use_cmax, pt_damping_width, pt_damping_factor       
1137       IF ( turbulent_inflow )  THEN
1138          IF ( .NOT. recycling_yshift ) THEN
1139             WRITE ( io, 319 )  recycling_width, recycling_plane, &
1140                                inflow_damping_height, inflow_damping_width
1141          ELSE
1142             WRITE ( io, 322 )  recycling_width, recycling_plane, &
1143                                inflow_damping_height, inflow_damping_width
1144          END IF
1145       ENDIF
1146    ENDIF
1147
1148!
1149!-- Initial Profiles
1150    WRITE ( io, 321 )
1151!
1152!-- Initial wind profiles
1153    IF ( u_profile(1) /= 9999999.9_wp )  WRITE ( io, 427 )
1154
1155!
1156!-- Initial temperature profile
1157!-- Building output strings, starting with surface temperature
1158    WRITE ( temperatures, '(F6.2)' )  pt_surface
1159    gradients = '------'
1160    slices = '     0'
1161    coordinates = '   0.0'
1162    i = 1
1163    DO  WHILE ( pt_vertical_gradient_level_ind(i) /= -9999 )
1164
1165       WRITE (coor_chr,'(F7.2)')  pt_init(pt_vertical_gradient_level_ind(i))
1166       temperatures = TRIM( temperatures ) // ' ' // TRIM( coor_chr )
1167
1168       WRITE (coor_chr,'(F7.2)')  pt_vertical_gradient(i)
1169       gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
1170
1171       WRITE (coor_chr,'(I7)')  pt_vertical_gradient_level_ind(i)
1172       slices = TRIM( slices ) // ' ' // TRIM( coor_chr )
1173
1174       WRITE (coor_chr,'(F7.1)')  pt_vertical_gradient_level(i)
1175       coordinates = TRIM( coordinates ) // ' '  // TRIM( coor_chr )
1176
1177       IF ( i == 10 )  THEN
1178          EXIT
1179       ELSE
1180          i = i + 1
1181       ENDIF
1182
1183    ENDDO
1184
1185    IF ( .NOT. nudging )  THEN
1186       WRITE ( io, 420 )  TRIM( coordinates ), TRIM( temperatures ), &
1187                          TRIM( gradients ), TRIM( slices )
1188    ELSE
1189       WRITE ( io, 428 ) 
1190    ENDIF
1191
1192!
1193!-- Initial humidity profile
1194!-- Building output strings, starting with surface humidity
1195    IF ( humidity  .OR.  passive_scalar )  THEN
1196       WRITE ( temperatures, '(E8.1)' )  q_surface
1197       gradients = '--------'
1198       slices = '       0'
1199       coordinates = '     0.0'
1200       i = 1
1201       DO  WHILE ( q_vertical_gradient_level_ind(i) /= -9999 )
1202         
1203          WRITE (coor_chr,'(E8.1,4X)')  q_init(q_vertical_gradient_level_ind(i))
1204          temperatures = TRIM( temperatures ) // '  ' // TRIM( coor_chr )
1205
1206          WRITE (coor_chr,'(E8.1,4X)')  q_vertical_gradient(i)
1207          gradients = TRIM( gradients ) // '  ' // TRIM( coor_chr )
1208         
1209          WRITE (coor_chr,'(I8,4X)')  q_vertical_gradient_level_ind(i)
1210          slices = TRIM( slices ) // '  ' // TRIM( coor_chr )
1211         
1212          WRITE (coor_chr,'(F8.1,4X)')  q_vertical_gradient_level(i)
1213          coordinates = TRIM( coordinates ) // '  '  // TRIM( coor_chr )
1214
1215          IF ( i == 10 )  THEN
1216             EXIT
1217          ELSE
1218             i = i + 1
1219          ENDIF
1220
1221       ENDDO
1222
1223       IF ( humidity )  THEN
1224          IF ( .NOT. nudging )  THEN
1225             WRITE ( io, 421 )  TRIM( coordinates ), TRIM( temperatures ), &
1226                                TRIM( gradients ), TRIM( slices )
1227          ENDIF
1228       ELSE
1229          WRITE ( io, 422 )  TRIM( coordinates ), TRIM( temperatures ), &
1230                             TRIM( gradients ), TRIM( slices )
1231       ENDIF
1232    ENDIF
1233
1234!
1235!-- Initial salinity profile
1236!-- Building output strings, starting with surface salinity
1237    IF ( ocean )  THEN
1238       WRITE ( temperatures, '(F6.2)' )  sa_surface
1239       gradients = '------'
1240       slices = '     0'
1241       coordinates = '   0.0'
1242       i = 1
1243       DO  WHILE ( sa_vertical_gradient_level_ind(i) /= -9999 )
1244
1245          WRITE (coor_chr,'(F7.2)')  sa_init(sa_vertical_gradient_level_ind(i))
1246          temperatures = TRIM( temperatures ) // ' ' // TRIM( coor_chr )
1247
1248          WRITE (coor_chr,'(F7.2)')  sa_vertical_gradient(i)
1249          gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
1250
1251          WRITE (coor_chr,'(I7)')  sa_vertical_gradient_level_ind(i)
1252          slices = TRIM( slices ) // ' ' // TRIM( coor_chr )
1253
1254          WRITE (coor_chr,'(F7.1)')  sa_vertical_gradient_level(i)
1255          coordinates = TRIM( coordinates ) // ' '  // TRIM( coor_chr )
1256
1257          IF ( i == 10 )  THEN
1258             EXIT
1259          ELSE
1260             i = i + 1
1261          ENDIF
1262
1263       ENDDO
1264
1265       WRITE ( io, 425 )  TRIM( coordinates ), TRIM( temperatures ), &
1266                          TRIM( gradients ), TRIM( slices )
1267    ENDIF
1268
1269
1270!
1271!-- Listing of 1D-profiles
1272    WRITE ( io, 325 )  dt_dopr_listing
1273    IF ( averaging_interval_pr /= 0.0_wp )  THEN
1274       WRITE ( io, 326 )  averaging_interval_pr, dt_averaging_input_pr
1275    ENDIF
1276
1277!
1278!-- DATA output
1279    WRITE ( io, 330 )
1280    IF ( averaging_interval_pr /= 0.0_wp )  THEN
1281       WRITE ( io, 326 )  averaging_interval_pr, dt_averaging_input_pr
1282    ENDIF
1283
1284!
1285!-- 1D-profiles
1286    dopr_chr = 'Profile:'
1287    IF ( dopr_n /= 0 )  THEN
1288       WRITE ( io, 331 )
1289
1290       output_format = ''
1291       output_format = output_format_netcdf
1292       WRITE ( io, 344 )  output_format
1293
1294       DO  i = 1, dopr_n
1295          dopr_chr = TRIM( dopr_chr ) // ' ' // TRIM( data_output_pr(i) ) // ','
1296          IF ( LEN_TRIM( dopr_chr ) >= 60 )  THEN
1297             WRITE ( io, 332 )  dopr_chr
1298             dopr_chr = '       :'
1299          ENDIF
1300       ENDDO
1301
1302       IF ( dopr_chr /= '' )  THEN
1303          WRITE ( io, 332 )  dopr_chr
1304       ENDIF
1305       WRITE ( io, 333 )  dt_dopr, averaging_interval_pr, dt_averaging_input_pr
1306       IF ( skip_time_dopr /= 0.0_wp )  WRITE ( io, 339 )  skip_time_dopr
1307    ENDIF
1308
1309!
1310!-- 2D-arrays
1311    DO  av = 0, 1
1312
1313       i = 1
1314       do2d_xy = ''
1315       do2d_xz = ''
1316       do2d_yz = ''
1317       DO  WHILE ( do2d(av,i) /= ' ' )
1318
1319          l = MAX( 2, LEN_TRIM( do2d(av,i) ) )
1320          do2d_mode = do2d(av,i)(l-1:l)
1321
1322          SELECT CASE ( do2d_mode )
1323             CASE ( 'xy' )
1324                ll = LEN_TRIM( do2d_xy )
1325                do2d_xy = do2d_xy(1:ll) // ' ' // do2d(av,i)(1:l-3) // ','
1326             CASE ( 'xz' )
1327                ll = LEN_TRIM( do2d_xz )
1328                do2d_xz = do2d_xz(1:ll) // ' ' // do2d(av,i)(1:l-3) // ','
1329             CASE ( 'yz' )
1330                ll = LEN_TRIM( do2d_yz )
1331                do2d_yz = do2d_yz(1:ll) // ' ' // do2d(av,i)(1:l-3) // ','
1332          END SELECT
1333
1334          i = i + 1
1335
1336       ENDDO
1337
1338       IF ( ( ( do2d_xy /= ''  .AND.  section(1,1) /= -9999 )  .OR.    &
1339              ( do2d_xz /= ''  .AND.  section(1,2) /= -9999 )  .OR.    &
1340              ( do2d_yz /= ''  .AND.  section(1,3) /= -9999 ) ) )  THEN
1341
1342          IF (  av == 0 )  THEN
1343             WRITE ( io, 334 )  ''
1344          ELSE
1345             WRITE ( io, 334 )  '(time-averaged)'
1346          ENDIF
1347
1348          IF ( do2d_at_begin )  THEN
1349             begin_chr = 'and at the start'
1350          ELSE
1351             begin_chr = ''
1352          ENDIF
1353
1354          output_format = ''
1355          output_format = output_format_netcdf
1356          WRITE ( io, 344 )  output_format
1357
1358          IF ( do2d_xy /= ''  .AND.  section(1,1) /= -9999 )  THEN
1359             i = 1
1360             slices = '/'
1361             coordinates = '/'
1362!
1363!--          Building strings with index and coordinate information of the
1364!--          slices
1365             DO  WHILE ( section(i,1) /= -9999 )
1366
1367                WRITE (section_chr,'(I5)')  section(i,1)
1368                section_chr = ADJUSTL( section_chr )
1369                slices = TRIM( slices ) // TRIM( section_chr ) // '/'
1370
1371                IF ( section(i,1) == -1 )  THEN
1372                   WRITE (coor_chr,'(F10.1)')  -1.0_wp
1373                ELSE
1374                   WRITE (coor_chr,'(F10.1)')  zu(section(i,1))
1375                ENDIF
1376                coor_chr = ADJUSTL( coor_chr )
1377                coordinates = TRIM( coordinates ) // TRIM( coor_chr ) // '/'
1378
1379                i = i + 1
1380             ENDDO
1381             IF ( av == 0 )  THEN
1382                WRITE ( io, 335 )  'XY', do2d_xy, dt_do2d_xy, &
1383                                   TRIM( begin_chr ), 'k', TRIM( slices ), &
1384                                   TRIM( coordinates )
1385                IF ( skip_time_do2d_xy /= 0.0_wp )  THEN
1386                   WRITE ( io, 339 )  skip_time_do2d_xy
1387                ENDIF
1388             ELSE
1389                WRITE ( io, 342 )  'XY', do2d_xy, dt_data_output_av, &
1390                                   TRIM( begin_chr ), averaging_interval, &
1391                                   dt_averaging_input, 'k', TRIM( slices ), &
1392                                   TRIM( coordinates )
1393                IF ( skip_time_data_output_av /= 0.0_wp )  THEN
1394                   WRITE ( io, 339 )  skip_time_data_output_av
1395                ENDIF
1396             ENDIF
1397             IF ( netcdf_data_format > 4 )  THEN
1398                WRITE ( io, 352 )  ntdim_2d_xy(av)
1399             ELSE
1400                WRITE ( io, 353 )
1401             ENDIF
1402          ENDIF
1403
1404          IF ( do2d_xz /= ''  .AND.  section(1,2) /= -9999 )  THEN
1405             i = 1
1406             slices = '/'
1407             coordinates = '/'
1408!
1409!--          Building strings with index and coordinate information of the
1410!--          slices
1411             DO  WHILE ( section(i,2) /= -9999 )
1412
1413                WRITE (section_chr,'(I5)')  section(i,2)
1414                section_chr = ADJUSTL( section_chr )
1415                slices = TRIM( slices ) // TRIM( section_chr ) // '/'
1416
1417                WRITE (coor_chr,'(F10.1)')  section(i,2) * dy
1418                coor_chr = ADJUSTL( coor_chr )
1419                coordinates = TRIM( coordinates ) // TRIM( coor_chr ) // '/'
1420
1421                i = i + 1
1422             ENDDO
1423             IF ( av == 0 )  THEN
1424                WRITE ( io, 335 )  'XZ', do2d_xz, dt_do2d_xz, &
1425                                   TRIM( begin_chr ), 'j', TRIM( slices ), &
1426                                   TRIM( coordinates )
1427                IF ( skip_time_do2d_xz /= 0.0_wp )  THEN
1428                   WRITE ( io, 339 )  skip_time_do2d_xz
1429                ENDIF
1430             ELSE
1431                WRITE ( io, 342 )  'XZ', do2d_xz, dt_data_output_av, &
1432                                   TRIM( begin_chr ), averaging_interval, &
1433                                   dt_averaging_input, 'j', TRIM( slices ), &
1434                                   TRIM( coordinates )
1435                IF ( skip_time_data_output_av /= 0.0_wp )  THEN
1436                   WRITE ( io, 339 )  skip_time_data_output_av
1437                ENDIF
1438             ENDIF
1439             IF ( netcdf_data_format > 4 )  THEN
1440                WRITE ( io, 352 )  ntdim_2d_xz(av)
1441             ELSE
1442                WRITE ( io, 353 )
1443             ENDIF
1444          ENDIF
1445
1446          IF ( do2d_yz /= ''  .AND.  section(1,3) /= -9999 )  THEN
1447             i = 1
1448             slices = '/'
1449             coordinates = '/'
1450!
1451!--          Building strings with index and coordinate information of the
1452!--          slices
1453             DO  WHILE ( section(i,3) /= -9999 )
1454
1455                WRITE (section_chr,'(I5)')  section(i,3)
1456                section_chr = ADJUSTL( section_chr )
1457                slices = TRIM( slices ) // TRIM( section_chr ) // '/'
1458
1459                WRITE (coor_chr,'(F10.1)')  section(i,3) * dx
1460                coor_chr = ADJUSTL( coor_chr )
1461                coordinates = TRIM( coordinates ) // TRIM( coor_chr ) // '/'
1462
1463                i = i + 1
1464             ENDDO
1465             IF ( av == 0 )  THEN
1466                WRITE ( io, 335 )  'YZ', do2d_yz, dt_do2d_yz, &
1467                                   TRIM( begin_chr ), 'i', TRIM( slices ), &
1468                                   TRIM( coordinates )
1469                IF ( skip_time_do2d_yz /= 0.0_wp )  THEN
1470                   WRITE ( io, 339 )  skip_time_do2d_yz
1471                ENDIF
1472             ELSE
1473                WRITE ( io, 342 )  'YZ', do2d_yz, dt_data_output_av, &
1474                                   TRIM( begin_chr ), averaging_interval, &
1475                                   dt_averaging_input, 'i', TRIM( slices ), &
1476                                   TRIM( coordinates )
1477                IF ( skip_time_data_output_av /= 0.0_wp )  THEN
1478                   WRITE ( io, 339 )  skip_time_data_output_av
1479                ENDIF
1480             ENDIF
1481             IF ( netcdf_data_format > 4 )  THEN
1482                WRITE ( io, 352 )  ntdim_2d_yz(av)
1483             ELSE
1484                WRITE ( io, 353 )
1485             ENDIF
1486          ENDIF
1487
1488       ENDIF
1489
1490    ENDDO
1491
1492!
1493!-- 3d-arrays
1494    DO  av = 0, 1
1495
1496       i = 1
1497       do3d_chr = ''
1498       DO  WHILE ( do3d(av,i) /= ' ' )
1499
1500          do3d_chr = TRIM( do3d_chr ) // ' ' // TRIM( do3d(av,i) ) // ','
1501          i = i + 1
1502
1503       ENDDO
1504
1505       IF ( do3d_chr /= '' )  THEN
1506          IF ( av == 0 )  THEN
1507             WRITE ( io, 336 )  ''
1508          ELSE
1509             WRITE ( io, 336 )  '(time-averaged)'
1510          ENDIF
1511
1512          output_format = output_format_netcdf
1513          WRITE ( io, 344 )  output_format
1514
1515          IF ( do3d_at_begin )  THEN
1516             begin_chr = 'and at the start'
1517          ELSE
1518             begin_chr = ''
1519          ENDIF
1520          IF ( av == 0 )  THEN
1521             WRITE ( io, 337 )  do3d_chr, dt_do3d, TRIM( begin_chr ), &
1522                                zu(nz_do3d), nz_do3d
1523          ELSE
1524             WRITE ( io, 343 )  do3d_chr, dt_data_output_av,           &
1525                                TRIM( begin_chr ), averaging_interval, &
1526                                dt_averaging_input, zu(nz_do3d), nz_do3d
1527          ENDIF
1528
1529          IF ( netcdf_data_format > 4 )  THEN
1530             WRITE ( io, 352 )  ntdim_3d(av)
1531          ELSE
1532             WRITE ( io, 353 )
1533          ENDIF
1534
1535          IF ( av == 0 )  THEN
1536             IF ( skip_time_do3d /= 0.0_wp )  THEN
1537                WRITE ( io, 339 )  skip_time_do3d
1538             ENDIF
1539          ELSE
1540             IF ( skip_time_data_output_av /= 0.0_wp )  THEN
1541                WRITE ( io, 339 )  skip_time_data_output_av
1542             ENDIF
1543          ENDIF
1544
1545       ENDIF
1546
1547    ENDDO
1548
1549!
1550!-- masked arrays
1551    IF ( masks > 0 )  WRITE ( io, 345 )  &
1552         mask_scale_x, mask_scale_y, mask_scale_z
1553    DO  mid = 1, masks
1554       DO  av = 0, 1
1555
1556          i = 1
1557          domask_chr = ''
1558          DO  WHILE ( domask(mid,av,i) /= ' ' )
1559             domask_chr = TRIM( domask_chr ) // ' ' //  &
1560                          TRIM( domask(mid,av,i) ) // ','
1561             i = i + 1
1562          ENDDO
1563
1564          IF ( domask_chr /= '' )  THEN
1565             IF ( av == 0 )  THEN
1566                WRITE ( io, 346 )  '', mid
1567             ELSE
1568                WRITE ( io, 346 )  ' (time-averaged)', mid
1569             ENDIF
1570
1571             output_format = output_format_netcdf
1572!--          Parallel output not implemented for mask data, hence
1573!--          output_format must be adjusted.
1574             IF ( netcdf_data_format == 5 ) output_format = 'netCDF4/HDF5'
1575             IF ( netcdf_data_format == 6 ) output_format = 'netCDF4/HDF5 classic'
1576             WRITE ( io, 344 )  output_format
1577
1578             IF ( av == 0 )  THEN
1579                WRITE ( io, 347 )  domask_chr, dt_domask(mid)
1580             ELSE
1581                WRITE ( io, 348 )  domask_chr, dt_data_output_av, &
1582                                   averaging_interval, dt_averaging_input
1583             ENDIF
1584
1585             IF ( av == 0 )  THEN
1586                IF ( skip_time_domask(mid) /= 0.0_wp )  THEN
1587                   WRITE ( io, 339 )  skip_time_domask(mid)
1588                ENDIF
1589             ELSE
1590                IF ( skip_time_data_output_av /= 0.0_wp )  THEN
1591                   WRITE ( io, 339 )  skip_time_data_output_av
1592                ENDIF
1593             ENDIF
1594!
1595!--          output locations
1596             DO  dim = 1, 3
1597                IF ( mask(mid,dim,1) >= 0.0_wp )  THEN
1598                   count = 0
1599                   DO  WHILE ( mask(mid,dim,count+1) >= 0.0_wp )
1600                      count = count + 1
1601                   ENDDO
1602                   WRITE ( io, 349 )  dir(dim), dir(dim), mid, dir(dim), &
1603                                      mask(mid,dim,:count)
1604                ELSEIF ( mask_loop(mid,dim,1) < 0.0_wp .AND.  &
1605                         mask_loop(mid,dim,2) < 0.0_wp .AND.  &
1606                         mask_loop(mid,dim,3) == 0.0_wp )  THEN
1607                   WRITE ( io, 350 )  dir(dim), dir(dim)
1608                ELSEIF ( mask_loop(mid,dim,3) == 0.0_wp )  THEN
1609                   WRITE ( io, 351 )  dir(dim), dir(dim), mid, dir(dim), &
1610                                      mask_loop(mid,dim,1:2)
1611                ELSE
1612                   WRITE ( io, 351 )  dir(dim), dir(dim), mid, dir(dim), &
1613                                      mask_loop(mid,dim,1:3)
1614                ENDIF
1615             ENDDO
1616          ENDIF
1617
1618       ENDDO
1619    ENDDO
1620
1621!
1622!-- Timeseries
1623    IF ( dt_dots /= 9999999.9_wp )  THEN
1624       WRITE ( io, 340 )
1625
1626       output_format = output_format_netcdf
1627       WRITE ( io, 344 )  output_format
1628       WRITE ( io, 341 )  dt_dots
1629    ENDIF
1630
1631#if defined( __dvrp_graphics )
1632!
1633!-- Dvrp-output
1634    IF ( dt_dvrp /= 9999999.9_wp )  THEN
1635       WRITE ( io, 360 )  dt_dvrp, TRIM( dvrp_output ), TRIM( dvrp_host ), &
1636                          TRIM( dvrp_username ), TRIM( dvrp_directory )
1637       i = 1
1638       l = 0
1639       m = 0
1640       DO WHILE ( mode_dvrp(i) /= ' ' )
1641          IF ( mode_dvrp(i)(1:10) == 'isosurface' )  THEN
1642             READ ( mode_dvrp(i), '(10X,I2)' )  j
1643             l = l + 1
1644             IF ( do3d(0,j) /= ' ' )  THEN
1645                WRITE ( io, 361 )  TRIM( do3d(0,j) ), threshold(l), &
1646                                   isosurface_color(:,l)
1647             ENDIF
1648          ELSEIF ( mode_dvrp(i)(1:6) == 'slicer' )  THEN
1649             READ ( mode_dvrp(i), '(6X,I2)' )  j
1650             m = m + 1
1651             IF ( do2d(0,j) /= ' ' )  THEN
1652                WRITE ( io, 362 )  TRIM( do2d(0,j) ), &
1653                                   slicer_range_limits_dvrp(:,m)
1654             ENDIF
1655          ELSEIF ( mode_dvrp(i)(1:9) == 'particles' )  THEN
1656             WRITE ( io, 363 )  dvrp_psize
1657             IF ( particle_dvrpsize /= 'none' )  THEN
1658                WRITE ( io, 364 )  'size', TRIM( particle_dvrpsize ), &
1659                                   dvrpsize_interval
1660             ENDIF
1661             IF ( particle_color /= 'none' )  THEN
1662                WRITE ( io, 364 )  'color', TRIM( particle_color ), &
1663                                   color_interval
1664             ENDIF
1665          ENDIF
1666          i = i + 1
1667       ENDDO
1668
1669       WRITE ( io, 365 )  groundplate_color, superelevation_x, &
1670                          superelevation_y, superelevation, clip_dvrp_l, &
1671                          clip_dvrp_r, clip_dvrp_s, clip_dvrp_n
1672
1673       IF ( TRIM( topography ) /= 'flat' )  THEN
1674          WRITE ( io, 366 )  topography_color
1675          IF ( cluster_size > 1 )  THEN
1676             WRITE ( io, 367 )  cluster_size
1677          ENDIF
1678       ENDIF
1679
1680    ENDIF
1681#endif
1682
1683#if defined( __spectra )
1684!
1685!-- Spectra output
1686    IF ( dt_dosp /= 9999999.9_wp )  THEN
1687       WRITE ( io, 370 )
1688
1689       output_format = output_format_netcdf
1690       WRITE ( io, 344 )  output_format
1691       WRITE ( io, 371 )  dt_dosp
1692       IF ( skip_time_dosp /= 0.0_wp )  WRITE ( io, 339 )  skip_time_dosp
1693       WRITE ( io, 372 )  ( data_output_sp(i), i = 1,10 ),     &
1694                          ( spectra_direction(i), i = 1,10 ),  &
1695                          ( comp_spectra_level(i), i = 1,100 ), &
1696                          ( plot_spectra_level(i), i = 1,100 ), &
1697                          averaging_interval_sp, dt_averaging_input_pr
1698    ENDIF
1699#endif
1700
1701    WRITE ( io, 99 )
1702
1703!
1704!-- Physical quantities
1705    WRITE ( io, 400 )
1706
1707!
1708!-- Geostrophic parameters
1709    IF ( radiation .AND. radiation_scheme /= 'constant' )  THEN
1710       WRITE ( io, 417 )  lambda
1711    ENDIF
1712    WRITE ( io, 410 )  phi, omega, f, fs
1713
1714!
1715!-- Other quantities
1716    WRITE ( io, 411 )  g
1717    IF ( radiation .AND. radiation_scheme /= 'constant' )  THEN
1718       WRITE ( io, 418 )  day_init, time_utc_init
1719    ENDIF
1720
1721    WRITE ( io, 412 )  TRIM( reference_state )
1722    IF ( use_single_reference_value )  THEN
1723       IF ( ocean )  THEN
1724          WRITE ( io, 413 )  prho_reference
1725       ELSE
1726          WRITE ( io, 414 )  pt_reference
1727       ENDIF
1728    ENDIF
1729
1730!
1731!-- Cloud physics parameters
1732    IF ( cloud_physics )  THEN
1733       WRITE ( io, 415 )
1734       WRITE ( io, 416 ) surface_pressure, r_d, rho_surface, cp, l_v
1735       IF ( icloud_scheme == 0 )  THEN
1736          WRITE ( io, 510 ) 1.0E-6_wp * nc_const
1737          IF ( precipitation )  WRITE ( io, 511 ) c_sedimentation
1738       ENDIF
1739    ENDIF
1740
1741!
1742!-- Cloud physcis parameters / quantities / numerical methods
1743    WRITE ( io, 430 )
1744    IF ( humidity .AND. .NOT. cloud_physics .AND. .NOT. cloud_droplets)  THEN
1745       WRITE ( io, 431 )
1746    ELSEIF ( humidity  .AND.  cloud_physics )  THEN
1747       WRITE ( io, 432 )
1748       IF ( cloud_top_radiation )  WRITE ( io, 132 )
1749       IF ( icloud_scheme == 1 )  THEN
1750          IF ( precipitation )  WRITE ( io, 133 )
1751       ELSEIF ( icloud_scheme == 0 )  THEN
1752          IF ( drizzle )  WRITE ( io, 506 )
1753          IF ( precipitation )  THEN
1754             WRITE ( io, 505 )
1755             IF ( turbulence )  WRITE ( io, 507 )
1756             IF ( ventilation_effect )  WRITE ( io, 508 )
1757             IF ( limiter_sedimentation )  WRITE ( io, 509 )
1758          ENDIF
1759       ENDIF
1760    ELSEIF ( humidity  .AND.  cloud_droplets )  THEN
1761       WRITE ( io, 433 )
1762       IF ( curvature_solution_effects )  WRITE ( io, 434 )
1763       IF ( collision_kernel /= 'none' )  THEN
1764          WRITE ( io, 435 )  TRIM( collision_kernel )
1765          IF ( collision_kernel(6:9) == 'fast' )  THEN
1766             WRITE ( io, 436 )  radius_classes, dissipation_classes
1767          ENDIF
1768       ELSE
1769          WRITE ( io, 437 )
1770       ENDIF
1771    ENDIF
1772
1773!
1774!-- LES / turbulence parameters
1775    WRITE ( io, 450 )
1776
1777!--
1778! ... LES-constants used must still be added here
1779!--
1780    IF ( constant_diffusion )  THEN
1781       WRITE ( io, 451 )  km_constant, km_constant/prandtl_number, &
1782                          prandtl_number
1783    ENDIF
1784    IF ( .NOT. constant_diffusion)  THEN
1785       IF ( e_init > 0.0_wp )  WRITE ( io, 455 )  e_init
1786       IF ( e_min > 0.0_wp )  WRITE ( io, 454 )  e_min
1787       IF ( wall_adjustment )  WRITE ( io, 453 )  wall_adjustment_factor
1788    ENDIF
1789
1790!
1791!-- Special actions during the run
1792    WRITE ( io, 470 )
1793    IF ( create_disturbances )  THEN
1794       WRITE ( io, 471 )  dt_disturb, disturbance_amplitude,                   &
1795                          zu(disturbance_level_ind_b), disturbance_level_ind_b,&
1796                          zu(disturbance_level_ind_t), disturbance_level_ind_t
1797       IF ( .NOT. bc_lr_cyc  .OR.  .NOT. bc_ns_cyc )  THEN
1798          WRITE ( io, 472 )  inflow_disturbance_begin, inflow_disturbance_end
1799       ELSE
1800          WRITE ( io, 473 )  disturbance_energy_limit
1801       ENDIF
1802       WRITE ( io, 474 )  TRIM( random_generator )
1803    ENDIF
1804    IF ( pt_surface_initial_change /= 0.0_wp )  THEN
1805       WRITE ( io, 475 )  pt_surface_initial_change
1806    ENDIF
1807    IF ( humidity  .AND.  q_surface_initial_change /= 0.0_wp )  THEN
1808       WRITE ( io, 476 )  q_surface_initial_change       
1809    ENDIF
1810    IF ( passive_scalar  .AND.  q_surface_initial_change /= 0.0_wp )  THEN
1811       WRITE ( io, 477 )  q_surface_initial_change       
1812    ENDIF
1813
1814    IF ( particle_advection )  THEN
1815!
1816!--    Particle attributes
1817       WRITE ( io, 480 )  particle_advection_start, dt_prel, bc_par_lr, &
1818                          bc_par_ns, bc_par_b, bc_par_t, particle_maximum_age, &
1819                          end_time_prel
1820       IF ( use_sgs_for_particles )  WRITE ( io, 488 )  dt_min_part
1821       IF ( random_start_position )  WRITE ( io, 481 )
1822       IF ( seed_follows_topography )  WRITE ( io, 496 )
1823       IF ( particles_per_point > 1 )  WRITE ( io, 489 )  particles_per_point
1824       WRITE ( io, 495 )  total_number_of_particles
1825       IF ( use_particle_tails  .AND.  maximum_number_of_tailpoints /= 0 )  THEN
1826          WRITE ( io, 483 )  maximum_number_of_tailpoints
1827          IF ( minimum_tailpoint_distance /= 0 )  THEN
1828             WRITE ( io, 484 )  total_number_of_tails,      &
1829                                minimum_tailpoint_distance, &
1830                                maximum_tailpoint_age
1831          ENDIF
1832       ENDIF
1833       IF ( dt_write_particle_data /= 9999999.9_wp )  THEN
1834          WRITE ( io, 485 )  dt_write_particle_data
1835          IF ( netcdf_data_format > 1 )  THEN
1836             output_format = 'netcdf (64 bit offset) and binary'
1837          ELSE
1838             output_format = 'netcdf and binary'
1839          ENDIF
1840          WRITE ( io, 344 )  output_format
1841       ENDIF
1842       IF ( dt_dopts /= 9999999.9_wp )  WRITE ( io, 494 )  dt_dopts
1843       IF ( write_particle_statistics )  WRITE ( io, 486 )
1844
1845       WRITE ( io, 487 )  number_of_particle_groups
1846
1847       DO  i = 1, number_of_particle_groups
1848          IF ( i == 1  .AND.  density_ratio(i) == 9999999.9_wp )  THEN
1849             WRITE ( io, 490 )  i, 0.0_wp
1850             WRITE ( io, 492 )
1851          ELSE
1852             WRITE ( io, 490 )  i, radius(i)
1853             IF ( density_ratio(i) /= 0.0_wp )  THEN
1854                WRITE ( io, 491 )  density_ratio(i)
1855             ELSE
1856                WRITE ( io, 492 )
1857             ENDIF
1858          ENDIF
1859          WRITE ( io, 493 )  psl(i), psr(i), pss(i), psn(i), psb(i), pst(i), &
1860                             pdx(i), pdy(i), pdz(i)
1861          IF ( .NOT. vertical_particle_advection(i) )  WRITE ( io, 482 )
1862       ENDDO
1863
1864    ENDIF
1865
1866
1867!
1868!-- Parameters of 1D-model
1869    IF ( INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 )  THEN
1870       WRITE ( io, 500 )  end_time_1d, dt_run_control_1d, dt_pr_1d, &
1871                          mixing_length_1d, dissipation_1d
1872       IF ( damp_level_ind_1d /= nzt+1 )  THEN
1873          WRITE ( io, 502 )  zu(damp_level_ind_1d), damp_level_ind_1d
1874       ENDIF
1875    ENDIF
1876
1877!
1878!-- User-defined information
1879    CALL user_header( io )
1880
1881    WRITE ( io, 99 )
1882
1883!
1884!-- Write buffer contents to disc immediately
1885    CALL local_flush( io )
1886
1887!
1888!-- Here the FORMATs start
1889
1890 99 FORMAT (1X,78('-'))
1891100 FORMAT (/1X,'******************************',4X,44('-')/        &
1892            1X,'* ',A,' *',4X,A/                               &
1893            1X,'******************************',4X,44('-'))
1894101 FORMAT (35X,'coupled run using MPI-',I1,': ',A/ &
1895            35X,42('-'))
1896102 FORMAT (/' Date:                 ',A8,4X,'Run:       ',A20/      &
1897            ' Time:                 ',A8,4X,'Run-No.:   ',I2.2/     &
1898            ' Run on host:        ',A10)
1899#if defined( __parallel )
1900103 FORMAT (' Number of PEs:',10X,I6,4X,'Processor grid (x,y): (',I4,',',I4, &
1901              ')',1X,A)
1902104 FORMAT (' Number of PEs:',10X,I6,4X,'Tasks:',I4,'   threads per task:',I4/ &
1903              35X,'Processor grid (x,y): (',I4,',',I4,')',1X,A)
1904105 FORMAT (35X,'One additional PE is used to handle'/37X,'the dvrp output!')
1905106 FORMAT (35X,'A 1d-decomposition along x is forced'/ &
1906            35X,'because the job is running on an SMP-cluster')
1907107 FORMAT (35X,'A 1d-decomposition along ',A,' is used')
1908108 FORMAT (35X,'Max. # of parallel I/O streams is ',I5)
1909109 FORMAT (35X,'Precursor run for coupled atmos-ocean run'/ &
1910            35X,42('-'))
1911114 FORMAT (35X,'Coupled atmosphere-ocean run following'/ &
1912            35X,'independent precursor runs'/             &
1913            35X,42('-'))
1914117 FORMAT (' Accelerator boards / node:  ',I2)
1915#endif
1916110 FORMAT (/' Numerical Schemes:'/ &
1917             ' -----------------'/)
1918111 FORMAT (' --> Solve perturbation pressure via FFT using ',A,' routines')
1919112 FORMAT (' --> Solve perturbation pressure via SOR-Red/Black-Schema'/ &
1920            '     Iterations (initial/other): ',I3,'/',I3,'  omega = ',F5.3)
1921113 FORMAT (' --> Momentum advection via Piascek-Williams-Scheme (Form C3)', &
1922                  ' or Upstream')
1923115 FORMAT ('     FFT and transpositions are overlapping')
1924116 FORMAT (' --> Scalar advection via Piascek-Williams-Scheme (Form C3)', &
1925                  ' or Upstream')
1926118 FORMAT (' --> Scalar advection via Bott-Chlond-Scheme')
1927119 FORMAT (' --> Galilei-Transform applied to horizontal advection:'/ &
1928            '     translation velocity = ',A/ &
1929            '     distance advected ',A,':  ',F8.3,' km(x)  ',F8.3,' km(y)')
1930120 FORMAT (' Accelerator boards: ',8X,I2)
1931122 FORMAT (' --> Time differencing scheme: ',A)
1932123 FORMAT (' --> Rayleigh-Damping active, starts ',A,' z = ',F8.2,' m'/ &
1933            '     maximum damping coefficient: ',F5.3, ' 1/s')
1934129 FORMAT (' --> Additional prognostic equation for the specific humidity')
1935130 FORMAT (' --> Additional prognostic equation for the total water content')
1936131 FORMAT (' --> No pt-equation solved. Neutral stratification with pt = ', &
1937                  F6.2, ' K assumed')
1938132 FORMAT ('     Parameterization of long-wave radiation processes via'/ &
1939            '     effective emissivity scheme')
1940133 FORMAT ('     Precipitation parameterization via Kessler-Scheme')
1941134 FORMAT (' --> Additional prognostic equation for a passive scalar')
1942135 FORMAT (' --> Solve perturbation pressure via ',A,' method (', &
1943                  A,'-cycle)'/ &
1944            '     number of grid levels:                   ',I2/ &
1945            '     Gauss-Seidel red/black iterations:       ',I2)
1946136 FORMAT ('     gridpoints of coarsest subdomain (x,y,z): (',I3,',',I3,',', &
1947                  I3,')')
1948137 FORMAT ('     level data gathered on PE0 at level:     ',I2/ &
1949            '     gridpoints of coarsest subdomain (x,y,z): (',I3,',',I3,',', &
1950                  I3,')'/ &
1951            '     gridpoints of coarsest domain (x,y,z):    (',I3,',',I3,',', &
1952                  I3,')')
1953139 FORMAT (' --> Loop optimization method: ',A)
1954140 FORMAT ('     maximum residual allowed:                ',E10.3)
1955141 FORMAT ('     fixed number of multigrid cycles:        ',I4)
1956142 FORMAT ('     perturbation pressure is calculated at every Runge-Kutta ', &
1957                  'step')
1958143 FORMAT ('     Euler/upstream scheme is used for the SGS turbulent ', &
1959                  'kinetic energy')
1960144 FORMAT ('     masking method is used')
1961150 FORMAT (' --> Volume flow at the right and north boundary will be ', &
1962                  'conserved'/ &
1963            '     using the ',A,' mode')
1964151 FORMAT ('     with u_bulk = ',F7.3,' m/s and v_bulk = ',F7.3,' m/s')
1965152 FORMAT (' --> External pressure gradient directly prescribed by the user:',&
1966           /'     ',2(1X,E12.5),'Pa/m in x/y direction', &
1967           /'     starting from dp_level_b =', F8.3, 'm', A /)
1968160 FORMAT (//' Large scale forcing and nudging:'/ &
1969              ' -------------------------------'/)
1970161 FORMAT (' --> No large scale forcing from external is used (default) ')
1971162 FORMAT (' --> Large scale forcing from external file LSF_DATA is used: ')
1972163 FORMAT ('     - large scale advection tendencies ')
1973164 FORMAT ('     - large scale subsidence velocity w_subs ')
1974165 FORMAT ('     - large scale subsidence tendencies ')
1975167 FORMAT ('     - and geostrophic wind components ug and vg')
1976168 FORMAT (' --> Large-scale vertical motion is used in the ', &
1977                  'prognostic equation(s) for')
1978169 FORMAT ('     the scalar(s) only')
1979170 FORMAT (' --> Nudging is used')
1980171 FORMAT (' --> No nudging is used (default) ')
1981180 FORMAT ('     - prescribed surface values for temperature')
1982181 FORMAT ('     - prescribed surface fluxes for temperature')
1983182 FORMAT ('     - prescribed surface values for humidity')
1984183 FORMAT ('     - prescribed surface fluxes for humidity')
1985200 FORMAT (//' Run time and time step information:'/ &
1986             ' ----------------------------------'/)
1987201 FORMAT ( ' Timestep:             variable     maximum value: ',F6.3,' s', &
1988             '    CFL-factor: ',F4.2)
1989202 FORMAT ( ' Timestep:          dt = ',F6.3,' s'/)
1990203 FORMAT ( ' Start time:          ',F9.3,' s'/ &
1991             ' End time:            ',F9.3,' s')
1992204 FORMAT ( A,F9.3,' s')
1993205 FORMAT ( A,F9.3,' s',5X,'restart every',17X,F9.3,' s')
1994206 FORMAT (/' Time reached:        ',F9.3,' s'/ &
1995             ' CPU-time used:       ',F9.3,' s     per timestep:               ', &
1996               '  ',F9.3,' s'/                                                    &
1997             '                                      per second of simulated tim', &
1998               'e: ',F9.3,' s')
1999207 FORMAT ( ' Coupling start time: ',F9.3,' s')
2000250 FORMAT (//' Computational grid and domain size:'/ &
2001              ' ----------------------------------'// &
2002              ' Grid length:      dx =    ',F7.3,' m    dy =    ',F7.3, &
2003              ' m    dz =    ',F7.3,' m'/ &
2004              ' Domain size:       x = ',F10.3,' m     y = ',F10.3, &
2005              ' m  z(u) = ',F10.3,' m'/)
2006252 FORMAT (' dz constant up to ',F10.3,' m (k=',I4,'), then stretched by', &
2007              ' factor: ',F5.3/ &
2008            ' maximum dz not to be exceeded is dz_max = ',F10.3,' m'/)
2009254 FORMAT (' Number of gridpoints (x,y,z):  (0:',I4,', 0:',I4,', 0:',I4,')'/ &
2010            ' Subdomain size (x,y,z):        (  ',I4,',   ',I4,',   ',I4,')'/)
2011260 FORMAT (/' The model has a slope in x-direction. Inclination angle: ',F6.2,&
2012             ' degrees')
2013270 FORMAT (//' Topography information:'/ &
2014              ' ----------------------'// &
2015              1X,'Topography: ',A)
2016271 FORMAT (  ' Building size (x/y/z) in m: ',F5.1,' / ',F5.1,' / ',F5.1/ &
2017              ' Horizontal index bounds (l/r/s/n): ',I4,' / ',I4,' / ',I4, &
2018                ' / ',I4)
2019272 FORMAT (  ' Single quasi-2D street canyon of infinite length in ',A, &
2020              ' direction' / &
2021              ' Canyon height: ', F6.2, 'm, ch = ', I4, '.'      / &
2022              ' Canyon position (',A,'-walls): cxl = ', I4,', cxr = ', I4, '.')
2023278 FORMAT (' Topography grid definition convention:'/ &
2024            ' cell edge (staggered grid points'/  &
2025            ' (u in x-direction, v in y-direction))' /)
2026279 FORMAT (' Topography grid definition convention:'/ &
2027            ' cell center (scalar grid points)' /)
2028280 FORMAT (//' Vegetation canopy (drag) model:'/ &
2029              ' ------------------------------'// &
2030              ' Canopy mode: ', A / &
2031              ' Canopy height: ',F6.2,'m (',I4,' grid points)' / &
2032              ' Leaf drag coefficient: ',F6.2 /)
2033281 FORMAT (/ ' Scalar exchange coefficient: ',F6.2 / &
2034              ' Scalar concentration at leaf surfaces in kg/m**3: ',F6.2 /)
2035282 FORMAT (' Predefined constant heatflux at the top of the vegetation: ',F6.2,' K m/s')
2036283 FORMAT (/ ' Characteristic levels of the leaf area density:'// &
2037              ' Height:              ',A,'  m'/ &
2038              ' Leaf area density:   ',A,'  m**2/m**3'/ &
2039              ' Gradient:            ',A,'  m**2/m**4'/ &
2040              ' Gridpoint:           ',A)
2041284 FORMAT (//' Characteristic levels of the leaf area density and coefficients:'// &
2042              ' Height:              ',A,'  m'/ &
2043              ' Leaf area density:   ',A,'  m**2/m**3'/ &
2044              ' Coefficient alpha: ',F6.2 / &
2045              ' Coefficient beta: ',F6.2 / &
2046              ' Leaf area index: ',F6.2,'  m**2/m**2' /)
2047               
2048300 FORMAT (//' Boundary conditions:'/ &
2049             ' -------------------'// &
2050             '                     p                    uv             ', &
2051             '                     pt'// &
2052             ' B. bound.: ',A/ &
2053             ' T. bound.: ',A)
2054301 FORMAT (/'                     ',A// &
2055             ' B. bound.: ',A/ &
2056             ' T. bound.: ',A)
2057303 FORMAT (/' Bottom surface fluxes are used in diffusion terms at k=1')
2058304 FORMAT (/' Top surface fluxes are used in diffusion terms at k=nzt')
2059305 FORMAT (//'    Prandtl-Layer between bottom surface and first ', &
2060               'computational u,v-level:'// &
2061             '       zp = ',F6.2,' m   z0 = ',F6.4,' m   z0h = ',F7.5,&
2062             ' m   kappa = ',F4.2/ &
2063             '       Rif value range:   ',F6.2,' <= rif <=',F6.2)
2064306 FORMAT ('       Predefined constant heatflux:   ',F9.6,' K m/s')
2065307 FORMAT ('       Heatflux has a random normal distribution')
2066308 FORMAT ('       Predefined surface temperature')
2067309 FORMAT ('       Predefined constant salinityflux:   ',F9.6,' psu m/s')
2068310 FORMAT (//'    1D-Model:'// &
2069             '       Rif value range:   ',F6.2,' <= rif <=',F6.2)
2070311 FORMAT ('       Predefined constant humidity flux: ',E10.3,' m/s')
2071312 FORMAT ('       Predefined surface humidity')
2072313 FORMAT ('       Predefined constant scalar flux: ',E10.3,' kg/(m**2 s)')
2073314 FORMAT ('       Predefined scalar value at the surface')
2074315 FORMAT ('       Humidity / scalar flux at top surface is 0.0')
2075316 FORMAT ('       Sensible heatflux and momentum flux from coupled ', &
2076                    'atmosphere model')
2077317 FORMAT (//' Lateral boundaries:'/ &
2078            '       left/right:  ',A/    &
2079            '       north/south: ',A)
2080318 FORMAT (/'       use_cmax: ',L1 / &
2081            '       pt damping layer width = ',F8.2,' m, pt ', &
2082                    'damping factor = ',F6.4)
2083319 FORMAT ('       turbulence recycling at inflow switched on'/ &
2084            '       width of recycling domain: ',F7.1,' m   grid index: ',I4/ &
2085            '       inflow damping height: ',F6.1,' m   width: ',F6.1,' m')
2086320 FORMAT ('       Predefined constant momentumflux:  u: ',F9.6,' m**2/s**2'/ &
2087            '                                          v: ',F9.6,' m**2/s**2')
2088321 FORMAT (//' Initial profiles:'/ &
2089              ' ----------------')
2090322 FORMAT ('       turbulence recycling at inflow switched on'/ &
2091            '       y shift of the recycled inflow turbulence switched on'/ &
2092            '       width of recycling domain: ',F7.1,' m   grid index: ',I4/ &
2093            '       inflow damping height: ',F6.1,' m   width: ',F6.1,' m'/)
2094325 FORMAT (//' List output:'/ &
2095             ' -----------'//  &
2096            '    1D-Profiles:'/    &
2097            '       Output every             ',F8.2,' s')
2098326 FORMAT ('       Time averaged over       ',F8.2,' s'/ &
2099            '       Averaging input every    ',F8.2,' s')
2100330 FORMAT (//' Data output:'/ &
2101             ' -----------'/)
2102331 FORMAT (/'    1D-Profiles:')
2103332 FORMAT (/'       ',A)
2104333 FORMAT ('       Output every             ',F8.2,' s',/ &
2105            '       Time averaged over       ',F8.2,' s'/ &
2106            '       Averaging input every    ',F8.2,' s')
2107334 FORMAT (/'    2D-Arrays',A,':')
2108335 FORMAT (/'       ',A2,'-cross-section  Arrays: ',A/ &
2109            '       Output every             ',F8.2,' s  ',A/ &
2110            '       Cross sections at ',A1,' = ',A/ &
2111            '       scalar-coordinates:   ',A,' m'/)
2112336 FORMAT (/'    3D-Arrays',A,':')
2113337 FORMAT (/'       Arrays: ',A/ &
2114            '       Output every             ',F8.2,' s  ',A/ &
2115            '       Upper output limit at    ',F8.2,' m  (GP ',I4,')'/)
2116339 FORMAT ('       No output during initial ',F8.2,' s')
2117340 FORMAT (/'    Time series:')
2118341 FORMAT ('       Output every             ',F8.2,' s'/)
2119342 FORMAT (/'       ',A2,'-cross-section  Arrays: ',A/ &
2120            '       Output every             ',F8.2,' s  ',A/ &
2121            '       Time averaged over       ',F8.2,' s'/ &
2122            '       Averaging input every    ',F8.2,' s'/ &
2123            '       Cross sections at ',A1,' = ',A/ &
2124            '       scalar-coordinates:   ',A,' m'/)
2125343 FORMAT (/'       Arrays: ',A/ &
2126            '       Output every             ',F8.2,' s  ',A/ &
2127            '       Time averaged over       ',F8.2,' s'/ &
2128            '       Averaging input every    ',F8.2,' s'/ &
2129            '       Upper output limit at    ',F8.2,' m  (GP ',I4,')'/)
2130344 FORMAT ('       Output format: ',A/)
2131345 FORMAT (/'    Scaling lengths for output locations of all subsequent mask IDs:',/ &
2132            '       mask_scale_x (in x-direction): ',F9.3, ' m',/ &
2133            '       mask_scale_y (in y-direction): ',F9.3, ' m',/ &
2134            '       mask_scale_z (in z-direction): ',F9.3, ' m' )
2135346 FORMAT (/'    Masked data output',A,' for mask ID ',I2, ':')
2136347 FORMAT ('       Variables: ',A/ &
2137            '       Output every             ',F8.2,' s')
2138348 FORMAT ('       Variables: ',A/ &
2139            '       Output every             ',F8.2,' s'/ &
2140            '       Time averaged over       ',F8.2,' s'/ &
2141            '       Averaging input every    ',F8.2,' s')
2142349 FORMAT (/'       Output locations in ',A,'-direction in multiples of ', &
2143            'mask_scale_',A,' predefined by array mask_',I2.2,'_',A,':'/ &
2144            13('       ',8(F8.2,',')/) )
2145350 FORMAT (/'       Output locations in ',A,'-direction: ', &
2146            'all gridpoints along ',A,'-direction (default).' )
2147351 FORMAT (/'       Output locations in ',A,'-direction in multiples of ', &
2148            'mask_scale_',A,' constructed from array mask_',I2.2,'_',A,'_loop:'/ &
2149            '          loop begin:',F8.2,', end:',F8.2,', stride:',F8.2 )
2150352 FORMAT  (/'       Number of output time levels allowed: ',I3 /)
2151353 FORMAT  (/'       Number of output time levels allowed: unlimited' /)
2152#if defined( __dvrp_graphics )
2153360 FORMAT ('    Plot-Sequence with dvrp-software:'/ &
2154            '       Output every      ',F7.1,' s'/ &
2155            '       Output mode:      ',A/ &
2156            '       Host / User:      ',A,' / ',A/ &
2157            '       Directory:        ',A// &
2158            '       The sequence contains:')
2159361 FORMAT (/'       Isosurface of "',A,'"    Threshold value: ', E12.3/ &
2160            '          Isosurface color: (',F4.2,',',F4.2,',',F4.2,') (R,G,B)')
2161362 FORMAT (/'       Slicer plane ',A/ &
2162            '       Slicer limits: [',F6.2,',',F6.2,']')
2163363 FORMAT (/'       Particles'/ &
2164            '          particle size:  ',F7.2,' m')
2165364 FORMAT ('          particle ',A,' controlled by "',A,'" with interval [', &
2166                       F6.2,',',F6.2,']')
2167365 FORMAT (/'       Groundplate color: (',F4.2,',',F4.2,',',F4.2,') (R,G,B)'/ &
2168            '       Superelevation along (x,y,z): (',F4.1,',',F4.1,',',F4.1, &
2169                     ')'/ &
2170            '       Clipping limits: from x = ',F9.1,' m to x = ',F9.1,' m'/ &
2171            '                        from y = ',F9.1,' m to y = ',F9.1,' m')
2172366 FORMAT (/'       Topography color: (',F4.2,',',F4.2,',',F4.2,') (R,G,B)')
2173367 FORMAT ('       Polygon reduction for topography: cluster_size = ', I1)
2174#endif
2175#if defined( __spectra )
2176370 FORMAT ('    Spectra:')
2177371 FORMAT ('       Output every ',F7.1,' s'/)
2178372 FORMAT ('       Arrays:     ', 10(A5,',')/                         &
2179            '       Directions: ', 10(A5,',')/                         &
2180            '       height levels  k = ', 20(I3,',')/                  &
2181            '                          ', 20(I3,',')/                  &
2182            '                          ', 20(I3,',')/                  &
2183            '                          ', 20(I3,',')/                  &
2184            '                          ', 19(I3,','),I3,'.'/           &
2185            '       height levels selected for standard plot:'/        &
2186            '                      k = ', 20(I3,',')/                  &
2187            '                          ', 20(I3,',')/                  &
2188            '                          ', 20(I3,',')/                  &
2189            '                          ', 20(I3,',')/                  &
2190            '                          ', 19(I3,','),I3,'.'/           &
2191            '       Time averaged over ', F7.1, ' s,' /                &
2192            '       Profiles for the time averaging are taken every ', &
2193                    F6.1,' s')
2194#endif
2195400 FORMAT (//' Physical quantities:'/ &
2196              ' -------------------'/)
2197410 FORMAT ('    Geograph. latitude  :   phi    = ',F4.1,' degr'/   &
2198            '    Angular velocity    :   omega  = ',E9.3,' rad/s'/  &
2199            '    Coriolis parameter  :   f      = ',F9.6,' 1/s'/    &
2200            '                            f*     = ',F9.6,' 1/s')
2201411 FORMAT (/'    Gravity             :   g      = ',F4.1,' m/s**2')
2202412 FORMAT (/'    Reference state used in buoyancy terms: ',A)
2203413 FORMAT ('       Reference density in buoyancy terms: ',F8.3,' kg/m**3')
2204414 FORMAT ('       Reference temperature in buoyancy terms: ',F8.4,' K')
2205415 FORMAT (/' Cloud physics parameters:'/ &
2206             ' ------------------------'/)
2207416 FORMAT ('    Surface pressure   :   p_0   = ',F7.2,' hPa'/      &
2208            '    Gas constant       :   R     = ',F5.1,' J/(kg K)'/ &
2209            '    Density of air     :   rho_0 = ',F5.3,' kg/m**3'/  &
2210            '    Specific heat cap. :   c_p   = ',F6.1,' J/(kg K)'/ &
2211            '    Vapourization heat :   L_v   = ',E8.2,' J/kg')
2212417 FORMAT ('    Geograph. longitude :   lambda = ',F4.1,' degr')
2213418 FORMAT (/'    Day of the year at model start :   day_init      =     ',I3 &
2214            /'    UTC time at model start        :   time_utc_init = ',F7.1' s')
2215419 FORMAT (//' Land surface model information:'/ &
2216              ' ------------------------------'/)
2217420 FORMAT (/'    Characteristic levels of the initial temperature profile:'// &
2218            '       Height:        ',A,'  m'/ &
2219            '       Temperature:   ',A,'  K'/ &
2220            '       Gradient:      ',A,'  K/100m'/ &
2221            '       Gridpoint:     ',A)
2222421 FORMAT (/'    Characteristic levels of the initial humidity profile:'// &
2223            '       Height:      ',A,'  m'/ &
2224            '       Humidity:    ',A,'  kg/kg'/ &
2225            '       Gradient:    ',A,'  (kg/kg)/100m'/ &
2226            '       Gridpoint:   ',A)
2227422 FORMAT (/'    Characteristic levels of the initial scalar profile:'// &
2228            '       Height:                  ',A,'  m'/ &
2229            '       Scalar concentration:    ',A,'  kg/m**3'/ &
2230            '       Gradient:                ',A,'  (kg/m**3)/100m'/ &
2231            '       Gridpoint:               ',A)
2232423 FORMAT (/'    Characteristic levels of the geo. wind component ug:'// &
2233            '       Height:      ',A,'  m'/ &
2234            '       ug:          ',A,'  m/s'/ &
2235            '       Gradient:    ',A,'  1/100s'/ &
2236            '       Gridpoint:   ',A)
2237424 FORMAT (/'    Characteristic levels of the geo. wind component vg:'// &
2238            '       Height:      ',A,'  m'/ &
2239            '       vg:          ',A,'  m/s'/ &
2240            '       Gradient:    ',A,'  1/100s'/ &
2241            '       Gridpoint:   ',A)
2242425 FORMAT (/'    Characteristic levels of the initial salinity profile:'// &
2243            '       Height:     ',A,'  m'/ &
2244            '       Salinity:   ',A,'  psu'/ &
2245            '       Gradient:   ',A,'  psu/100m'/ &
2246            '       Gridpoint:  ',A)
2247426 FORMAT (/'    Characteristic levels of the subsidence/ascent profile:'// &
2248            '       Height:      ',A,'  m'/ &
2249            '       w_subs:      ',A,'  m/s'/ &
2250            '       Gradient:    ',A,'  (m/s)/100m'/ &
2251            '       Gridpoint:   ',A)
2252427 FORMAT (/'    Initial wind profiles (u,v) are interpolated from given'// &
2253                  ' profiles')
2254428 FORMAT (/'    Initial profiles (u, v, pt, q) are taken from file '/ &
2255             '    NUDGING_DATA')
2256430 FORMAT (//' Cloud physics quantities / methods:'/ &
2257              ' ----------------------------------'/)
2258431 FORMAT ('    Humidity is treated as purely passive scalar (no condensati', &
2259                 'on)')
2260432 FORMAT ('    Bulk scheme with liquid water potential temperature and'/ &
2261            '    total water content is used.'/ &
2262            '    Condensation is parameterized via 0% - or 100% scheme.')
2263433 FORMAT ('    Cloud droplets treated explicitly using the Lagrangian part', &
2264                 'icle model')
2265434 FORMAT ('    Curvature and solution effecs are considered for growth of', &
2266                 ' droplets < 1.0E-6 m')
2267435 FORMAT ('    Droplet collision is handled by ',A,'-kernel')
2268436 FORMAT ('       Fast kernel with fixed radius- and dissipation classes ', &
2269                    'are used'/ &
2270            '          number of radius classes:       ',I3,'    interval ', &
2271                       '[1.0E-6,2.0E-4] m'/ &
2272            '          number of dissipation classes:   ',I2,'    interval ', &
2273                       '[0,1000] cm**2/s**3')
2274437 FORMAT ('    Droplet collision is switched off')
2275438 FORMAT (' --> Land surface type  : ',A,/ &
2276            ' --> Soil porosity type : ',A)
2277439 FORMAT (/'    Initial soil temperature and moisture profile:'// &
2278            '       Height:        ',A,'  m'/ &
2279            '       Temperature:   ',A,'  K'/ &
2280            '       Moisture:      ',A,'  m**3/m**3'/ &
2281            '       Root fraction: ',A,'  '/ &
2282            '       Gridpoint:     ',A)
2283440 FORMAT (/' --> Dewfall is allowed (default)')
2284441 FORMAT (' --> Dewfall is inhibited')
2285442 FORMAT (' --> Soil bottom is closed (water content is conserved, default)')
2286443 FORMAT (' --> Soil bottom is open (water content is not conserved)')
2287444 FORMAT (//' Radiation model information:'/                                 &
2288              ' ----------------------------'/)
2289445 FORMAT (' --> Using constant net radiation: net_radiation = ', F6.2, '  W/m**2')
2290446 FORMAT (' --> Simple radiation scheme for clear sky is used (no clouds,',  &
2291                   ' default)')
2292447 FORMAT (' --> RRTMG scheme is used')
2293448 FORMAT (/'     User-specific surface albedo: albedo = ', F5.3)
2294449 FORMAT  ('     Timestep: dt_radiation = ', F5.2, '  s')
2295
2296450 FORMAT (//' LES / Turbulence quantities:'/ &
2297              ' ---------------------------'/)
2298451 FORMAT ('    Diffusion coefficients are constant:'/ &
2299            '    Km = ',F6.2,' m**2/s   Kh = ',F6.2,' m**2/s   Pr = ',F5.2)
2300453 FORMAT ('    Mixing length is limited to ',F4.2,' * z')
2301454 FORMAT ('    TKE is not allowed to fall below ',E9.2,' (m/s)**2')
2302455 FORMAT ('    initial TKE is prescribed as ',E9.2,' (m/s)**2')
2303456 FORMAT (/'    Albedo is set for land surface type: ', A)
2304457 FORMAT (/'    --> Albedo is fixed during the run')
2305458 FORMAT (/'    --> Longwave radiation is disabled')
2306459 FORMAT (/'    --> Shortwave radiation is disabled.')
2307470 FORMAT (//' Actions during the simulation:'/ &
2308              ' -----------------------------'/)
2309471 FORMAT ('    Disturbance impulse (u,v) every :   ',F6.2,' s'/            &
2310            '    Disturbance amplitude           :     ',F4.2, ' m/s'/       &
2311            '    Lower disturbance level         : ',F8.2,' m (GP ',I4,')'/  &
2312            '    Upper disturbance level         : ',F8.2,' m (GP ',I4,')')
2313472 FORMAT ('    Disturbances continued during the run from i/j =',I4, &
2314                 ' to i/j =',I4)
2315473 FORMAT ('    Disturbances cease as soon as the disturbance energy exceeds',&
2316                 1X,F5.3, ' m**2/s**2')
2317474 FORMAT ('    Random number generator used    : ',A/)
2318475 FORMAT ('    The surface temperature is increased (or decreased, ', &
2319                 'respectively, if'/ &
2320            '    the value is negative) by ',F5.2,' K at the beginning of the',&
2321                 ' 3D-simulation'/)
2322476 FORMAT ('    The surface humidity is increased (or decreased, ',&
2323                 'respectively, if the'/ &
2324            '    value is negative) by ',E8.1,' kg/kg at the beginning of', &
2325                 ' the 3D-simulation'/)
2326477 FORMAT ('    The scalar value is increased at the surface (or decreased, ',&
2327                 'respectively, if the'/ &
2328            '    value is negative) by ',E8.1,' kg/m**3 at the beginning of', &
2329                 ' the 3D-simulation'/)
2330480 FORMAT ('    Particles:'/ &
2331            '    ---------'// &
2332            '       Particle advection is active (switched on at t = ', F7.1, &
2333                    ' s)'/ &
2334            '       Start of new particle generations every  ',F6.1,' s'/ &
2335            '       Boundary conditions: left/right: ', A, ' north/south: ', A/&
2336            '                            bottom:     ', A, ' top:         ', A/&
2337            '       Maximum particle age:                 ',F9.1,' s'/ &
2338            '       Advection stopped at t = ',F9.1,' s'/)
2339481 FORMAT ('       Particles have random start positions'/)
2340482 FORMAT ('          Particles are advected only horizontally'/)
2341483 FORMAT ('       Particles have tails with a maximum of ',I3,' points')
2342484 FORMAT ('            Number of tails of the total domain: ',I10/ &
2343            '            Minimum distance between tailpoints: ',F8.2,' m'/ &
2344            '            Maximum age of the end of the tail:  ',F8.2,' s')
2345485 FORMAT ('       Particle data are written on file every ', F9.1, ' s')
2346486 FORMAT ('       Particle statistics are written on file'/)
2347487 FORMAT ('       Number of particle groups: ',I2/)
2348488 FORMAT ('       SGS velocity components are used for particle advection'/ &
2349            '          minimum timestep for advection: ', F7.5/)
2350489 FORMAT ('       Number of particles simultaneously released at each ', &
2351                    'point: ', I5/)
2352490 FORMAT ('       Particle group ',I2,':'/ &
2353            '          Particle radius: ',E10.3, 'm')
2354491 FORMAT ('          Particle inertia is activated'/ &
2355            '             density_ratio (rho_fluid/rho_particle) = ',F5.3/)
2356492 FORMAT ('          Particles are advected only passively (no inertia)'/)
2357493 FORMAT ('          Boundaries of particle source: x:',F8.1,' - ',F8.1,' m'/&
2358            '                                         y:',F8.1,' - ',F8.1,' m'/&
2359            '                                         z:',F8.1,' - ',F8.1,' m'/&
2360            '          Particle distances:  dx = ',F8.1,' m  dy = ',F8.1, &
2361                       ' m  dz = ',F8.1,' m'/)
2362494 FORMAT ('       Output of particle time series in NetCDF format every ', &
2363                    F8.2,' s'/)
2364495 FORMAT ('       Number of particles in total domain: ',I10/)
2365496 FORMAT ('       Initial vertical particle positions are interpreted ', &
2366                    'as relative to the given topography')
2367500 FORMAT (//' 1D-Model parameters:'/                           &
2368              ' -------------------'//                           &
2369            '    Simulation time:                   ',F8.1,' s'/ &
2370            '    Run-controll output every:         ',F8.1,' s'/ &
2371            '    Vertical profile output every:     ',F8.1,' s'/ &
2372            '    Mixing length calculation:         ',A/         &
2373            '    Dissipation calculation:           ',A/)
2374502 FORMAT ('    Damping layer starts from ',F7.1,' m (GP ',I4,')'/)
2375503 FORMAT (' --> Momentum advection via Wicker-Skamarock-Scheme 5th order')
2376504 FORMAT (' --> Scalar advection via Wicker-Skamarock-Scheme 5th order')
2377505 FORMAT ('    Precipitation parameterization via Seifert-Beheng-Scheme')
2378506 FORMAT ('    Drizzle parameterization via Stokes law')
2379507 FORMAT ('    Turbulence effects on precipitation process')
2380508 FORMAT ('    Ventilation effects on evaporation of rain drops')
2381509 FORMAT ('    Slope limiter used for sedimentation process')
2382510 FORMAT ('    Droplet density    :   N_c   = ',F6.1,' 1/cm**3')
2383511 FORMAT ('    Sedimentation Courant number:                  '/&
2384            '                               C_s   = ',F3.1,'        ')
2385512 FORMAT (/' Date:                 ',A8,6X,'Run:       ',A20/      &
2386            ' Time:                 ',A8,6X,'Run-No.:   ',I2.2/     &
2387            ' Run on host:        ',A10,6X,'En-No.:    ',I2.2)
2388513 FORMAT (' --> Scalar advection via Wicker-Skamarock-Scheme 5th order ' // & 
2389            '+ monotonic adjustment')
2390
2391
2392 END SUBROUTINE header
Note: See TracBrowser for help on using the repository browser.