Home

Context Navigation

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

Last change on this file since 1558 was 1558, checked in by suehring, 10 years ago
last commit documented
Property svn:keywords set to `Id`
File size: 89.7 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |