Home

Context Navigation

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

Last change on this file since 481 was 449, checked in by raasch, 14 years ago
branch revision comments from Marcus (rev 410) replaced by normal revision comments
Property svn:keywords set to `Id`
File size: 75.7 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |