Home

Context Navigation

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

Last change on this file since 1322 was 1322, checked in by raasch, 10 years ago
REAL functions and a lot of REAL constants provided with KIND-attribute, some small bugfixes
Property svn:keywords set to `Id`
File size: 82.2 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |