Home

Context Navigation

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

Last change on this file since 2768 was 2746, checked in by suehring, 6 years ago
Read information from statitic driver for resolved vegetation independently from land- or urban-surface model
Property svn:keywords set to `Id`
File size: 91.0 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |