Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |