Home

Context Navigation

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

Last change on this file since 2322 was 2320, checked in by suehring, 7 years ago
large-scale forcing and nudging modularized
Property svn:keywords set to `Id`
File size: 88.8 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |