Home

Context Navigation

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

Last change on this file since 1248 was 1242, checked in by heinze, 11 years ago
Last commmit documented
Property svn:keywords set to `Id`
File size: 80.1 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |