Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |