Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |