Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |