Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |