Home

Context Navigation

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

Last change on this file since 4331 was 4309, checked in by suehring, 5 years ago
Synthetic turbulence generator: Computation of velocity seeds optimized. This implies that random numbers are computed now using the parallel random number generator. Random number are now only computed and normalized locally, while distributed over all mpi ranks afterwards, instead of computing random numbers on a global array. urther, the number of calls for the time-consuming velocity-seed generation is reduced - now the left and right, as well as the north and south boundary share the same velocity-seed matrices.
Property svn:keywords set to `Id`
File size: 76.6 KB

Line
1	! !> @file header.f90
2	!------------------------------------------------------------------------------!
3	! This file is part of the PALM model system.
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-2019 Leibniz Universitaet Hannover
18	!------------------------------------------------------------------------------!
19	!
20	! Current revisions:
21	! -----------------
22	!
23	!
24	! Former revisions:
25	! -----------------
26	! $Id: header.f90 4309 2019-11-26 18:49:59Z suehring $
27	! replaced recycling_yshift by y_shift
28	!
29	! 4301 2019-11-22 12:09:09Z oliver.maas
30	!
31	! 4297 2019-11-21 10:37:50Z oliver.maas
32	! Adjusted format for simulated time and related quantities
33	!
34	! 4297 2019-11-21 10:37:50Z oliver.maas
35	! adjusted message to the changed parameter recycling_yshift
36	!
37	! 4227 2019-09-10 18:04:34Z gronemeier
38	! implement new palm_date_time_mod
39	!
40	! 4223 2019-09-10 09:20:47Z gronemeier
41	! Write information about rotation angle
42	!
43	! 4182 2019-08-22 15:20:23Z scharf
44	! Corrected "Former revisions" section
45	!
46	! 4168 2019-08-16 13:50:17Z suehring
47	! Replace function get_topography_top_index by topo_top_ind
48	!
49	! 4069 2019-07-01 14:05:51Z Giersch
50	! Masked output running index mid has been introduced as a local variable to
51	! avoid runtime error (Loop variable has been modified) in time_integration
52	!
53	! 4023 2019-06-12 13:20:01Z maronga
54	! Renamed "coupling start time" to "spinup time"
55	!
56	! 4017 2019-06-06 12:16:46Z schwenkel
57	! unused variable removed
58	!
59	! 3655 2019-01-07 16:51:22Z knoop
60	! Implementation of the PALM module interface
61	!
62	! Revision 1.1 1997/08/11 06:17:20 raasch
63	! Initial revision
64	!
65	!
66	! Description:
67	! ------------
68	!> Writing a header with all important information about the current run.
69	!> This subroutine is called three times, two times at the beginning
70	!> (writing information on files RUN_CONTROL and HEADER) and one time at the
71	!> end of the run, then writing additional information about CPU-usage on file
72	!> header.
73	!-----------------------------------------------------------------------------!
74	SUBROUTINE header
75
76
77	USE arrays_3d, &
78	ONLY: pt_init, q_init, s_init, sa_init, ug, vg, w_subs, zu, zw
79
80	USE basic_constants_and_equations_mod, &
81	ONLY: g, kappa
82
83	USE bulk_cloud_model_mod, &
84	ONLY: bulk_cloud_model
85
86	USE control_parameters
87
88	USE cpulog, &
89	ONLY: log_point_s
90
91	USE grid_variables, &
92	ONLY: dx, dy
93
94	USE indices, &
95	ONLY: mg_loc_ind, nnx, nny, nnz, nx, ny, nxl_mg, nxr_mg, nyn_mg, &
96	nys_mg, nzt, nzt_mg, topo_top_ind
97
98	USE kinds
99
100	USE model_1d_mod, &
101	ONLY: damp_level_ind_1d, dt_pr_1d, dt_run_control_1d, end_time_1d
102
103	USE module_interface, &
104	ONLY: module_interface_header
105
106	USE netcdf_interface, &
107	ONLY: netcdf_data_format, netcdf_data_format_string, netcdf_deflate
108
109	USE ocean_mod, &
110	ONLY: ibc_sa_t, prho_reference, sa_surface, &
111	sa_vertical_gradient, sa_vertical_gradient_level, &
112	sa_vertical_gradient_level_ind
113
114	USE palm_date_time_mod, &
115	ONLY: get_date_time
116
117	USE pegrid
118
119	#if defined( __parallel )
120	USE pmc_handle_communicator, &
121	ONLY: pmc_get_model_info
122	#endif
123
124	USE pmc_interface, &
125	ONLY: nested_run, nesting_datatransfer_mode, nesting_mode
126
127	USE surface_mod, &
128	ONLY: surf_def_h
129
130	USE turbulence_closure_mod, &
131	ONLY: rans_const_c, rans_const_sigma
132
133	IMPLICIT NONE
134
135
136	CHARACTER (LEN=2) :: do2d_mode !< mode of 2D data output (xy, xz, yz)
137
138	CHARACTER (LEN=5) :: section_chr !< string indicating grid information where to output 2D slices
139
140	CHARACTER (LEN=10) :: coor_chr !< string for subsidence velocities in large-scale forcing
141	CHARACTER (LEN=10) :: host_chr !< string for hostname
142
143	CHARACTER (LEN=16) :: begin_chr !< string indication start time for the data output
144
145	CHARACTER (LEN=26) :: ver_rev !< string for run identification
146
147	CHARACTER (LEN=32) :: cpl_name !< name of child domain (nesting mode only)
148
149	CHARACTER (LEN=40) :: output_format !< netcdf format
150
151	CHARACTER (LEN=70) :: char1 !< dummy varialbe used for various strings
152	CHARACTER (LEN=70) :: char2 !< string containing informating about the advected distance in case of Galilei transformation
153	CHARACTER (LEN=23) :: date_time_str !< string for date and time information
154	CHARACTER (LEN=70) :: dopr_chr !< string indicating profile output variables
155	CHARACTER (LEN=70) :: do2d_xy !< string indicating 2D-xy output variables
156	CHARACTER (LEN=70) :: do2d_xz !< string indicating 2D-xz output variables
157	CHARACTER (LEN=70) :: do2d_yz !< string indicating 2D-yz output variables
158	CHARACTER (LEN=70) :: do3d_chr !< string indicating 3D output variables
159	CHARACTER (LEN=70) :: domask_chr !< string indicating masked output variables
160	CHARACTER (LEN=70) :: run_classification !< string classifying type of run, e.g. nested, coupled, etc.
161
162	CHARACTER (LEN=85) :: r_upper !< string indicating model top boundary condition for various quantities
163	CHARACTER (LEN=85) :: r_lower !< string indicating bottom boundary condition for various quantities
164
165	CHARACTER (LEN=86) :: coordinates !< string indicating height coordinates for profile-prescribed variables
166	CHARACTER (LEN=86) :: gradients !< string indicating gradients of profile-prescribed variables between the prescribed height coordinates
167	CHARACTER (LEN=86) :: slices !< string indicating grid coordinates of profile-prescribed subsidence velocity
168	CHARACTER (LEN=86) :: temperatures !< string indicating profile-prescribed subsidence velocities
169	CHARACTER (LEN=86) :: ugcomponent !< string indicating profile-prescribed geostrophic u-component
170	CHARACTER (LEN=86) :: vgcomponent !< string indicating profile-prescribed geostrophic v-component
171
172	CHARACTER (LEN=1), DIMENSION(1:3) :: dir = (/ 'x', 'y', 'z' /) !< string indicating masking steps along certain direction
173
174	INTEGER(iwp) :: av !< index indicating average output quantities
175	INTEGER(iwp) :: bh !< building height in generic single-building setup
176	INTEGER(iwp) :: blx !< building width in grid points along x in generic single-building setup
177	INTEGER(iwp) :: bly !< building width in grid points along y in generic single-building setup
178	INTEGER(iwp) :: bxl !< index for left building wall in generic single-building setup
179	INTEGER(iwp) :: bxr !< index for right building wall in generic single-building setup
180	INTEGER(iwp) :: byn !< index for north building wall in generic single-building setup
181	INTEGER(iwp) :: bys !< index for south building wall in generic single-building setup
182	INTEGER(iwp) :: ch !< canyon depth in generic street-canyon setup
183	INTEGER(iwp) :: count !< number of masked output locations
184	INTEGER(iwp) :: cpl_parent_id !< parent ID for the respective child model
185	INTEGER(iwp) :: cwx !< canyon width along x in generic street-canyon setup
186	INTEGER(iwp) :: cwy !< canyon width along y in generic street-canyon setup
187	INTEGER(iwp) :: cxl !< index for left canyon wall in generic street-canyon setup
188	INTEGER(iwp) :: cxr !< index for right canyon wall in generic street-canyon setup
189	INTEGER(iwp) :: cyn !< index for north canyon wall in generic street-canyon setup
190	INTEGER(iwp) :: cys !< index for south canyon wall in generic street-canyon setup
191	INTEGER(iwp) :: dim !< running index for masking output locations
192	INTEGER(iwp) :: i !< running index for various loops
193	INTEGER(iwp) :: io !< file unit of HEADER file
194	INTEGER(iwp) :: l !< substring length
195	INTEGER(iwp) :: ll !< substring length
196	INTEGER(iwp) :: mid !< masked output running index
197	INTEGER(iwp) :: my_cpl_id !< run id in a nested model setup
198	INTEGER(iwp) :: n !< running index over number of couplers in a nested model setup
199	INTEGER(iwp) :: ncpl !< number of coupler in a nested model setup
200	INTEGER(iwp) :: npe_total !< number of total PEs in a coupler (parent + child)
201
202
203	REAL(wp) :: cpuseconds_per_simulated_second !< CPU time (in s) per simulated second
204	REAL(wp) :: lower_left_coord_x !< x-coordinate of nest domain
205	REAL(wp) :: lower_left_coord_y !< y-coordinate of nest domain
206
207	!
208	!-- Open the output file. At the end of the simulation, output is directed
209	!-- to unit 19.
210	IF ( ( runnr == 0 .OR. force_print_header ) .AND. &
211	.NOT. simulated_time_at_begin /= simulated_time ) THEN
212	io = 15 ! header output on file RUN_CONTROL
213	ELSE
214	io = 19 ! header output on file HEADER
215	ENDIF
216	CALL check_open( io )
217
218	!
219	!-- At the end of the run, output file (HEADER) will be rewritten with
220	!-- new information
221	IF ( io == 19 .AND. simulated_time_at_begin /= simulated_time ) REWIND( 19 )
222
223	!
224	!-- Determine kind of model run
225	IF ( TRIM( initializing_actions ) == 'read_restart_data' ) THEN
226	run_classification = 'restart run'
227	ELSEIF ( TRIM( initializing_actions ) == 'cyclic_fill' ) THEN
228	run_classification = 'run with cyclic fill of 3D - prerun data'
229	ELSEIF ( INDEX( initializing_actions, 'set_constant_profiles' ) /= 0 ) THEN
230	run_classification = 'run without 1D - prerun'
231	ELSEIF ( INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 ) THEN
232	run_classification = 'run with 1D - prerun'
233	ELSEIF ( INDEX( initializing_actions, 'inifor' ) /= 0 ) THEN
234	run_classification = 'run initialized with COSMO data'
235	ELSEIF ( INDEX( initializing_actions, 'by_user' ) /=0 ) THEN
236	run_classification = 'run initialized by user'
237	ELSEIF ( INDEX( initializing_actions, 'initialize_vortex' ) /=0 ) THEN
238	run_classification = 'run additionally initialized by a Rankine-vortex'
239	ELSEIF ( INDEX( initializing_actions, 'initialize_ptanom' ) /=0 ) THEN
240	run_classification = 'run additionally initialized by temperature anomaly'
241	ELSE
242	message_string = ' unknown action(s): ' // TRIM( initializing_actions )
243	CALL message( 'header', 'PA0191', 0, 0, 0, 6, 0 )
244	ENDIF
245	IF ( nested_run ) run_classification = 'nested ' // run_classification(1:63)
246	IF ( ocean_mode ) THEN
247	run_classification = 'ocean - ' // run_classification(1:61)
248	ELSE
249	run_classification = 'atmosphere - ' // run_classification(1:57)
250	ENDIF
251
252	!
253	!-- Run-identification, date, time, host
254	host_chr = host(1:10)
255	ver_rev = TRIM( version ) // ' ' // TRIM( revision )
256	WRITE ( io, 100 ) ver_rev, TRIM( run_classification )
257	IF ( TRIM( coupling_mode ) /= 'uncoupled' ) THEN
258	WRITE ( io, 101 ) coupling_mode
259	ENDIF
260	#if defined( __parallel )
261	IF ( coupling_start_time /= 0.0_wp .AND. .NOT. spinup ) THEN
262	IF ( coupling_start_time > simulated_time_at_begin ) THEN
263	WRITE ( io, 109 )
264	ELSE
265	WRITE ( io, 114 )
266	ENDIF
267	ENDIF
268	#endif
269	IF ( ensemble_member_nr /= 0 ) THEN
270	WRITE ( io, 512 ) run_date, run_identifier, run_time, runnr, &
271	ADJUSTR( host_chr ), ensemble_member_nr
272	ELSE
273	WRITE ( io, 102 ) run_date, run_identifier, run_time, runnr, &
274	ADJUSTR( host_chr )
275	ENDIF
276	#if defined( __parallel )
277	IF ( npex == -1 .AND. npey == -1 ) THEN
278	char1 = 'calculated'
279	ELSE
280	char1 = 'predefined'
281	ENDIF
282	IF ( threads_per_task == 1 ) THEN
283	WRITE ( io, 103 ) numprocs, pdims(1), pdims(2), TRIM( char1 )
284	ELSE
285	WRITE ( io, 104 ) numprocs*threads_per_task, numprocs, &
286	threads_per_task, pdims(1), pdims(2), TRIM( char1 )
287	ENDIF
288
289	IF ( pdims(2) == 1 ) THEN
290	WRITE ( io, 107 ) 'x'
291	ELSEIF ( pdims(1) == 1 ) THEN
292	WRITE ( io, 107 ) 'y'
293	ENDIF
294	IF ( numprocs /= maximum_parallel_io_streams ) THEN
295	WRITE ( io, 108 ) maximum_parallel_io_streams
296	ENDIF
297	#endif
298
299	!
300	!-- Nesting informations
301	IF ( nested_run ) THEN
302
303	#if defined( __parallel )
304	WRITE ( io, 600 ) TRIM( nesting_mode ), &
305	TRIM( nesting_datatransfer_mode )
306	CALL pmc_get_model_info( ncpl = ncpl, cpl_id = my_cpl_id )
307
308	DO n = 1, ncpl
309	CALL pmc_get_model_info( request_for_cpl_id = n, cpl_name = cpl_name,&
310	cpl_parent_id = cpl_parent_id, &
311	lower_left_x = lower_left_coord_x, &
312	lower_left_y = lower_left_coord_y, &
313	npe_total = npe_total )
314	IF ( n == my_cpl_id ) THEN
315	char1 = '*'
316	ELSE
317	char1 = ' '
318	ENDIF
319	WRITE ( io, 601 ) TRIM( char1 ), n, cpl_parent_id, npe_total, &
320	lower_left_coord_x, lower_left_coord_y, &
321	TRIM( cpl_name )
322	ENDDO
323	#endif
324
325	ENDIF
326	WRITE ( io, 99 )
327
328	!
329	!-- Numerical schemes
330	WRITE ( io, 110 )
331	IF ( rans_mode ) THEN
332	WRITE ( io, 124 ) TRIM( turbulence_closure ), 'RANS'
333	ELSE
334	WRITE ( io, 124 ) TRIM( turbulence_closure ), 'LES'
335	ENDIF
336	WRITE ( io, 121 ) TRIM( approximation )
337	IF ( psolver(1:7) == 'poisfft' ) THEN
338	WRITE ( io, 111 ) TRIM( fft_method )
339	IF ( transpose_compute_overlap ) WRITE( io, 115 )
340	ELSEIF ( psolver == 'sor' ) THEN
341	WRITE ( io, 112 ) nsor_ini, nsor, omega_sor
342	ELSEIF ( psolver(1:9) == 'multigrid' ) THEN
343	WRITE ( io, 135 ) TRIM(psolver), cycle_mg, maximum_grid_level, ngsrb
344	IF ( mg_cycles == -1 ) THEN
345	WRITE ( io, 140 ) residual_limit
346	ELSE
347	WRITE ( io, 141 ) mg_cycles
348	ENDIF
349	IF ( mg_switch_to_pe0_level == 0 ) THEN
350	WRITE ( io, 136 ) nxr_mg(1)-nxl_mg(1)+1, nyn_mg(1)-nys_mg(1)+1, &
351	nzt_mg(1)
352	ELSEIF ( mg_switch_to_pe0_level /= -1 ) THEN
353	WRITE ( io, 137 ) mg_switch_to_pe0_level, &
354	mg_loc_ind(2,0)-mg_loc_ind(1,0)+1, &
355	mg_loc_ind(4,0)-mg_loc_ind(3,0)+1, &
356	nzt_mg(mg_switch_to_pe0_level), &
357	nxr_mg(1)-nxl_mg(1)+1, nyn_mg(1)-nys_mg(1)+1, &
358	nzt_mg(1)
359	ENDIF
360	IF ( psolver == 'multigrid_noopt' .AND. masking_method ) WRITE ( io, 144 )
361	ENDIF
362	IF ( call_psolver_at_all_substeps .AND. timestep_scheme(1:5) == 'runge' ) &
363	THEN
364	WRITE ( io, 142 )
365	ENDIF
366
367	IF ( momentum_advec == 'pw-scheme' ) THEN
368	WRITE ( io, 113 )
369	ELSEIF (momentum_advec == 'ws-scheme' ) THEN
370	WRITE ( io, 503 )
371	ENDIF
372	IF ( scalar_advec == 'pw-scheme' ) THEN
373	WRITE ( io, 116 )
374	ELSEIF ( scalar_advec == 'ws-scheme' ) THEN
375	WRITE ( io, 504 )
376	ELSE
377	WRITE ( io, 118 )
378	ENDIF
379
380	WRITE ( io, 139 ) TRIM( loop_optimization )
381
382	IF ( galilei_transformation ) THEN
383	IF ( use_ug_for_galilei_tr ) THEN
384	char1 = '0.6 * geostrophic wind'
385	ELSE
386	char1 = 'mean wind in model domain'
387	ENDIF
388	IF ( simulated_time_at_begin == simulated_time ) THEN
389	char2 = 'at the start of the run'
390	ELSE
391	char2 = 'at the end of the run'
392	ENDIF
393	WRITE ( io, 119 ) TRIM( char1 ), TRIM( char2 ), &
394	advected_distance_x/1000.0_wp, &
395	advected_distance_y/1000.0_wp
396	ENDIF
397	WRITE ( io, 122 ) timestep_scheme
398	IF ( use_upstream_for_tke ) WRITE ( io, 143 )
399	IF ( rayleigh_damping_factor /= 0.0_wp ) THEN
400	IF ( .NOT. ocean_mode ) THEN
401	WRITE ( io, 123 ) 'above', rayleigh_damping_height, &
402	rayleigh_damping_factor
403	ELSE
404	WRITE ( io, 123 ) 'below', rayleigh_damping_height, &
405	rayleigh_damping_factor
406	ENDIF
407	ENDIF
408	IF ( neutral ) WRITE ( io, 131 ) pt_surface
409	IF ( humidity ) THEN
410	IF ( .NOT. bulk_cloud_model ) THEN
411	WRITE ( io, 129 )
412	ELSE
413	WRITE ( io, 130 )
414	ENDIF
415	ENDIF
416	IF ( passive_scalar ) WRITE ( io, 134 )
417	IF ( conserve_volume_flow ) THEN
418	WRITE ( io, 150 ) conserve_volume_flow_mode
419	IF ( TRIM( conserve_volume_flow_mode ) == 'bulk_velocity' ) THEN
420	WRITE ( io, 151 ) u_bulk, v_bulk
421	ENDIF
422	ELSEIF ( dp_external ) THEN
423	IF ( dp_smooth ) THEN
424	WRITE ( io, 152 ) dpdxy, dp_level_b, ', vertically smoothed.'
425	ELSE
426	WRITE ( io, 152 ) dpdxy, dp_level_b, '.'
427	ENDIF
428	ENDIF
429	WRITE ( io, 99 )
430
431	!
432	!-- Runtime and timestep information
433	WRITE ( io, 200 )
434	IF ( .NOT. dt_fixed ) THEN
435	WRITE ( io, 201 ) dt_max, cfl_factor
436	ELSE
437	WRITE ( io, 202 ) dt
438	ENDIF
439	WRITE ( io, 203 ) simulated_time_at_begin, end_time
440
441	IF ( time_restart /= 9999999.9_wp .AND. &
442	simulated_time_at_begin == simulated_time ) THEN
443	IF ( dt_restart == 9999999.9_wp ) THEN
444	WRITE ( io, 204 ) ' Restart at: ',time_restart
445	ELSE
446	WRITE ( io, 205 ) ' Restart at: ',time_restart, dt_restart
447	ENDIF
448	ENDIF
449
450	IF ( simulated_time_at_begin /= simulated_time ) THEN
451	i = MAX ( log_point_s(10)%counts, 1 )
452	IF ( ( simulated_time - simulated_time_at_begin ) == 0.0_wp ) THEN
453	cpuseconds_per_simulated_second = 0.0_wp
454	ELSE
455	cpuseconds_per_simulated_second = log_point_s(10)%sum / &
456	( simulated_time - &
457	simulated_time_at_begin )
458	ENDIF
459	WRITE ( io, 206 ) simulated_time, log_point_s(10)%sum, &
460	log_point_s(10)%sum / REAL( i, KIND=wp ), &
461	cpuseconds_per_simulated_second
462	IF ( time_restart /= 9999999.9_wp .AND. time_restart < end_time ) THEN
463	IF ( dt_restart == 9999999.9_wp ) THEN
464	WRITE ( io, 204 ) ' Next restart at: ',time_restart
465	ELSE
466	WRITE ( io, 205 ) ' Next restart at: ',time_restart, dt_restart
467	ENDIF
468	ENDIF
469	ENDIF
470
471
472	!
473	!-- Start time for coupled runs, if independent precursor runs for atmosphere
474	!-- and ocean are used or have been used. In this case, coupling_start_time
475	!-- defines the time when the coupling is switched on.
476	IF ( coupling_start_time /= 0.0_wp ) THEN
477	WRITE ( io, 207 ) coupling_start_time
478	ENDIF
479
480	!
481	!-- Computational grid
482	IF ( .NOT. ocean_mode ) THEN
483	WRITE ( io, 250 ) dx, dy
484
485	DO i = 1, number_stretch_level_start+1
486	WRITE ( io, 253 ) i, dz(i)
487	ENDDO
488
489	WRITE( io, 251 ) (nx+1)dx, (ny+1)dy, zu(nzt+1)
490
491	IF ( ANY( dz_stretch_level_start_index < nzt+1 ) ) THEN
492	WRITE( io, '(A)', advance='no') ' Vertical stretching starts at height:'
493	DO i = 1, number_stretch_level_start
494	WRITE ( io, '(F10.1,A3)', advance='no' ) dz_stretch_level_start(i), ' m,'
495	ENDDO
496	WRITE( io, '(/,A)', advance='no') ' Vertical stretching starts at index: '
497	DO i = 1, number_stretch_level_start
498	WRITE ( io, '(I12,A1)', advance='no' ) dz_stretch_level_start_index(i), ','
499	ENDDO
500	WRITE( io, '(/,A)', advance='no') ' Vertical stretching ends at height: '
501	DO i = 1, number_stretch_level_start
502	WRITE ( io, '(F10.1,A3)', advance='no' ) dz_stretch_level_end(i), ' m,'
503	ENDDO
504	WRITE( io, '(/,A)', advance='no') ' Vertical stretching ends at index: '
505	DO i = 1, number_stretch_level_start
506	WRITE ( io, '(I12,A1)', advance='no' ) dz_stretch_level_end_index(i), ','
507	ENDDO
508	WRITE( io, '(/,A)', advance='no') ' Factor used for stretching: '
509	DO i = 1, number_stretch_level_start
510	WRITE ( io, '(F12.3,A1)', advance='no' ) dz_stretch_factor_array(i), ','
511	ENDDO
512	ENDIF
513
514	ELSE
515	WRITE ( io, 250 ) dx, dy
516	DO i = 1, number_stretch_level_start+1
517	WRITE ( io, 253 ) i, dz(i)
518	ENDDO
519
520	WRITE ( io, 251 ) (nx+1)dx, (ny+1)dy, zu(0)
521
522	IF ( ANY( dz_stretch_level_start_index > 0 ) ) THEN
523	WRITE( io, '(A)', advance='no') ' Vertical stretching starts at height:'
524	DO i = 1, number_stretch_level_start
525	WRITE ( io, '(F10.1,A3)', advance='no' ) dz_stretch_level_start(i), ' m,'
526	ENDDO
527	WRITE( io, '(/,A)', advance='no') ' Vertical stretching starts at index: '
528	DO i = 1, number_stretch_level_start
529	WRITE ( io, '(I12,A1)', advance='no' ) dz_stretch_level_start_index(i), ','
530	ENDDO
531	WRITE( io, '(/,A)', advance='no') ' Vertical stretching ends at height: '
532	DO i = 1, number_stretch_level_start
533	WRITE ( io, '(F10.1,A3)', advance='no' ) dz_stretch_level_end(i), ' m,'
534	ENDDO
535	WRITE( io, '(/,A)', advance='no') ' Vertical stretching ends at index: '
536	DO i = 1, number_stretch_level_start
537	WRITE ( io, '(I12,A1)', advance='no' ) dz_stretch_level_end_index(i), ','
538	ENDDO
539	WRITE( io, '(/,A)', advance='no') ' Factor used for stretching: '
540	DO i = 1, number_stretch_level_start
541	WRITE ( io, '(F12.3,A1)', advance='no' ) dz_stretch_factor_array(i), ','
542	ENDDO
543	ENDIF
544	ENDIF
545	WRITE ( io, 254 ) nx, ny, nzt+1, MIN( nnx, nx+1 ), MIN( nny, ny+1 ), &
546	MIN( nnz+2, nzt+2 )
547	IF ( sloping_surface ) WRITE ( io, 260 ) alpha_surface
548
549	!
550	!-- Profile for the large scale vertial velocity
551	!-- Building output strings, starting with surface value
552	IF ( large_scale_subsidence ) THEN
553	temperatures = ' 0.0'
554	gradients = '------'
555	slices = ' 0'
556	coordinates = ' 0.0'
557	i = 1
558	DO WHILE ( subs_vertical_gradient_level_i(i) /= -9999 )
559
560	WRITE (coor_chr,'(E10.2,7X)') &
561	w_subs(subs_vertical_gradient_level_i(i))
562	temperatures = TRIM( temperatures ) // ' ' // TRIM( coor_chr )
563
564	WRITE (coor_chr,'(E10.2,7X)') subs_vertical_gradient(i)
565	gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
566
567	WRITE (coor_chr,'(I10,7X)') subs_vertical_gradient_level_i(i)
568	slices = TRIM( slices ) // ' ' // TRIM( coor_chr )
569
570	WRITE (coor_chr,'(F10.2,7X)') subs_vertical_gradient_level(i)
571	coordinates = TRIM( coordinates ) // ' ' // TRIM( coor_chr )
572
573	IF ( i == 10 ) THEN
574	EXIT
575	ELSE
576	i = i + 1
577	ENDIF
578
579	ENDDO
580
581
582	IF ( .NOT. large_scale_forcing ) THEN
583	WRITE ( io, 426 ) TRIM( coordinates ), TRIM( temperatures ), &
584	TRIM( gradients ), TRIM( slices )
585	ENDIF
586
587
588	ENDIF
589
590	!-- Profile of the geostrophic wind (component ug)
591	!-- Building output strings
592	WRITE ( ugcomponent, '(F6.2)' ) ug_surface
593	gradients = '------'
594	slices = ' 0'
595	coordinates = ' 0.0'
596	i = 1
597	DO WHILE ( ug_vertical_gradient_level_ind(i) /= -9999 )
598
599	WRITE (coor_chr,'(F6.2,1X)') ug(ug_vertical_gradient_level_ind(i))
600	ugcomponent = TRIM( ugcomponent ) // ' ' // TRIM( coor_chr )
601
602	WRITE (coor_chr,'(F6.2,1X)') ug_vertical_gradient(i)
603	gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
604
605	WRITE (coor_chr,'(I6,1X)') ug_vertical_gradient_level_ind(i)
606	slices = TRIM( slices ) // ' ' // TRIM( coor_chr )
607
608	WRITE (coor_chr,'(F6.1,1X)') ug_vertical_gradient_level(i)
609	coordinates = TRIM( coordinates ) // ' ' // TRIM( coor_chr )
610
611	IF ( i == 10 ) THEN
612	EXIT
613	ELSE
614	i = i + 1
615	ENDIF
616
617	ENDDO
618
619	IF ( .NOT. large_scale_forcing ) THEN
620	WRITE ( io, 423 ) TRIM( coordinates ), TRIM( ugcomponent ), &
621	TRIM( gradients ), TRIM( slices )
622	ENDIF
623
624	!-- Profile of the geostrophic wind (component vg)
625	!-- Building output strings
626	WRITE ( vgcomponent, '(F6.2)' ) vg_surface
627	gradients = '------'
628	slices = ' 0'
629	coordinates = ' 0.0'
630	i = 1
631	DO WHILE ( vg_vertical_gradient_level_ind(i) /= -9999 )
632
633	WRITE (coor_chr,'(F6.2,1X)') vg(vg_vertical_gradient_level_ind(i))
634	vgcomponent = TRIM( vgcomponent ) // ' ' // TRIM( coor_chr )
635
636	WRITE (coor_chr,'(F6.2,1X)') vg_vertical_gradient(i)
637	gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
638
639	WRITE (coor_chr,'(I6,1X)') vg_vertical_gradient_level_ind(i)
640	slices = TRIM( slices ) // ' ' // TRIM( coor_chr )
641
642	WRITE (coor_chr,'(F6.1,1X)') vg_vertical_gradient_level(i)
643	coordinates = TRIM( coordinates ) // ' ' // TRIM( coor_chr )
644
645	IF ( i == 10 ) THEN
646	EXIT
647	ELSE
648	i = i + 1
649	ENDIF
650
651	ENDDO
652
653	IF ( .NOT. large_scale_forcing ) THEN
654	WRITE ( io, 424 ) TRIM( coordinates ), TRIM( vgcomponent ), &
655	TRIM( gradients ), TRIM( slices )
656	ENDIF
657
658	!
659	!-- Topography
660	WRITE ( io, 270 ) topography
661	SELECT CASE ( TRIM( topography ) )
662
663	CASE ( 'flat' )
664	! no actions necessary
665
666	CASE ( 'single_building' )
667	blx = INT( building_length_x / dx )
668	bly = INT( building_length_y / dy )
669	bh = MINLOC( ABS( zw - building_height ), 1 ) - 1
670	IF ( ABS( zw(bh ) - building_height ) == &
671	ABS( zw(bh+1) - building_height ) ) bh = bh + 1
672
673	IF ( building_wall_left == 9999999.9_wp ) THEN
674	building_wall_left = ( nx + 1 - blx ) / 2 * dx
675	ENDIF
676	bxl = INT ( building_wall_left / dx + 0.5_wp )
677	bxr = bxl + blx
678
679	IF ( building_wall_south == 9999999.9_wp ) THEN
680	building_wall_south = ( ny + 1 - bly ) / 2 * dy
681	ENDIF
682	bys = INT ( building_wall_south / dy + 0.5_wp )
683	byn = bys + bly
684
685	WRITE ( io, 271 ) building_length_x, building_length_y, &
686	building_height, bxl, bxr, bys, byn
687
688	CASE ( 'single_street_canyon' )
689	ch = MINLOC( ABS( zw - canyon_height ), 1 ) - 1
690	IF ( ABS( zw(ch ) - canyon_height ) == &
691	ABS( zw(ch+1) - canyon_height ) ) ch = ch + 1
692	IF ( canyon_width_x /= 9999999.9_wp ) THEN
693	!
694	!-- Street canyon in y direction
695	cwx = NINT( canyon_width_x / dx )
696	IF ( canyon_wall_left == 9999999.9_wp ) THEN
697	canyon_wall_left = ( nx + 1 - cwx ) / 2 * dx
698	ENDIF
699	cxl = NINT( canyon_wall_left / dx )
700	cxr = cxl + cwx
701	WRITE ( io, 272 ) 'y', canyon_height, ch, 'u', cxl, cxr
702
703	ELSEIF ( canyon_width_y /= 9999999.9_wp ) THEN
704	!
705	!-- Street canyon in x direction
706	cwy = NINT( canyon_width_y / dy )
707	IF ( canyon_wall_south == 9999999.9_wp ) THEN
708	canyon_wall_south = ( ny + 1 - cwy ) / 2 * dy
709	ENDIF
710	cys = NINT( canyon_wall_south / dy )
711	cyn = cys + cwy
712	WRITE ( io, 272 ) 'x', canyon_height, ch, 'v', cys, cyn
713	ENDIF
714
715	CASE ( 'tunnel' )
716	IF ( tunnel_width_x /= 9999999.9_wp ) THEN
717	!
718	!-- Tunnel axis in y direction
719	IF ( tunnel_length == 9999999.9_wp .OR. &
720	tunnel_length >= ( nx + 1 ) * dx ) THEN
721	WRITE ( io, 273 ) 'y', tunnel_height, tunnel_wall_depth, &
722	tunnel_width_x
723	ELSE
724	WRITE ( io, 274 ) 'y', tunnel_height, tunnel_wall_depth, &
725	tunnel_width_x, tunnel_length
726	ENDIF
727
728	ELSEIF ( tunnel_width_y /= 9999999.9_wp ) THEN
729	!
730	!-- Tunnel axis in x direction
731	IF ( tunnel_length == 9999999.9_wp .OR. &
732	tunnel_length >= ( ny + 1 ) * dy ) THEN
733	WRITE ( io, 273 ) 'x', tunnel_height, tunnel_wall_depth, &
734	tunnel_width_y
735	ELSE
736	WRITE ( io, 274 ) 'x', tunnel_height, tunnel_wall_depth, &
737	tunnel_width_y, tunnel_length
738	ENDIF
739	ENDIF
740
741	END SELECT
742
743	IF ( TRIM( topography ) /= 'flat' ) THEN
744	IF ( TRIM( topography_grid_convention ) == ' ' ) THEN
745	IF ( TRIM( topography ) == 'single_building' .OR. &
746	TRIM( topography ) == 'single_street_canyon' ) THEN
747	WRITE ( io, 278 )
748	ELSEIF ( TRIM( topography ) == 'read_from_file' ) THEN
749	WRITE ( io, 279 )
750	ENDIF
751	ELSEIF ( TRIM( topography_grid_convention ) == 'cell_edge' ) THEN
752	WRITE ( io, 278 )
753	ELSEIF ( TRIM( topography_grid_convention ) == 'cell_center' ) THEN
754	WRITE ( io, 279 )
755	ENDIF
756	ENDIF
757
758	!-- Complex terrain
759	IF ( complex_terrain ) THEN
760	WRITE( io, 280 )
761	IF ( turbulent_inflow ) THEN
762	WRITE( io, 281 ) zu(topo_top_ind(0,0,0))
763	ENDIF
764	IF ( TRIM( initializing_actions ) == 'cyclic_fill' ) THEN
765	WRITE( io, 282 )
766	ENDIF
767	ENDIF
768	!
769	!-- Boundary conditions
770	IF ( ibc_p_b == 0 ) THEN
771	r_lower = 'p(0) = 0 \|'
772	ELSEIF ( ibc_p_b == 1 ) THEN
773	r_lower = 'p(0) = p(1) \|'
774	ENDIF
775	IF ( ibc_p_t == 0 ) THEN
776	r_upper = 'p(nzt+1) = 0 \|'
777	ELSE
778	r_upper = 'p(nzt+1) = p(nzt) \|'
779	ENDIF
780
781	IF ( ibc_uv_b == 0 ) THEN
782	r_lower = TRIM( r_lower ) // ' uv(0) = -uv(1) \|'
783	ELSE
784	r_lower = TRIM( r_lower ) // ' uv(0) = uv(1) \|'
785	ENDIF
786	IF ( TRIM( bc_uv_t ) == 'dirichlet_0' ) THEN
787	r_upper = TRIM( r_upper ) // ' uv(nzt+1) = 0 \|'
788	ELSEIF ( ibc_uv_t == 0 ) THEN
789	r_upper = TRIM( r_upper ) // ' uv(nzt+1) = ug(nzt+1), vg(nzt+1) \|'
790	ELSE
791	r_upper = TRIM( r_upper ) // ' uv(nzt+1) = uv(nzt) \|'
792	ENDIF
793
794	IF ( ibc_pt_b == 0 ) THEN
795	IF ( land_surface ) THEN
796	r_lower = TRIM( r_lower ) // ' pt(0) = from soil model'
797	ELSE
798	r_lower = TRIM( r_lower ) // ' pt(0) = pt_surface'
799	ENDIF
800	ELSEIF ( ibc_pt_b == 1 ) THEN
801	r_lower = TRIM( r_lower ) // ' pt(0) = pt(1)'
802	ELSEIF ( ibc_pt_b == 2 ) THEN
803	r_lower = TRIM( r_lower ) // ' pt(0) = from coupled model'
804	ENDIF
805	IF ( ibc_pt_t == 0 ) THEN
806	r_upper = TRIM( r_upper ) // ' pt(nzt+1) = pt_top'
807	ELSEIF( ibc_pt_t == 1 ) THEN
808	r_upper = TRIM( r_upper ) // ' pt(nzt+1) = pt(nzt)'
809	ELSEIF( ibc_pt_t == 2 ) THEN
810	r_upper = TRIM( r_upper ) // ' pt(nzt+1) = pt(nzt) + dpt/dz_ini'
811
812	ENDIF
813
814	WRITE ( io, 300 ) r_lower, r_upper
815
816	IF ( .NOT. constant_diffusion ) THEN
817	IF ( ibc_e_b == 1 ) THEN
818	r_lower = 'e(0) = e(1)'
819	ELSE
820	r_lower = 'e(0) = e(1) = (u/0.1)*2'
821	ENDIF
822	r_upper = 'e(nzt+1) = e(nzt) = e(nzt-1)'
823
824	WRITE ( io, 301 ) 'e', r_lower, r_upper
825
826	ENDIF
827
828	IF ( ocean_mode ) THEN
829	r_lower = 'sa(0) = sa(1)'
830	IF ( ibc_sa_t == 0 ) THEN
831	r_upper = 'sa(nzt+1) = sa_surface'
832	ELSE
833	r_upper = 'sa(nzt+1) = sa(nzt)'
834	ENDIF
835	WRITE ( io, 301 ) 'sa', r_lower, r_upper
836	ENDIF
837
838	IF ( humidity ) THEN
839	IF ( ibc_q_b == 0 ) THEN
840	IF ( land_surface ) THEN
841	r_lower = 'q(0) = from soil model'
842	ELSE
843	r_lower = 'q(0) = q_surface'
844	ENDIF
845
846	ELSE
847	r_lower = 'q(0) = q(1)'
848	ENDIF
849	IF ( ibc_q_t == 0 ) THEN
850	r_upper = 'q(nzt+1) = q_top'
851	ELSE
852	r_upper = 'q(nzt+1) = q(nzt) + dq/dz'
853	ENDIF
854	WRITE ( io, 301 ) 'q', r_lower, r_upper
855	ENDIF
856
857	IF ( passive_scalar ) THEN
858	IF ( ibc_s_b == 0 ) THEN
859	r_lower = 's(0) = s_surface'
860	ELSE
861	r_lower = 's(0) = s(1)'
862	ENDIF
863	IF ( ibc_s_t == 0 ) THEN
864	r_upper = 's(nzt+1) = s_top'
865	ELSEIF ( ibc_s_t == 1 ) THEN
866	r_upper = 's(nzt+1) = s(nzt)'
867	ELSEIF ( ibc_s_t == 2 ) THEN
868	r_upper = 's(nzt+1) = s(nzt) + ds/dz'
869	ENDIF
870	WRITE ( io, 301 ) 's', r_lower, r_upper
871	ENDIF
872
873	IF ( use_surface_fluxes ) THEN
874	WRITE ( io, 303 )
875	IF ( constant_heatflux ) THEN
876	IF ( large_scale_forcing .AND. lsf_surf ) THEN
877	IF ( surf_def_h(0)%ns >= 1 ) WRITE ( io, 306 ) surf_def_h(0)%shf(1)
878	ELSE
879	WRITE ( io, 306 ) surface_heatflux
880	ENDIF
881	IF ( random_heatflux ) WRITE ( io, 307 )
882	ENDIF
883	IF ( humidity .AND. constant_waterflux ) THEN
884	IF ( large_scale_forcing .AND. lsf_surf ) THEN
885	WRITE ( io, 311 ) surf_def_h(0)%qsws(1)
886	ELSE
887	WRITE ( io, 311 ) surface_waterflux
888	ENDIF
889	ENDIF
890	IF ( passive_scalar .AND. constant_scalarflux ) THEN
891	WRITE ( io, 313 ) surface_scalarflux
892	ENDIF
893	ENDIF
894
895	IF ( use_top_fluxes ) THEN
896	WRITE ( io, 304 )
897	IF ( coupling_mode == 'uncoupled' ) THEN
898	WRITE ( io, 320 ) top_momentumflux_u, top_momentumflux_v
899	IF ( constant_top_heatflux ) THEN
900	WRITE ( io, 306 ) top_heatflux
901	ENDIF
902	ELSEIF ( coupling_mode == 'ocean_to_atmosphere' ) THEN
903	WRITE ( io, 316 )
904	ENDIF
905	IF ( ocean_mode .AND. constant_top_salinityflux ) &
906	WRITE ( io, 309 ) top_salinityflux
907	IF ( humidity ) WRITE ( io, 315 )
908	IF ( passive_scalar .AND. constant_top_scalarflux ) &
909	WRITE ( io, 302 ) top_scalarflux
910	ENDIF
911
912	IF ( constant_flux_layer ) THEN
913	WRITE ( io, 305 ) (zu(1)-zu(0)), roughness_length, &
914	z0h_factor*roughness_length, kappa, &
915	zeta_min, zeta_max
916	IF ( .NOT. constant_heatflux ) WRITE ( io, 308 )
917	IF ( humidity .AND. .NOT. constant_waterflux ) THEN
918	WRITE ( io, 312 )
919	ENDIF
920	IF ( passive_scalar .AND. .NOT. constant_scalarflux ) THEN
921	WRITE ( io, 314 )
922	ENDIF
923	ELSE
924	IF ( INDEX(initializing_actions, 'set_1d-model_profiles') /= 0 ) THEN
925	WRITE ( io, 310 ) zeta_min, zeta_max
926	ENDIF
927	ENDIF
928
929	WRITE ( io, 317 ) bc_lr, bc_ns
930	IF ( .NOT. bc_lr_cyc .OR. .NOT. bc_ns_cyc ) THEN
931	WRITE ( io, 318 ) use_cmax, pt_damping_width, pt_damping_factor
932	IF ( turbulent_inflow ) THEN
933	IF ( y_shift == 0 ) THEN
934	WRITE ( io, 319 ) recycling_width, recycling_plane, &
935	inflow_damping_height, inflow_damping_width
936	ELSE
937	WRITE ( io, 322 ) y_shift, recycling_width, recycling_plane, &
938	inflow_damping_height, inflow_damping_width
939	END IF
940	ENDIF
941	IF ( turbulent_outflow ) THEN
942	WRITE ( io, 323 ) outflow_source_plane, INT(outflow_source_plane/dx)
943	ENDIF
944	ENDIF
945
946	!
947	!-- Initial Profiles
948	WRITE ( io, 321 )
949	!
950	!-- Initial wind profiles
951	IF ( u_profile(1) /= 9999999.9_wp ) WRITE ( io, 427 )
952
953	!
954	!-- Initial temperature profile
955	!-- Building output strings, starting with surface temperature
956	WRITE ( temperatures, '(F6.2)' ) pt_surface
957	gradients = '------'
958	slices = ' 0'
959	coordinates = ' 0.0'
960	i = 1
961	DO WHILE ( pt_vertical_gradient_level_ind(i) /= -9999 )
962
963	WRITE (coor_chr,'(F7.2)') pt_init(pt_vertical_gradient_level_ind(i))
964	temperatures = TRIM( temperatures ) // ' ' // TRIM( coor_chr )
965
966	WRITE (coor_chr,'(F7.2)') pt_vertical_gradient(i)
967	gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
968
969	WRITE (coor_chr,'(I7)') pt_vertical_gradient_level_ind(i)
970	slices = TRIM( slices ) // ' ' // TRIM( coor_chr )
971
972	WRITE (coor_chr,'(F7.1)') pt_vertical_gradient_level(i)
973	coordinates = TRIM( coordinates ) // ' ' // TRIM( coor_chr )
974
975	IF ( i == 10 ) THEN
976	EXIT
977	ELSE
978	i = i + 1
979	ENDIF
980
981	ENDDO
982
983	IF ( .NOT. nudging ) THEN
984	WRITE ( io, 420 ) TRIM( coordinates ), TRIM( temperatures ), &
985	TRIM( gradients ), TRIM( slices )
986	ELSE
987	WRITE ( io, 428 )
988	ENDIF
989
990	!
991	!-- Initial humidity profile
992	!-- Building output strings, starting with surface humidity
993	IF ( humidity ) THEN
994	WRITE ( temperatures, '(E8.1)' ) q_surface
995	gradients = '--------'
996	slices = ' 0'
997	coordinates = ' 0.0'
998	i = 1
999	DO WHILE ( q_vertical_gradient_level_ind(i) /= -9999 )
1000
1001	WRITE (coor_chr,'(E8.1,4X)') q_init(q_vertical_gradient_level_ind(i))
1002	temperatures = TRIM( temperatures ) // ' ' // TRIM( coor_chr )
1003
1004	WRITE (coor_chr,'(E8.1,4X)') q_vertical_gradient(i)
1005	gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
1006
1007	WRITE (coor_chr,'(I8,4X)') q_vertical_gradient_level_ind(i)
1008	slices = TRIM( slices ) // ' ' // TRIM( coor_chr )
1009
1010	WRITE (coor_chr,'(F8.1,4X)') q_vertical_gradient_level(i)
1011	coordinates = TRIM( coordinates ) // ' ' // TRIM( coor_chr )
1012
1013	IF ( i == 10 ) THEN
1014	EXIT
1015	ELSE
1016	i = i + 1
1017	ENDIF
1018
1019	ENDDO
1020
1021	IF ( .NOT. nudging ) THEN
1022	WRITE ( io, 421 ) TRIM( coordinates ), TRIM( temperatures ), &
1023	TRIM( gradients ), TRIM( slices )
1024	ENDIF
1025	ENDIF
1026	!
1027	!-- Initial scalar profile
1028	!-- Building output strings, starting with surface humidity
1029	IF ( passive_scalar ) THEN
1030	WRITE ( temperatures, '(E8.1)' ) s_surface
1031	gradients = '--------'
1032	slices = ' 0'
1033	coordinates = ' 0.0'
1034	i = 1
1035	DO WHILE ( s_vertical_gradient_level_ind(i) /= -9999 )
1036
1037	WRITE (coor_chr,'(E8.1,4X)') s_init(s_vertical_gradient_level_ind(i))
1038	temperatures = TRIM( temperatures ) // ' ' // TRIM( coor_chr )
1039
1040	WRITE (coor_chr,'(E8.1,4X)') s_vertical_gradient(i)
1041	gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
1042
1043	WRITE (coor_chr,'(I8,4X)') s_vertical_gradient_level_ind(i)
1044	slices = TRIM( slices ) // ' ' // TRIM( coor_chr )
1045
1046	WRITE (coor_chr,'(F8.1,4X)') s_vertical_gradient_level(i)
1047	coordinates = TRIM( coordinates ) // ' ' // TRIM( coor_chr )
1048
1049	IF ( i == 10 ) THEN
1050	EXIT
1051	ELSE
1052	i = i + 1
1053	ENDIF
1054
1055	ENDDO
1056
1057	WRITE ( io, 422 ) TRIM( coordinates ), TRIM( temperatures ), &
1058	TRIM( gradients ), TRIM( slices )
1059	ENDIF
1060
1061	!
1062	!-- Initial salinity profile
1063	!-- Building output strings, starting with surface salinity
1064	IF ( ocean_mode ) THEN
1065	WRITE ( temperatures, '(F6.2)' ) sa_surface
1066	gradients = '------'
1067	slices = ' 0'
1068	coordinates = ' 0.0'
1069	i = 1
1070	DO WHILE ( sa_vertical_gradient_level_ind(i) /= -9999 )
1071
1072	WRITE (coor_chr,'(F7.2)') sa_init(sa_vertical_gradient_level_ind(i))
1073	temperatures = TRIM( temperatures ) // ' ' // TRIM( coor_chr )
1074
1075	WRITE (coor_chr,'(F7.2)') sa_vertical_gradient(i)
1076	gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
1077
1078	WRITE (coor_chr,'(I7)') sa_vertical_gradient_level_ind(i)
1079	slices = TRIM( slices ) // ' ' // TRIM( coor_chr )
1080
1081	WRITE (coor_chr,'(F7.1)') sa_vertical_gradient_level(i)
1082	coordinates = TRIM( coordinates ) // ' ' // TRIM( coor_chr )
1083
1084	IF ( i == 10 ) THEN
1085	EXIT
1086	ELSE
1087	i = i + 1
1088	ENDIF
1089
1090	ENDDO
1091
1092	WRITE ( io, 425 ) TRIM( coordinates ), TRIM( temperatures ), &
1093	TRIM( gradients ), TRIM( slices )
1094	ENDIF
1095
1096
1097	!
1098	!-- Listing of 1D-profiles
1099	WRITE ( io, 325 ) dt_dopr_listing
1100	IF ( averaging_interval_pr /= 0.0_wp ) THEN
1101	WRITE ( io, 326 ) averaging_interval_pr, dt_averaging_input_pr
1102	ENDIF
1103
1104	!
1105	!-- DATA output
1106	WRITE ( io, 330 )
1107	IF ( averaging_interval_pr /= 0.0_wp ) THEN
1108	WRITE ( io, 326 ) averaging_interval_pr, dt_averaging_input_pr
1109	ENDIF
1110
1111	!
1112	!-- 1D-profiles
1113	dopr_chr = 'Profile:'
1114	IF ( dopr_n /= 0 ) THEN
1115	WRITE ( io, 331 )
1116
1117	output_format = ''
1118	output_format = netcdf_data_format_string
1119	IF ( netcdf_deflate == 0 ) THEN
1120	WRITE ( io, 344 ) output_format
1121	ELSE
1122	WRITE ( io, 354 ) TRIM( output_format ), netcdf_deflate
1123	ENDIF
1124
1125	DO i = 1, dopr_n
1126	dopr_chr = TRIM( dopr_chr ) // ' ' // TRIM( data_output_pr(i) ) // ','
1127	IF ( LEN_TRIM( dopr_chr ) >= 60 ) THEN
1128	WRITE ( io, 332 ) dopr_chr
1129	dopr_chr = ' :'
1130	ENDIF
1131	ENDDO
1132
1133	IF ( dopr_chr /= '' ) THEN
1134	WRITE ( io, 332 ) dopr_chr
1135	ENDIF
1136	WRITE ( io, 333 ) dt_dopr, averaging_interval_pr, dt_averaging_input_pr
1137	IF ( skip_time_dopr /= 0.0_wp ) WRITE ( io, 339 ) skip_time_dopr
1138	ENDIF
1139
1140	!
1141	!-- 2D-arrays
1142	DO av = 0, 1
1143
1144	i = 1
1145	do2d_xy = ''
1146	do2d_xz = ''
1147	do2d_yz = ''
1148	DO WHILE ( do2d(av,i) /= ' ' )
1149
1150	l = MAX( 2, LEN_TRIM( do2d(av,i) ) )
1151	do2d_mode = do2d(av,i)(l-1:l)
1152
1153	SELECT CASE ( do2d_mode )
1154	CASE ( 'xy' )
1155	ll = LEN_TRIM( do2d_xy )
1156	do2d_xy = do2d_xy(1:ll) // ' ' // do2d(av,i)(1:l-3) // ','
1157	CASE ( 'xz' )
1158	ll = LEN_TRIM( do2d_xz )
1159	do2d_xz = do2d_xz(1:ll) // ' ' // do2d(av,i)(1:l-3) // ','
1160	CASE ( 'yz' )
1161	ll = LEN_TRIM( do2d_yz )
1162	do2d_yz = do2d_yz(1:ll) // ' ' // do2d(av,i)(1:l-3) // ','
1163	END SELECT
1164
1165	i = i + 1
1166
1167	ENDDO
1168
1169	IF ( ( ( do2d_xy /= '' .AND. section(1,1) /= -9999 ) .OR. &
1170	( do2d_xz /= '' .AND. section(1,2) /= -9999 ) .OR. &
1171	( do2d_yz /= '' .AND. section(1,3) /= -9999 ) ) ) THEN
1172
1173	IF ( av == 0 ) THEN
1174	WRITE ( io, 334 ) ''
1175	ELSE
1176	WRITE ( io, 334 ) '(time-averaged)'
1177	ENDIF
1178
1179	IF ( do2d_at_begin ) THEN
1180	begin_chr = 'and at the start'
1181	ELSE
1182	begin_chr = ''
1183	ENDIF
1184
1185	output_format = ''
1186	output_format = netcdf_data_format_string
1187	IF ( netcdf_deflate == 0 ) THEN
1188	WRITE ( io, 344 ) output_format
1189	ELSE
1190	WRITE ( io, 354 ) TRIM( output_format ), netcdf_deflate
1191	ENDIF
1192
1193	IF ( do2d_xy /= '' .AND. section(1,1) /= -9999 ) THEN
1194	i = 1
1195	slices = '/'
1196	coordinates = '/'
1197	!
1198	!-- Building strings with index and coordinate information of the
1199	!-- slices
1200	DO WHILE ( section(i,1) /= -9999 )
1201
1202	WRITE (section_chr,'(I5)') section(i,1)
1203	section_chr = ADJUSTL( section_chr )
1204	slices = TRIM( slices ) // TRIM( section_chr ) // '/'
1205
1206	IF ( section(i,1) == -1 ) THEN
1207	WRITE (coor_chr,'(F10.1)') -1.0_wp
1208	ELSE
1209	WRITE (coor_chr,'(F10.1)') zu(section(i,1))
1210	ENDIF
1211	coor_chr = ADJUSTL( coor_chr )
1212	coordinates = TRIM( coordinates ) // TRIM( coor_chr ) // '/'
1213
1214	i = i + 1
1215	ENDDO
1216	IF ( av == 0 ) THEN
1217	WRITE ( io, 335 ) 'XY', do2d_xy, dt_do2d_xy, &
1218	TRIM( begin_chr ), 'k', TRIM( slices ), &
1219	TRIM( coordinates )
1220	IF ( skip_time_do2d_xy /= 0.0_wp ) THEN
1221	WRITE ( io, 339 ) skip_time_do2d_xy
1222	ENDIF
1223	ELSE
1224	WRITE ( io, 342 ) 'XY', do2d_xy, dt_data_output_av, &
1225	TRIM( begin_chr ), averaging_interval, &
1226	dt_averaging_input, 'k', TRIM( slices ), &
1227	TRIM( coordinates )
1228	IF ( skip_time_data_output_av /= 0.0_wp ) THEN
1229	WRITE ( io, 339 ) skip_time_data_output_av
1230	ENDIF
1231	ENDIF
1232	IF ( netcdf_data_format > 4 ) THEN
1233	WRITE ( io, 352 ) ntdim_2d_xy(av)
1234	ELSE
1235	WRITE ( io, 353 )
1236	ENDIF
1237	ENDIF
1238
1239	IF ( do2d_xz /= '' .AND. section(1,2) /= -9999 ) THEN
1240	i = 1
1241	slices = '/'
1242	coordinates = '/'
1243	!
1244	!-- Building strings with index and coordinate information of the
1245	!-- slices
1246	DO WHILE ( section(i,2) /= -9999 )
1247
1248	WRITE (section_chr,'(I5)') section(i,2)
1249	section_chr = ADJUSTL( section_chr )
1250	slices = TRIM( slices ) // TRIM( section_chr ) // '/'
1251
1252	WRITE (coor_chr,'(F10.1)') section(i,2) * dy
1253	coor_chr = ADJUSTL( coor_chr )
1254	coordinates = TRIM( coordinates ) // TRIM( coor_chr ) // '/'
1255
1256	i = i + 1
1257	ENDDO
1258	IF ( av == 0 ) THEN
1259	WRITE ( io, 335 ) 'XZ', do2d_xz, dt_do2d_xz, &
1260	TRIM( begin_chr ), 'j', TRIM( slices ), &
1261	TRIM( coordinates )
1262	IF ( skip_time_do2d_xz /= 0.0_wp ) THEN
1263	WRITE ( io, 339 ) skip_time_do2d_xz
1264	ENDIF
1265	ELSE
1266	WRITE ( io, 342 ) 'XZ', do2d_xz, dt_data_output_av, &
1267	TRIM( begin_chr ), averaging_interval, &
1268	dt_averaging_input, 'j', TRIM( slices ), &
1269	TRIM( coordinates )
1270	IF ( skip_time_data_output_av /= 0.0_wp ) THEN
1271	WRITE ( io, 339 ) skip_time_data_output_av
1272	ENDIF
1273	ENDIF
1274	IF ( netcdf_data_format > 4 ) THEN
1275	WRITE ( io, 352 ) ntdim_2d_xz(av)
1276	ELSE
1277	WRITE ( io, 353 )
1278	ENDIF
1279	ENDIF
1280
1281	IF ( do2d_yz /= '' .AND. section(1,3) /= -9999 ) THEN
1282	i = 1
1283	slices = '/'
1284	coordinates = '/'
1285	!
1286	!-- Building strings with index and coordinate information of the
1287	!-- slices
1288	DO WHILE ( section(i,3) /= -9999 )
1289
1290	WRITE (section_chr,'(I5)') section(i,3)
1291	section_chr = ADJUSTL( section_chr )
1292	slices = TRIM( slices ) // TRIM( section_chr ) // '/'
1293
1294	WRITE (coor_chr,'(F10.1)') section(i,3) * dx
1295	coor_chr = ADJUSTL( coor_chr )
1296	coordinates = TRIM( coordinates ) // TRIM( coor_chr ) // '/'
1297
1298	i = i + 1
1299	ENDDO
1300	IF ( av == 0 ) THEN
1301	WRITE ( io, 335 ) 'YZ', do2d_yz, dt_do2d_yz, &
1302	TRIM( begin_chr ), 'i', TRIM( slices ), &
1303	TRIM( coordinates )
1304	IF ( skip_time_do2d_yz /= 0.0_wp ) THEN
1305	WRITE ( io, 339 ) skip_time_do2d_yz
1306	ENDIF
1307	ELSE
1308	WRITE ( io, 342 ) 'YZ', do2d_yz, dt_data_output_av, &
1309	TRIM( begin_chr ), averaging_interval, &
1310	dt_averaging_input, 'i', TRIM( slices ), &
1311	TRIM( coordinates )
1312	IF ( skip_time_data_output_av /= 0.0_wp ) THEN
1313	WRITE ( io, 339 ) skip_time_data_output_av
1314	ENDIF
1315	ENDIF
1316	IF ( netcdf_data_format > 4 ) THEN
1317	WRITE ( io, 352 ) ntdim_2d_yz(av)
1318	ELSE
1319	WRITE ( io, 353 )
1320	ENDIF
1321	ENDIF
1322
1323	ENDIF
1324
1325	ENDDO
1326
1327	!
1328	!-- 3d-arrays
1329	DO av = 0, 1
1330
1331	i = 1
1332	do3d_chr = ''
1333	DO WHILE ( do3d(av,i) /= ' ' )
1334
1335	do3d_chr = TRIM( do3d_chr ) // ' ' // TRIM( do3d(av,i) ) // ','
1336	i = i + 1
1337
1338	ENDDO
1339
1340	IF ( do3d_chr /= '' ) THEN
1341	IF ( av == 0 ) THEN
1342	WRITE ( io, 336 ) ''
1343	ELSE
1344	WRITE ( io, 336 ) '(time-averaged)'
1345	ENDIF
1346
1347	output_format = netcdf_data_format_string
1348	IF ( netcdf_deflate == 0 ) THEN
1349	WRITE ( io, 344 ) output_format
1350	ELSE
1351	WRITE ( io, 354 ) TRIM( output_format ), netcdf_deflate
1352	ENDIF
1353
1354	IF ( do3d_at_begin ) THEN
1355	begin_chr = 'and at the start'
1356	ELSE
1357	begin_chr = ''
1358	ENDIF
1359	IF ( av == 0 ) THEN
1360	WRITE ( io, 337 ) do3d_chr, dt_do3d, TRIM( begin_chr ), &
1361	zu(nz_do3d), nz_do3d
1362	ELSE
1363	WRITE ( io, 343 ) do3d_chr, dt_data_output_av, &
1364	TRIM( begin_chr ), averaging_interval, &
1365	dt_averaging_input, zu(nz_do3d), nz_do3d
1366	ENDIF
1367
1368	IF ( netcdf_data_format > 4 ) THEN
1369	WRITE ( io, 352 ) ntdim_3d(av)
1370	ELSE
1371	WRITE ( io, 353 )
1372	ENDIF
1373
1374	IF ( av == 0 ) THEN
1375	IF ( skip_time_do3d /= 0.0_wp ) THEN
1376	WRITE ( io, 339 ) skip_time_do3d
1377	ENDIF
1378	ELSE
1379	IF ( skip_time_data_output_av /= 0.0_wp ) THEN
1380	WRITE ( io, 339 ) skip_time_data_output_av
1381	ENDIF
1382	ENDIF
1383
1384	ENDIF
1385
1386	ENDDO
1387
1388	!
1389	!-- masked arrays
1390	IF ( masks > 0 ) WRITE ( io, 345 ) &
1391	mask_scale_x, mask_scale_y, mask_scale_z
1392	DO mid = 1, masks
1393	DO av = 0, 1
1394
1395	i = 1
1396	domask_chr = ''
1397	DO WHILE ( domask(mid,av,i) /= ' ' )
1398	domask_chr = TRIM( domask_chr ) // ' ' // &
1399	TRIM( domask(mid,av,i) ) // ','
1400	i = i + 1
1401	ENDDO
1402
1403	IF ( domask_chr /= '' ) THEN
1404	IF ( av == 0 ) THEN
1405	WRITE ( io, 346 ) '', mid
1406	ELSE
1407	WRITE ( io, 346 ) ' (time-averaged)', mid
1408	ENDIF
1409
1410	output_format = netcdf_data_format_string
1411	!-- Parallel output not implemented for mask data, hence
1412	!-- output_format must be adjusted.
1413	IF ( netcdf_data_format == 5 ) output_format = 'netCDF4/HDF5'
1414	IF ( netcdf_data_format == 6 ) output_format = 'netCDF4/HDF5 classic'
1415	IF ( netcdf_deflate == 0 ) THEN
1416	WRITE ( io, 344 ) output_format
1417	ELSE
1418	WRITE ( io, 354 ) TRIM( output_format ), netcdf_deflate
1419	ENDIF
1420
1421	IF ( av == 0 ) THEN
1422	WRITE ( io, 347 ) domask_chr, dt_domask(mid)
1423	ELSE
1424	WRITE ( io, 348 ) domask_chr, dt_data_output_av, &
1425	averaging_interval, dt_averaging_input
1426	ENDIF
1427
1428	IF ( av == 0 ) THEN
1429	IF ( skip_time_domask(mid) /= 0.0_wp ) THEN
1430	WRITE ( io, 339 ) skip_time_domask(mid)
1431	ENDIF
1432	ELSE
1433	IF ( skip_time_data_output_av /= 0.0_wp ) THEN
1434	WRITE ( io, 339 ) skip_time_data_output_av
1435	ENDIF
1436	ENDIF
1437	!
1438	!-- output locations
1439	DO dim = 1, 3
1440	IF ( mask(mid,dim,1) >= 0.0_wp ) THEN
1441	count = 0
1442	DO WHILE ( mask(mid,dim,count+1) >= 0.0_wp )
1443	count = count + 1
1444	ENDDO
1445	WRITE ( io, 349 ) dir(dim), dir(dim), mid, dir(dim), &
1446	mask(mid,dim,:count)
1447	ELSEIF ( mask_loop(mid,dim,1) < 0.0_wp .AND. &
1448	mask_loop(mid,dim,2) < 0.0_wp .AND. &
1449	mask_loop(mid,dim,3) == 0.0_wp ) THEN
1450	WRITE ( io, 350 ) dir(dim), dir(dim)
1451	ELSEIF ( mask_loop(mid,dim,3) == 0.0_wp ) THEN
1452	WRITE ( io, 351 ) dir(dim), dir(dim), mid, dir(dim), &
1453	mask_loop(mid,dim,1:2)
1454	ELSE
1455	WRITE ( io, 351 ) dir(dim), dir(dim), mid, dir(dim), &
1456	mask_loop(mid,dim,1:3)
1457	ENDIF
1458	ENDDO
1459	ENDIF
1460
1461	ENDDO
1462	ENDDO
1463
1464	!
1465	!-- Timeseries
1466	IF ( dt_dots /= 9999999.9_wp ) THEN
1467	WRITE ( io, 340 )
1468
1469	output_format = netcdf_data_format_string
1470	IF ( netcdf_deflate == 0 ) THEN
1471	WRITE ( io, 344 ) output_format
1472	ELSE
1473	WRITE ( io, 354 ) TRIM( output_format ), netcdf_deflate
1474	ENDIF
1475	WRITE ( io, 341 ) dt_dots
1476	ENDIF
1477
1478	WRITE ( io, 99 )
1479
1480	!
1481	!-- Physical quantities
1482	WRITE ( io, 400 )
1483
1484	!
1485	!-- Geostrophic parameters
1486	WRITE ( io, 410 ) latitude, longitude, rotation_angle, omega, f, fs
1487
1488	!
1489	!-- Day and time during model start
1490	CALL get_date_time( 0.0_wp, date_time_str=date_time_str )
1491	WRITE ( io, 456 ) TRIM( date_time_str )
1492
1493	!
1494	!-- Other quantities
1495	WRITE ( io, 411 ) g
1496
1497	WRITE ( io, 412 ) TRIM( reference_state )
1498	IF ( use_single_reference_value ) THEN
1499	IF ( ocean_mode ) THEN
1500	WRITE ( io, 413 ) prho_reference
1501	ELSE
1502	WRITE ( io, 414 ) pt_reference
1503	ENDIF
1504	ENDIF
1505
1506	!
1507	!-- Cloud physcis parameters / quantities / numerical methods
1508	WRITE ( io, 430 )
1509	IF ( humidity .AND. .NOT. bulk_cloud_model .AND. .NOT. cloud_droplets) THEN
1510	WRITE ( io, 431 )
1511	ENDIF
1512	!
1513	!-- LES / turbulence parameters
1514	WRITE ( io, 450 )
1515
1516	!--
1517	! ... LES-constants used must still be added here
1518	!--
1519	IF ( constant_diffusion ) THEN
1520	WRITE ( io, 451 ) km_constant, km_constant/prandtl_number, &
1521	prandtl_number
1522	ENDIF
1523	IF ( .NOT. constant_diffusion) THEN
1524	IF ( e_init > 0.0_wp ) WRITE ( io, 455 ) e_init
1525	IF ( e_min > 0.0_wp ) WRITE ( io, 454 ) e_min
1526	IF ( wall_adjustment ) WRITE ( io, 453 ) wall_adjustment_factor
1527	ENDIF
1528	IF ( rans_mode ) THEN
1529	WRITE ( io, 457 ) rans_const_c, rans_const_sigma
1530	ENDIF
1531	!
1532	!-- Special actions during the run
1533	WRITE ( io, 470 )
1534	IF ( create_disturbances ) THEN
1535	WRITE ( io, 471 ) dt_disturb, disturbance_amplitude, &
1536	zu(disturbance_level_ind_b), disturbance_level_ind_b,&
1537	zu(disturbance_level_ind_t), disturbance_level_ind_t
1538	IF ( .NOT. bc_lr_cyc .OR. .NOT. bc_ns_cyc ) THEN
1539	WRITE ( io, 472 ) inflow_disturbance_begin, inflow_disturbance_end
1540	ELSE
1541	WRITE ( io, 473 ) disturbance_energy_limit
1542	ENDIF
1543	WRITE ( io, 474 ) TRIM( random_generator )
1544	ENDIF
1545	IF ( pt_surface_initial_change /= 0.0_wp ) THEN
1546	WRITE ( io, 475 ) pt_surface_initial_change
1547	ENDIF
1548	IF ( humidity .AND. q_surface_initial_change /= 0.0_wp ) THEN
1549	WRITE ( io, 476 ) q_surface_initial_change
1550	ENDIF
1551	IF ( passive_scalar .AND. q_surface_initial_change /= 0.0_wp ) THEN
1552	WRITE ( io, 477 ) q_surface_initial_change
1553	ENDIF
1554
1555	!
1556	!-- Parameters of 1D-model
1557	IF ( INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 ) THEN
1558	WRITE ( io, 500 ) end_time_1d, dt_run_control_1d, dt_pr_1d, &
1559	mixing_length_1d, dissipation_1d
1560	IF ( damp_level_ind_1d /= nzt+1 ) THEN
1561	WRITE ( io, 502 ) zu(damp_level_ind_1d), damp_level_ind_1d
1562	ENDIF
1563	ENDIF
1564
1565	!
1566	!-- Header information from other modules
1567	CALL module_interface_header( io )
1568
1569
1570	WRITE ( io, 99 )
1571
1572	!
1573	!-- Write buffer contents to disc immediately
1574	FLUSH( io )
1575
1576	!
1577	!-- Here the FORMATs start
1578
1579	99 FORMAT (1X,78('-'))
1580	100 FORMAT (/1X,'******************************',4X,44('-')/ &
1581	1X,'* ',A,' *',4X,A/ &
1582	1X,'******************************',4X,44('-'))
1583	101 FORMAT (35X,'coupled run: ',A/ &
1584	35X,42('-'))
1585	102 FORMAT (/' Date: ',A10,4X,'Run: ',A34/ &
1586	' Time: ',A8,4X,'Run-No.: ',I2.2/ &
1587	' Run on host: ',A10)
1588	#if defined( __parallel )
1589	103 FORMAT (' Number of PEs:',10X,I6,4X,'Processor grid (x,y): (',I4,',',I4, &
1590	')',1X,A)
1591	104 FORMAT (' Number of PEs:',10X,I6,4X,'Tasks:',I4,' threads per task:',I4/ &
1592	35X,'Processor grid (x,y): (',I4,',',I4,')',1X,A)
1593	107 FORMAT (35X,'A 1d-decomposition along ',A,' is used')
1594	108 FORMAT (35X,'Max. # of parallel I/O streams is ',I5)
1595	109 FORMAT (35X,'Precursor run for coupled atmos-ocean run'/ &
1596	35X,42('-'))
1597	114 FORMAT (35X,'Coupled atmosphere-ocean run following'/ &
1598	35X,'independent precursor runs'/ &
1599	35X,42('-'))
1600	#endif
1601	110 FORMAT (/' Numerical Schemes:'/ &
1602	' -----------------'/)
1603	124 FORMAT (' --> Use the ',A,' turbulence closure (',A,' mode).')
1604	121 FORMAT (' --> Use the ',A,' approximation for the model equations.')
1605	111 FORMAT (' --> Solve perturbation pressure via FFT using ',A,' routines')
1606	112 FORMAT (' --> Solve perturbation pressure via SOR-Red/Black-Schema'/ &
1607	' Iterations (initial/other): ',I3,'/',I3,' omega =',F6.3)
1608	113 FORMAT (' --> Momentum advection via Piascek-Williams-Scheme (Form C3)', &
1609	' or Upstream')
1610	115 FORMAT (' FFT and transpositions are overlapping')
1611	116 FORMAT (' --> Scalar advection via Piascek-Williams-Scheme (Form C3)', &
1612	' or Upstream')
1613	118 FORMAT (' --> Scalar advection via Bott-Chlond-Scheme')
1614	119 FORMAT (' --> Galilei-Transform applied to horizontal advection:'/ &
1615	' translation velocity = ',A/ &
1616	' distance advected ',A,': ',F8.3,' km(x) ',F8.3,' km(y)')
1617	122 FORMAT (' --> Time differencing scheme: ',A)
1618	123 FORMAT (' --> Rayleigh-Damping active, starts ',A,' z = ',F8.2,' m'/ &
1619	' maximum damping coefficient:',F6.3, ' 1/s')
1620	129 FORMAT (' --> Additional prognostic equation for the specific humidity')
1621	130 FORMAT (' --> Additional prognostic equation for the total water content')
1622	131 FORMAT (' --> No pt-equation solved. Neutral stratification with pt = ', &
1623	F6.2, ' K assumed')
1624	134 FORMAT (' --> Additional prognostic equation for a passive scalar')
1625	135 FORMAT (' --> Solve perturbation pressure via ',A,' method (', &
1626	A,'-cycle)'/ &
1627	' number of grid levels: ',I2/ &
1628	' Gauss-Seidel red/black iterations: ',I2)
1629	136 FORMAT (' gridpoints of coarsest subdomain (x,y,z): (',I3,',',I3,',', &
1630	I3,')')
1631	137 FORMAT (' level data gathered on PE0 at level: ',I2/ &
1632	' gridpoints of coarsest subdomain (x,y,z): (',I3,',',I3,',', &
1633	I3,')'/ &
1634	' gridpoints of coarsest domain (x,y,z): (',I3,',',I3,',', &
1635	I3,')')
1636	139 FORMAT (' --> Loop optimization method: ',A)
1637	140 FORMAT (' maximum residual allowed: ',E10.3)
1638	141 FORMAT (' fixed number of multigrid cycles: ',I4)
1639	142 FORMAT (' perturbation pressure is calculated at every Runge-Kutta ', &
1640	'step')
1641	143 FORMAT (' Euler/upstream scheme is used for the SGS turbulent ', &
1642	'kinetic energy')
1643	144 FORMAT (' masking method is used')
1644	150 FORMAT (' --> Volume flow at the right and north boundary will be ', &
1645	'conserved'/ &
1646	' using the ',A,' mode')
1647	151 FORMAT (' with u_bulk = ',F7.3,' m/s and v_bulk = ',F7.3,' m/s')
1648	152 FORMAT (' --> External pressure gradient directly prescribed by the user:',&
1649	/' ',2(1X,E12.5),'Pa/m in x/y direction', &
1650	/' starting from dp_level_b =', F8.3, 'm', A /)
1651	200 FORMAT (//' Run time and time step information:'/ &
1652	' ----------------------------------'/)
1653	201 FORMAT ( ' Timestep: variable maximum value: ',F6.3,' s', &
1654	' CFL-factor:',F5.2)
1655	202 FORMAT ( ' Timestep: dt = ',F6.3,' s'/)
1656	203 FORMAT ( ' Start time: ',F11.3,' s'/ &
1657	' End time: ',F11.3,' s')
1658	204 FORMAT ( A,F11.3,' s')
1659	205 FORMAT ( A,F11.3,' s',5X,'restart every',17X,F11.3,' s')
1660	206 FORMAT (/' Time reached: ',F11.3,' s'/ &
1661	' CPU-time used: ',F9.3,' s per timestep: ',F9.3,' s'/ &
1662	' per second of simulated time: ',F9.3,' s')
1663	207 FORMAT ( ' Spinup time: ',F11.3,' s')
1664	250 FORMAT (//' Computational grid and domain size:'/ &
1665	' ----------------------------------'// &
1666	' Grid length: dx = ',F8.3,' m dy = ',F8.3, ' m')
1667	251 FORMAT ( /' Domain size: x = ',F10.3,' m y = ',F10.3, &
1668	' m z(u) = ',F10.3,' m'/)
1669	253 FORMAT ( ' dz(',I1,') = ', F8.3, ' m')
1670	254 FORMAT (//' Number of gridpoints (x,y,z): (0:',I4,', 0:',I4,', 0:',I4,')'/ &
1671	' Subdomain size (x,y,z): ( ',I4,', ',I4,', ',I4,')'/)
1672	260 FORMAT (/' The model has a slope in x-direction. Inclination angle: ',F6.2,&
1673	' degrees')
1674	270 FORMAT (//' Topography information:'/ &
1675	' ----------------------'// &
1676	1X,'Topography: ',A)
1677	271 FORMAT ( ' Building size (x/y/z) in m: ',F5.1,' / ',F5.1,' / ',F5.1/ &
1678	' Horizontal index bounds (l/r/s/n): ',I4,' / ',I4,' / ',I4, &
1679	' / ',I4)
1680	272 FORMAT ( ' Single quasi-2D street canyon of infinite length in ',A, &
1681	' direction' / &
1682	' Canyon height: ', F6.2, 'm, ch = ', I4, '.' / &
1683	' Canyon position (',A,'-walls): cxl = ', I4,', cxr = ', I4, '.')
1684	273 FORMAT ( ' Tunnel of infinite length in ',A, &
1685	' direction' / &
1686	' Tunnel height: ', F6.2, / &
1687	' Tunnel-wall depth: ', F6.2 / &
1688	' Tunnel width: ', F6.2 )
1689	274 FORMAT ( ' Tunnel in ', A, ' direction.' / &
1690	' Tunnel height: ', F6.2, / &
1691	' Tunnel-wall depth: ', F6.2 / &
1692	' Tunnel width: ', F6.2, / &
1693	' Tunnel length: ', F6.2 )
1694	278 FORMAT (' Topography grid definition convention:'/ &
1695	' cell edge (staggered grid points'/ &
1696	' (u in x-direction, v in y-direction))' /)
1697	279 FORMAT (' Topography grid definition convention:'/ &
1698	' cell center (scalar grid points)' /)
1699	280 FORMAT (' Complex terrain simulation is activated.')
1700	281 FORMAT (' --> Mean inflow profiles are adjusted.' / &
1701	' --> Elevation of inflow boundary: ', F7.1, ' m' )
1702	282 FORMAT (' --> Initial data from 3D-precursor run is shifted' / &
1703	' vertically depending on local surface height.')
1704	300 FORMAT (//' Boundary conditions:'/ &
1705	' -------------------'// &
1706	' p uv ', &
1707	' pt'// &
1708	' B. bound.: ',A/ &
1709	' T. bound.: ',A)
1710	301 FORMAT (/' ',A// &
1711	' B. bound.: ',A/ &
1712	' T. bound.: ',A)
1713	303 FORMAT (/' Bottom surface fluxes are used in diffusion terms at k=1')
1714	304 FORMAT (/' Top surface fluxes are used in diffusion terms at k=nzt')
1715	305 FORMAT (//' Constant flux layer between bottom surface and first ', &
1716	'computational u,v-level:'// &
1717	' z_mo = ',F6.2,' m z0 =',F7.4,' m z0h =',F8.5,&
1718	' m kappa =',F5.2/ &
1719	' Rif value range: ',F8.2,' <= rif <=',F6.2)
1720	306 FORMAT (' Predefined constant heatflux: ',F9.6,' K m/s')
1721	307 FORMAT (' Heatflux has a random normal distribution')
1722	308 FORMAT (' Predefined surface temperature')
1723	309 FORMAT (' Predefined constant salinityflux: ',F9.6,' psu m/s')
1724	310 FORMAT (//' 1D-Model:'// &
1725	' Rif value range: ',F6.2,' <= rif <=',F6.2)
1726	311 FORMAT (' Predefined constant humidity flux: ',E10.3,' kg/kg m/s')
1727	312 FORMAT (' Predefined surface humidity')
1728	313 FORMAT (' Predefined constant scalar flux: ',E10.3,' kg/(m**2 s)')
1729	314 FORMAT (' Predefined scalar value at the surface')
1730	302 FORMAT (' Predefined constant scalarflux: ',F9.6,' kg/(m**2 s)')
1731	315 FORMAT (' Humidity flux at top surface is 0.0')
1732	316 FORMAT (' Sensible heatflux and momentum flux from coupled ', &
1733	'atmosphere model')
1734	317 FORMAT (//' Lateral boundaries:'/ &
1735	' left/right: ',A/ &
1736	' north/south: ',A)
1737	318 FORMAT (/' use_cmax: ',L1 / &
1738	' pt damping layer width = ',F8.2,' m, pt ', &
1739	'damping factor =',F7.4)
1740	319 FORMAT (' turbulence recycling at inflow switched on'/ &
1741	' width of recycling domain: ',F7.1,' m grid index: ',I4/ &
1742	' inflow damping height: ',F6.1,' m width: ',F6.1,' m')
1743	320 FORMAT (' Predefined constant momentumflux: u: ',F9.6,' m2/s2'/ &
1744	' v: ',F9.6,' m2/s2')
1745	321 FORMAT (//' Initial profiles:'/ &
1746	' ----------------')
1747	322 FORMAT (' turbulence recycling at inflow switched on'/ &
1748	' y-shift of the recycled inflow turbulence is',I3,' PE'/ &
1749	' width of recycling domain: ',F7.1,' m grid index: ',I4/ &
1750	' inflow damping height: ',F6.1,' m width: ',F6.1,' m'/)
1751	323 FORMAT (' turbulent outflow conditon switched on'/ &
1752	' position of outflow source plane: ',F7.1,' m ', &
1753	'grid index: ', I4)
1754	325 FORMAT (//' List output:'/ &
1755	' -----------'// &
1756	' 1D-Profiles:'/ &
1757	' Output every ',F10.2,' s')
1758	326 FORMAT (' Time averaged over ',F8.2,' s'/ &
1759	' Averaging input every ',F8.2,' s')
1760	330 FORMAT (//' Data output:'/ &
1761	' -----------'/)
1762	331 FORMAT (/' 1D-Profiles:')
1763	332 FORMAT (/' ',A)
1764	333 FORMAT (' Output every ',F8.2,' s',/ &
1765	' Time averaged over ',F8.2,' s'/ &
1766	' Averaging input every ',F8.2,' s')
1767	334 FORMAT (/' 2D-Arrays',A,':')
1768	335 FORMAT (/' ',A2,'-cross-section Arrays: ',A/ &
1769	' Output every ',F8.2,' s ',A/ &
1770	' Cross sections at ',A1,' = ',A/ &
1771	' scalar-coordinates: ',A,' m'/)
1772	336 FORMAT (/' 3D-Arrays',A,':')
1773	337 FORMAT (/' Arrays: ',A/ &
1774	' Output every ',F8.2,' s ',A/ &
1775	' Upper output limit at ',F8.2,' m (GP ',I4,')'/)
1776	339 FORMAT (' No output during initial ',F8.2,' s')
1777	340 FORMAT (/' Time series:')
1778	341 FORMAT (' Output every ',F8.2,' s'/)
1779	342 FORMAT (/' ',A2,'-cross-section Arrays: ',A/ &
1780	' Output every ',F8.2,' s ',A/ &
1781	' Time averaged over ',F8.2,' s'/ &
1782	' Averaging input every ',F8.2,' s'/ &
1783	' Cross sections at ',A1,' = ',A/ &
1784	' scalar-coordinates: ',A,' m'/)
1785	343 FORMAT (/' Arrays: ',A/ &
1786	' Output every ',F8.2,' s ',A/ &
1787	' Time averaged over ',F8.2,' s'/ &
1788	' Averaging input every ',F8.2,' s'/ &
1789	' Upper output limit at ',F8.2,' m (GP ',I4,')'/)
1790	344 FORMAT (' Output format: ',A/)
1791	345 FORMAT (/' Scaling lengths for output locations of all subsequent mask IDs:',/ &
1792	' mask_scale_x (in x-direction): ',F9.3, ' m',/ &
1793	' mask_scale_y (in y-direction): ',F9.3, ' m',/ &
1794	' mask_scale_z (in z-direction): ',F9.3, ' m' )
1795	346 FORMAT (/' Masked data output',A,' for mask ID ',I2, ':')
1796	347 FORMAT (' Variables: ',A/ &
1797	' Output every ',F8.2,' s')
1798	348 FORMAT (' Variables: ',A/ &
1799	' Output every ',F8.2,' s'/ &
1800	' Time averaged over ',F8.2,' s'/ &
1801	' Averaging input every ',F8.2,' s')
1802	349 FORMAT (/' Output locations in ',A,'-direction in multiples of ', &
1803	'mask_scale_',A,' predefined by array mask_',I2.2,'_',A,':'/ &
1804	13(' ',8(F8.2,',')/) )
1805	350 FORMAT (/' Output locations in ',A,'-direction: ', &
1806	'all gridpoints along ',A,'-direction (default).' )
1807	351 FORMAT (/' Output locations in ',A,'-direction in multiples of ', &
1808	'mask_scale_',A,' constructed from array mask_',I2.2,'_',A,'_loop:'/ &
1809	' loop begin:',F8.2,', end:',F8.2,', stride:',F8.2 )
1810	352 FORMAT (/' Number of output time levels allowed: ',I3 /)
1811	353 FORMAT (/' Number of output time levels allowed: unlimited' /)
1812	354 FORMAT (' Output format: ',A, ' compressed with level: ',I1/)
1813	400 FORMAT (//' Physical quantities:'/ &
1814	' -------------------'/)
1815	410 FORMAT (' Geograph. latitude : latitude = ',F5.1,' degr'/ &
1816	' Geograph. longitude : longitude = ',F5.1,' degr'/ &
1817	' Rotation angle : rotation_angle = ',F5.1,' degr'/ &
1818	' Angular velocity : omega =',E10.3,' rad/s'/ &
1819	' Coriolis parameter : f = ',F9.6,' 1/s'/ &
1820	' f* = ',F9.6,' 1/s')
1821	411 FORMAT (/' Gravity : g = ',F4.1,' m/s**2')
1822	412 FORMAT (/' Reference state used in buoyancy terms: ',A)
1823	413 FORMAT (' Reference density in buoyancy terms: ',F8.3,' kg/m**3')
1824	414 FORMAT (' Reference temperature in buoyancy terms: ',F8.4,' K')
1825	420 FORMAT (/' Characteristic levels of the initial temperature profile:'// &
1826	' Height: ',A,' m'/ &
1827	' Temperature: ',A,' K'/ &
1828	' Gradient: ',A,' K/100m'/ &
1829	' Gridpoint: ',A)
1830	421 FORMAT (/' Characteristic levels of the initial humidity profile:'// &
1831	' Height: ',A,' m'/ &
1832	' Humidity: ',A,' kg/kg'/ &
1833	' Gradient: ',A,' (kg/kg)/100m'/ &
1834	' Gridpoint: ',A)
1835	422 FORMAT (/' Characteristic levels of the initial scalar profile:'// &
1836	' Height: ',A,' m'/ &
1837	' Scalar concentration: ',A,' kg/m**3'/ &
1838	' Gradient: ',A,' (kg/m**3)/100m'/ &
1839	' Gridpoint: ',A)
1840	423 FORMAT (/' Characteristic levels of the geo. wind component ug:'// &
1841	' Height: ',A,' m'/ &
1842	' ug: ',A,' m/s'/ &
1843	' Gradient: ',A,' 1/100s'/ &
1844	' Gridpoint: ',A)
1845	424 FORMAT (/' Characteristic levels of the geo. wind component vg:'// &
1846	' Height: ',A,' m'/ &
1847	' vg: ',A,' m/s'/ &
1848	' Gradient: ',A,' 1/100s'/ &
1849	' Gridpoint: ',A)
1850	425 FORMAT (/' Characteristic levels of the initial salinity profile:'// &
1851	' Height: ',A,' m'/ &
1852	' Salinity: ',A,' psu'/ &
1853	' Gradient: ',A,' psu/100m'/ &
1854	' Gridpoint: ',A)
1855	426 FORMAT (/' Characteristic levels of the subsidence/ascent profile:'// &
1856	' Height: ',A,' m'/ &
1857	' w_subs: ',A,' m/s'/ &
1858	' Gradient: ',A,' (m/s)/100m'/ &
1859	' Gridpoint: ',A)
1860	427 FORMAT (/' Initial wind profiles (u,v) are interpolated from given'// &
1861	' profiles')
1862	428 FORMAT (/' Initial profiles (u, v, pt, q) are taken from file '/ &
1863	' NUDGING_DATA')
1864	430 FORMAT (//' Cloud physics quantities / methods:'/ &
1865	' ----------------------------------'/)
1866	431 FORMAT (' Humidity is considered, bu no condensation')
1867	450 FORMAT (//' LES / Turbulence quantities:'/ &
1868	' ---------------------------'/)
1869	451 FORMAT (' Diffusion coefficients are constant:'/ &
1870	' Km = ',F6.2,' m2/s Kh = ',F6.2,' m2/s Pr = ',F5.2)
1871	453 FORMAT (' Mixing length is limited to',F5.2,' * z')
1872	454 FORMAT (' TKE is not allowed to fall below ',E9.2,' (m/s)**2')
1873	455 FORMAT (' initial TKE is prescribed as ',E9.2,' (m/s)**2')
1874	456 FORMAT (/' Date and time at model start : ',A)
1875	457 FORMAT (' RANS-mode constants: c_0 = ',F9.5/ &
1876	' c_1 = ',F9.5/ &
1877	' c_2 = ',F9.5/ &
1878	' c_3 = ',F9.5/ &
1879	' c_4 = ',F9.5/ &
1880	' sigma_e = ',F9.5/ &
1881	' sigma_diss = ',F9.5)
1882	470 FORMAT (//' Actions during the simulation:'/ &
1883	' -----------------------------'/)
1884	471 FORMAT (' Disturbance impulse (u,v) every : ',F6.2,' s'/ &
1885	' Disturbance amplitude : ',F5.2, ' m/s'/ &
1886	' Lower disturbance level : ',F8.2,' m (GP ',I4,')'/ &
1887	' Upper disturbance level : ',F8.2,' m (GP ',I4,')')
1888	472 FORMAT (' Disturbances continued during the run from i/j =',I4, &
1889	' to i/j =',I4)
1890	473 FORMAT (' Disturbances cease as soon as the disturbance energy exceeds',&
1891	F6.3, ' m2/s2')
1892	474 FORMAT (' Random number generator used : ',A/)
1893	475 FORMAT (' The surface temperature is increased (or decreased, ', &
1894	'respectively, if'/ &
1895	' the value is negative) by ',F5.2,' K at the beginning of the',&
1896	' 3D-simulation'/)
1897	476 FORMAT (' The surface humidity is increased (or decreased, ',&
1898	'respectively, if the'/ &
1899	' value is negative) by ',E8.1,' kg/kg at the beginning of', &
1900	' the 3D-simulation'/)
1901	477 FORMAT (' The scalar value is increased at the surface (or decreased, ',&
1902	'respectively, if the'/ &
1903	' value is negative) by ',E8.1,' kg/m**3 at the beginning of', &
1904	' the 3D-simulation'/)
1905	500 FORMAT (//' 1D-Model parameters:'/ &
1906	' -------------------'// &
1907	' Simulation time: ',F8.1,' s'/ &
1908	' Run-controll output every: ',F8.1,' s'/ &
1909	' Vertical profile output every: ',F8.1,' s'/ &
1910	' Mixing length calculation: ',A/ &
1911	' Dissipation calculation: ',A/)
1912	502 FORMAT (' Damping layer starts from ',F7.1,' m (GP ',I4,')'/)
1913	503 FORMAT (' --> Momentum advection via Wicker-Skamarock-Scheme 5th order')
1914	504 FORMAT (' --> Scalar advection via Wicker-Skamarock-Scheme 5th order')
1915	512 FORMAT (/' Date: ',A10,6X,'Run: ',A34/ &
1916	' Time: ',A8,6X,'Run-No.: ',I2.2/ &
1917	' Run on host: ',A10,6X,'En-No.: ',I2.2)
1918	600 FORMAT (/' Nesting informations:'/ &
1919	' --------------------'/ &
1920	' Nesting mode: ',A/ &
1921	' Nesting-datatransfer mode: ',A// &
1922	' Nest id parent number lower left coordinates name'/ &
1923	' (*=me) id of PEs x (m) y (m)' )
1924	601 FORMAT (2X,A1,1X,I2.2,6X,I2.2,5X,I5,5X,F8.2,2X,F8.2,5X,A)
1925
1926	END SUBROUTINE header

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |