Home

Context Navigation

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

Last change on this file since 2284 was 2270, checked in by maronga, 8 years ago
major revisions in land surface model
Property svn:keywords set to `Id`
File size: 90.7 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |