Home

Context Navigation

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

Last change on this file since 2338 was 2338, checked in by gronemeier, 7 years ago
modularized 1d model
Property svn:keywords set to `Id`
File size: 88.8 KB

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

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |