Home

Context Navigation

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

Last change on this file since 2118 was 2118, checked in by raasch, 7 years ago
all OpenACC directives and related parts removed from the code
Property svn:keywords set to `Id`
File size: 88.5 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |