Home

Context Navigation

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

Last change on this file since 1590 was 1590, checked in by maronga, 9 years ago
bugfix in definition of character arrays within the LSM and radiation code
Property svn:keywords set to `Id`
File size: 91.5 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	! Bugfix: Added TRIM statements for character strings for LSM and radiation code
23	!
24	! Former revisions:
25	! -----------------
26	! $Id: header.f90 1590 2015-05-08 13:56:27Z maronga $
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 ) TRIM( veg_type_name(veg_type) ), &
925	TRIM (soil_type_name(soil_type) )
926	WRITE( io, 439 ) TRIM( coordinates ), TRIM( temperatures ), &
927	TRIM( gradients ), TRIM( roots ), TRIM( slices )
928
929
930	ENDIF
931
932	IF ( radiation ) THEN
933	!
934	!-- Write radiation model header
935	WRITE( io, 444 )
936
937	IF ( radiation_scheme == "constant" ) THEN
938	WRITE( io, 445 ) net_radiation
939	ELSEIF ( radiation_scheme == "clear-sky" ) THEN
940	WRITE( io, 446 )
941	ELSEIF ( radiation_scheme == "rrtmg" ) THEN
942	WRITE( io, 447 )
943	IF ( .NOT. lw_radiation ) WRITE( io, 458 )
944	IF ( .NOT. sw_radiation ) WRITE( io, 459 )
945	ENDIF
946
947	IF ( albedo_type == 0 ) THEN
948	WRITE( io, 448 ) albedo
949	ELSE
950	WRITE( io, 456 ) TRIM( albedo_type_name(albedo_type) )
951	ENDIF
952	IF ( constant_albedo ) THEN
953	WRITE( io, 457 )
954	ENDIF
955	WRITE( io, 449 ) dt_radiation
956	ENDIF
957
958
959	!
960	!-- Boundary conditions
961	IF ( ibc_p_b == 0 ) THEN
962	runten = 'p(0) = 0 \|'
963	ELSEIF ( ibc_p_b == 1 ) THEN
964	runten = 'p(0) = p(1) \|'
965	ENDIF
966	IF ( ibc_p_t == 0 ) THEN
967	roben = 'p(nzt+1) = 0 \|'
968	ELSE
969	roben = 'p(nzt+1) = p(nzt) \|'
970	ENDIF
971
972	IF ( ibc_uv_b == 0 ) THEN
973	runten = TRIM( runten ) // ' uv(0) = -uv(1) \|'
974	ELSE
975	runten = TRIM( runten ) // ' uv(0) = uv(1) \|'
976	ENDIF
977	IF ( TRIM( bc_uv_t ) == 'dirichlet_0' ) THEN
978	roben = TRIM( roben ) // ' uv(nzt+1) = 0 \|'
979	ELSEIF ( ibc_uv_t == 0 ) THEN
980	roben = TRIM( roben ) // ' uv(nzt+1) = ug(nzt+1), vg(nzt+1) \|'
981	ELSE
982	roben = TRIM( roben ) // ' uv(nzt+1) = uv(nzt) \|'
983	ENDIF
984
985	IF ( ibc_pt_b == 0 ) THEN
986	IF ( land_surface ) THEN
987	runten = TRIM( runten ) // ' pt(0) = from soil model'
988	ELSE
989	runten = TRIM( runten ) // ' pt(0) = pt_surface'
990	ENDIF
991	ELSEIF ( ibc_pt_b == 1 ) THEN
992	runten = TRIM( runten ) // ' pt(0) = pt(1)'
993	ELSEIF ( ibc_pt_b == 2 ) THEN
994	runten = TRIM( runten ) // ' pt(0) = from coupled model'
995	ENDIF
996	IF ( ibc_pt_t == 0 ) THEN
997	roben = TRIM( roben ) // ' pt(nzt+1) = pt_top'
998	ELSEIF( ibc_pt_t == 1 ) THEN
999	roben = TRIM( roben ) // ' pt(nzt+1) = pt(nzt)'
1000	ELSEIF( ibc_pt_t == 2 ) THEN
1001	roben = TRIM( roben ) // ' pt(nzt+1) = pt(nzt) + dpt/dz_ini'
1002
1003	ENDIF
1004
1005	WRITE ( io, 300 ) runten, roben
1006
1007	IF ( .NOT. constant_diffusion ) THEN
1008	IF ( ibc_e_b == 1 ) THEN
1009	runten = 'e(0) = e(1)'
1010	ELSE
1011	runten = 'e(0) = e(1) = (u/0.1)*2'
1012	ENDIF
1013	roben = 'e(nzt+1) = e(nzt) = e(nzt-1)'
1014
1015	WRITE ( io, 301 ) 'e', runten, roben
1016
1017	ENDIF
1018
1019	IF ( ocean ) THEN
1020	runten = 'sa(0) = sa(1)'
1021	IF ( ibc_sa_t == 0 ) THEN
1022	roben = 'sa(nzt+1) = sa_surface'
1023	ELSE
1024	roben = 'sa(nzt+1) = sa(nzt)'
1025	ENDIF
1026	WRITE ( io, 301 ) 'sa', runten, roben
1027	ENDIF
1028
1029	IF ( humidity ) THEN
1030	IF ( ibc_q_b == 0 ) THEN
1031	IF ( land_surface ) THEN
1032	runten = 'q(0) = from soil model'
1033	ELSE
1034	runten = 'q(0) = q_surface'
1035	ENDIF
1036
1037	ELSE
1038	runten = 'q(0) = q(1)'
1039	ENDIF
1040	IF ( ibc_q_t == 0 ) THEN
1041	roben = 'q(nzt) = q_top'
1042	ELSE
1043	roben = 'q(nzt) = q(nzt-1) + dq/dz'
1044	ENDIF
1045	WRITE ( io, 301 ) 'q', runten, roben
1046	ENDIF
1047
1048	IF ( passive_scalar ) THEN
1049	IF ( ibc_q_b == 0 ) THEN
1050	runten = 's(0) = s_surface'
1051	ELSE
1052	runten = 's(0) = s(1)'
1053	ENDIF
1054	IF ( ibc_q_t == 0 ) THEN
1055	roben = 's(nzt) = s_top'
1056	ELSE
1057	roben = 's(nzt) = s(nzt-1) + ds/dz'
1058	ENDIF
1059	WRITE ( io, 301 ) 's', runten, roben
1060	ENDIF
1061
1062	IF ( use_surface_fluxes ) THEN
1063	WRITE ( io, 303 )
1064	IF ( constant_heatflux ) THEN
1065	IF ( large_scale_forcing .AND. lsf_surf ) THEN
1066	WRITE ( io, 306 ) shf(0,0)
1067	ELSE
1068	WRITE ( io, 306 ) surface_heatflux
1069	ENDIF
1070	IF ( random_heatflux ) WRITE ( io, 307 )
1071	ENDIF
1072	IF ( humidity .AND. constant_waterflux ) THEN
1073	IF ( large_scale_forcing .AND. lsf_surf ) THEN
1074	WRITE ( io, 311 ) qsws(0,0)
1075	ELSE
1076	WRITE ( io, 311 ) surface_waterflux
1077	ENDIF
1078	ENDIF
1079	IF ( passive_scalar .AND. constant_waterflux ) THEN
1080	WRITE ( io, 313 ) surface_waterflux
1081	ENDIF
1082	ENDIF
1083
1084	IF ( use_top_fluxes ) THEN
1085	WRITE ( io, 304 )
1086	IF ( coupling_mode == 'uncoupled' ) THEN
1087	WRITE ( io, 320 ) top_momentumflux_u, top_momentumflux_v
1088	IF ( constant_top_heatflux ) THEN
1089	WRITE ( io, 306 ) top_heatflux
1090	ENDIF
1091	ELSEIF ( coupling_mode == 'ocean_to_atmosphere' ) THEN
1092	WRITE ( io, 316 )
1093	ENDIF
1094	IF ( ocean .AND. constant_top_salinityflux ) THEN
1095	WRITE ( io, 309 ) top_salinityflux
1096	ENDIF
1097	IF ( humidity .OR. passive_scalar ) THEN
1098	WRITE ( io, 315 )
1099	ENDIF
1100	ENDIF
1101
1102	IF ( prandtl_layer ) THEN
1103	WRITE ( io, 305 ) (zu(1)-zu(0)), roughness_length, &
1104	z0h_factor*roughness_length, kappa, &
1105	rif_min, rif_max
1106	IF ( .NOT. constant_heatflux ) WRITE ( io, 308 )
1107	IF ( humidity .AND. .NOT. constant_waterflux ) THEN
1108	WRITE ( io, 312 )
1109	ENDIF
1110	IF ( passive_scalar .AND. .NOT. constant_waterflux ) THEN
1111	WRITE ( io, 314 )
1112	ENDIF
1113	ELSE
1114	IF ( INDEX(initializing_actions, 'set_1d-model_profiles') /= 0 ) THEN
1115	WRITE ( io, 310 ) rif_min, rif_max
1116	ENDIF
1117	ENDIF
1118
1119	WRITE ( io, 317 ) bc_lr, bc_ns
1120	IF ( .NOT. bc_lr_cyc .OR. .NOT. bc_ns_cyc ) THEN
1121	WRITE ( io, 318 ) use_cmax, pt_damping_width, pt_damping_factor
1122	IF ( turbulent_inflow ) THEN
1123	IF ( .NOT. recycling_yshift ) THEN
1124	WRITE ( io, 319 ) recycling_width, recycling_plane, &
1125	inflow_damping_height, inflow_damping_width
1126	ELSE
1127	WRITE ( io, 322 ) recycling_width, recycling_plane, &
1128	inflow_damping_height, inflow_damping_width
1129	END IF
1130	ENDIF
1131	ENDIF
1132
1133	!
1134	!-- Initial Profiles
1135	WRITE ( io, 321 )
1136	!
1137	!-- Initial wind profiles
1138	IF ( u_profile(1) /= 9999999.9_wp ) WRITE ( io, 427 )
1139
1140	!
1141	!-- Initial temperature profile
1142	!-- Building output strings, starting with surface temperature
1143	WRITE ( temperatures, '(F6.2)' ) pt_surface
1144	gradients = '------'
1145	slices = ' 0'
1146	coordinates = ' 0.0'
1147	i = 1
1148	DO WHILE ( pt_vertical_gradient_level_ind(i) /= -9999 )
1149
1150	WRITE (coor_chr,'(F7.2)') pt_init(pt_vertical_gradient_level_ind(i))
1151	temperatures = TRIM( temperatures ) // ' ' // TRIM( coor_chr )
1152
1153	WRITE (coor_chr,'(F7.2)') pt_vertical_gradient(i)
1154	gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
1155
1156	WRITE (coor_chr,'(I7)') pt_vertical_gradient_level_ind(i)
1157	slices = TRIM( slices ) // ' ' // TRIM( coor_chr )
1158
1159	WRITE (coor_chr,'(F7.1)') pt_vertical_gradient_level(i)
1160	coordinates = TRIM( coordinates ) // ' ' // TRIM( coor_chr )
1161
1162	IF ( i == 10 ) THEN
1163	EXIT
1164	ELSE
1165	i = i + 1
1166	ENDIF
1167
1168	ENDDO
1169
1170	IF ( .NOT. nudging ) THEN
1171	WRITE ( io, 420 ) TRIM( coordinates ), TRIM( temperatures ), &
1172	TRIM( gradients ), TRIM( slices )
1173	ELSE
1174	WRITE ( io, 428 )
1175	ENDIF
1176
1177	!
1178	!-- Initial humidity profile
1179	!-- Building output strings, starting with surface humidity
1180	IF ( humidity .OR. passive_scalar ) THEN
1181	WRITE ( temperatures, '(E8.1)' ) q_surface
1182	gradients = '--------'
1183	slices = ' 0'
1184	coordinates = ' 0.0'
1185	i = 1
1186	DO WHILE ( q_vertical_gradient_level_ind(i) /= -9999 )
1187
1188	WRITE (coor_chr,'(E8.1,4X)') q_init(q_vertical_gradient_level_ind(i))
1189	temperatures = TRIM( temperatures ) // ' ' // TRIM( coor_chr )
1190
1191	WRITE (coor_chr,'(E8.1,4X)') q_vertical_gradient(i)
1192	gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
1193
1194	WRITE (coor_chr,'(I8,4X)') q_vertical_gradient_level_ind(i)
1195	slices = TRIM( slices ) // ' ' // TRIM( coor_chr )
1196
1197	WRITE (coor_chr,'(F8.1,4X)') q_vertical_gradient_level(i)
1198	coordinates = TRIM( coordinates ) // ' ' // TRIM( coor_chr )
1199
1200	IF ( i == 10 ) THEN
1201	EXIT
1202	ELSE
1203	i = i + 1
1204	ENDIF
1205
1206	ENDDO
1207
1208	IF ( humidity ) THEN
1209	IF ( .NOT. nudging ) THEN
1210	WRITE ( io, 421 ) TRIM( coordinates ), TRIM( temperatures ), &
1211	TRIM( gradients ), TRIM( slices )
1212	ENDIF
1213	ELSE
1214	WRITE ( io, 422 ) TRIM( coordinates ), TRIM( temperatures ), &
1215	TRIM( gradients ), TRIM( slices )
1216	ENDIF
1217	ENDIF
1218
1219	!
1220	!-- Initial salinity profile
1221	!-- Building output strings, starting with surface salinity
1222	IF ( ocean ) THEN
1223	WRITE ( temperatures, '(F6.2)' ) sa_surface
1224	gradients = '------'
1225	slices = ' 0'
1226	coordinates = ' 0.0'
1227	i = 1
1228	DO WHILE ( sa_vertical_gradient_level_ind(i) /= -9999 )
1229
1230	WRITE (coor_chr,'(F7.2)') sa_init(sa_vertical_gradient_level_ind(i))
1231	temperatures = TRIM( temperatures ) // ' ' // TRIM( coor_chr )
1232
1233	WRITE (coor_chr,'(F7.2)') sa_vertical_gradient(i)
1234	gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
1235
1236	WRITE (coor_chr,'(I7)') sa_vertical_gradient_level_ind(i)
1237	slices = TRIM( slices ) // ' ' // TRIM( coor_chr )
1238
1239	WRITE (coor_chr,'(F7.1)') sa_vertical_gradient_level(i)
1240	coordinates = TRIM( coordinates ) // ' ' // TRIM( coor_chr )
1241
1242	IF ( i == 10 ) THEN
1243	EXIT
1244	ELSE
1245	i = i + 1
1246	ENDIF
1247
1248	ENDDO
1249
1250	WRITE ( io, 425 ) TRIM( coordinates ), TRIM( temperatures ), &
1251	TRIM( gradients ), TRIM( slices )
1252	ENDIF
1253
1254
1255	!
1256	!-- Listing of 1D-profiles
1257	WRITE ( io, 325 ) dt_dopr_listing
1258	IF ( averaging_interval_pr /= 0.0_wp ) THEN
1259	WRITE ( io, 326 ) averaging_interval_pr, dt_averaging_input_pr
1260	ENDIF
1261
1262	!
1263	!-- DATA output
1264	WRITE ( io, 330 )
1265	IF ( averaging_interval_pr /= 0.0_wp ) THEN
1266	WRITE ( io, 326 ) averaging_interval_pr, dt_averaging_input_pr
1267	ENDIF
1268
1269	!
1270	!-- 1D-profiles
1271	dopr_chr = 'Profile:'
1272	IF ( dopr_n /= 0 ) THEN
1273	WRITE ( io, 331 )
1274
1275	output_format = ''
1276	output_format = output_format_netcdf
1277	WRITE ( io, 344 ) output_format
1278
1279	DO i = 1, dopr_n
1280	dopr_chr = TRIM( dopr_chr ) // ' ' // TRIM( data_output_pr(i) ) // ','
1281	IF ( LEN_TRIM( dopr_chr ) >= 60 ) THEN
1282	WRITE ( io, 332 ) dopr_chr
1283	dopr_chr = ' :'
1284	ENDIF
1285	ENDDO
1286
1287	IF ( dopr_chr /= '' ) THEN
1288	WRITE ( io, 332 ) dopr_chr
1289	ENDIF
1290	WRITE ( io, 333 ) dt_dopr, averaging_interval_pr, dt_averaging_input_pr
1291	IF ( skip_time_dopr /= 0.0_wp ) WRITE ( io, 339 ) skip_time_dopr
1292	ENDIF
1293
1294	!
1295	!-- 2D-arrays
1296	DO av = 0, 1
1297
1298	i = 1
1299	do2d_xy = ''
1300	do2d_xz = ''
1301	do2d_yz = ''
1302	DO WHILE ( do2d(av,i) /= ' ' )
1303
1304	l = MAX( 2, LEN_TRIM( do2d(av,i) ) )
1305	do2d_mode = do2d(av,i)(l-1:l)
1306
1307	SELECT CASE ( do2d_mode )
1308	CASE ( 'xy' )
1309	ll = LEN_TRIM( do2d_xy )
1310	do2d_xy = do2d_xy(1:ll) // ' ' // do2d(av,i)(1:l-3) // ','
1311	CASE ( 'xz' )
1312	ll = LEN_TRIM( do2d_xz )
1313	do2d_xz = do2d_xz(1:ll) // ' ' // do2d(av,i)(1:l-3) // ','
1314	CASE ( 'yz' )
1315	ll = LEN_TRIM( do2d_yz )
1316	do2d_yz = do2d_yz(1:ll) // ' ' // do2d(av,i)(1:l-3) // ','
1317	END SELECT
1318
1319	i = i + 1
1320
1321	ENDDO
1322
1323	IF ( ( ( do2d_xy /= '' .AND. section(1,1) /= -9999 ) .OR. &
1324	( do2d_xz /= '' .AND. section(1,2) /= -9999 ) .OR. &
1325	( do2d_yz /= '' .AND. section(1,3) /= -9999 ) ) ) THEN
1326
1327	IF ( av == 0 ) THEN
1328	WRITE ( io, 334 ) ''
1329	ELSE
1330	WRITE ( io, 334 ) '(time-averaged)'
1331	ENDIF
1332
1333	IF ( do2d_at_begin ) THEN
1334	begin_chr = 'and at the start'
1335	ELSE
1336	begin_chr = ''
1337	ENDIF
1338
1339	output_format = ''
1340	output_format = output_format_netcdf
1341	WRITE ( io, 344 ) output_format
1342
1343	IF ( do2d_xy /= '' .AND. section(1,1) /= -9999 ) THEN
1344	i = 1
1345	slices = '/'
1346	coordinates = '/'
1347	!
1348	!-- Building strings with index and coordinate information of the
1349	!-- slices
1350	DO WHILE ( section(i,1) /= -9999 )
1351
1352	WRITE (section_chr,'(I5)') section(i,1)
1353	section_chr = ADJUSTL( section_chr )
1354	slices = TRIM( slices ) // TRIM( section_chr ) // '/'
1355
1356	IF ( section(i,1) == -1 ) THEN
1357	WRITE (coor_chr,'(F10.1)') -1.0_wp
1358	ELSE
1359	WRITE (coor_chr,'(F10.1)') zu(section(i,1))
1360	ENDIF
1361	coor_chr = ADJUSTL( coor_chr )
1362	coordinates = TRIM( coordinates ) // TRIM( coor_chr ) // '/'
1363
1364	i = i + 1
1365	ENDDO
1366	IF ( av == 0 ) THEN
1367	WRITE ( io, 335 ) 'XY', do2d_xy, dt_do2d_xy, &
1368	TRIM( begin_chr ), 'k', TRIM( slices ), &
1369	TRIM( coordinates )
1370	IF ( skip_time_do2d_xy /= 0.0_wp ) THEN
1371	WRITE ( io, 339 ) skip_time_do2d_xy
1372	ENDIF
1373	ELSE
1374	WRITE ( io, 342 ) 'XY', do2d_xy, dt_data_output_av, &
1375	TRIM( begin_chr ), averaging_interval, &
1376	dt_averaging_input, 'k', TRIM( slices ), &
1377	TRIM( coordinates )
1378	IF ( skip_time_data_output_av /= 0.0_wp ) THEN
1379	WRITE ( io, 339 ) skip_time_data_output_av
1380	ENDIF
1381	ENDIF
1382	IF ( netcdf_data_format > 4 ) THEN
1383	WRITE ( io, 352 ) ntdim_2d_xy(av)
1384	ELSE
1385	WRITE ( io, 353 )
1386	ENDIF
1387	ENDIF
1388
1389	IF ( do2d_xz /= '' .AND. section(1,2) /= -9999 ) THEN
1390	i = 1
1391	slices = '/'
1392	coordinates = '/'
1393	!
1394	!-- Building strings with index and coordinate information of the
1395	!-- slices
1396	DO WHILE ( section(i,2) /= -9999 )
1397
1398	WRITE (section_chr,'(I5)') section(i,2)
1399	section_chr = ADJUSTL( section_chr )
1400	slices = TRIM( slices ) // TRIM( section_chr ) // '/'
1401
1402	WRITE (coor_chr,'(F10.1)') section(i,2) * dy
1403	coor_chr = ADJUSTL( coor_chr )
1404	coordinates = TRIM( coordinates ) // TRIM( coor_chr ) // '/'
1405
1406	i = i + 1
1407	ENDDO
1408	IF ( av == 0 ) THEN
1409	WRITE ( io, 335 ) 'XZ', do2d_xz, dt_do2d_xz, &
1410	TRIM( begin_chr ), 'j', TRIM( slices ), &
1411	TRIM( coordinates )
1412	IF ( skip_time_do2d_xz /= 0.0_wp ) THEN
1413	WRITE ( io, 339 ) skip_time_do2d_xz
1414	ENDIF
1415	ELSE
1416	WRITE ( io, 342 ) 'XZ', do2d_xz, dt_data_output_av, &
1417	TRIM( begin_chr ), averaging_interval, &
1418	dt_averaging_input, 'j', TRIM( slices ), &
1419	TRIM( coordinates )
1420	IF ( skip_time_data_output_av /= 0.0_wp ) THEN
1421	WRITE ( io, 339 ) skip_time_data_output_av
1422	ENDIF
1423	ENDIF
1424	IF ( netcdf_data_format > 4 ) THEN
1425	WRITE ( io, 352 ) ntdim_2d_xz(av)
1426	ELSE
1427	WRITE ( io, 353 )
1428	ENDIF
1429	ENDIF
1430
1431	IF ( do2d_yz /= '' .AND. section(1,3) /= -9999 ) THEN
1432	i = 1
1433	slices = '/'
1434	coordinates = '/'
1435	!
1436	!-- Building strings with index and coordinate information of the
1437	!-- slices
1438	DO WHILE ( section(i,3) /= -9999 )
1439
1440	WRITE (section_chr,'(I5)') section(i,3)
1441	section_chr = ADJUSTL( section_chr )
1442	slices = TRIM( slices ) // TRIM( section_chr ) // '/'
1443
1444	WRITE (coor_chr,'(F10.1)') section(i,3) * dx
1445	coor_chr = ADJUSTL( coor_chr )
1446	coordinates = TRIM( coordinates ) // TRIM( coor_chr ) // '/'
1447
1448	i = i + 1
1449	ENDDO
1450	IF ( av == 0 ) THEN
1451	WRITE ( io, 335 ) 'YZ', do2d_yz, dt_do2d_yz, &
1452	TRIM( begin_chr ), 'i', TRIM( slices ), &
1453	TRIM( coordinates )
1454	IF ( skip_time_do2d_yz /= 0.0_wp ) THEN
1455	WRITE ( io, 339 ) skip_time_do2d_yz
1456	ENDIF
1457	ELSE
1458	WRITE ( io, 342 ) 'YZ', do2d_yz, dt_data_output_av, &
1459	TRIM( begin_chr ), averaging_interval, &
1460	dt_averaging_input, 'i', TRIM( slices ), &
1461	TRIM( coordinates )
1462	IF ( skip_time_data_output_av /= 0.0_wp ) THEN
1463	WRITE ( io, 339 ) skip_time_data_output_av
1464	ENDIF
1465	ENDIF
1466	IF ( netcdf_data_format > 4 ) THEN
1467	WRITE ( io, 352 ) ntdim_2d_yz(av)
1468	ELSE
1469	WRITE ( io, 353 )
1470	ENDIF
1471	ENDIF
1472
1473	ENDIF
1474
1475	ENDDO
1476
1477	!
1478	!-- 3d-arrays
1479	DO av = 0, 1
1480
1481	i = 1
1482	do3d_chr = ''
1483	DO WHILE ( do3d(av,i) /= ' ' )
1484
1485	do3d_chr = TRIM( do3d_chr ) // ' ' // TRIM( do3d(av,i) ) // ','
1486	i = i + 1
1487
1488	ENDDO
1489
1490	IF ( do3d_chr /= '' ) THEN
1491	IF ( av == 0 ) THEN
1492	WRITE ( io, 336 ) ''
1493	ELSE
1494	WRITE ( io, 336 ) '(time-averaged)'
1495	ENDIF
1496
1497	output_format = output_format_netcdf
1498	WRITE ( io, 344 ) output_format
1499
1500	IF ( do3d_at_begin ) THEN
1501	begin_chr = 'and at the start'
1502	ELSE
1503	begin_chr = ''
1504	ENDIF
1505	IF ( av == 0 ) THEN
1506	WRITE ( io, 337 ) do3d_chr, dt_do3d, TRIM( begin_chr ), &
1507	zu(nz_do3d), nz_do3d
1508	ELSE
1509	WRITE ( io, 343 ) do3d_chr, dt_data_output_av, &
1510	TRIM( begin_chr ), averaging_interval, &
1511	dt_averaging_input, zu(nz_do3d), nz_do3d
1512	ENDIF
1513
1514	IF ( netcdf_data_format > 4 ) THEN
1515	WRITE ( io, 352 ) ntdim_3d(av)
1516	ELSE
1517	WRITE ( io, 353 )
1518	ENDIF
1519
1520	IF ( av == 0 ) THEN
1521	IF ( skip_time_do3d /= 0.0_wp ) THEN
1522	WRITE ( io, 339 ) skip_time_do3d
1523	ENDIF
1524	ELSE
1525	IF ( skip_time_data_output_av /= 0.0_wp ) THEN
1526	WRITE ( io, 339 ) skip_time_data_output_av
1527	ENDIF
1528	ENDIF
1529
1530	ENDIF
1531
1532	ENDDO
1533
1534	!
1535	!-- masked arrays
1536	IF ( masks > 0 ) WRITE ( io, 345 ) &
1537	mask_scale_x, mask_scale_y, mask_scale_z
1538	DO mid = 1, masks
1539	DO av = 0, 1
1540
1541	i = 1
1542	domask_chr = ''
1543	DO WHILE ( domask(mid,av,i) /= ' ' )
1544	domask_chr = TRIM( domask_chr ) // ' ' // &
1545	TRIM( domask(mid,av,i) ) // ','
1546	i = i + 1
1547	ENDDO
1548
1549	IF ( domask_chr /= '' ) THEN
1550	IF ( av == 0 ) THEN
1551	WRITE ( io, 346 ) '', mid
1552	ELSE
1553	WRITE ( io, 346 ) ' (time-averaged)', mid
1554	ENDIF
1555
1556	output_format = output_format_netcdf
1557	!-- Parallel output not implemented for mask data, hence
1558	!-- output_format must be adjusted.
1559	IF ( netcdf_data_format == 5 ) output_format = 'netCDF4/HDF5'
1560	IF ( netcdf_data_format == 6 ) output_format = 'netCDF4/HDF5 classic'
1561	WRITE ( io, 344 ) output_format
1562
1563	IF ( av == 0 ) THEN
1564	WRITE ( io, 347 ) domask_chr, dt_domask(mid)
1565	ELSE
1566	WRITE ( io, 348 ) domask_chr, dt_data_output_av, &
1567	averaging_interval, dt_averaging_input
1568	ENDIF
1569
1570	IF ( av == 0 ) THEN
1571	IF ( skip_time_domask(mid) /= 0.0_wp ) THEN
1572	WRITE ( io, 339 ) skip_time_domask(mid)
1573	ENDIF
1574	ELSE
1575	IF ( skip_time_data_output_av /= 0.0_wp ) THEN
1576	WRITE ( io, 339 ) skip_time_data_output_av
1577	ENDIF
1578	ENDIF
1579	!
1580	!-- output locations
1581	DO dim = 1, 3
1582	IF ( mask(mid,dim,1) >= 0.0_wp ) THEN
1583	count = 0
1584	DO WHILE ( mask(mid,dim,count+1) >= 0.0_wp )
1585	count = count + 1
1586	ENDDO
1587	WRITE ( io, 349 ) dir(dim), dir(dim), mid, dir(dim), &
1588	mask(mid,dim,:count)
1589	ELSEIF ( mask_loop(mid,dim,1) < 0.0_wp .AND. &
1590	mask_loop(mid,dim,2) < 0.0_wp .AND. &
1591	mask_loop(mid,dim,3) == 0.0_wp ) THEN
1592	WRITE ( io, 350 ) dir(dim), dir(dim)
1593	ELSEIF ( mask_loop(mid,dim,3) == 0.0_wp ) THEN
1594	WRITE ( io, 351 ) dir(dim), dir(dim), mid, dir(dim), &
1595	mask_loop(mid,dim,1:2)
1596	ELSE
1597	WRITE ( io, 351 ) dir(dim), dir(dim), mid, dir(dim), &
1598	mask_loop(mid,dim,1:3)
1599	ENDIF
1600	ENDDO
1601	ENDIF
1602
1603	ENDDO
1604	ENDDO
1605
1606	!
1607	!-- Timeseries
1608	IF ( dt_dots /= 9999999.9_wp ) THEN
1609	WRITE ( io, 340 )
1610
1611	output_format = output_format_netcdf
1612	WRITE ( io, 344 ) output_format
1613	WRITE ( io, 341 ) dt_dots
1614	ENDIF
1615
1616	#if defined( __dvrp_graphics )
1617	!
1618	!-- Dvrp-output
1619	IF ( dt_dvrp /= 9999999.9_wp ) THEN
1620	WRITE ( io, 360 ) dt_dvrp, TRIM( dvrp_output ), TRIM( dvrp_host ), &
1621	TRIM( dvrp_username ), TRIM( dvrp_directory )
1622	i = 1
1623	l = 0
1624	m = 0
1625	DO WHILE ( mode_dvrp(i) /= ' ' )
1626	IF ( mode_dvrp(i)(1:10) == 'isosurface' ) THEN
1627	READ ( mode_dvrp(i), '(10X,I2)' ) j
1628	l = l + 1
1629	IF ( do3d(0,j) /= ' ' ) THEN
1630	WRITE ( io, 361 ) TRIM( do3d(0,j) ), threshold(l), &
1631	isosurface_color(:,l)
1632	ENDIF
1633	ELSEIF ( mode_dvrp(i)(1:6) == 'slicer' ) THEN
1634	READ ( mode_dvrp(i), '(6X,I2)' ) j
1635	m = m + 1
1636	IF ( do2d(0,j) /= ' ' ) THEN
1637	WRITE ( io, 362 ) TRIM( do2d(0,j) ), &
1638	slicer_range_limits_dvrp(:,m)
1639	ENDIF
1640	ELSEIF ( mode_dvrp(i)(1:9) == 'particles' ) THEN
1641	WRITE ( io, 363 ) dvrp_psize
1642	IF ( particle_dvrpsize /= 'none' ) THEN
1643	WRITE ( io, 364 ) 'size', TRIM( particle_dvrpsize ), &
1644	dvrpsize_interval
1645	ENDIF
1646	IF ( particle_color /= 'none' ) THEN
1647	WRITE ( io, 364 ) 'color', TRIM( particle_color ), &
1648	color_interval
1649	ENDIF
1650	ENDIF
1651	i = i + 1
1652	ENDDO
1653
1654	WRITE ( io, 365 ) groundplate_color, superelevation_x, &
1655	superelevation_y, superelevation, clip_dvrp_l, &
1656	clip_dvrp_r, clip_dvrp_s, clip_dvrp_n
1657
1658	IF ( TRIM( topography ) /= 'flat' ) THEN
1659	WRITE ( io, 366 ) topography_color
1660	IF ( cluster_size > 1 ) THEN
1661	WRITE ( io, 367 ) cluster_size
1662	ENDIF
1663	ENDIF
1664
1665	ENDIF
1666	#endif
1667
1668	#if defined( __spectra )
1669	!
1670	!-- Spectra output
1671	IF ( dt_dosp /= 9999999.9_wp ) THEN
1672	WRITE ( io, 370 )
1673
1674	output_format = output_format_netcdf
1675	WRITE ( io, 344 ) output_format
1676	WRITE ( io, 371 ) dt_dosp
1677	IF ( skip_time_dosp /= 0.0_wp ) WRITE ( io, 339 ) skip_time_dosp
1678	WRITE ( io, 372 ) ( data_output_sp(i), i = 1,10 ), &
1679	( spectra_direction(i), i = 1,10 ), &
1680	( comp_spectra_level(i), i = 1,100 ), &
1681	( plot_spectra_level(i), i = 1,100 ), &
1682	averaging_interval_sp, dt_averaging_input_pr
1683	ENDIF
1684	#endif
1685
1686	WRITE ( io, 99 )
1687
1688	!
1689	!-- Physical quantities
1690	WRITE ( io, 400 )
1691
1692	!
1693	!-- Geostrophic parameters
1694	IF ( radiation .AND. radiation_scheme /= 'constant' ) THEN
1695	WRITE ( io, 417 ) lambda
1696	ENDIF
1697	WRITE ( io, 410 ) phi, omega, f, fs
1698
1699	!
1700	!-- Other quantities
1701	WRITE ( io, 411 ) g
1702	IF ( radiation .AND. radiation_scheme /= 'constant' ) THEN
1703	WRITE ( io, 418 ) day_init, time_utc_init
1704	ENDIF
1705
1706	WRITE ( io, 412 ) TRIM( reference_state )
1707	IF ( use_single_reference_value ) THEN
1708	IF ( ocean ) THEN
1709	WRITE ( io, 413 ) prho_reference
1710	ELSE
1711	WRITE ( io, 414 ) pt_reference
1712	ENDIF
1713	ENDIF
1714
1715	!
1716	!-- Cloud physics parameters
1717	IF ( cloud_physics ) THEN
1718	WRITE ( io, 415 )
1719	WRITE ( io, 416 ) surface_pressure, r_d, rho_surface, cp, l_v
1720	IF ( icloud_scheme == 0 ) THEN
1721	WRITE ( io, 510 ) 1.0E-6_wp * nc_const
1722	IF ( precipitation ) WRITE ( io, 511 ) c_sedimentation
1723	ENDIF
1724	ENDIF
1725
1726	!
1727	!-- Cloud physcis parameters / quantities / numerical methods
1728	WRITE ( io, 430 )
1729	IF ( humidity .AND. .NOT. cloud_physics .AND. .NOT. cloud_droplets) THEN
1730	WRITE ( io, 431 )
1731	ELSEIF ( humidity .AND. cloud_physics ) THEN
1732	WRITE ( io, 432 )
1733	IF ( cloud_top_radiation ) WRITE ( io, 132 )
1734	IF ( icloud_scheme == 1 ) THEN
1735	IF ( precipitation ) WRITE ( io, 133 )
1736	ELSEIF ( icloud_scheme == 0 ) THEN
1737	IF ( drizzle ) WRITE ( io, 506 )
1738	IF ( precipitation ) THEN
1739	WRITE ( io, 505 )
1740	IF ( turbulence ) WRITE ( io, 507 )
1741	IF ( ventilation_effect ) WRITE ( io, 508 )
1742	IF ( limiter_sedimentation ) WRITE ( io, 509 )
1743	ENDIF
1744	ENDIF
1745	ELSEIF ( humidity .AND. cloud_droplets ) THEN
1746	WRITE ( io, 433 )
1747	IF ( curvature_solution_effects ) WRITE ( io, 434 )
1748	IF ( collision_kernel /= 'none' ) THEN
1749	WRITE ( io, 435 ) TRIM( collision_kernel )
1750	IF ( collision_kernel(6:9) == 'fast' ) THEN
1751	WRITE ( io, 436 ) radius_classes, dissipation_classes
1752	ENDIF
1753	ELSE
1754	WRITE ( io, 437 )
1755	ENDIF
1756	ENDIF
1757
1758	!
1759	!-- LES / turbulence parameters
1760	WRITE ( io, 450 )
1761
1762	!--
1763	! ... LES-constants used must still be added here
1764	!--
1765	IF ( constant_diffusion ) THEN
1766	WRITE ( io, 451 ) km_constant, km_constant/prandtl_number, &
1767	prandtl_number
1768	ENDIF
1769	IF ( .NOT. constant_diffusion) THEN
1770	IF ( e_init > 0.0_wp ) WRITE ( io, 455 ) e_init
1771	IF ( e_min > 0.0_wp ) WRITE ( io, 454 ) e_min
1772	IF ( wall_adjustment ) WRITE ( io, 453 ) wall_adjustment_factor
1773	ENDIF
1774
1775	!
1776	!-- Special actions during the run
1777	WRITE ( io, 470 )
1778	IF ( create_disturbances ) THEN
1779	WRITE ( io, 471 ) dt_disturb, disturbance_amplitude, &
1780	zu(disturbance_level_ind_b), disturbance_level_ind_b,&
1781	zu(disturbance_level_ind_t), disturbance_level_ind_t
1782	IF ( .NOT. bc_lr_cyc .OR. .NOT. bc_ns_cyc ) THEN
1783	WRITE ( io, 472 ) inflow_disturbance_begin, inflow_disturbance_end
1784	ELSE
1785	WRITE ( io, 473 ) disturbance_energy_limit
1786	ENDIF
1787	WRITE ( io, 474 ) TRIM( random_generator )
1788	ENDIF
1789	IF ( pt_surface_initial_change /= 0.0_wp ) THEN
1790	WRITE ( io, 475 ) pt_surface_initial_change
1791	ENDIF
1792	IF ( humidity .AND. q_surface_initial_change /= 0.0_wp ) THEN
1793	WRITE ( io, 476 ) q_surface_initial_change
1794	ENDIF
1795	IF ( passive_scalar .AND. q_surface_initial_change /= 0.0_wp ) THEN
1796	WRITE ( io, 477 ) q_surface_initial_change
1797	ENDIF
1798
1799	IF ( particle_advection ) THEN
1800	!
1801	!-- Particle attributes
1802	WRITE ( io, 480 ) particle_advection_start, dt_prel, bc_par_lr, &
1803	bc_par_ns, bc_par_b, bc_par_t, particle_maximum_age, &
1804	end_time_prel
1805	IF ( use_sgs_for_particles ) WRITE ( io, 488 ) dt_min_part
1806	IF ( random_start_position ) WRITE ( io, 481 )
1807	IF ( seed_follows_topography ) WRITE ( io, 496 )
1808	IF ( particles_per_point > 1 ) WRITE ( io, 489 ) particles_per_point
1809	WRITE ( io, 495 ) total_number_of_particles
1810	IF ( use_particle_tails .AND. maximum_number_of_tailpoints /= 0 ) THEN
1811	WRITE ( io, 483 ) maximum_number_of_tailpoints
1812	IF ( minimum_tailpoint_distance /= 0 ) THEN
1813	WRITE ( io, 484 ) total_number_of_tails, &
1814	minimum_tailpoint_distance, &
1815	maximum_tailpoint_age
1816	ENDIF
1817	ENDIF
1818	IF ( dt_write_particle_data /= 9999999.9_wp ) THEN
1819	WRITE ( io, 485 ) dt_write_particle_data
1820	IF ( netcdf_data_format > 1 ) THEN
1821	output_format = 'netcdf (64 bit offset) and binary'
1822	ELSE
1823	output_format = 'netcdf and binary'
1824	ENDIF
1825	WRITE ( io, 344 ) output_format
1826	ENDIF
1827	IF ( dt_dopts /= 9999999.9_wp ) WRITE ( io, 494 ) dt_dopts
1828	IF ( write_particle_statistics ) WRITE ( io, 486 )
1829
1830	WRITE ( io, 487 ) number_of_particle_groups
1831
1832	DO i = 1, number_of_particle_groups
1833	IF ( i == 1 .AND. density_ratio(i) == 9999999.9_wp ) THEN
1834	WRITE ( io, 490 ) i, 0.0_wp
1835	WRITE ( io, 492 )
1836	ELSE
1837	WRITE ( io, 490 ) i, radius(i)
1838	IF ( density_ratio(i) /= 0.0_wp ) THEN
1839	WRITE ( io, 491 ) density_ratio(i)
1840	ELSE
1841	WRITE ( io, 492 )
1842	ENDIF
1843	ENDIF
1844	WRITE ( io, 493 ) psl(i), psr(i), pss(i), psn(i), psb(i), pst(i), &
1845	pdx(i), pdy(i), pdz(i)
1846	IF ( .NOT. vertical_particle_advection(i) ) WRITE ( io, 482 )
1847	ENDDO
1848
1849	ENDIF
1850
1851
1852	!
1853	!-- Parameters of 1D-model
1854	IF ( INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 ) THEN
1855	WRITE ( io, 500 ) end_time_1d, dt_run_control_1d, dt_pr_1d, &
1856	mixing_length_1d, dissipation_1d
1857	IF ( damp_level_ind_1d /= nzt+1 ) THEN
1858	WRITE ( io, 502 ) zu(damp_level_ind_1d), damp_level_ind_1d
1859	ENDIF
1860	ENDIF
1861
1862	!
1863	!-- User-defined information
1864	CALL user_header( io )
1865
1866	WRITE ( io, 99 )
1867
1868	!
1869	!-- Write buffer contents to disc immediately
1870	CALL local_flush( io )
1871
1872	!
1873	!-- Here the FORMATs start
1874
1875	99 FORMAT (1X,78('-'))
1876	100 FORMAT (/1X,'******************************',4X,44('-')/ &
1877	1X,'* ',A,' *',4X,A/ &
1878	1X,'******************************',4X,44('-'))
1879	101 FORMAT (35X,'coupled run using MPI-',I1,': ',A/ &
1880	35X,42('-'))
1881	102 FORMAT (/' Date: ',A8,4X,'Run: ',A20/ &
1882	' Time: ',A8,4X,'Run-No.: ',I2.2/ &
1883	' Run on host: ',A10)
1884	#if defined( __parallel )
1885	103 FORMAT (' Number of PEs:',10X,I6,4X,'Processor grid (x,y): (',I4,',',I4, &
1886	')',1X,A)
1887	104 FORMAT (' Number of PEs:',10X,I6,4X,'Tasks:',I4,' threads per task:',I4/ &
1888	35X,'Processor grid (x,y): (',I4,',',I4,')',1X,A)
1889	105 FORMAT (35X,'One additional PE is used to handle'/37X,'the dvrp output!')
1890	106 FORMAT (35X,'A 1d-decomposition along x is forced'/ &
1891	35X,'because the job is running on an SMP-cluster')
1892	107 FORMAT (35X,'A 1d-decomposition along ',A,' is used')
1893	108 FORMAT (35X,'Max. # of parallel I/O streams is ',I5)
1894	109 FORMAT (35X,'Precursor run for coupled atmos-ocean run'/ &
1895	35X,42('-'))
1896	114 FORMAT (35X,'Coupled atmosphere-ocean run following'/ &
1897	35X,'independent precursor runs'/ &
1898	35X,42('-'))
1899	117 FORMAT (' Accelerator boards / node: ',I2)
1900	#endif
1901	110 FORMAT (/' Numerical Schemes:'/ &
1902	' -----------------'/)
1903	111 FORMAT (' --> Solve perturbation pressure via FFT using ',A,' routines')
1904	112 FORMAT (' --> Solve perturbation pressure via SOR-Red/Black-Schema'/ &
1905	' Iterations (initial/other): ',I3,'/',I3,' omega = ',F5.3)
1906	113 FORMAT (' --> Momentum advection via Piascek-Williams-Scheme (Form C3)', &
1907	' or Upstream')
1908	115 FORMAT (' FFT and transpositions are overlapping')
1909	116 FORMAT (' --> Scalar advection via Piascek-Williams-Scheme (Form C3)', &
1910	' or Upstream')
1911	118 FORMAT (' --> Scalar advection via Bott-Chlond-Scheme')
1912	119 FORMAT (' --> Galilei-Transform applied to horizontal advection:'/ &
1913	' translation velocity = ',A/ &
1914	' distance advected ',A,': ',F8.3,' km(x) ',F8.3,' km(y)')
1915	120 FORMAT (' Accelerator boards: ',8X,I2)
1916	122 FORMAT (' --> Time differencing scheme: ',A)
1917	123 FORMAT (' --> Rayleigh-Damping active, starts ',A,' z = ',F8.2,' m'/ &
1918	' maximum damping coefficient: ',F5.3, ' 1/s')
1919	129 FORMAT (' --> Additional prognostic equation for the specific humidity')
1920	130 FORMAT (' --> Additional prognostic equation for the total water content')
1921	131 FORMAT (' --> No pt-equation solved. Neutral stratification with pt = ', &
1922	F6.2, ' K assumed')
1923	132 FORMAT (' Parameterization of long-wave radiation processes via'/ &
1924	' effective emissivity scheme')
1925	133 FORMAT (' Precipitation parameterization via Kessler-Scheme')
1926	134 FORMAT (' --> Additional prognostic equation for a passive scalar')
1927	135 FORMAT (' --> Solve perturbation pressure via ',A,' method (', &
1928	A,'-cycle)'/ &
1929	' number of grid levels: ',I2/ &
1930	' Gauss-Seidel red/black iterations: ',I2)
1931	136 FORMAT (' gridpoints of coarsest subdomain (x,y,z): (',I3,',',I3,',', &
1932	I3,')')
1933	137 FORMAT (' level data gathered on PE0 at level: ',I2/ &
1934	' gridpoints of coarsest subdomain (x,y,z): (',I3,',',I3,',', &
1935	I3,')'/ &
1936	' gridpoints of coarsest domain (x,y,z): (',I3,',',I3,',', &
1937	I3,')')
1938	139 FORMAT (' --> Loop optimization method: ',A)
1939	140 FORMAT (' maximum residual allowed: ',E10.3)
1940	141 FORMAT (' fixed number of multigrid cycles: ',I4)
1941	142 FORMAT (' perturbation pressure is calculated at every Runge-Kutta ', &
1942	'step')
1943	143 FORMAT (' Euler/upstream scheme is used for the SGS turbulent ', &
1944	'kinetic energy')
1945	144 FORMAT (' masking method is used')
1946	150 FORMAT (' --> Volume flow at the right and north boundary will be ', &
1947	'conserved'/ &
1948	' using the ',A,' mode')
1949	151 FORMAT (' with u_bulk = ',F7.3,' m/s and v_bulk = ',F7.3,' m/s')
1950	152 FORMAT (' --> External pressure gradient directly prescribed by the user:',&
1951	/' ',2(1X,E12.5),'Pa/m in x/y direction', &
1952	/' starting from dp_level_b =', F8.3, 'm', A /)
1953	160 FORMAT (//' Large scale forcing and nudging:'/ &
1954	' -------------------------------'/)
1955	161 FORMAT (' --> No large scale forcing from external is used (default) ')
1956	162 FORMAT (' --> Large scale forcing from external file LSF_DATA is used: ')
1957	163 FORMAT (' - large scale advection tendencies ')
1958	164 FORMAT (' - large scale subsidence velocity w_subs ')
1959	165 FORMAT (' - large scale subsidence tendencies ')
1960	167 FORMAT (' - and geostrophic wind components ug and vg')
1961	168 FORMAT (' --> Large-scale vertical motion is used in the ', &
1962	'prognostic equation(s) for')
1963	169 FORMAT (' the scalar(s) only')
1964	170 FORMAT (' --> Nudging is used')
1965	171 FORMAT (' --> No nudging is used (default) ')
1966	180 FORMAT (' - prescribed surface values for temperature')
1967	181 FORMAT (' - prescribed surface fluxes for temperature')
1968	182 FORMAT (' - prescribed surface values for humidity')
1969	183 FORMAT (' - prescribed surface fluxes for humidity')
1970	200 FORMAT (//' Run time and time step information:'/ &
1971	' ----------------------------------'/)
1972	201 FORMAT ( ' Timestep: variable maximum value: ',F6.3,' s', &
1973	' CFL-factor: ',F4.2)
1974	202 FORMAT ( ' Timestep: dt = ',F6.3,' s'/)
1975	203 FORMAT ( ' Start time: ',F9.3,' s'/ &
1976	' End time: ',F9.3,' s')
1977	204 FORMAT ( A,F9.3,' s')
1978	205 FORMAT ( A,F9.3,' s',5X,'restart every',17X,F9.3,' s')
1979	206 FORMAT (/' Time reached: ',F9.3,' s'/ &
1980	' CPU-time used: ',F9.3,' s per timestep: ', &
1981	' ',F9.3,' s'/ &
1982	' per second of simulated tim', &
1983	'e: ',F9.3,' s')
1984	207 FORMAT ( ' Coupling start time: ',F9.3,' s')
1985	250 FORMAT (//' Computational grid and domain size:'/ &
1986	' ----------------------------------'// &
1987	' Grid length: dx = ',F7.3,' m dy = ',F7.3, &
1988	' m dz = ',F7.3,' m'/ &
1989	' Domain size: x = ',F10.3,' m y = ',F10.3, &
1990	' m z(u) = ',F10.3,' m'/)
1991	252 FORMAT (' dz constant up to ',F10.3,' m (k=',I4,'), then stretched by', &
1992	' factor: ',F5.3/ &
1993	' maximum dz not to be exceeded is dz_max = ',F10.3,' m'/)
1994	254 FORMAT (' Number of gridpoints (x,y,z): (0:',I4,', 0:',I4,', 0:',I4,')'/ &
1995	' Subdomain size (x,y,z): ( ',I4,', ',I4,', ',I4,')'/)
1996	260 FORMAT (/' The model has a slope in x-direction. Inclination angle: ',F6.2,&
1997	' degrees')
1998	270 FORMAT (//' Topography information:'/ &
1999	' ----------------------'// &
2000	1X,'Topography: ',A)
2001	271 FORMAT ( ' Building size (x/y/z) in m: ',F5.1,' / ',F5.1,' / ',F5.1/ &
2002	' Horizontal index bounds (l/r/s/n): ',I4,' / ',I4,' / ',I4, &
2003	' / ',I4)
2004	272 FORMAT ( ' Single quasi-2D street canyon of infinite length in ',A, &
2005	' direction' / &
2006	' Canyon height: ', F6.2, 'm, ch = ', I4, '.' / &
2007	' Canyon position (',A,'-walls): cxl = ', I4,', cxr = ', I4, '.')
2008	278 FORMAT (' Topography grid definition convention:'/ &
2009	' cell edge (staggered grid points'/ &
2010	' (u in x-direction, v in y-direction))' /)
2011	279 FORMAT (' Topography grid definition convention:'/ &
2012	' cell center (scalar grid points)' /)
2013	280 FORMAT (//' Vegetation canopy (drag) model:'/ &
2014	' ------------------------------'// &
2015	' Canopy mode: ', A / &
2016	' Canopy height: ',F6.2,'m (',I4,' grid points)' / &
2017	' Leaf drag coefficient: ',F6.2 /)
2018	281 FORMAT (/ ' Scalar exchange coefficient: ',F6.2 / &
2019	' Scalar concentration at leaf surfaces in kg/m**3: ',F6.2 /)
2020	282 FORMAT (' Predefined constant heatflux at the top of the vegetation: ',F6.2,' K m/s')
2021	283 FORMAT (/ ' Characteristic levels of the leaf area density:'// &
2022	' Height: ',A,' m'/ &
2023	' Leaf area density: ',A,' m2/m3'/ &
2024	' Gradient: ',A,' m2/m4'/ &
2025	' Gridpoint: ',A)
2026	284 FORMAT (//' Characteristic levels of the leaf area density and coefficients:'// &
2027	' Height: ',A,' m'/ &
2028	' Leaf area density: ',A,' m2/m3'/ &
2029	' Coefficient alpha: ',F6.2 / &
2030	' Coefficient beta: ',F6.2 / &
2031	' Leaf area index: ',F6.2,' m2/m2' /)
2032
2033	300 FORMAT (//' Boundary conditions:'/ &
2034	' -------------------'// &
2035	' p uv ', &
2036	' pt'// &
2037	' B. bound.: ',A/ &
2038	' T. bound.: ',A)
2039	301 FORMAT (/' ',A// &
2040	' B. bound.: ',A/ &
2041	' T. bound.: ',A)
2042	303 FORMAT (/' Bottom surface fluxes are used in diffusion terms at k=1')
2043	304 FORMAT (/' Top surface fluxes are used in diffusion terms at k=nzt')
2044	305 FORMAT (//' Prandtl-Layer between bottom surface and first ', &
2045	'computational u,v-level:'// &
2046	' zp = ',F6.2,' m z0 = ',F6.4,' m z0h = ',F7.5,&
2047	' m kappa = ',F4.2/ &
2048	' Rif value range: ',F6.2,' <= rif <=',F6.2)
2049	306 FORMAT (' Predefined constant heatflux: ',F9.6,' K m/s')
2050	307 FORMAT (' Heatflux has a random normal distribution')
2051	308 FORMAT (' Predefined surface temperature')
2052	309 FORMAT (' Predefined constant salinityflux: ',F9.6,' psu m/s')
2053	310 FORMAT (//' 1D-Model:'// &
2054	' Rif value range: ',F6.2,' <= rif <=',F6.2)
2055	311 FORMAT (' Predefined constant humidity flux: ',E10.3,' m/s')
2056	312 FORMAT (' Predefined surface humidity')
2057	313 FORMAT (' Predefined constant scalar flux: ',E10.3,' kg/(m**2 s)')
2058	314 FORMAT (' Predefined scalar value at the surface')
2059	315 FORMAT (' Humidity / scalar flux at top surface is 0.0')
2060	316 FORMAT (' Sensible heatflux and momentum flux from coupled ', &
2061	'atmosphere model')
2062	317 FORMAT (//' Lateral boundaries:'/ &
2063	' left/right: ',A/ &
2064	' north/south: ',A)
2065	318 FORMAT (/' use_cmax: ',L1 / &
2066	' pt damping layer width = ',F8.2,' m, pt ', &
2067	'damping factor = ',F6.4)
2068	319 FORMAT (' turbulence recycling at inflow switched on'/ &
2069	' width of recycling domain: ',F7.1,' m grid index: ',I4/ &
2070	' inflow damping height: ',F6.1,' m width: ',F6.1,' m')
2071	320 FORMAT (' Predefined constant momentumflux: u: ',F9.6,' m2/s2'/ &
2072	' v: ',F9.6,' m2/s2')
2073	321 FORMAT (//' Initial profiles:'/ &
2074	' ----------------')
2075	322 FORMAT (' turbulence recycling at inflow switched on'/ &
2076	' y shift of the recycled inflow turbulence switched on'/ &
2077	' width of recycling domain: ',F7.1,' m grid index: ',I4/ &
2078	' inflow damping height: ',F6.1,' m width: ',F6.1,' m'/q)
2079	325 FORMAT (//' List output:'/ &
2080	' -----------'// &
2081	' 1D-Profiles:'/ &
2082	' Output every ',F8.2,' s')
2083	326 FORMAT (' Time averaged over ',F8.2,' s'/ &
2084	' Averaging input every ',F8.2,' s')
2085	330 FORMAT (//' Data output:'/ &
2086	' -----------'/)
2087	331 FORMAT (/' 1D-Profiles:')
2088	332 FORMAT (/' ',A)
2089	333 FORMAT (' Output every ',F8.2,' s',/ &
2090	' Time averaged over ',F8.2,' s'/ &
2091	' Averaging input every ',F8.2,' s')
2092	334 FORMAT (/' 2D-Arrays',A,':')
2093	335 FORMAT (/' ',A2,'-cross-section Arrays: ',A/ &
2094	' Output every ',F8.2,' s ',A/ &
2095	' Cross sections at ',A1,' = ',A/ &
2096	' scalar-coordinates: ',A,' m'/)
2097	336 FORMAT (/' 3D-Arrays',A,':')
2098	337 FORMAT (/' Arrays: ',A/ &
2099	' Output every ',F8.2,' s ',A/ &
2100	' Upper output limit at ',F8.2,' m (GP ',I4,')'/)
2101	339 FORMAT (' No output during initial ',F8.2,' s')
2102	340 FORMAT (/' Time series:')
2103	341 FORMAT (' Output every ',F8.2,' s'/)
2104	342 FORMAT (/' ',A2,'-cross-section Arrays: ',A/ &
2105	' Output every ',F8.2,' s ',A/ &
2106	' Time averaged over ',F8.2,' s'/ &
2107	' Averaging input every ',F8.2,' s'/ &
2108	' Cross sections at ',A1,' = ',A/ &
2109	' scalar-coordinates: ',A,' m'/)
2110	343 FORMAT (/' Arrays: ',A/ &
2111	' Output every ',F8.2,' s ',A/ &
2112	' Time averaged over ',F8.2,' s'/ &
2113	' Averaging input every ',F8.2,' s'/ &
2114	' Upper output limit at ',F8.2,' m (GP ',I4,')'/)
2115	344 FORMAT (' Output format: ',A/)
2116	345 FORMAT (/' Scaling lengths for output locations of all subsequent mask IDs:',/ &
2117	' mask_scale_x (in x-direction): ',F9.3, ' m',/ &
2118	' mask_scale_y (in y-direction): ',F9.3, ' m',/ &
2119	' mask_scale_z (in z-direction): ',F9.3, ' m' )
2120	346 FORMAT (/' Masked data output',A,' for mask ID ',I2, ':')
2121	347 FORMAT (' Variables: ',A/ &
2122	' Output every ',F8.2,' s')
2123	348 FORMAT (' Variables: ',A/ &
2124	' Output every ',F8.2,' s'/ &
2125	' Time averaged over ',F8.2,' s'/ &
2126	' Averaging input every ',F8.2,' s')
2127	349 FORMAT (/' Output locations in ',A,'-direction in multiples of ', &
2128	'mask_scale_',A,' predefined by array mask_',I2.2,'_',A,':'/ &
2129	13(' ',8(F8.2,',')/) )
2130	350 FORMAT (/' Output locations in ',A,'-direction: ', &
2131	'all gridpoints along ',A,'-direction (default).' )
2132	351 FORMAT (/' Output locations in ',A,'-direction in multiples of ', &
2133	'mask_scale_',A,' constructed from array mask_',I2.2,'_',A,'_loop:'/ &
2134	' loop begin:',F8.2,', end:',F8.2,', stride:',F8.2 )
2135	352 FORMAT (/' Number of output time levels allowed: ',I3 /)
2136	353 FORMAT (/' Number of output time levels allowed: unlimited' /)
2137	#if defined( __dvrp_graphics )
2138	360 FORMAT (' Plot-Sequence with dvrp-software:'/ &
2139	' Output every ',F7.1,' s'/ &
2140	' Output mode: ',A/ &
2141	' Host / User: ',A,' / ',A/ &
2142	' Directory: ',A// &
2143	' The sequence contains:')
2144	361 FORMAT (/' Isosurface of "',A,'" Threshold value: ', E12.3/ &
2145	' Isosurface color: (',F4.2,',',F4.2,',',F4.2,') (R,G,B)')
2146	362 FORMAT (/' Slicer plane ',A/ &
2147	' Slicer limits: [',F6.2,',',F6.2,']')
2148	363 FORMAT (/' Particles'/ &
2149	' particle size: ',F7.2,' m')
2150	364 FORMAT (' particle ',A,' controlled by "',A,'" with interval [', &
2151	F6.2,',',F6.2,']')
2152	365 FORMAT (/' Groundplate color: (',F4.2,',',F4.2,',',F4.2,') (R,G,B)'/ &
2153	' Superelevation along (x,y,z): (',F4.1,',',F4.1,',',F4.1, &
2154	')'/ &
2155	' Clipping limits: from x = ',F9.1,' m to x = ',F9.1,' m'/ &
2156	' from y = ',F9.1,' m to y = ',F9.1,' m')
2157	366 FORMAT (/' Topography color: (',F4.2,',',F4.2,',',F4.2,') (R,G,B)')
2158	367 FORMAT (' Polygon reduction for topography: cluster_size = ', I1)
2159	#endif
2160	#if defined( __spectra )
2161	370 FORMAT (' Spectra:')
2162	371 FORMAT (' Output every ',F7.1,' s'/)
2163	372 FORMAT (' Arrays: ', 10(A5,',')/ &
2164	' Directions: ', 10(A5,',')/ &
2165	' height levels k = ', 20(I3,',')/ &
2166	' ', 20(I3,',')/ &
2167	' ', 20(I3,',')/ &
2168	' ', 20(I3,',')/ &
2169	' ', 19(I3,','),I3,'.'/ &
2170	' height levels selected for standard plot:'/ &
2171	' k = ', 20(I3,',')/ &
2172	' ', 20(I3,',')/ &
2173	' ', 20(I3,',')/ &
2174	' ', 20(I3,',')/ &
2175	' ', 19(I3,','),I3,'.'/ &
2176	' Time averaged over ', F7.1, ' s,' / &
2177	' Profiles for the time averaging are taken every ', &
2178	F6.1,' s')
2179	#endif
2180	400 FORMAT (//' Physical quantities:'/ &
2181	' -------------------'/)
2182	410 FORMAT (' Geograph. latitude : phi = ',F4.1,' degr'/ &
2183	' Angular velocity : omega = ',E9.3,' rad/s'/ &
2184	' Coriolis parameter : f = ',F9.6,' 1/s'/ &
2185	' f* = ',F9.6,' 1/s')
2186	411 FORMAT (/' Gravity : g = ',F4.1,' m/s**2')
2187	412 FORMAT (/' Reference state used in buoyancy terms: ',A)
2188	413 FORMAT (' Reference density in buoyancy terms: ',F8.3,' kg/m**3')
2189	414 FORMAT (' Reference temperature in buoyancy terms: ',F8.4,' K')
2190	415 FORMAT (/' Cloud physics parameters:'/ &
2191	' ------------------------'/)
2192	416 FORMAT (' Surface pressure : p_0 = ',F7.2,' hPa'/ &
2193	' Gas constant : R = ',F5.1,' J/(kg K)'/ &
2194	' Density of air : rho_0 = ',F5.3,' kg/m**3'/ &
2195	' Specific heat cap. : c_p = ',F6.1,' J/(kg K)'/ &
2196	' Vapourization heat : L_v = ',E8.2,' J/kg')
2197	417 FORMAT (' Geograph. longitude : lambda = ',F4.1,' degr')
2198	418 FORMAT (/' Day of the year at model start : day_init = ',I3 &
2199	/' UTC time at model start : time_utc_init = ',F7.1' s')
2200	419 FORMAT (//' Land surface model information:'/ &
2201	' ------------------------------'/)
2202	420 FORMAT (/' Characteristic levels of the initial temperature profile:'// &
2203	' Height: ',A,' m'/ &
2204	' Temperature: ',A,' K'/ &
2205	' Gradient: ',A,' K/100m'/ &
2206	' Gridpoint: ',A)
2207	421 FORMAT (/' Characteristic levels of the initial humidity profile:'// &
2208	' Height: ',A,' m'/ &
2209	' Humidity: ',A,' kg/kg'/ &
2210	' Gradient: ',A,' (kg/kg)/100m'/ &
2211	' Gridpoint: ',A)
2212	422 FORMAT (/' Characteristic levels of the initial scalar profile:'// &
2213	' Height: ',A,' m'/ &
2214	' Scalar concentration: ',A,' kg/m**3'/ &
2215	' Gradient: ',A,' (kg/m**3)/100m'/ &
2216	' Gridpoint: ',A)
2217	423 FORMAT (/' Characteristic levels of the geo. wind component ug:'// &
2218	' Height: ',A,' m'/ &
2219	' ug: ',A,' m/s'/ &
2220	' Gradient: ',A,' 1/100s'/ &
2221	' Gridpoint: ',A)
2222	424 FORMAT (/' Characteristic levels of the geo. wind component vg:'// &
2223	' Height: ',A,' m'/ &
2224	' vg: ',A,' m/s'/ &
2225	' Gradient: ',A,' 1/100s'/ &
2226	' Gridpoint: ',A)
2227	425 FORMAT (/' Characteristic levels of the initial salinity profile:'// &
2228	' Height: ',A,' m'/ &
2229	' Salinity: ',A,' psu'/ &
2230	' Gradient: ',A,' psu/100m'/ &
2231	' Gridpoint: ',A)
2232	426 FORMAT (/' Characteristic levels of the subsidence/ascent profile:'// &
2233	' Height: ',A,' m'/ &
2234	' w_subs: ',A,' m/s'/ &
2235	' Gradient: ',A,' (m/s)/100m'/ &
2236	' Gridpoint: ',A)
2237	427 FORMAT (/' Initial wind profiles (u,v) are interpolated from given'// &
2238	' profiles')
2239	428 FORMAT (/' Initial profiles (u, v, pt, q) are taken from file '/ &
2240	' NUDGING_DATA')
2241	430 FORMAT (//' Cloud physics quantities / methods:'/ &
2242	' ----------------------------------'/)
2243	431 FORMAT (' Humidity is treated as purely passive scalar (no condensati', &
2244	'on)')
2245	432 FORMAT (' Bulk scheme with liquid water potential temperature and'/ &
2246	' total water content is used.'/ &
2247	' Condensation is parameterized via 0% - or 100% scheme.')
2248	433 FORMAT (' Cloud droplets treated explicitly using the Lagrangian part', &
2249	'icle model')
2250	434 FORMAT (' Curvature and solution effecs are considered for growth of', &
2251	' droplets < 1.0E-6 m')
2252	435 FORMAT (' Droplet collision is handled by ',A,'-kernel')
2253	436 FORMAT (' Fast kernel with fixed radius- and dissipation classes ', &
2254	'are used'/ &
2255	' number of radius classes: ',I3,' interval ', &
2256	'[1.0E-6,2.0E-4] m'/ &
2257	' number of dissipation classes: ',I2,' interval ', &
2258	'[0,1000] cm2/s3')
2259	437 FORMAT (' Droplet collision is switched off')
2260	438 FORMAT (' --> Land surface type : ',A,/ &
2261	' --> Soil porosity type : ',A)
2262	439 FORMAT (/' Initial soil temperature and moisture profile:'// &
2263	' Height: ',A,' m'/ &
2264	' Temperature: ',A,' K'/ &
2265	' Moisture: ',A,' m3/m3'/ &
2266	' Root fraction: ',A,' '/ &
2267	' Gridpoint: ',A)
2268	440 FORMAT (/' --> Dewfall is allowed (default)')
2269	441 FORMAT (' --> Dewfall is inhibited')
2270	442 FORMAT (' --> Soil bottom is closed (water content is conserved, default)')
2271	443 FORMAT (' --> Soil bottom is open (water content is not conserved)')
2272	444 FORMAT (//' Radiation model information:'/ &
2273	' ----------------------------'/)
2274	445 FORMAT (' --> Using constant net radiation: net_radiation = ', F6.2, ' W/m**2')
2275	446 FORMAT (' --> Simple radiation scheme for clear sky is used (no clouds,', &
2276	' default)')
2277	447 FORMAT (' --> RRTMG scheme is used')
2278	448 FORMAT (/' User-specific surface albedo: albedo = ', F5.3)
2279	449 FORMAT (' Timestep: dt_radiation = ', F5.2, ' s')
2280
2281	450 FORMAT (//' LES / Turbulence quantities:'/ &
2282	' ---------------------------'/)
2283	451 FORMAT (' Diffusion coefficients are constant:'/ &
2284	' Km = ',F6.2,' m2/s Kh = ',F6.2,' m2/s Pr = ',F5.2)
2285	453 FORMAT (' Mixing length is limited to ',F4.2,' * z')
2286	454 FORMAT (' TKE is not allowed to fall below ',E9.2,' (m/s)**2')
2287	455 FORMAT (' initial TKE is prescribed as ',E9.2,' (m/s)**2')
2288	456 FORMAT (/' Albedo is set for land surface type: ', A)
2289	457 FORMAT (/' --> Albedo is fixed during the run')
2290	458 FORMAT (/' --> Longwave radiation is disabled')
2291	459 FORMAT (/' --> Shortwave radiation is disabled.')
2292	470 FORMAT (//' Actions during the simulation:'/ &
2293	' -----------------------------'/)
2294	471 FORMAT (' Disturbance impulse (u,v) every : ',F6.2,' s'/ &
2295	' Disturbance amplitude : ',F4.2, ' m/s'/ &
2296	' Lower disturbance level : ',F8.2,' m (GP ',I4,')'/ &
2297	' Upper disturbance level : ',F8.2,' m (GP ',I4,')')
2298	472 FORMAT (' Disturbances continued during the run from i/j =',I4, &
2299	' to i/j =',I4)
2300	473 FORMAT (' Disturbances cease as soon as the disturbance energy exceeds',&
2301	1X,F5.3, ' m2/s2')
2302	474 FORMAT (' Random number generator used : ',A/)
2303	475 FORMAT (' The surface temperature is increased (or decreased, ', &
2304	'respectively, if'/ &
2305	' the value is negative) by ',F5.2,' K at the beginning of the',&
2306	' 3D-simulation'/)
2307	476 FORMAT (' The surface humidity is increased (or decreased, ',&
2308	'respectively, if the'/ &
2309	' value is negative) by ',E8.1,' kg/kg at the beginning of', &
2310	' the 3D-simulation'/)
2311	477 FORMAT (' The scalar value is increased at the surface (or decreased, ',&
2312	'respectively, if the'/ &
2313	' value is negative) by ',E8.1,' kg/m**3 at the beginning of', &
2314	' the 3D-simulation'/)
2315	480 FORMAT (' Particles:'/ &
2316	' ---------'// &
2317	' Particle advection is active (switched on at t = ', F7.1, &
2318	' s)'/ &
2319	' Start of new particle generations every ',F6.1,' s'/ &
2320	' Boundary conditions: left/right: ', A, ' north/south: ', A/&
2321	' bottom: ', A, ' top: ', A/&
2322	' Maximum particle age: ',F9.1,' s'/ &
2323	' Advection stopped at t = ',F9.1,' s'/)
2324	481 FORMAT (' Particles have random start positions'/)
2325	482 FORMAT (' Particles are advected only horizontally'/)
2326	483 FORMAT (' Particles have tails with a maximum of ',I3,' points')
2327	484 FORMAT (' Number of tails of the total domain: ',I10/ &
2328	' Minimum distance between tailpoints: ',F8.2,' m'/ &
2329	' Maximum age of the end of the tail: ',F8.2,' s')
2330	485 FORMAT (' Particle data are written on file every ', F9.1, ' s')
2331	486 FORMAT (' Particle statistics are written on file'/)
2332	487 FORMAT (' Number of particle groups: ',I2/)
2333	488 FORMAT (' SGS velocity components are used for particle advection'/ &
2334	' minimum timestep for advection: ', F7.5/)
2335	489 FORMAT (' Number of particles simultaneously released at each ', &
2336	'point: ', I5/)
2337	490 FORMAT (' Particle group ',I2,':'/ &
2338	' Particle radius: ',E10.3, 'm')
2339	491 FORMAT (' Particle inertia is activated'/ &
2340	' density_ratio (rho_fluid/rho_particle) = ',F5.3/)
2341	492 FORMAT (' Particles are advected only passively (no inertia)'/)
2342	493 FORMAT (' Boundaries of particle source: x:',F8.1,' - ',F8.1,' m'/&
2343	' y:',F8.1,' - ',F8.1,' m'/&
2344	' z:',F8.1,' - ',F8.1,' m'/&
2345	' Particle distances: dx = ',F8.1,' m dy = ',F8.1, &
2346	' m dz = ',F8.1,' m'/)
2347	494 FORMAT (' Output of particle time series in NetCDF format every ', &
2348	F8.2,' s'/)
2349	495 FORMAT (' Number of particles in total domain: ',I10/)
2350	496 FORMAT (' Initial vertical particle positions are interpreted ', &
2351	'as relative to the given topography')
2352	500 FORMAT (//' 1D-Model parameters:'/ &
2353	' -------------------'// &
2354	' Simulation time: ',F8.1,' s'/ &
2355	' Run-controll output every: ',F8.1,' s'/ &
2356	' Vertical profile output every: ',F8.1,' s'/ &
2357	' Mixing length calculation: ',A/ &
2358	' Dissipation calculation: ',A/)
2359	502 FORMAT (' Damping layer starts from ',F7.1,' m (GP ',I4,')'/)
2360	503 FORMAT (' --> Momentum advection via Wicker-Skamarock-Scheme 5th order')
2361	504 FORMAT (' --> Scalar advection via Wicker-Skamarock-Scheme 5th order')
2362	505 FORMAT (' Precipitation parameterization via Seifert-Beheng-Scheme')
2363	506 FORMAT (' Drizzle parameterization via Stokes law')
2364	507 FORMAT (' Turbulence effects on precipitation process')
2365	508 FORMAT (' Ventilation effects on evaporation of rain drops')
2366	509 FORMAT (' Slope limiter used for sedimentation process')
2367	510 FORMAT (' Droplet density : N_c = ',F6.1,' 1/cm**3')
2368	511 FORMAT (' Sedimentation Courant number: '/&
2369	' C_s = ',F3.1,' ')
2370	512 FORMAT (/' Date: ',A8,6X,'Run: ',A20/ &
2371	' Time: ',A8,6X,'Run-No.: ',I2.2/ &
2372	' Run on host: ',A10,6X,'En-No.: ',I2.2)
2373	513 FORMAT (' --> Scalar advection via Wicker-Skamarock-Scheme 5th order ' // &
2374	'+ monotonic adjustment')
2375
2376
2377	END SUBROUTINE header

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |