Home

Context Navigation

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

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |