Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |