Home

Context Navigation

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

Last change on this file since 1832 was 1832, checked in by hoffmann, 8 years ago
last commit documented
Property svn:keywords set to `Id`
File size: 87.6 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	!
22	!
23	! Former revisions:
24	! -----------------
25	! $Id: header.f90 1832 2016-04-07 13:28:15Z hoffmann $
26	!
27	! 1831 2016-04-07 13:15:51Z hoffmann
28	! turbulence renamed collision_turbulence,
29	! drizzle renamed cloud_water_sedimentation
30	!
31	! 1826 2016-04-07 12:01:39Z maronga
32	! Moved radiation model header output to the respective module.
33	! Moved canopy model header output to the respective module.
34	!
35	! 1822 2016-04-07 07:49:42Z hoffmann
36	! Tails removed. icloud_scheme replaced by microphysics_*
37	!
38	! 1817 2016-04-06 15:44:20Z maronga
39	! Moved land_surface_model header output to the respective module.
40	!
41	! 1808 2016-04-05 19:44:00Z raasch
42	! routine local_flush replaced by FORTRAN statement
43	!
44	! 1797 2016-03-21 16:50:28Z raasch
45	! output of nesting datatransfer mode
46	!
47	! 1791 2016-03-11 10:41:25Z raasch
48	! output of nesting informations of all domains
49	!
50	! 1788 2016-03-10 11:01:04Z maronga
51	! Parameter dewfall removed
52	!
53	! 1786 2016-03-08 05:49:27Z raasch
54	! cpp-direktives for spectra removed
55	!
56	! 1783 2016-03-06 18:36:17Z raasch
57	! netcdf module and variable names changed, output of netcdf_deflate
58	!
59	! 1764 2016-02-28 12:45:19Z raasch
60	! output of nesting informations
61	!
62	! 1697 2015-10-28 17:14:10Z raasch
63	! small E- and F-FORMAT changes to avoid informative compiler messages about
64	! insufficient field width
65	!
66	! 1691 2015-10-26 16:17:44Z maronga
67	! Renamed prandtl_layer to constant_flux_layer, renames rif_min/rif_max to
68	! zeta_min/zeta_max.
69	!
70	! 1682 2015-10-07 23:56:08Z knoop
71	! Code annotations made doxygen readable
72	!
73	! 1675 2015-10-02 08:28:59Z gronemeier
74	! Bugfix: Definition of topography grid levels
75	!
76	! 1660 2015-09-21 08:15:16Z gronemeier
77	! Bugfix: Definition of building/street canyon height if vertical grid stretching
78	! starts below the maximum topography height.
79	!
80	! 1590 2015-05-08 13:56:27Z maronga
81	! Bugfix: Added TRIM statements for character strings for LSM and radiation code
82	!
83	! 1585 2015-04-30 07:05:52Z maronga
84	! Further output for radiation model(s).
85	!
86	! 1575 2015-03-27 09:56:27Z raasch
87	! adjustments for psolver-queries, output of seed_follows_topography
88	!
89	! 1560 2015-03-06 10:48:54Z keck
90	! output for recycling y shift
91	!
92	! 1557 2015-03-05 16:43:04Z suehring
93	! output for monotonic limiter
94	!
95	! 1551 2015-03-03 14:18:16Z maronga
96	! Added informal output for land surface model and radiation model. Removed typo.
97	!
98	! 1496 2014-12-02 17:25:50Z maronga
99	! Renamed: "radiation -> "cloud_top_radiation"
100	!
101	! 1484 2014-10-21 10:53:05Z kanani
102	! Changes due to new module structure of the plant canopy model:
103	! module plant_canopy_model_mod and output for new canopy model parameters
104	! (alpha_lad, beta_lad, lai_beta,...) added,
105	! drag_coefficient, leaf_surface_concentration and scalar_exchange_coefficient
106	! renamed to canopy_drag_coeff, leaf_surface_conc and leaf_scalar_exch_coeff,
107	! learde renamed leaf_area_density.
108	! Bugfix: DO-WHILE-loop for lad header information additionally restricted
109	! by maximum number of gradient levels (currently 10)
110	!
111	! 1482 2014-10-18 12:34:45Z raasch
112	! information about calculated or predefined virtual processor topology adjusted
113	!
114	! 1468 2014-09-24 14:06:57Z maronga
115	! Adapted for use on up to 6-digit processor cores
116	!
117	! 1429 2014-07-15 12:53:45Z knoop
118	! header exended to provide ensemble_member_nr if specified
119	!
120	! 1376 2014-04-26 11:21:22Z boeske
121	! Correction of typos
122	!
123	! 1365 2014-04-22 15:03:56Z boeske
124	! New section 'Large scale forcing and nudging':
125	! output of large scale forcing and nudging information,
126	! new section for initial profiles created
127	!
128	! 1359 2014-04-11 17:15:14Z hoffmann
129	! dt_sort_particles removed
130	!
131	! 1353 2014-04-08 15:21:23Z heinze
132	! REAL constants provided with KIND-attribute
133	!
134	! 1327 2014-03-21 11:00:16Z raasch
135	! parts concerning iso2d and avs output removed,
136	! -netcdf output queries
137	!
138	! 1324 2014-03-21 09:13:16Z suehring
139	! Bugfix: module spectrum added
140	!
141	! 1322 2014-03-20 16:38:49Z raasch
142	! REAL functions provided with KIND-attribute,
143	! some REAL constants defined as wp-kind
144	!
145	! 1320 2014-03-20 08:40:49Z raasch
146	! ONLY-attribute added to USE-statements,
147	! kind-parameters added to all INTEGER and REAL declaration statements,
148	! kinds are defined in new module kinds,
149	! revision history before 2012 removed,
150	! comment fields (!:) to be used for variable explanations added to
151	! all variable declaration statements
152	!
153	! 1308 2014-03-13 14:58:42Z fricke
154	! output of the fixed number of output time levels
155	! output_format adjusted for masked data if netcdf_data_format > 5
156	!
157	! 1299 2014-03-06 13:15:21Z heinze
158	! output for using large_scale subsidence in combination
159	! with large_scale_forcing
160	! reformatting, more detailed explanations
161	!
162	! 1241 2013-10-30 11:36:58Z heinze
163	! output for nudging + large scale forcing from external file
164	!
165	! 1216 2013-08-26 09:31:42Z raasch
166	! output for transpose_compute_overlap
167	!
168	! 1212 2013-08-15 08:46:27Z raasch
169	! output for poisfft_hybrid removed
170	!
171	! 1179 2013-06-14 05:57:58Z raasch
172	! output of reference_state, use_reference renamed use_single_reference_value
173	!
174	! 1159 2013-05-21 11:58:22Z fricke
175	! +use_cmax
176	!
177	! 1115 2013-03-26 18:16:16Z hoffmann
178	! descriptions for Seifert-Beheng-cloud-physics-scheme added
179	!
180	! 1111 2013-03-08 23:54:10Z raasch
181	! output of accelerator board information
182	! ibc_p_b = 2 removed
183	!
184	! 1108 2013-03-05 07:03:32Z raasch
185	! bugfix for r1106
186	!
187	! 1106 2013-03-04 05:31:38Z raasch
188	! some format changes for coupled runs
189	!
190	! 1092 2013-02-02 11:24:22Z raasch
191	! unused variables removed
192	!
193	! 1036 2012-10-22 13:43:42Z raasch
194	! code put under GPL (PALM 3.9)
195	!
196	! 1031 2012-10-19 14:35:30Z raasch
197	! output of netCDF data format modified
198	!
199	! 1015 2012-09-27 09:23:24Z raasch
200	! output of Adjustment of mixing length to the Prandtl mixing length at first
201	! grid point above ground removed
202	!
203	! 1003 2012-09-14 14:35:53Z raasch
204	! output of information about equal/unequal subdomain size removed
205	!
206	! 1001 2012-09-13 14:08:46Z raasch
207	! all actions concerning leapfrog- and upstream-spline-scheme removed
208	!
209	! 978 2012-08-09 08:28:32Z fricke
210	! -km_damp_max, outflow_damping_width
211	! +pt_damping_factor, pt_damping_width
212	! +z0h
213	!
214	! 964 2012-07-26 09:14:24Z raasch
215	! output of profil-related quantities removed
216	!
217	! 940 2012-07-09 14:31:00Z raasch
218	! Output in case of simulations for pure neutral stratification (no pt-equation
219	! solved)
220	!
221	! 927 2012-06-06 19:15:04Z raasch
222	! output of masking_method for mg-solver
223	!
224	! 868 2012-03-28 12:21:07Z raasch
225	! translation velocity in Galilean transformation changed to 0.6 * ug
226	!
227	! 833 2012-02-22 08:55:55Z maronga
228	! Adjusted format for leaf area density
229	!
230	! 828 2012-02-21 12:00:36Z raasch
231	! output of dissipation_classes + radius_classes
232	!
233	! 825 2012-02-19 03:03:44Z raasch
234	! Output of cloud physics parameters/quantities complemented and restructured
235	!
236	! Revision 1.1 1997/08/11 06:17:20 raasch
237	! Initial revision
238	!
239	!
240	! Description:
241	! ------------
242	!> Writing a header with all important information about the current run.
243	!> This subroutine is called three times, two times at the beginning
244	!> (writing information on files RUN_CONTROL and HEADER) and one time at the
245	!> end of the run, then writing additional information about CPU-usage on file
246	!> header.
247	!-----------------------------------------------------------------------------!
248	SUBROUTINE header
249
250
251	USE arrays_3d, &
252	ONLY: pt_init, qsws, q_init, sa_init, shf, ug, vg, w_subs, zu, zw
253
254	USE control_parameters
255
256	USE cloud_parameters, &
257	ONLY: cloud_water_sedimentation, collision_turbulence, cp, &
258	c_sedimentation, limiter_sedimentation, l_v, nc_const, &
259	r_d, ventilation_effect
260
261	USE cpulog, &
262	ONLY: log_point_s
263
264	USE dvrp_variables, &
265	ONLY: use_seperate_pe_for_dvrp_output
266
267	USE grid_variables, &
268	ONLY: dx, dy
269
270	USE indices, &
271	ONLY: mg_loc_ind, nnx, nny, nnz, nx, ny, nxl_mg, nxr_mg, nyn_mg, &
272	nys_mg, nzt, nzt_mg
273
274	USE kinds
275
276	USE land_surface_model_mod, &
277	ONLY: land_surface, lsm_header
278
279	USE model_1d, &
280	ONLY: damp_level_ind_1d, dt_pr_1d, dt_run_control_1d, end_time_1d
281
282	USE netcdf_interface, &
283	ONLY: netcdf_data_format, netcdf_data_format_string, netcdf_deflate
284
285	USE particle_attributes, &
286	ONLY: bc_par_b, bc_par_lr, bc_par_ns, bc_par_t, collision_kernel, &
287	curvature_solution_effects, &
288	density_ratio, dissipation_classes, dt_min_part, dt_prel, &
289	dt_write_particle_data, end_time_prel, &
290	number_of_particle_groups, particle_advection, &
291	particle_advection_start, &
292	particles_per_point, pdx, pdy, pdz, psb, psl, psn, psr, pss, &
293	pst, radius, radius_classes, random_start_position, &
294	seed_follows_topography, &
295	total_number_of_particles, use_sgs_for_particles, &
296	vertical_particle_advection, write_particle_statistics
297
298	USE pegrid
299
300	USE plant_canopy_model_mod, &
301	ONLY: pcm_header, plant_canopy
302
303	USE pmc_handle_communicator, &
304	ONLY: pmc_get_model_info
305
306	USE pmc_interface, &
307	ONLY: nested_run, nesting_datatransfer_mode, nesting_mode
308
309	USE radiation_model_mod, &
310	ONLY: radiation, radiation_header
311
312	USE spectrum, &
313	ONLY: comp_spectra_level, data_output_sp, plot_spectra_level, &
314	spectra_direction
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	!-- Spectra output
1634	IF ( dt_dosp /= 9999999.9_wp ) THEN
1635	WRITE ( io, 370 )
1636
1637	output_format = netcdf_data_format_string
1638	IF ( netcdf_deflate == 0 ) THEN
1639	WRITE ( io, 344 ) output_format
1640	ELSE
1641	WRITE ( io, 354 ) TRIM( output_format ), netcdf_deflate
1642	ENDIF
1643	WRITE ( io, 371 ) dt_dosp
1644	IF ( skip_time_dosp /= 0.0_wp ) WRITE ( io, 339 ) skip_time_dosp
1645	WRITE ( io, 372 ) ( data_output_sp(i), i = 1,10 ), &
1646	( spectra_direction(i), i = 1,10 ), &
1647	( comp_spectra_level(i), i = 1,100 ), &
1648	( plot_spectra_level(i), i = 1,100 ), &
1649	averaging_interval_sp, dt_averaging_input_pr
1650	ENDIF
1651
1652	WRITE ( io, 99 )
1653
1654	!
1655	!-- Physical quantities
1656	WRITE ( io, 400 )
1657
1658	!
1659	!-- Geostrophic parameters
1660	WRITE ( io, 410 ) phi, omega, f, fs
1661
1662	!
1663	!-- Other quantities
1664	WRITE ( io, 411 ) g
1665
1666	WRITE ( io, 412 ) TRIM( reference_state )
1667	IF ( use_single_reference_value ) THEN
1668	IF ( ocean ) THEN
1669	WRITE ( io, 413 ) prho_reference
1670	ELSE
1671	WRITE ( io, 414 ) pt_reference
1672	ENDIF
1673	ENDIF
1674
1675	!
1676	!-- Cloud physics parameters
1677	IF ( cloud_physics ) THEN
1678	WRITE ( io, 415 )
1679	WRITE ( io, 416 ) surface_pressure, r_d, rho_surface, cp, l_v
1680	IF ( microphysics_seifert ) THEN
1681	WRITE ( io, 510 ) 1.0E-6_wp * nc_const
1682	WRITE ( io, 511 ) c_sedimentation
1683	ENDIF
1684	ENDIF
1685
1686	!
1687	!-- Cloud physcis parameters / quantities / numerical methods
1688	WRITE ( io, 430 )
1689	IF ( humidity .AND. .NOT. cloud_physics .AND. .NOT. cloud_droplets) THEN
1690	WRITE ( io, 431 )
1691	ELSEIF ( humidity .AND. cloud_physics ) THEN
1692	WRITE ( io, 432 )
1693	IF ( cloud_top_radiation ) WRITE ( io, 132 )
1694	IF ( microphysics_kessler ) THEN
1695	WRITE ( io, 133 )
1696	ELSEIF ( microphysics_seifert ) THEN
1697	IF ( cloud_water_sedimentation ) WRITE ( io, 506 )
1698	WRITE ( io, 505 )
1699	IF ( collision_turbulence ) WRITE ( io, 507 )
1700	IF ( ventilation_effect ) WRITE ( io, 508 )
1701	IF ( limiter_sedimentation ) WRITE ( io, 509 )
1702	ENDIF
1703	ELSEIF ( humidity .AND. cloud_droplets ) THEN
1704	WRITE ( io, 433 )
1705	IF ( curvature_solution_effects ) WRITE ( io, 434 )
1706	IF ( collision_kernel /= 'none' ) THEN
1707	WRITE ( io, 435 ) TRIM( collision_kernel )
1708	IF ( collision_kernel(6:9) == 'fast' ) THEN
1709	WRITE ( io, 436 ) radius_classes, dissipation_classes
1710	ENDIF
1711	ELSE
1712	WRITE ( io, 437 )
1713	ENDIF
1714	ENDIF
1715
1716	!
1717	!-- LES / turbulence parameters
1718	WRITE ( io, 450 )
1719
1720	!--
1721	! ... LES-constants used must still be added here
1722	!--
1723	IF ( constant_diffusion ) THEN
1724	WRITE ( io, 451 ) km_constant, km_constant/prandtl_number, &
1725	prandtl_number
1726	ENDIF
1727	IF ( .NOT. constant_diffusion) THEN
1728	IF ( e_init > 0.0_wp ) WRITE ( io, 455 ) e_init
1729	IF ( e_min > 0.0_wp ) WRITE ( io, 454 ) e_min
1730	IF ( wall_adjustment ) WRITE ( io, 453 ) wall_adjustment_factor
1731	ENDIF
1732
1733	!
1734	!-- Special actions during the run
1735	WRITE ( io, 470 )
1736	IF ( create_disturbances ) THEN
1737	WRITE ( io, 471 ) dt_disturb, disturbance_amplitude, &
1738	zu(disturbance_level_ind_b), disturbance_level_ind_b,&
1739	zu(disturbance_level_ind_t), disturbance_level_ind_t
1740	IF ( .NOT. bc_lr_cyc .OR. .NOT. bc_ns_cyc ) THEN
1741	WRITE ( io, 472 ) inflow_disturbance_begin, inflow_disturbance_end
1742	ELSE
1743	WRITE ( io, 473 ) disturbance_energy_limit
1744	ENDIF
1745	WRITE ( io, 474 ) TRIM( random_generator )
1746	ENDIF
1747	IF ( pt_surface_initial_change /= 0.0_wp ) THEN
1748	WRITE ( io, 475 ) pt_surface_initial_change
1749	ENDIF
1750	IF ( humidity .AND. q_surface_initial_change /= 0.0_wp ) THEN
1751	WRITE ( io, 476 ) q_surface_initial_change
1752	ENDIF
1753	IF ( passive_scalar .AND. q_surface_initial_change /= 0.0_wp ) THEN
1754	WRITE ( io, 477 ) q_surface_initial_change
1755	ENDIF
1756
1757	IF ( particle_advection ) THEN
1758	!
1759	!-- Particle attributes
1760	WRITE ( io, 480 ) particle_advection_start, dt_prel, bc_par_lr, &
1761	bc_par_ns, bc_par_b, bc_par_t, particle_maximum_age, &
1762	end_time_prel
1763	IF ( use_sgs_for_particles ) WRITE ( io, 488 ) dt_min_part
1764	IF ( random_start_position ) WRITE ( io, 481 )
1765	IF ( seed_follows_topography ) WRITE ( io, 496 )
1766	IF ( particles_per_point > 1 ) WRITE ( io, 489 ) particles_per_point
1767	WRITE ( io, 495 ) total_number_of_particles
1768	IF ( dt_write_particle_data /= 9999999.9_wp ) THEN
1769	WRITE ( io, 485 ) dt_write_particle_data
1770	IF ( netcdf_data_format > 1 ) THEN
1771	output_format = 'netcdf (64 bit offset) and binary'
1772	ELSE
1773	output_format = 'netcdf and binary'
1774	ENDIF
1775	IF ( netcdf_deflate == 0 ) THEN
1776	WRITE ( io, 344 ) output_format
1777	ELSE
1778	WRITE ( io, 354 ) TRIM( output_format ), netcdf_deflate
1779	ENDIF
1780	ENDIF
1781	IF ( dt_dopts /= 9999999.9_wp ) WRITE ( io, 494 ) dt_dopts
1782	IF ( write_particle_statistics ) WRITE ( io, 486 )
1783
1784	WRITE ( io, 487 ) number_of_particle_groups
1785
1786	DO i = 1, number_of_particle_groups
1787	IF ( i == 1 .AND. density_ratio(i) == 9999999.9_wp ) THEN
1788	WRITE ( io, 490 ) i, 0.0_wp
1789	WRITE ( io, 492 )
1790	ELSE
1791	WRITE ( io, 490 ) i, radius(i)
1792	IF ( density_ratio(i) /= 0.0_wp ) THEN
1793	WRITE ( io, 491 ) density_ratio(i)
1794	ELSE
1795	WRITE ( io, 492 )
1796	ENDIF
1797	ENDIF
1798	WRITE ( io, 493 ) psl(i), psr(i), pss(i), psn(i), psb(i), pst(i), &
1799	pdx(i), pdy(i), pdz(i)
1800	IF ( .NOT. vertical_particle_advection(i) ) WRITE ( io, 482 )
1801	ENDDO
1802
1803	ENDIF
1804
1805
1806	!
1807	!-- Parameters of 1D-model
1808	IF ( INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 ) THEN
1809	WRITE ( io, 500 ) end_time_1d, dt_run_control_1d, dt_pr_1d, &
1810	mixing_length_1d, dissipation_1d
1811	IF ( damp_level_ind_1d /= nzt+1 ) THEN
1812	WRITE ( io, 502 ) zu(damp_level_ind_1d), damp_level_ind_1d
1813	ENDIF
1814	ENDIF
1815
1816	!
1817	!-- User-defined information
1818	CALL user_header( io )
1819
1820	WRITE ( io, 99 )
1821
1822	!
1823	!-- Write buffer contents to disc immediately
1824	FLUSH( io )
1825
1826	!
1827	!-- Here the FORMATs start
1828
1829	99 FORMAT (1X,78('-'))
1830	100 FORMAT (/1X,'******************************',4X,44('-')/ &
1831	1X,'* ',A,' *',4X,A/ &
1832	1X,'******************************',4X,44('-'))
1833	101 FORMAT (35X,'coupled run using MPI-',I1,': ',A/ &
1834	35X,42('-'))
1835	102 FORMAT (/' Date: ',A8,4X,'Run: ',A20/ &
1836	' Time: ',A8,4X,'Run-No.: ',I2.2/ &
1837	' Run on host: ',A10)
1838	#if defined( __parallel )
1839	103 FORMAT (' Number of PEs:',10X,I6,4X,'Processor grid (x,y): (',I4,',',I4, &
1840	')',1X,A)
1841	104 FORMAT (' Number of PEs:',10X,I6,4X,'Tasks:',I4,' threads per task:',I4/ &
1842	35X,'Processor grid (x,y): (',I4,',',I4,')',1X,A)
1843	105 FORMAT (35X,'One additional PE is used to handle'/37X,'the dvrp output!')
1844	106 FORMAT (35X,'A 1d-decomposition along x is forced'/ &
1845	35X,'because the job is running on an SMP-cluster')
1846	107 FORMAT (35X,'A 1d-decomposition along ',A,' is used')
1847	108 FORMAT (35X,'Max. # of parallel I/O streams is ',I5)
1848	109 FORMAT (35X,'Precursor run for coupled atmos-ocean run'/ &
1849	35X,42('-'))
1850	114 FORMAT (35X,'Coupled atmosphere-ocean run following'/ &
1851	35X,'independent precursor runs'/ &
1852	35X,42('-'))
1853	117 FORMAT (' Accelerator boards / node: ',I2)
1854	#endif
1855	110 FORMAT (/' Numerical Schemes:'/ &
1856	' -----------------'/)
1857	111 FORMAT (' --> Solve perturbation pressure via FFT using ',A,' routines')
1858	112 FORMAT (' --> Solve perturbation pressure via SOR-Red/Black-Schema'/ &
1859	' Iterations (initial/other): ',I3,'/',I3,' omega =',F6.3)
1860	113 FORMAT (' --> Momentum advection via Piascek-Williams-Scheme (Form C3)', &
1861	' or Upstream')
1862	115 FORMAT (' FFT and transpositions are overlapping')
1863	116 FORMAT (' --> Scalar advection via Piascek-Williams-Scheme (Form C3)', &
1864	' or Upstream')
1865	118 FORMAT (' --> Scalar advection via Bott-Chlond-Scheme')
1866	119 FORMAT (' --> Galilei-Transform applied to horizontal advection:'/ &
1867	' translation velocity = ',A/ &
1868	' distance advected ',A,': ',F8.3,' km(x) ',F8.3,' km(y)')
1869	120 FORMAT (' Accelerator boards: ',8X,I2)
1870	122 FORMAT (' --> Time differencing scheme: ',A)
1871	123 FORMAT (' --> Rayleigh-Damping active, starts ',A,' z = ',F8.2,' m'/ &
1872	' maximum damping coefficient:',F6.3, ' 1/s')
1873	129 FORMAT (' --> Additional prognostic equation for the specific humidity')
1874	130 FORMAT (' --> Additional prognostic equation for the total water content')
1875	131 FORMAT (' --> No pt-equation solved. Neutral stratification with pt = ', &
1876	F6.2, ' K assumed')
1877	132 FORMAT (' Parameterization of long-wave radiation processes via'/ &
1878	' effective emissivity scheme')
1879	133 FORMAT (' Precipitation parameterization via Kessler-Scheme')
1880	134 FORMAT (' --> Additional prognostic equation for a passive scalar')
1881	135 FORMAT (' --> Solve perturbation pressure via ',A,' method (', &
1882	A,'-cycle)'/ &
1883	' number of grid levels: ',I2/ &
1884	' Gauss-Seidel red/black iterations: ',I2)
1885	136 FORMAT (' gridpoints of coarsest subdomain (x,y,z): (',I3,',',I3,',', &
1886	I3,')')
1887	137 FORMAT (' level data gathered on PE0 at level: ',I2/ &
1888	' gridpoints of coarsest subdomain (x,y,z): (',I3,',',I3,',', &
1889	I3,')'/ &
1890	' gridpoints of coarsest domain (x,y,z): (',I3,',',I3,',', &
1891	I3,')')
1892	139 FORMAT (' --> Loop optimization method: ',A)
1893	140 FORMAT (' maximum residual allowed: ',E10.3)
1894	141 FORMAT (' fixed number of multigrid cycles: ',I4)
1895	142 FORMAT (' perturbation pressure is calculated at every Runge-Kutta ', &
1896	'step')
1897	143 FORMAT (' Euler/upstream scheme is used for the SGS turbulent ', &
1898	'kinetic energy')
1899	144 FORMAT (' masking method is used')
1900	150 FORMAT (' --> Volume flow at the right and north boundary will be ', &
1901	'conserved'/ &
1902	' using the ',A,' mode')
1903	151 FORMAT (' with u_bulk = ',F7.3,' m/s and v_bulk = ',F7.3,' m/s')
1904	152 FORMAT (' --> External pressure gradient directly prescribed by the user:',&
1905	/' ',2(1X,E12.5),'Pa/m in x/y direction', &
1906	/' starting from dp_level_b =', F8.3, 'm', A /)
1907	160 FORMAT (//' Large scale forcing and nudging:'/ &
1908	' -------------------------------'/)
1909	161 FORMAT (' --> No large scale forcing from external is used (default) ')
1910	162 FORMAT (' --> Large scale forcing from external file LSF_DATA is used: ')
1911	163 FORMAT (' - large scale advection tendencies ')
1912	164 FORMAT (' - large scale subsidence velocity w_subs ')
1913	165 FORMAT (' - large scale subsidence tendencies ')
1914	167 FORMAT (' - and geostrophic wind components ug and vg')
1915	168 FORMAT (' --> Large-scale vertical motion is used in the ', &
1916	'prognostic equation(s) for')
1917	169 FORMAT (' the scalar(s) only')
1918	170 FORMAT (' --> Nudging is used')
1919	171 FORMAT (' --> No nudging is used (default) ')
1920	180 FORMAT (' - prescribed surface values for temperature')
1921	181 FORMAT (' - prescribed surface fluxes for temperature')
1922	182 FORMAT (' - prescribed surface values for humidity')
1923	183 FORMAT (' - prescribed surface fluxes for humidity')
1924	200 FORMAT (//' Run time and time step information:'/ &
1925	' ----------------------------------'/)
1926	201 FORMAT ( ' Timestep: variable maximum value: ',F6.3,' s', &
1927	' CFL-factor:',F5.2)
1928	202 FORMAT ( ' Timestep: dt = ',F6.3,' s'/)
1929	203 FORMAT ( ' Start time: ',F9.3,' s'/ &
1930	' End time: ',F9.3,' s')
1931	204 FORMAT ( A,F9.3,' s')
1932	205 FORMAT ( A,F9.3,' s',5X,'restart every',17X,F9.3,' s')
1933	206 FORMAT (/' Time reached: ',F9.3,' s'/ &
1934	' CPU-time used: ',F9.3,' s per timestep: ', &
1935	' ',F9.3,' s'/ &
1936	' per second of simulated tim', &
1937	'e: ',F9.3,' s')
1938	207 FORMAT ( ' Coupling start time: ',F9.3,' s')
1939	250 FORMAT (//' Computational grid and domain size:'/ &
1940	' ----------------------------------'// &
1941	' Grid length: dx = ',F7.3,' m dy = ',F7.3, &
1942	' m dz = ',F7.3,' m'/ &
1943	' Domain size: x = ',F10.3,' m y = ',F10.3, &
1944	' m z(u) = ',F10.3,' m'/)
1945	252 FORMAT (' dz constant up to ',F10.3,' m (k=',I4,'), then stretched by', &
1946	' factor:',F6.3/ &
1947	' maximum dz not to be exceeded is dz_max = ',F10.3,' m'/)
1948	254 FORMAT (' Number of gridpoints (x,y,z): (0:',I4,', 0:',I4,', 0:',I4,')'/ &
1949	' Subdomain size (x,y,z): ( ',I4,', ',I4,', ',I4,')'/)
1950	260 FORMAT (/' The model has a slope in x-direction. Inclination angle: ',F6.2,&
1951	' degrees')
1952	270 FORMAT (//' Topography information:'/ &
1953	' ----------------------'// &
1954	1X,'Topography: ',A)
1955	271 FORMAT ( ' Building size (x/y/z) in m: ',F5.1,' / ',F5.1,' / ',F5.1/ &
1956	' Horizontal index bounds (l/r/s/n): ',I4,' / ',I4,' / ',I4, &
1957	' / ',I4)
1958	272 FORMAT ( ' Single quasi-2D street canyon of infinite length in ',A, &
1959	' direction' / &
1960	' Canyon height: ', F6.2, 'm, ch = ', I4, '.' / &
1961	' Canyon position (',A,'-walls): cxl = ', I4,', cxr = ', I4, '.')
1962	278 FORMAT (' Topography grid definition convention:'/ &
1963	' cell edge (staggered grid points'/ &
1964	' (u in x-direction, v in y-direction))' /)
1965	279 FORMAT (' Topography grid definition convention:'/ &
1966	' cell center (scalar grid points)' /)
1967	300 FORMAT (//' Boundary conditions:'/ &
1968	' -------------------'// &
1969	' p uv ', &
1970	' pt'// &
1971	' B. bound.: ',A/ &
1972	' T. bound.: ',A)
1973	301 FORMAT (/' ',A// &
1974	' B. bound.: ',A/ &
1975	' T. bound.: ',A)
1976	303 FORMAT (/' Bottom surface fluxes are used in diffusion terms at k=1')
1977	304 FORMAT (/' Top surface fluxes are used in diffusion terms at k=nzt')
1978	305 FORMAT (//' Prandtl-Layer between bottom surface and first ', &
1979	'computational u,v-level:'// &
1980	' zp = ',F6.2,' m z0 =',F7.4,' m z0h =',F8.5,&
1981	' m kappa =',F5.2/ &
1982	' Rif value range: ',F8.2,' <= rif <=',F6.2)
1983	306 FORMAT (' Predefined constant heatflux: ',F9.6,' K m/s')
1984	307 FORMAT (' Heatflux has a random normal distribution')
1985	308 FORMAT (' Predefined surface temperature')
1986	309 FORMAT (' Predefined constant salinityflux: ',F9.6,' psu m/s')
1987	310 FORMAT (//' 1D-Model:'// &
1988	' Rif value range: ',F6.2,' <= rif <=',F6.2)
1989	311 FORMAT (' Predefined constant humidity flux: ',E10.3,' m/s')
1990	312 FORMAT (' Predefined surface humidity')
1991	313 FORMAT (' Predefined constant scalar flux: ',E10.3,' kg/(m**2 s)')
1992	314 FORMAT (' Predefined scalar value at the surface')
1993	315 FORMAT (' Humidity / scalar flux at top surface is 0.0')
1994	316 FORMAT (' Sensible heatflux and momentum flux from coupled ', &
1995	'atmosphere model')
1996	317 FORMAT (//' Lateral boundaries:'/ &
1997	' left/right: ',A/ &
1998	' north/south: ',A)
1999	318 FORMAT (/' use_cmax: ',L1 / &
2000	' pt damping layer width = ',F8.2,' m, pt ', &
2001	'damping factor =',F7.4)
2002	319 FORMAT (' turbulence recycling at inflow switched on'/ &
2003	' width of recycling domain: ',F7.1,' m grid index: ',I4/ &
2004	' inflow damping height: ',F6.1,' m width: ',F6.1,' m')
2005	320 FORMAT (' Predefined constant momentumflux: u: ',F9.6,' m2/s2'/ &
2006	' v: ',F9.6,' m2/s2')
2007	321 FORMAT (//' Initial profiles:'/ &
2008	' ----------------')
2009	322 FORMAT (' turbulence recycling at inflow switched on'/ &
2010	' y shift of the recycled inflow turbulence switched on'/ &
2011	' width of recycling domain: ',F7.1,' m grid index: ',I4/ &
2012	' inflow damping height: ',F6.1,' m width: ',F6.1,' m'/)
2013	325 FORMAT (//' List output:'/ &
2014	' -----------'// &
2015	' 1D-Profiles:'/ &
2016	' Output every ',F8.2,' s')
2017	326 FORMAT (' Time averaged over ',F8.2,' s'/ &
2018	' Averaging input every ',F8.2,' s')
2019	330 FORMAT (//' Data output:'/ &
2020	' -----------'/)
2021	331 FORMAT (/' 1D-Profiles:')
2022	332 FORMAT (/' ',A)
2023	333 FORMAT (' Output every ',F8.2,' s',/ &
2024	' Time averaged over ',F8.2,' s'/ &
2025	' Averaging input every ',F8.2,' s')
2026	334 FORMAT (/' 2D-Arrays',A,':')
2027	335 FORMAT (/' ',A2,'-cross-section Arrays: ',A/ &
2028	' Output every ',F8.2,' s ',A/ &
2029	' Cross sections at ',A1,' = ',A/ &
2030	' scalar-coordinates: ',A,' m'/)
2031	336 FORMAT (/' 3D-Arrays',A,':')
2032	337 FORMAT (/' Arrays: ',A/ &
2033	' Output every ',F8.2,' s ',A/ &
2034	' Upper output limit at ',F8.2,' m (GP ',I4,')'/)
2035	339 FORMAT (' No output during initial ',F8.2,' s')
2036	340 FORMAT (/' Time series:')
2037	341 FORMAT (' Output every ',F8.2,' s'/)
2038	342 FORMAT (/' ',A2,'-cross-section Arrays: ',A/ &
2039	' Output every ',F8.2,' s ',A/ &
2040	' Time averaged over ',F8.2,' s'/ &
2041	' Averaging input every ',F8.2,' s'/ &
2042	' Cross sections at ',A1,' = ',A/ &
2043	' scalar-coordinates: ',A,' m'/)
2044	343 FORMAT (/' Arrays: ',A/ &
2045	' Output every ',F8.2,' s ',A/ &
2046	' Time averaged over ',F8.2,' s'/ &
2047	' Averaging input every ',F8.2,' s'/ &
2048	' Upper output limit at ',F8.2,' m (GP ',I4,')'/)
2049	344 FORMAT (' Output format: ',A/)
2050	345 FORMAT (/' Scaling lengths for output locations of all subsequent mask IDs:',/ &
2051	' mask_scale_x (in x-direction): ',F9.3, ' m',/ &
2052	' mask_scale_y (in y-direction): ',F9.3, ' m',/ &
2053	' mask_scale_z (in z-direction): ',F9.3, ' m' )
2054	346 FORMAT (/' Masked data output',A,' for mask ID ',I2, ':')
2055	347 FORMAT (' Variables: ',A/ &
2056	' Output every ',F8.2,' s')
2057	348 FORMAT (' Variables: ',A/ &
2058	' Output every ',F8.2,' s'/ &
2059	' Time averaged over ',F8.2,' s'/ &
2060	' Averaging input every ',F8.2,' s')
2061	349 FORMAT (/' Output locations in ',A,'-direction in multiples of ', &
2062	'mask_scale_',A,' predefined by array mask_',I2.2,'_',A,':'/ &
2063	13(' ',8(F8.2,',')/) )
2064	350 FORMAT (/' Output locations in ',A,'-direction: ', &
2065	'all gridpoints along ',A,'-direction (default).' )
2066	351 FORMAT (/' Output locations in ',A,'-direction in multiples of ', &
2067	'mask_scale_',A,' constructed from array mask_',I2.2,'_',A,'_loop:'/ &
2068	' loop begin:',F8.2,', end:',F8.2,', stride:',F8.2 )
2069	352 FORMAT (/' Number of output time levels allowed: ',I3 /)
2070	353 FORMAT (/' Number of output time levels allowed: unlimited' /)
2071	354 FORMAT (' Output format: ',A, ' compressed with level: ',I1/)
2072	#if defined( __dvrp_graphics )
2073	360 FORMAT (' Plot-Sequence with dvrp-software:'/ &
2074	' Output every ',F7.1,' s'/ &
2075	' Output mode: ',A/ &
2076	' Host / User: ',A,' / ',A/ &
2077	' Directory: ',A// &
2078	' The sequence contains:')
2079	361 FORMAT (/' Isosurface of "',A,'" Threshold value: ', E12.3/ &
2080	' Isosurface color: (',F4.2,',',F4.2,',',F4.2,') (R,G,B)')
2081	362 FORMAT (/' Slicer plane ',A/ &
2082	' Slicer limits: [',F6.2,',',F6.2,']')
2083	365 FORMAT (/' Groundplate color: (',F4.2,',',F4.2,',',F4.2,') (R,G,B)'/ &
2084	' Superelevation along (x,y,z): (',F4.1,',',F4.1,',',F4.1, &
2085	')'/ &
2086	' Clipping limits: from x = ',F9.1,' m to x = ',F9.1,' m'/ &
2087	' from y = ',F9.1,' m to y = ',F9.1,' m')
2088	366 FORMAT (/' Topography color: (',F4.2,',',F4.2,',',F4.2,') (R,G,B)')
2089	367 FORMAT (' Polygon reduction for topography: cluster_size = ', I1)
2090	#endif
2091	370 FORMAT (' Spectra:')
2092	371 FORMAT (' Output every ',F7.1,' s'/)
2093	372 FORMAT (' Arrays: ', 10(A5,',')/ &
2094	' Directions: ', 10(A5,',')/ &
2095	' height levels k = ', 20(I3,',')/ &
2096	' ', 20(I3,',')/ &
2097	' ', 20(I3,',')/ &
2098	' ', 20(I3,',')/ &
2099	' ', 19(I3,','),I3,'.'/ &
2100	' height levels selected for standard plot:'/ &
2101	' k = ', 20(I3,',')/ &
2102	' ', 20(I3,',')/ &
2103	' ', 20(I3,',')/ &
2104	' ', 20(I3,',')/ &
2105	' ', 19(I3,','),I3,'.'/ &
2106	' Time averaged over ', F7.1, ' s,' / &
2107	' Profiles for the time averaging are taken every ', &
2108	F6.1,' s')
2109	400 FORMAT (//' Physical quantities:'/ &
2110	' -------------------'/)
2111	410 FORMAT (' Geograph. latitude : phi = ',F4.1,' degr'/ &
2112	' Angular velocity : omega =',E10.3,' rad/s'/ &
2113	' Coriolis parameter : f = ',F9.6,' 1/s'/ &
2114	' f* = ',F9.6,' 1/s')
2115	411 FORMAT (/' Gravity : g = ',F4.1,' m/s**2')
2116	412 FORMAT (/' Reference state used in buoyancy terms: ',A)
2117	413 FORMAT (' Reference density in buoyancy terms: ',F8.3,' kg/m**3')
2118	414 FORMAT (' Reference temperature in buoyancy terms: ',F8.4,' K')
2119	415 FORMAT (/' Cloud physics parameters:'/ &
2120	' ------------------------'/)
2121	416 FORMAT (' Surface pressure : p_0 = ',F7.2,' hPa'/ &
2122	' Gas constant : R = ',F5.1,' J/(kg K)'/ &
2123	' Density of air : rho_0 =',F6.3,' kg/m**3'/ &
2124	' Specific heat cap. : c_p = ',F6.1,' J/(kg K)'/ &
2125	' Vapourization heat : L_v =',E9.2,' J/kg')
2126	417 FORMAT (' Geograph. longitude : lambda = ',F4.1,' degr')
2127	418 FORMAT (/' Day of the year at model start : day_init = ',I3 &
2128	/' UTC time at model start : time_utc_init = ',F7.1' s')
2129	420 FORMAT (/' Characteristic levels of the initial temperature profile:'// &
2130	' Height: ',A,' m'/ &
2131	' Temperature: ',A,' K'/ &
2132	' Gradient: ',A,' K/100m'/ &
2133	' Gridpoint: ',A)
2134	421 FORMAT (/' Characteristic levels of the initial humidity profile:'// &
2135	' Height: ',A,' m'/ &
2136	' Humidity: ',A,' kg/kg'/ &
2137	' Gradient: ',A,' (kg/kg)/100m'/ &
2138	' Gridpoint: ',A)
2139	422 FORMAT (/' Characteristic levels of the initial scalar profile:'// &
2140	' Height: ',A,' m'/ &
2141	' Scalar concentration: ',A,' kg/m**3'/ &
2142	' Gradient: ',A,' (kg/m**3)/100m'/ &
2143	' Gridpoint: ',A)
2144	423 FORMAT (/' Characteristic levels of the geo. wind component ug:'// &
2145	' Height: ',A,' m'/ &
2146	' ug: ',A,' m/s'/ &
2147	' Gradient: ',A,' 1/100s'/ &
2148	' Gridpoint: ',A)
2149	424 FORMAT (/' Characteristic levels of the geo. wind component vg:'// &
2150	' Height: ',A,' m'/ &
2151	' vg: ',A,' m/s'/ &
2152	' Gradient: ',A,' 1/100s'/ &
2153	' Gridpoint: ',A)
2154	425 FORMAT (/' Characteristic levels of the initial salinity profile:'// &
2155	' Height: ',A,' m'/ &
2156	' Salinity: ',A,' psu'/ &
2157	' Gradient: ',A,' psu/100m'/ &
2158	' Gridpoint: ',A)
2159	426 FORMAT (/' Characteristic levels of the subsidence/ascent profile:'// &
2160	' Height: ',A,' m'/ &
2161	' w_subs: ',A,' m/s'/ &
2162	' Gradient: ',A,' (m/s)/100m'/ &
2163	' Gridpoint: ',A)
2164	427 FORMAT (/' Initial wind profiles (u,v) are interpolated from given'// &
2165	' profiles')
2166	428 FORMAT (/' Initial profiles (u, v, pt, q) are taken from file '/ &
2167	' NUDGING_DATA')
2168	430 FORMAT (//' Cloud physics quantities / methods:'/ &
2169	' ----------------------------------'/)
2170	431 FORMAT (' Humidity is treated as purely passive scalar (no condensati', &
2171	'on)')
2172	432 FORMAT (' Bulk scheme with liquid water potential temperature and'/ &
2173	' total water content is used.'/ &
2174	' Condensation is parameterized via 0% - or 100% scheme.')
2175	433 FORMAT (' Cloud droplets treated explicitly using the Lagrangian part', &
2176	'icle model')
2177	434 FORMAT (' Curvature and solution effecs are considered for growth of', &
2178	' droplets < 1.0E-6 m')
2179	435 FORMAT (' Droplet collision is handled by ',A,'-kernel')
2180	436 FORMAT (' Fast kernel with fixed radius- and dissipation classes ', &
2181	'are used'/ &
2182	' number of radius classes: ',I3,' interval ', &
2183	'[1.0E-6,2.0E-4] m'/ &
2184	' number of dissipation classes: ',I2,' interval ', &
2185	'[0,1000] cm2/s3')
2186	437 FORMAT (' Droplet collision is switched off')
2187	450 FORMAT (//' LES / Turbulence quantities:'/ &
2188	' ---------------------------'/)
2189	451 FORMAT (' Diffusion coefficients are constant:'/ &
2190	' Km = ',F6.2,' m2/s Kh = ',F6.2,' m2/s Pr = ',F5.2)
2191	453 FORMAT (' Mixing length is limited to',F5.2,' * z')
2192	454 FORMAT (' TKE is not allowed to fall below ',E9.2,' (m/s)**2')
2193	455 FORMAT (' initial TKE is prescribed as ',E9.2,' (m/s)**2')
2194	470 FORMAT (//' Actions during the simulation:'/ &
2195	' -----------------------------'/)
2196	471 FORMAT (' Disturbance impulse (u,v) every : ',F6.2,' s'/ &
2197	' Disturbance amplitude : ',F5.2, ' m/s'/ &
2198	' Lower disturbance level : ',F8.2,' m (GP ',I4,')'/ &
2199	' Upper disturbance level : ',F8.2,' m (GP ',I4,')')
2200	472 FORMAT (' Disturbances continued during the run from i/j =',I4, &
2201	' to i/j =',I4)
2202	473 FORMAT (' Disturbances cease as soon as the disturbance energy exceeds',&
2203	F6.3, ' m2/s2')
2204	474 FORMAT (' Random number generator used : ',A/)
2205	475 FORMAT (' The surface temperature is increased (or decreased, ', &
2206	'respectively, if'/ &
2207	' the value is negative) by ',F5.2,' K at the beginning of the',&
2208	' 3D-simulation'/)
2209	476 FORMAT (' The surface humidity is increased (or decreased, ',&
2210	'respectively, if the'/ &
2211	' value is negative) by ',E8.1,' kg/kg at the beginning of', &
2212	' the 3D-simulation'/)
2213	477 FORMAT (' The scalar value is increased at the surface (or decreased, ',&
2214	'respectively, if the'/ &
2215	' value is negative) by ',E8.1,' kg/m**3 at the beginning of', &
2216	' the 3D-simulation'/)
2217	480 FORMAT (' Particles:'/ &
2218	' ---------'// &
2219	' Particle advection is active (switched on at t = ', F7.1, &
2220	' s)'/ &
2221	' Start of new particle generations every ',F6.1,' s'/ &
2222	' Boundary conditions: left/right: ', A, ' north/south: ', A/&
2223	' bottom: ', A, ' top: ', A/&
2224	' Maximum particle age: ',F9.1,' s'/ &
2225	' Advection stopped at t = ',F9.1,' s'/)
2226	481 FORMAT (' Particles have random start positions'/)
2227	482 FORMAT (' Particles are advected only horizontally'/)
2228	485 FORMAT (' Particle data are written on file every ', F9.1, ' s')
2229	486 FORMAT (' Particle statistics are written on file'/)
2230	487 FORMAT (' Number of particle groups: ',I2/)
2231	488 FORMAT (' SGS velocity components are used for particle advection'/ &
2232	' minimum timestep for advection:', F8.5/)
2233	489 FORMAT (' Number of particles simultaneously released at each ', &
2234	'point: ', I5/)
2235	490 FORMAT (' Particle group ',I2,':'/ &
2236	' Particle radius: ',E10.3, 'm')
2237	491 FORMAT (' Particle inertia is activated'/ &
2238	' density_ratio (rho_fluid/rho_particle) =',F6.3/)
2239	492 FORMAT (' Particles are advected only passively (no inertia)'/)
2240	493 FORMAT (' Boundaries of particle source: x:',F8.1,' - ',F8.1,' m'/&
2241	' y:',F8.1,' - ',F8.1,' m'/&
2242	' z:',F8.1,' - ',F8.1,' m'/&
2243	' Particle distances: dx = ',F8.1,' m dy = ',F8.1, &
2244	' m dz = ',F8.1,' m'/)
2245	494 FORMAT (' Output of particle time series in NetCDF format every ', &
2246	F8.2,' s'/)
2247	495 FORMAT (' Number of particles in total domain: ',I10/)
2248	496 FORMAT (' Initial vertical particle positions are interpreted ', &
2249	'as relative to the given topography')
2250	500 FORMAT (//' 1D-Model parameters:'/ &
2251	' -------------------'// &
2252	' Simulation time: ',F8.1,' s'/ &
2253	' Run-controll output every: ',F8.1,' s'/ &
2254	' Vertical profile output every: ',F8.1,' s'/ &
2255	' Mixing length calculation: ',A/ &
2256	' Dissipation calculation: ',A/)
2257	502 FORMAT (' Damping layer starts from ',F7.1,' m (GP ',I4,')'/)
2258	503 FORMAT (' --> Momentum advection via Wicker-Skamarock-Scheme 5th order')
2259	504 FORMAT (' --> Scalar advection via Wicker-Skamarock-Scheme 5th order')
2260	505 FORMAT (' Precipitation parameterization via Seifert-Beheng-Scheme')
2261	506 FORMAT (' Cloud water sedimentation parameterization via Stokes law')
2262	507 FORMAT (' Turbulence effects on precipitation process')
2263	508 FORMAT (' Ventilation effects on evaporation of rain drops')
2264	509 FORMAT (' Slope limiter used for sedimentation process')
2265	510 FORMAT (' Droplet density : N_c = ',F6.1,' 1/cm**3')
2266	511 FORMAT (' Sedimentation Courant number: '/&
2267	' C_s =',F4.1,' ')
2268	512 FORMAT (/' Date: ',A8,6X,'Run: ',A20/ &
2269	' Time: ',A8,6X,'Run-No.: ',I2.2/ &
2270	' Run on host: ',A10,6X,'En-No.: ',I2.2)
2271	513 FORMAT (' --> Scalar advection via Wicker-Skamarock-Scheme 5th order ' // &
2272	'+ monotonic adjustment')
2273	600 FORMAT (/' Nesting informations:'/ &
2274	' --------------------'/ &
2275	' Nesting mode: ',A/ &
2276	' Nesting-datatransfer mode: ',A// &
2277	' Nest id parent number lower left coordinates name'/ &
2278	' (*=me) id of PEs x (m) y (m)' )
2279	601 FORMAT (2X,A1,1X,I2.2,6X,I2.2,5X,I5,5X,F8.2,2X,F8.2,5X,A)
2280
2281	END SUBROUTINE header

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |