Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |