Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |