Home

Context Navigation

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

Last change on this file since 1661 was 1660, checked in by gronemeier, 10 years ago
Bugfix: Definition of topography height in case of grid stretching
Property svn:keywords set to `Id`
File size: 91.9 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |