Home

Context Navigation

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

Last change on this file since 2186 was 2119, checked in by raasch, 8 years ago
last commit documented
Property svn:keywords set to `Id`
File size: 88.5 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |