Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |