Home

Context Navigation

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

Last change on this file since 1833 was 1833, checked in by raasch, 8 years ago
spectrum renamed spactra_par and further modularized, POINTER-attributes added in coupler-routines to avoid gfortran error messages
Property svn:keywords set to `Id`
File size: 85.8 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |