Home

Context Navigation

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

Last change on this file since 1323 was 1323, checked in by raasch, 10 years ago
last commit documented
Property svn:keywords set to `Id`
File size: 82.2 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |