Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |