Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |