Home

Context Navigation

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

Last change on this file since 1717 was 1698, checked in by raasch, 9 years ago
last commit documented example_cbl_rc updated
Property svn:keywords set to `Id`
File size: 92.5 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |