Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |