Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |