Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |