Home

Context Navigation

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

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |