Home

Context Navigation

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

Last change on this file since 1355 was 1354, checked in by heinze, 11 years ago
last commit documented
Property svn:keywords set to `Id`
File size: 80.9 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |