Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |