Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |