Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |