Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |