Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |