Home

Context Navigation

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

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |