Home

Context Navigation

check_parameters.f90 @ 707

Last change on this file since 707 was 707, checked in by raasch, 13 years ago

New:
---

In case of multigrid method, on coarse grid levels, gathered data are
identically processed on all PEs (before, on PE0 only), so that the subsequent
scattering of data is not neccessary any more. (modules, init_pegrid, poismg)

Changed:

Calculation of weighted average of p is now handled in the same way
regardless of the number of ghost layers (advection scheme). (pres)

multigrid and sor method are using p_loc for iterative
advancements of pressure. p_sub removed. (init_3d_model, modules, poismg, pres, sor)

bc_lr and bc_ns replaced by bc_lr_dirrad, bc_lr_raddir, bc_ns_dirrad, bc_ns_raddir
for speed optimization. (calc_spectra, check_parameters, exchange_horiz,
exchange_horiz_2d, header, init_3d_model, init_grid, init_pegrid, modules,
poismg, pres, sor, time_integration, timestep)

grid_level directly used as index for MPI data type arrays. (exchange_horiz,
poismg)

initial assignments of zero to array p for iterative solvers only (init_3d_model)

Errors:

localsum calculation modified for proper OpenMP reduction. (pres)

Bugfix: bottom (nzb) and top (nzt+1) boundary conditions set in routines
resid and restrict. They were missed before, which may have led to
unpredictable results. (poismg)

Property svn:keywords set to Id

File size: 127.7 KB

Rev	Line
[1]	1	SUBROUTINE check_parameters
	2
	3	!------------------------------------------------------------------------------!
[484]	4	! Current revisions:
[1]	5	! -----------------
[707]	6	! setting of bc_lr/ns_dirrad/raddir
[667]	7	!
[668]	8	! Former revisions:
	9	! -----------------
[687]	10	!
[690]	11	! 689 2011-02-20 19:31:12z gryschka
	12	! Bugfix for some logical expressions
	13	! (syntax was not compatible with all compilers)
[687]	14	!
[690]	15	! 680 2011-02-04 23:16:06Z gryschka
[681]	16	! init_vortex is not allowed with volume_flow_control
[668]	17	!
[674]	18	! 673 2011-01-18 16:19:48Z suehring
	19	! Declaration of ws_scheme_sca and ws_scheme_mom added (moved from advec_ws).
	20	!
[668]	21	! 667 2010-12-23 12:06:00Z suehring/gryschka
[667]	22	! Exchange of parameters between ocean and atmosphere via PE0
	23	! Check for illegal combination of ws-scheme and timestep scheme.
	24	! Check for topography and ws-scheme.
	25	! Check for not cyclic boundary conditions in combination with ws-scheme and
	26	! loop_optimization = 'vector'.
	27	! Check for call_psolver_at_all_substeps and ws-scheme for momentum_advec.
	28	! Different processor/grid topology in atmosphere and ocean is now allowed!
	29	! Bugfixes in checking for conserve_volume_flow_mode
[601]	30	! 600 2010-11-24 16:10:51Z raasch
	31	! change due to new default value of surface_waterflux
[581]	32	! 580 2010-10-05 13:59:11Z heinze
[601]	33	! renaming of ws_vertical_gradient_level to subs_vertical_gradient_level
[581]	34	!
[568]	35	! 567 2010-10-01 10:46:30Z helmke
	36	! calculating masks changed
	37	!
[565]	38	! 564 2010-09-30 13:18:59Z helmke
	39	! palm message identifiers of masked output changed, 20 replaced by max_masks
	40	!
[554]	41	! 553 2010-09-01 14:09:06Z weinreis
	42	! masks is calculated and removed from inipar
[532]	43	!
	44	! 531 2010-04-21 06:47:21Z heinze
	45	! Bugfix: unit of hyp changed to dbar
[392]	46	!
[525]	47	! 524 2010-03-30 02:04:51Z raasch
	48	! Bugfix: "/" in netcdf profile variable names replaced by ":"
	49	!
[494]	50	! 493 2010-03-01 08:30:24Z raasch
	51	! netcdf_data_format is checked
	52	!
[449]	53	! 411 2009-12-11 14:15:58Z heinze
	54	! Enabled passive scalar/humidity wall fluxes for non-flat topography
	55	! Initialization of large scale vertical motion (subsidence/ascent)
	56	!
	57	! 410 2009-12-04 17:05:40Z letzel
	58	! masked data output
	59	!
[392]	60	! 388 2009-09-23 09:40:33Z raasch
[388]	61	! Check profiles fpr prho and hyp.
[376]	62	! Bugfix: output of averaged 2d/3d quantities requires that an avaraging
	63	! interval has been set, respective error message is included
[366]	64	! bc_lr_cyc and bc_ns_cyc are set,
[328]	65	! initializing_actions='read_data_for_recycling' renamed to 'cyclic_fill'
[308]	66	! Check for illegal entries in section_xy\|xz\|yz that exceed nz+1\|ny+1\|nx+1
[291]	67	! Coupling with independent precursor runs.
[264]	68	! Check particle_color, particle_dvrpsize, color_interval, dvrpsize_interval
[232]	69	! Bugfix: pressure included for profile output
[240]	70	! Check pressure gradient conditions
[256]	71	! topography_grid_convention moved from user_check_parameters
	72	! 'single_street_canyon'
[367]	73	! Added shf* and qsws* to the list of available output data
[226]	74	!
	75	! 222 2009-01-12 16:04:16Z letzel
[217]	76	! +user_check_parameters
[213]	77	! Output of messages replaced by message handling routine.
[206]	78	! Implementation of an MPI-1 coupling: replaced myid with target_id,
	79	! deleted __mpi2 directives
[222]	80	! Check that PALM is called with mrun -K parallel for coupling
[198]	81	!
	82	! 197 2008-09-16 15:29:03Z raasch
[177]	83	! Bug fix: Construction of vertical profiles when 10 gradients have been
[181]	84	! specified in the parameter list (ug, vg, pt, q, sa, lad)
	85	!
	86	! Strict grid matching along z is not needed for mg-solver.
[153]	87	! Leaf area density (LAD) explicitly set to its surface value at k=0
[151]	88	! Case of reading data for recycling included in initializing_actions,
[197]	89	! check of turbulent_inflow and calculation of recycling_plane.
	90	! q*2 profile added
[139]	91	!
	92	! 138 2007-11-28 10:03:58Z letzel
[138]	93	! Plant canopy added
[132]	94	! Allow new case bc_uv_t = 'dirichlet_0' for channel flow.
[116]	95	! Multigrid solver allows topography, checking of dt_sort_particles
[121]	96	! Bugfix: initializing u_init and v_init in case of ocean runs
[110]	97	!
	98	! 109 2007-08-28 15:26:47Z letzel
[102]	99	! Check coupling_mode and set default (obligatory) values (like boundary
[106]	100	! conditions for temperature and fluxes) in case of coupled runs.
	101	! +profiles for wp and w"e
[104]	102	! Bugfix: Error message concerning output of particle concentration (pc)
	103	! modified
[108]	104	! More checks and more default values for coupled runs
	105	! allow data_output_pr= q, wq, w"q", wq for humidity = .T. (instead of
	106	! cloud_physics = .T.)
	107	! Rayleigh damping for ocean fixed.
[109]	108	! Check and, if necessary, set default value for dt_coupling
[98]	109	!
	110	! 97 2007-06-21 08:23:15Z raasch
[97]	111	! Initial salinity profile is calculated, salinity boundary conditions are
	112	! checked,
[94]	113	! z_max_do1d is checked only in case of ocean = .f.,
[97]	114	! +initial temperature and geostrophic velocity profiles for the ocean version,
	115	! use_pt_reference renamed use_reference
[77]	116	!
[90]	117	! 89 2007-05-25 12:08:31Z raasch
	118	! Check for user-defined profiles
	119	!
[77]	120	! 75 2007-03-22 09:54:05Z raasch
[51]	121	! "by_user" allowed as initializing action, -data_output_ts,
[63]	122	! leapfrog with non-flat topography not allowed any more, loop_optimization
[75]	123	! and pt_reference are checked, moisture renamed humidity,
[72]	124	! output of precipitation amount/rate and roughnes length + check
[73]	125	! possible negative humidities are avoided in initial profile,
[75]	126	! dirichlet/neumann changed to dirichlet/radiation, etc.,
	127	! revision added to run_description_header
[1]	128	!
[39]	129	! 20 2007-02-26 00:12:32Z raasch
	130	! Temperature and humidity gradients at top are now calculated for nzt+1,
	131	! top_heatflux and respective boundary condition bc_pt_t is checked
	132	!
[3]	133	! RCS Log replace by Id keyword, revision history cleaned up
	134	!
[1]	135	! Revision 1.61 2006/08/04 14:20:25 raasch
	136	! do2d_unit and do3d_unit now defined as 2d-arrays, check of
	137	! use_upstream_for_tke, default value for dt_dopts,
	138	! generation of file header moved from routines palm and header to here
	139	!
	140	! Revision 1.1 1997/08/26 06:29:23 raasch
	141	! Initial revision
	142	!
	143	!
	144	! Description:
	145	! ------------
	146	! Check control parameters and deduce further quantities.
	147	!------------------------------------------------------------------------------!
	148
	149	USE arrays_3d
	150	USE constants
	151	USE control_parameters
[264]	152	USE dvrp_variables
[1]	153	USE grid_variables
	154	USE indices
	155	USE model_1d
	156	USE netcdf_control
	157	USE particle_attributes
	158	USE pegrid
	159	USE profil_parameter
[411]	160	USE subsidence_mod
[1]	161	USE statistics
	162	USE transpose_indices
	163
	164	IMPLICIT NONE
	165
	166	CHARACTER (LEN=1) :: sq
	167	CHARACTER (LEN=6) :: var
	168	CHARACTER (LEN=7) :: unit
	169	CHARACTER (LEN=8) :: date
	170	CHARACTER (LEN=10) :: time
[102]	171	CHARACTER (LEN=40) :: coupling_string
[1]	172	CHARACTER (LEN=100) :: action
	173
[108]	174	INTEGER :: i, ilen, intervals, iremote = 0, iter, j, k, nnxh, nnyh, &
	175	position, prec
[1]	176	LOGICAL :: found, ldum
[291]	177	REAL :: gradient, maxn, maxp, remote = 0.0, &
	178	simulation_time_since_reference
[1]	179
	180	!
	181	!-- Warning, if host is not set
	182	IF ( host(1:1) == ' ' ) THEN
[213]	183	message_string = '"host" is not set. Please check that environment ' // &
	184	'variable "localhost" & is set before running PALM'
[226]	185	CALL message( 'check_parameters', 'PA0001', 0, 0, 0, 6, 0 )
[1]	186	ENDIF
	187
	188	!
[102]	189	!-- Check the coupling mode
	190	IF ( coupling_mode /= 'uncoupled' .AND. &
	191	coupling_mode /= 'atmosphere_to_ocean' .AND. &
	192	coupling_mode /= 'ocean_to_atmosphere' ) THEN
[213]	193	message_string = 'illegal coupling mode: ' // TRIM( coupling_mode )
[226]	194	CALL message( 'check_parameters', 'PA0002', 1, 2, 0, 6, 0 )
[102]	195	ENDIF
	196
	197	!
[108]	198	!-- Check dt_coupling, restart_time, dt_restart, end_time, dx, dy, nx and ny
[667]	199	IF ( coupling_mode /= 'uncoupled') THEN
[213]	200
[108]	201	IF ( dt_coupling == 9999999.9 ) THEN
[213]	202	message_string = 'dt_coupling is not set but required for coup' // &
	203	'ling mode "' // TRIM( coupling_mode ) // '"'
[226]	204	CALL message( 'check_parameters', 'PA0003', 1, 2, 0, 6, 0 )
[108]	205	ENDIF
[213]	206
[206]	207	#if defined( __parallel )
[667]	208	IF ( myid == 0 ) THEN
	209	CALL MPI_SEND( dt_coupling, 1, MPI_REAL, target_id, 11, comm_inter, &
	210	ierr )
	211	CALL MPI_RECV( remote, 1, MPI_REAL, target_id, 11, comm_inter, &
	212	status, ierr )
	213	ENDIF
	214	CALL MPI_BCAST( remote, 1, MPI_REAL, 0, comm2d, ierr)
	215
[108]	216	IF ( dt_coupling /= remote ) THEN
[213]	217	WRITE( message_string, * ) 'coupling mode "', TRIM( coupling_mode ), &
	218	'": dt_coupling = ', dt_coupling, '& is not equal to ', &
	219	'dt_coupling_remote = ', remote
[226]	220	CALL message( 'check_parameters', 'PA0004', 1, 2, 0, 6, 0 )
[108]	221	ENDIF
[109]	222	IF ( dt_coupling <= 0.0 ) THEN
[667]	223	IF ( myid == 0 ) THEN
	224	CALL MPI_SEND( dt_max, 1, MPI_REAL, target_id, 19, comm_inter, ierr )
	225	CALL MPI_RECV( remote, 1, MPI_REAL, target_id, 19, comm_inter, &
	226	status, ierr )
	227	ENDIF
	228	CALL MPI_BCAST( remote, 1, MPI_REAL, 0, comm2d, ierr)
	229
[109]	230	dt_coupling = MAX( dt_max, remote )
[213]	231	WRITE( message_string, * ) 'coupling mode "', TRIM( coupling_mode ), &
	232	'": dt_coupling <= 0.0 & is not allowed and is reset to ', &
	233	'MAX(dt_max(A,O)) = ', dt_coupling
[226]	234	CALL message( 'check_parameters', 'PA0005', 0, 1, 0, 6, 0 )
[109]	235	ENDIF
[667]	236	IF ( myid == 0 ) THEN
	237	CALL MPI_SEND( restart_time, 1, MPI_REAL, target_id, 12, comm_inter, &
	238	ierr )
	239	CALL MPI_RECV( remote, 1, MPI_REAL, target_id, 12, comm_inter, &
	240	status, ierr )
	241	ENDIF
	242	CALL MPI_BCAST( remote, 1, MPI_REAL, 0, comm2d, ierr)
	243
[108]	244	IF ( restart_time /= remote ) THEN
[213]	245	WRITE( message_string, * ) 'coupling mode "', TRIM( coupling_mode ), &
	246	'": restart_time = ', restart_time, '& is not equal to ', &
	247	'restart_time_remote = ', remote
[226]	248	CALL message( 'check_parameters', 'PA0006', 1, 2, 0, 6, 0 )
[108]	249	ENDIF
[667]	250	IF ( myid == 0 ) THEN
	251	CALL MPI_SEND( dt_restart, 1, MPI_REAL, target_id, 13, comm_inter, &
	252	ierr )
	253	CALL MPI_RECV( remote, 1, MPI_REAL, target_id, 13, comm_inter, &
	254	status, ierr )
	255	ENDIF
	256	CALL MPI_BCAST( remote, 1, MPI_REAL, 0, comm2d, ierr)
	257
[108]	258	IF ( dt_restart /= remote ) THEN
[213]	259	WRITE( message_string, * ) 'coupling mode "', TRIM( coupling_mode ), &
	260	'": dt_restart = ', dt_restart, '& is not equal to ', &
	261	'dt_restart_remote = ', remote
[226]	262	CALL message( 'check_parameters', 'PA0007', 1, 2, 0, 6, 0 )
[108]	263	ENDIF
[213]	264
[291]	265	simulation_time_since_reference = end_time - coupling_start_time
[667]	266	IF ( myid == 0 ) THEN
	267	CALL MPI_SEND( simulation_time_since_reference, 1, MPI_REAL, target_id, &
	268	14, comm_inter, ierr )
	269	CALL MPI_RECV( remote, 1, MPI_REAL, target_id, 14, comm_inter, &
	270	status, ierr )
	271	ENDIF
	272	CALL MPI_BCAST( remote, 1, MPI_REAL, 0, comm2d, ierr)
	273
[291]	274	IF ( simulation_time_since_reference /= remote ) THEN
[213]	275	WRITE( message_string, * ) 'coupling mode "', TRIM( coupling_mode ), &
[291]	276	'": simulation_time_since_reference = ', &
	277	simulation_time_since_reference, '& is not equal to ', &
	278	'simulation_time_since_reference_remote = ', remote
[226]	279	CALL message( 'check_parameters', 'PA0008', 1, 2, 0, 6, 0 )
[108]	280	ENDIF
[213]	281
[667]	282
	283	IF ( myid == 0 ) THEN
	284	CALL MPI_SEND( dx, 1, MPI_REAL, target_id, 15, comm_inter, ierr )
	285	CALL MPI_RECV( remote, 1, MPI_REAL, target_id, 15, comm_inter, &
	286	status, ierr )
[108]	287	ENDIF
[667]	288	CALL MPI_BCAST( remote, 1, MPI_REAL, 0, comm2d, ierr)
[213]	289
[667]	290
	291	IF ( coupling_mode == 'atmosphere_to_ocean') THEN
	292
	293	IF ( dx < remote ) THEN
	294	WRITE( message_string, * ) 'coupling mode "', &
	295	TRIM( coupling_mode ), &
	296	'": dx in Atmosphere is not equal to or not larger then dx in ocean'
	297	CALL message( 'check_parameters', 'PA0009', 1, 2, 0, 6, 0 )
	298	ENDIF
	299
	300	IF ( (nx_a+1)dx /= (nx_o+1)remote ) THEN
	301	WRITE( message_string, * ) 'coupling mode "', &
	302	TRIM( coupling_mode ), &
	303	'": Domain size in x-direction is not equal in ocean and atmosphere'
	304	CALL message( 'check_parameters', 'PA0010', 1, 2, 0, 6, 0 )
	305	ENDIF
	306
[108]	307	ENDIF
[213]	308
[667]	309	IF ( myid == 0) THEN
	310	CALL MPI_SEND( dy, 1, MPI_REAL, target_id, 16, comm_inter, ierr )
	311	CALL MPI_RECV( remote, 1, MPI_REAL, target_id, 16, comm_inter, &
	312	status, ierr )
[108]	313	ENDIF
[667]	314	CALL MPI_BCAST( remote, 1, MPI_REAL, 0, comm2d, ierr)
[213]	315
[667]	316	IF ( coupling_mode == 'atmosphere_to_ocean') THEN
	317
	318	IF ( dy < remote ) THEN
	319	WRITE( message_string, * ) 'coupling mode "', &
	320	TRIM( coupling_mode ), &
	321	'": dy in Atmosphere is not equal to or not larger then dy in ocean'
	322	CALL message( 'check_parameters', 'PA0011', 1, 2, 0, 6, 0 )
	323	ENDIF
	324
	325	IF ( (ny_a+1)dy /= (ny_o+1)remote ) THEN
	326	WRITE( message_string, * ) 'coupling mode "', &
	327	TRIM( coupling_mode ), &
	328	'": Domain size in y-direction is not equal in ocean and atmosphere'
	329	CALL message( 'check_parameters', 'PA0012', 1, 2, 0, 6, 0 )
	330	ENDIF
	331
	332	IF ( MOD(nx_o+1,nx_a+1) /= 0 ) THEN
	333	WRITE( message_string, * ) 'coupling mode "', &
	334	TRIM( coupling_mode ), &
	335	'": nx+1 in ocean is not divisible without remainder with nx+1 in', &
	336	' atmosphere'
	337	CALL message( 'check_parameters', 'PA0339', 1, 2, 0, 6, 0 )
	338	ENDIF
	339
	340	IF ( MOD(ny_o+1,ny_a+1) /= 0 ) THEN
	341	WRITE( message_string, * ) 'coupling mode "', &
	342	TRIM( coupling_mode ), &
	343	'": ny+1 in ocean is not divisible without remainder with ny+1 in', &
	344	' atmosphere'
	345	CALL message( 'check_parameters', 'PA0340', 1, 2, 0, 6, 0 )
	346	ENDIF
	347
[108]	348	ENDIF
[222]	349	#else
	350	WRITE( message_string, * ) 'coupling requires PALM to be called with', &
	351	' ''mrun -K parallel'''
[226]	352	CALL message( 'check_parameters', 'PA0141', 1, 2, 0, 6, 0 )
[108]	353	#endif
	354	ENDIF
	355
[206]	356	#if defined( __parallel )
[108]	357	!
	358	!-- Exchange via intercommunicator
[667]	359	IF ( coupling_mode == 'atmosphere_to_ocean' .AND. myid == 0 ) THEN
[206]	360	CALL MPI_SEND( humidity, 1, MPI_LOGICAL, target_id, 19, comm_inter, &
	361	ierr )
[667]	362	ELSEIF ( coupling_mode == 'ocean_to_atmosphere' .AND. myid == 0) THEN
[206]	363	CALL MPI_RECV( humidity_remote, 1, MPI_LOGICAL, target_id, 19, &
	364	comm_inter, status, ierr )
[108]	365	ENDIF
[667]	366	CALL MPI_BCAST( humidity_remote, 1, MPI_LOGICAL, 0, comm2d, ierr)
	367
[108]	368	#endif
	369
	370
	371	!
[1]	372	!-- Generate the file header which is used as a header for most of PALM's
	373	!-- output files
	374	CALL DATE_AND_TIME( date, time )
	375	run_date = date(7:8)//'-'//date(5:6)//'-'//date(3:4)
	376	run_time = time(1:2)//':'//time(3:4)//':'//time(5:6)
[102]	377	IF ( coupling_mode == 'uncoupled' ) THEN
	378	coupling_string = ''
	379	ELSEIF ( coupling_mode == 'atmosphere_to_ocean' ) THEN
	380	coupling_string = ' coupled (atmosphere)'
	381	ELSEIF ( coupling_mode == 'ocean_to_atmosphere' ) THEN
	382	coupling_string = ' coupled (ocean)'
	383	ENDIF
[1]	384
[102]	385	WRITE ( run_description_header, &
	386	'(A,2X,A,2X,A,A,A,I2.2,A,2X,A,A,2X,A,1X,A)' ) &
	387	TRIM( version ), TRIM( revision ), 'run: ', &
	388	TRIM( run_identifier ), '.', runnr, TRIM( coupling_string ), &
	389	'host: ', TRIM( host ), run_date, run_time
[1]	390
	391	!
[63]	392	!-- Check the general loop optimization method
	393	IF ( loop_optimization == 'default' ) THEN
	394	IF ( host(1:3) == 'nec' ) THEN
	395	loop_optimization = 'vector'
	396	ELSE
	397	loop_optimization = 'cache'
	398	ENDIF
	399	ENDIF
	400	IF ( loop_optimization /= 'noopt' .AND. loop_optimization /= 'cache' &
	401	.AND. loop_optimization /= 'vector' ) THEN
[213]	402	message_string = 'illegal value given for loop_optimization: "' // &
	403	TRIM( loop_optimization ) // '"'
[226]	404	CALL message( 'check_parameters', 'PA0013', 1, 2, 0, 6, 0 )
[63]	405	ENDIF
	406
	407	!
[1]	408	!-- Check topography setting (check for illegal parameter combinations)
	409	IF ( topography /= 'flat' ) THEN
	410	action = ' '
	411	IF ( scalar_advec /= 'pw-scheme' ) THEN
	412	WRITE( action, '(A,A)' ) 'scalar_advec = ', scalar_advec
	413	ENDIF
	414	IF ( momentum_advec /= 'pw-scheme' ) THEN
	415	WRITE( action, '(A,A)' ) 'momentum_advec = ', momentum_advec
	416	ENDIF
[51]	417	IF ( timestep_scheme(1:8) == 'leapfrog' ) THEN
	418	WRITE( action, '(A,A)' ) 'timestep_scheme = ', timestep_scheme
	419	ENDIF
[114]	420	IF ( psolver == 'sor' ) THEN
[1]	421	WRITE( action, '(A,A)' ) 'psolver = ', psolver
	422	ENDIF
	423	IF ( sloping_surface ) THEN
	424	WRITE( action, '(A)' ) 'sloping surface = .TRUE.'
	425	ENDIF
	426	IF ( galilei_transformation ) THEN
	427	WRITE( action, '(A)' ) 'galilei_transformation = .TRUE.'
	428	ENDIF
	429	IF ( cloud_physics ) THEN
	430	WRITE( action, '(A)' ) 'cloud_physics = .TRUE.'
	431	ENDIF
	432	IF ( cloud_droplets ) THEN
	433	WRITE( action, '(A)' ) 'cloud_droplets = .TRUE.'
	434	ENDIF
	435	IF ( .NOT. prandtl_layer ) THEN
	436	WRITE( action, '(A)' ) 'prandtl_layer = .FALSE.'
	437	ENDIF
	438	IF ( action /= ' ' ) THEN
[213]	439	message_string = 'a non-flat topography does not allow ' // &
	440	TRIM( action )
[226]	441	CALL message( 'check_parameters', 'PA0014', 1, 2, 0, 6, 0 )
[1]	442	ENDIF
[667]	443	IF ( momentum_advec == 'ws-scheme' .OR. scalar_advec == 'ws-scheme' ) &
	444	THEN
	445	message_string = 'topography is still not allowed with ' // &
	446	'momentum_advec = "' // TRIM( momentum_advec ) // &
	447	'"or scalar_advec = "' // TRIM( scalar_advec ) //'"'
	448	! message number still needs modification
	449	CALL message( 'check_parameters', 'PA0341', 1, 2, 0, 6, 0 )
	450	END IF
	451
[256]	452	!
	453	!-- In case of non-flat topography, check whether the convention how to
	454	!-- define the topography grid has been set correctly, or whether the default
	455	!-- is applicable. If this is not possible, abort.
	456	IF ( TRIM( topography_grid_convention ) == ' ' ) THEN
	457	IF ( TRIM( topography ) /= 'single_building' .AND. &
	458	TRIM( topography ) /= 'single_street_canyon' .AND. &
	459	TRIM( topography ) /= 'read_from_file' ) THEN
	460	!-- The default value is not applicable here, because it is only valid
	461	!-- for the two standard cases 'single_building' and 'read_from_file'
	462	!-- defined in init_grid.
	463	WRITE( message_string, * ) &
	464	'The value for "topography_grid_convention" ', &
	465	'is not set. Its default value is & only valid for ', &
	466	'"topography" = ''single_building'', ', &
	467	'''single_street_canyon'' & or ''read_from_file''.', &
	468	' & Choose ''cell_edge'' or ''cell_center''.'
	469	CALL message( 'user_check_parameters', 'PA0239', 1, 2, 0, 6, 0 )
	470	ELSE
	471	!-- The default value is applicable here.
	472	!-- Set convention according to topography.
	473	IF ( TRIM( topography ) == 'single_building' .OR. &
	474	TRIM( topography ) == 'single_street_canyon' ) THEN
	475	topography_grid_convention = 'cell_edge'
	476	ELSEIF ( TRIM( topography ) == 'read_from_file' ) THEN
	477	topography_grid_convention = 'cell_center'
	478	ENDIF
	479	ENDIF
	480	ELSEIF ( TRIM( topography_grid_convention ) /= 'cell_edge' .AND. &
	481	TRIM( topography_grid_convention ) /= 'cell_center' ) THEN
	482	WRITE( message_string, * ) &
	483	'The value for "topography_grid_convention" is ', &
	484	'not recognized. & Choose ''cell_edge'' or ''cell_center''.'
	485	CALL message( 'user_check_parameters', 'PA0240', 1, 2, 0, 6, 0 )
	486	ENDIF
	487
[1]	488	ENDIF
[94]	489
[1]	490	!
[94]	491	!-- Check ocean setting
	492	IF ( ocean ) THEN
[332]	493
[94]	494	action = ' '
	495	IF ( timestep_scheme(1:8) == 'leapfrog' ) THEN
	496	WRITE( action, '(A,A)' ) 'timestep_scheme = ', timestep_scheme
	497	ENDIF
[97]	498	IF ( momentum_advec == 'ups-scheme' ) THEN
	499	WRITE( action, '(A,A)' ) 'momentum_advec = ', momentum_advec
	500	ENDIF
[94]	501	IF ( action /= ' ' ) THEN
[213]	502	message_string = 'ocean = .T. does not allow ' // TRIM( action )
[226]	503	CALL message( 'check_parameters', 'PA0015', 1, 2, 0, 6, 0 )
[94]	504	ENDIF
	505
[332]	506	ELSEIF ( TRIM( coupling_mode ) == 'uncoupled' .AND. &
	507	TRIM( coupling_char ) == '_O' ) THEN
[343]	508
[94]	509	!
[332]	510	!-- Check whether an (uncoupled) atmospheric run has been declared as an
	511	!-- ocean run (this setting is done via mrun-option -y)
[343]	512
[291]	513	message_string = 'ocean = .F. does not allow coupling_char = "' // &
	514	TRIM( coupling_char ) // '" set by mrun-option "-y"'
	515	CALL message( 'check_parameters', 'PA0317', 1, 2, 0, 6, 0 )
[332]	516
[291]	517	ENDIF
	518
	519	!
[1]	520	!-- Check whether there are any illegal values
	521	!-- Pressure solver:
	522	IF ( psolver /= 'poisfft' .AND. psolver /= 'poisfft_hybrid' .AND. &
	523	psolver /= 'sor' .AND. psolver /= 'multigrid' ) THEN
[213]	524	message_string = 'unknown solver for perturbation pressure: psolver' // &
	525	' = "' // TRIM( psolver ) // '"'
[226]	526	CALL message( 'check_parameters', 'PA0016', 1, 2, 0, 6, 0 )
[1]	527	ENDIF
	528
	529	#if defined( __parallel )
	530	IF ( psolver == 'poisfft_hybrid' .AND. pdims(2) /= 1 ) THEN
[213]	531	message_string = 'psolver = "' // TRIM( psolver ) // '" only works ' // &
	532	'for a 1d domain-decomposition along x & please do' // &
	533	' not set npey/=1 in the parameter file'
[226]	534	CALL message( 'check_parameters', 'PA0017', 1, 2, 0, 6, 0 )
[1]	535	ENDIF
[181]	536	IF ( psolver == 'poisfft_hybrid' .AND. &
	537	( nxra > nxr .OR. nyna > nyn .OR. nza > nz ) .OR. &
	538	psolver == 'multigrid' .AND. &
	539	( nxra > nxr .OR. nyna > nyn ) ) THEN
[213]	540	message_string = 'psolver = "' // TRIM( psolver ) // '" does not ' // &
	541	'work for subdomains with unequal size & please ' // &
	542	'set grid_matching = ''strict'' in the parameter file'
[226]	543	CALL message( 'check_parameters', 'PA0018', 1, 2, 0, 6, 0 )
[1]	544	ENDIF
	545	#else
	546	IF ( psolver == 'poisfft_hybrid' ) THEN
[213]	547	message_string = 'psolver = "' // TRIM( psolver ) // '" only works' // &
	548	' for a parallel environment'
[226]	549	CALL message( 'check_parameters', 'PA0019', 1, 2, 0, 6, 0 )
[1]	550	ENDIF
	551	#endif
	552
	553	IF ( psolver == 'multigrid' ) THEN
	554	IF ( cycle_mg == 'w' ) THEN
	555	gamma_mg = 2
	556	ELSEIF ( cycle_mg == 'v' ) THEN
	557	gamma_mg = 1
	558	ELSE
[213]	559	message_string = 'unknown multigrid cycle: cycle_mg = "' // &
	560	TRIM( cycle_mg ) // '"'
[226]	561	CALL message( 'check_parameters', 'PA0020', 1, 2, 0, 6, 0 )
[1]	562	ENDIF
	563	ENDIF
	564
	565	IF ( fft_method /= 'singleton-algorithm' .AND. &
	566	fft_method /= 'temperton-algorithm' .AND. &
	567	fft_method /= 'system-specific' ) THEN
[213]	568	message_string = 'unknown fft-algorithm: fft_method = "' // &
	569	TRIM( fft_method ) // '"'
[226]	570	CALL message( 'check_parameters', 'PA0021', 1, 2, 0, 6, 0 )
[1]	571	ENDIF
[667]	572
	573	IF( momentum_advec == 'ws-scheme' .AND. &
[688]	574	.NOT. call_psolver_at_all_substeps ) THEN
[667]	575	message_string = 'psolver must be called at each RK3 substep when "'//&
	576	TRIM(momentum_advec) // ' "is used for momentum_advec'
[685]	577	CALL message( 'check_parameters', 'PA0344', 1, 2, 0, 6, 0 )
[667]	578	END IF
[1]	579	!
	580	!-- Advection schemes:
[673]	581	!
	582	!-- Set the LOGICALS to enhance the performance.
	583	IF ( momentum_advec == 'ws-scheme' ) ws_scheme_mom = .TRUE.
	584	IF ( scalar_advec == 'ws-scheme' ) ws_scheme_sca = .TRUE.
	585
[667]	586	IF ( momentum_advec /= 'pw-scheme' .AND. momentum_advec /= 'ws-scheme' .AND. &
	587	momentum_advec /= 'ups-scheme' ) THEN
[214]	588	message_string = 'unknown advection scheme: momentum_advec = "' // &
	589	TRIM( momentum_advec ) // '"'
[226]	590	CALL message( 'check_parameters', 'PA0022', 1, 2, 0, 6, 0 )
[1]	591	ENDIF
[667]	592	IF ((( momentum_advec == 'ups-scheme' .OR. scalar_advec == 'ups-scheme' )&
	593	.AND. timestep_scheme /= 'euler' ) .OR. (( momentum_advec == 'ws-scheme'&
	594	.OR. scalar_advec == 'ws-scheme') .AND. (timestep_scheme == 'euler' .OR. &
	595	timestep_scheme == 'leapfrog+euler' .OR. timestep_scheme == 'leapfrog' &
	596	.OR. timestep_scheme == 'runge-kutta-2'))) THEN
	597	message_string = 'momentum_advec or scalar_advec = "' &
	598	// TRIM( momentum_advec ) // '" is not allowed with timestep_scheme = "' // &
	599	TRIM( timestep_scheme ) // '"'
[226]	600	CALL message( 'check_parameters', 'PA0023', 1, 2, 0, 6, 0 )
[1]	601	ENDIF
[667]	602	IF ( scalar_advec /= 'pw-scheme' .AND. scalar_advec /= 'ws-scheme' .AND. &
	603	scalar_advec /= 'bc-scheme' .AND. scalar_advec /= 'ups-scheme' ) THEN
[214]	604	message_string = 'unknown advection scheme: scalar_advec = "' // &
	605	TRIM( scalar_advec ) // '"'
[226]	606	CALL message( 'check_parameters', 'PA0024', 1, 2, 0, 6, 0 )
[1]	607	ENDIF
	608
	609	IF ( use_sgs_for_particles .AND. .NOT. use_upstream_for_tke ) THEN
	610	use_upstream_for_tke = .TRUE.
[214]	611	message_string = 'use_upstream_for_tke set .TRUE. because ' // &
	612	'use_sgs_for_particles = .TRUE.'
[226]	613	CALL message( 'check_parameters', 'PA0025', 0, 1, 0, 6, 0 )
[1]	614	ENDIF
	615
	616	IF ( use_upstream_for_tke .AND. timestep_scheme(1:8) == 'leapfrog' ) THEN
[214]	617	message_string = 'use_upstream_for_tke = .TRUE. not allowed with ' // &
	618	'timestep_scheme = "' // TRIM( timestep_scheme ) // '"'
[226]	619	CALL message( 'check_parameters', 'PA0026', 1, 2, 0, 6, 0 )
[1]	620	ENDIF
	621
	622	!
	623	!-- Timestep schemes:
	624	SELECT CASE ( TRIM( timestep_scheme ) )
	625
	626	CASE ( 'euler' )
	627	intermediate_timestep_count_max = 1
	628	asselin_filter_factor = 0.0
	629
	630	CASE ( 'leapfrog', 'leapfrog+euler' )
	631	intermediate_timestep_count_max = 1
	632
	633	CASE ( 'runge-kutta-2' )
	634	intermediate_timestep_count_max = 2
	635	asselin_filter_factor = 0.0
	636
	637	CASE ( 'runge-kutta-3' )
	638	intermediate_timestep_count_max = 3
	639	asselin_filter_factor = 0.0
	640
	641	CASE DEFAULT
[214]	642	message_string = 'unknown timestep scheme: timestep_scheme = "' // &
	643	TRIM( timestep_scheme ) // '"'
[226]	644	CALL message( 'check_parameters', 'PA0027', 1, 2, 0, 6, 0 )
[1]	645
	646	END SELECT
	647
[63]	648	IF ( scalar_advec == 'ups-scheme' .AND. timestep_scheme(1:5) == 'runge' )&
[1]	649	THEN
[214]	650	message_string = 'scalar advection scheme "' // TRIM( scalar_advec ) // &
	651	'" & does not work with timestep_scheme "' // &
	652	TRIM( timestep_scheme ) // '"'
[226]	653	CALL message( 'check_parameters', 'PA0028', 1, 2, 0, 6, 0 )
[1]	654	ENDIF
	655
[667]	656	IF ( (momentum_advec /= 'pw-scheme' .AND. momentum_advec /= 'ws-scheme') &
	657	.AND. timestep_scheme(1:5) == 'runge' ) THEN
[214]	658	message_string = 'momentum advection scheme "' // &
	659	TRIM( momentum_advec ) // '" & does not work with ' // &
	660	'timestep_scheme "' // TRIM( timestep_scheme ) // '"'
[226]	661	CALL message( 'check_parameters', 'PA0029', 1, 2, 0, 6, 0 )
[1]	662	ENDIF
	663
[147]	664	IF ( TRIM( initializing_actions ) /= 'read_restart_data' .AND. &
[328]	665	TRIM( initializing_actions ) /= 'cyclic_fill' ) THEN
[1]	666	!
[214]	667	!-- No restart run: several initialising actions are possible
[1]	668	action = initializing_actions
	669	DO WHILE ( TRIM( action ) /= '' )
	670	position = INDEX( action, ' ' )
	671	SELECT CASE ( action(1:position-1) )
	672
	673	CASE ( 'set_constant_profiles', 'set_1d-model_profiles', &
[46]	674	'by_user', 'initialize_vortex', 'initialize_ptanom' )
[1]	675	action = action(position+1:)
	676
	677	CASE DEFAULT
[214]	678	message_string = 'initializing_action = "' // &
	679	TRIM( action ) // '" unkown or not allowed'
[226]	680	CALL message( 'check_parameters', 'PA0030', 1, 2, 0, 6, 0 )
[1]	681
	682	END SELECT
	683	ENDDO
	684	ENDIF
[214]	685
[680]	686	IF ( TRIM( initializing_actions ) == 'initialize_vortex' .AND. &
	687	conserve_volume_flow ) THEN
	688	message_string = 'initializing_actions = "initialize_vortex"' // &
	689	' ist not allowed with conserve_volume_flow = .T.'
	690	CALL message( 'check_parameters', 'PA0343', 1, 2, 0, 6, 0 )
	691	ENDIF
	692
	693
[1]	694	IF ( INDEX( initializing_actions, 'set_constant_profiles' ) /= 0 .AND. &
	695	INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 ) THEN
[214]	696	message_string = 'initializing_actions = "set_constant_profiles"' // &
	697	' and "set_1d-model_profiles" are not allowed ' // &
	698	'simultaneously'
[226]	699	CALL message( 'check_parameters', 'PA0031', 1, 2, 0, 6, 0 )
[1]	700	ENDIF
[214]	701
[46]	702	IF ( INDEX( initializing_actions, 'set_constant_profiles' ) /= 0 .AND. &
	703	INDEX( initializing_actions, 'by_user' ) /= 0 ) THEN
[214]	704	message_string = 'initializing_actions = "set_constant_profiles"' // &
	705	' and "by_user" are not allowed simultaneously'
[226]	706	CALL message( 'check_parameters', 'PA0032', 1, 2, 0, 6, 0 )
[46]	707	ENDIF
[214]	708
[46]	709	IF ( INDEX( initializing_actions, 'by_user' ) /= 0 .AND. &
	710	INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 ) THEN
[214]	711	message_string = 'initializing_actions = "by_user" and ' // &
	712	'"set_1d-model_profiles" are not allowed simultaneously'
[226]	713	CALL message( 'check_parameters', 'PA0033', 1, 2, 0, 6, 0 )
[46]	714	ENDIF
[1]	715
[75]	716	IF ( cloud_physics .AND. .NOT. humidity ) THEN
[214]	717	WRITE( message_string, * ) 'cloud_physics = ', cloud_physics, ' is ', &
	718	'not allowed with humidity = ', humidity
[226]	719	CALL message( 'check_parameters', 'PA0034', 1, 2, 0, 6, 0 )
[1]	720	ENDIF
	721
[72]	722	IF ( precipitation .AND. .NOT. cloud_physics ) THEN
[214]	723	WRITE( message_string, * ) 'precipitation = ', precipitation, ' is ', &
	724	'not allowed with cloud_physics = ', cloud_physics
[226]	725	CALL message( 'check_parameters', 'PA0035', 1, 2, 0, 6, 0 )
[72]	726	ENDIF
	727
[75]	728	IF ( humidity .AND. sloping_surface ) THEN
[214]	729	message_string = 'humidity = .TRUE. and sloping_surface = .TRUE. ' // &
	730	'are not allowed simultaneously'
[226]	731	CALL message( 'check_parameters', 'PA0036', 1, 2, 0, 6, 0 )
[1]	732	ENDIF
	733
[75]	734	IF ( humidity .AND. scalar_advec == 'ups-scheme' ) THEN
[214]	735	message_string = 'UPS-scheme is not implemented for humidity = .TRUE.'
[226]	736	CALL message( 'check_parameters', 'PA0037', 1, 2, 0, 6, 0 )
[1]	737	ENDIF
	738
[75]	739	IF ( passive_scalar .AND. humidity ) THEN
[214]	740	message_string = 'humidity = .TRUE. and passive_scalar = .TRUE. ' // &
	741	'is not allowed simultaneously'
[226]	742	CALL message( 'check_parameters', 'PA0038', 1, 2, 0, 6, 0 )
[1]	743	ENDIF
	744
	745	IF ( passive_scalar .AND. scalar_advec == 'ups-scheme' ) THEN
[214]	746	message_string = 'UPS-scheme is not implemented for passive_scalar' // &
	747	' = .TRUE.'
[226]	748	CALL message( 'check_parameters', 'PA0039', 1, 2, 0, 6, 0 )
[1]	749	ENDIF
	750
	751	IF ( grid_matching /= 'strict' .AND. grid_matching /= 'match' ) THEN
[214]	752	message_string = 'illegal value "' // TRIM( grid_matching ) // &
	753	'" found for parameter grid_matching'
[226]	754	CALL message( 'check_parameters', 'PA0040', 1, 2, 0, 6, 0 )
[1]	755	ENDIF
	756
[138]	757	IF ( plant_canopy .AND. ( drag_coefficient == 0.0 ) ) THEN
[214]	758	message_string = 'plant_canopy = .TRUE. requires a non-zero drag ' // &
	759	'coefficient & given value is drag_coefficient = 0.0'
[226]	760	CALL message( 'check_parameters', 'PA0041', 1, 2, 0, 6, 0 )
[138]	761	ENDIF
	762
[1]	763	!
	764	!-- In case of no model continuation run, check initialising parameters and
	765	!-- deduce further quantities
	766	IF ( TRIM( initializing_actions ) /= 'read_restart_data' ) THEN
	767
	768	!
	769	!-- Initial profiles for 1D and 3D model, respectively
	770	u_init = ug_surface
	771	v_init = vg_surface
	772	pt_init = pt_surface
[94]	773	IF ( humidity ) q_init = q_surface
	774	IF ( ocean ) sa_init = sa_surface
	775	IF ( passive_scalar ) q_init = s_surface
[138]	776	IF ( plant_canopy ) lad = 0.0
[1]	777
	778	!
	779	!--
	780	!-- If required, compute initial profile of the geostrophic wind
	781	!-- (component ug)
	782	i = 1
	783	gradient = 0.0
[97]	784
	785	IF ( .NOT. ocean ) THEN
	786
	787	ug_vertical_gradient_level_ind(1) = 0
	788	ug(0) = ug_surface
	789	DO k = 1, nzt+1
[177]	790	IF ( i < 11 ) THEN
	791	IF ( ug_vertical_gradient_level(i) < zu(k) .AND. &
	792	ug_vertical_gradient_level(i) >= 0.0 ) THEN
	793	gradient = ug_vertical_gradient(i) / 100.0
	794	ug_vertical_gradient_level_ind(i) = k - 1
	795	i = i + 1
[1]	796	ENDIF
[177]	797	ENDIF
[97]	798	IF ( gradient /= 0.0 ) THEN
	799	IF ( k /= 1 ) THEN
	800	ug(k) = ug(k-1) + dzu(k) * gradient
	801	ELSE
	802	ug(k) = ug_surface + 0.5 * dzu(k) * gradient
	803	ENDIF
[1]	804	ELSE
[97]	805	ug(k) = ug(k-1)
[1]	806	ENDIF
[97]	807	ENDDO
[1]	808
[97]	809	ELSE
	810
	811	ug_vertical_gradient_level_ind(1) = nzt+1
[121]	812	ug(nzt+1) = ug_surface
[667]	813	DO k = nzt, nzb, -1
[177]	814	IF ( i < 11 ) THEN
	815	IF ( ug_vertical_gradient_level(i) > zu(k) .AND. &
	816	ug_vertical_gradient_level(i) <= 0.0 ) THEN
	817	gradient = ug_vertical_gradient(i) / 100.0
	818	ug_vertical_gradient_level_ind(i) = k + 1
	819	i = i + 1
[97]	820	ENDIF
	821	ENDIF
	822	IF ( gradient /= 0.0 ) THEN
	823	IF ( k /= nzt ) THEN
	824	ug(k) = ug(k+1) - dzu(k+1) * gradient
	825	ELSE
	826	ug(k) = ug_surface - 0.5 * dzu(k+1) * gradient
	827	ug(k+1) = ug_surface + 0.5 * dzu(k+1) * gradient
	828	ENDIF
	829	ELSE
	830	ug(k) = ug(k+1)
	831	ENDIF
	832	ENDDO
	833
	834	ENDIF
	835
[1]	836	u_init = ug
	837
	838	!
	839	!-- In case of no given gradients for ug, choose a vanishing gradient
[97]	840	IF ( ug_vertical_gradient_level(1) == -9999999.9 ) THEN
[1]	841	ug_vertical_gradient_level(1) = 0.0
	842	ENDIF
	843
	844	!
	845	!--
	846	!-- If required, compute initial profile of the geostrophic wind
	847	!-- (component vg)
	848	i = 1
	849	gradient = 0.0
[97]	850
	851	IF ( .NOT. ocean ) THEN
	852
	853	vg_vertical_gradient_level_ind(1) = 0
	854	vg(0) = vg_surface
	855	DO k = 1, nzt+1
[177]	856	IF ( i < 11 ) THEN
	857	IF ( vg_vertical_gradient_level(i) < zu(k) .AND. &
	858	vg_vertical_gradient_level(i) >= 0.0 ) THEN
	859	gradient = vg_vertical_gradient(i) / 100.0
	860	vg_vertical_gradient_level_ind(i) = k - 1
	861	i = i + 1
[1]	862	ENDIF
	863	ENDIF
[97]	864	IF ( gradient /= 0.0 ) THEN
	865	IF ( k /= 1 ) THEN
	866	vg(k) = vg(k-1) + dzu(k) * gradient
	867	ELSE
	868	vg(k) = vg_surface + 0.5 * dzu(k) * gradient
	869	ENDIF
[1]	870	ELSE
[97]	871	vg(k) = vg(k-1)
[1]	872	ENDIF
[97]	873	ENDDO
[1]	874
[97]	875	ELSE
	876
[121]	877	vg_vertical_gradient_level_ind(1) = nzt+1
	878	vg(nzt+1) = vg_surface
[667]	879	DO k = nzt, nzb, -1
[177]	880	IF ( i < 11 ) THEN
	881	IF ( vg_vertical_gradient_level(i) > zu(k) .AND. &
	882	vg_vertical_gradient_level(i) <= 0.0 ) THEN
	883	gradient = vg_vertical_gradient(i) / 100.0
	884	vg_vertical_gradient_level_ind(i) = k + 1
	885	i = i + 1
[97]	886	ENDIF
	887	ENDIF
	888	IF ( gradient /= 0.0 ) THEN
	889	IF ( k /= nzt ) THEN
	890	vg(k) = vg(k+1) - dzu(k+1) * gradient
	891	ELSE
	892	vg(k) = vg_surface - 0.5 * dzu(k+1) * gradient
	893	vg(k+1) = vg_surface + 0.5 * dzu(k+1) * gradient
	894	ENDIF
	895	ELSE
	896	vg(k) = vg(k+1)
	897	ENDIF
	898	ENDDO
	899
	900	ENDIF
	901
[1]	902	v_init = vg
	903
	904	!
	905	!-- In case of no given gradients for vg, choose a vanishing gradient
[97]	906	IF ( vg_vertical_gradient_level(1) == -9999999.9 ) THEN
[1]	907	vg_vertical_gradient_level(1) = 0.0
	908	ENDIF
	909
	910	!
[94]	911	!-- Compute initial temperature profile using the given temperature gradients
[1]	912	i = 1
	913	gradient = 0.0
[94]	914
	915	IF ( .NOT. ocean ) THEN
	916
	917	pt_vertical_gradient_level_ind(1) = 0
	918	DO k = 1, nzt+1
[177]	919	IF ( i < 11 ) THEN
	920	IF ( pt_vertical_gradient_level(i) < zu(k) .AND. &
	921	pt_vertical_gradient_level(i) >= 0.0 ) THEN
	922	gradient = pt_vertical_gradient(i) / 100.0
	923	pt_vertical_gradient_level_ind(i) = k - 1
	924	i = i + 1
[1]	925	ENDIF
	926	ENDIF
[94]	927	IF ( gradient /= 0.0 ) THEN
	928	IF ( k /= 1 ) THEN
	929	pt_init(k) = pt_init(k-1) + dzu(k) * gradient
	930	ELSE
	931	pt_init(k) = pt_surface + 0.5 * dzu(k) * gradient
	932	ENDIF
[1]	933	ELSE
[94]	934	pt_init(k) = pt_init(k-1)
[1]	935	ENDIF
[94]	936	ENDDO
[1]	937
[94]	938	ELSE
	939
	940	pt_vertical_gradient_level_ind(1) = nzt+1
	941	DO k = nzt, 0, -1
[177]	942	IF ( i < 11 ) THEN
	943	IF ( pt_vertical_gradient_level(i) > zu(k) .AND. &
	944	pt_vertical_gradient_level(i) <= 0.0 ) THEN
	945	gradient = pt_vertical_gradient(i) / 100.0
	946	pt_vertical_gradient_level_ind(i) = k + 1
	947	i = i + 1
[94]	948	ENDIF
	949	ENDIF
	950	IF ( gradient /= 0.0 ) THEN
	951	IF ( k /= nzt ) THEN
	952	pt_init(k) = pt_init(k+1) - dzu(k+1) * gradient
	953	ELSE
	954	pt_init(k) = pt_surface - 0.5 * dzu(k+1) * gradient
	955	pt_init(k+1) = pt_surface + 0.5 * dzu(k+1) * gradient
	956	ENDIF
	957	ELSE
	958	pt_init(k) = pt_init(k+1)
	959	ENDIF
	960	ENDDO
	961
	962	ENDIF
	963
[1]	964	!
	965	!-- In case of no given temperature gradients, choose gradient of neutral
	966	!-- stratification
[94]	967	IF ( pt_vertical_gradient_level(1) == -9999999.9 ) THEN
[1]	968	pt_vertical_gradient_level(1) = 0.0
	969	ENDIF
	970
	971	!
[94]	972	!-- Store temperature gradient at the top boundary for possible Neumann
[1]	973	!-- boundary condition
[19]	974	bc_pt_t_val = ( pt_init(nzt+1) - pt_init(nzt) ) / dzu(nzt+1)
[1]	975
	976	!
	977	!-- If required, compute initial humidity or scalar profile using the given
	978	!-- humidity/scalar gradient. In case of scalar transport, initially store
	979	!-- values of the scalar parameters on humidity parameters
	980	IF ( passive_scalar ) THEN
	981	bc_q_b = bc_s_b
	982	bc_q_t = bc_s_t
	983	q_surface = s_surface
	984	q_surface_initial_change = s_surface_initial_change
	985	q_vertical_gradient = s_vertical_gradient
	986	q_vertical_gradient_level = s_vertical_gradient_level
	987	surface_waterflux = surface_scalarflux
[407]	988	wall_humidityflux = wall_scalarflux
[1]	989	ENDIF
	990
[75]	991	IF ( humidity .OR. passive_scalar ) THEN
[1]	992
	993	i = 1
	994	gradient = 0.0
	995	q_vertical_gradient_level_ind(1) = 0
	996	DO k = 1, nzt+1
[177]	997	IF ( i < 11 ) THEN
	998	IF ( q_vertical_gradient_level(i) < zu(k) .AND. &
	999	q_vertical_gradient_level(i) >= 0.0 ) THEN
	1000	gradient = q_vertical_gradient(i) / 100.0
	1001	q_vertical_gradient_level_ind(i) = k - 1
	1002	i = i + 1
[1]	1003	ENDIF
	1004	ENDIF
	1005	IF ( gradient /= 0.0 ) THEN
	1006	IF ( k /= 1 ) THEN
	1007	q_init(k) = q_init(k-1) + dzu(k) * gradient
	1008	ELSE
	1009	q_init(k) = q_init(k-1) + 0.5 * dzu(k) * gradient
	1010	ENDIF
	1011	ELSE
	1012	q_init(k) = q_init(k-1)
	1013	ENDIF
[72]	1014	!
	1015	!-- Avoid negative humidities
	1016	IF ( q_init(k) < 0.0 ) THEN
	1017	q_init(k) = 0.0
	1018	ENDIF
[1]	1019	ENDDO
	1020
	1021	!
	1022	!-- In case of no given humidity gradients, choose zero gradient
	1023	!-- conditions
	1024	IF ( q_vertical_gradient_level(1) == -1.0 ) THEN
	1025	q_vertical_gradient_level(1) = 0.0
	1026	ENDIF
	1027
	1028	!
	1029	!-- Store humidity gradient at the top boundary for possile Neumann
	1030	!-- boundary condition
[19]	1031	bc_q_t_val = ( q_init(nzt+1) - q_init(nzt) ) / dzu(nzt+1)
[1]	1032
	1033	ENDIF
	1034
[94]	1035	!
	1036	!-- If required, compute initial salinity profile using the given salinity
	1037	!-- gradients
	1038	IF ( ocean ) THEN
	1039
	1040	i = 1
	1041	gradient = 0.0
	1042
	1043	sa_vertical_gradient_level_ind(1) = nzt+1
	1044	DO k = nzt, 0, -1
[177]	1045	IF ( i < 11 ) THEN
	1046	IF ( sa_vertical_gradient_level(i) > zu(k) .AND. &
	1047	sa_vertical_gradient_level(i) <= 0.0 ) THEN
	1048	gradient = sa_vertical_gradient(i) / 100.0
	1049	sa_vertical_gradient_level_ind(i) = k + 1
	1050	i = i + 1
[94]	1051	ENDIF
	1052	ENDIF
	1053	IF ( gradient /= 0.0 ) THEN
	1054	IF ( k /= nzt ) THEN
	1055	sa_init(k) = sa_init(k+1) - dzu(k+1) * gradient
	1056	ELSE
	1057	sa_init(k) = sa_surface - 0.5 * dzu(k+1) * gradient
	1058	sa_init(k+1) = sa_surface + 0.5 * dzu(k+1) * gradient
	1059	ENDIF
	1060	ELSE
	1061	sa_init(k) = sa_init(k+1)
	1062	ENDIF
	1063	ENDDO
	1064
	1065	ENDIF
	1066
[138]	1067	!
[388]	1068	!-- If required compute the profile of leaf area density used in the plant
	1069	!-- canopy model
[138]	1070	IF ( plant_canopy ) THEN
	1071
	1072	i = 1
	1073	gradient = 0.0
[1]	1074
[138]	1075	IF ( .NOT. ocean ) THEN
[153]	1076
	1077	lad(0) = lad_surface
[138]	1078
	1079	lad_vertical_gradient_level_ind(1) = 0
	1080	DO k = 1, pch_index
[177]	1081	IF ( i < 11 ) THEN
	1082	IF ( lad_vertical_gradient_level(i) < zu(k) .AND. &
	1083	lad_vertical_gradient_level(i) >= 0.0 ) THEN
	1084	gradient = lad_vertical_gradient(i)
	1085	lad_vertical_gradient_level_ind(i) = k - 1
	1086	i = i + 1
[138]	1087	ENDIF
	1088	ENDIF
	1089	IF ( gradient /= 0.0 ) THEN
	1090	IF ( k /= 1 ) THEN
	1091	lad(k) = lad(k-1) + dzu(k) * gradient
	1092	ELSE
	1093	lad(k) = lad_surface + 0.5 * dzu(k) *gradient
	1094	ENDIF
	1095	ELSE
	1096	lad(k) = lad(k-1)
	1097	ENDIF
	1098	ENDDO
	1099
	1100	ENDIF
	1101
[1]	1102	!
[388]	1103	!-- In case of no given leaf area density gradients, choose a vanishing
	1104	!-- gradient
[138]	1105	IF ( lad_vertical_gradient_level(1) == -9999999.9 ) THEN
	1106	lad_vertical_gradient_level(1) = 0.0
	1107	ENDIF
	1108
	1109	ENDIF
	1110
	1111	ENDIF
[411]	1112
	1113	!
	1114	!-- Initialize large scale subsidence if required
[580]	1115	IF ( subs_vertical_gradient_level(1) /= -9999999.9 ) THEN
[411]	1116	large_scale_subsidence = .TRUE.
	1117	CALL init_w_subsidence
	1118	END IF
	1119
[138]	1120
[667]	1121
[138]	1122	!
[1]	1123	!-- Compute Coriolis parameter
	1124	f = 2.0 * omega * SIN( phi / 180.0 * pi )
	1125	fs = 2.0 * omega * COS( phi / 180.0 * pi )
	1126
	1127	!
[97]	1128	!-- Ocean runs always use reference values in the buoyancy term. Therefore
	1129	!-- set the reference temperature equal to the surface temperature.
	1130	IF ( ocean .AND. pt_reference == 9999999.9 ) pt_reference = pt_surface
[57]	1131
	1132	!
[97]	1133	!-- Reference value has to be used in buoyancy terms
	1134	IF ( pt_reference /= 9999999.9 ) use_reference = .TRUE.
	1135
	1136	!
	1137	!-- Sign of buoyancy/stability terms
	1138	IF ( ocean ) atmos_ocean_sign = -1.0
	1139
	1140	!
[108]	1141	!-- Ocean version must use flux boundary conditions at the top
	1142	IF ( ocean .AND. .NOT. use_top_fluxes ) THEN
[215]	1143	message_string = 'use_top_fluxes must be .TRUE. in ocean version'
[226]	1144	CALL message( 'check_parameters', 'PA0042', 1, 2, 0, 6, 0 )
[108]	1145	ENDIF
[97]	1146
	1147	!
[1]	1148	!-- In case of a given slope, compute the relevant quantities
	1149	IF ( alpha_surface /= 0.0 ) THEN
	1150	IF ( ABS( alpha_surface ) > 90.0 ) THEN
[215]	1151	WRITE( message_string, * ) 'ABS( alpha_surface = ', alpha_surface, &
	1152	' ) must be < 90.0'
[226]	1153	CALL message( 'check_parameters', 'PA0043', 1, 2, 0, 6, 0 )
[1]	1154	ENDIF
	1155	sloping_surface = .TRUE.
	1156	cos_alpha_surface = COS( alpha_surface / 180.0 * pi )
	1157	sin_alpha_surface = SIN( alpha_surface / 180.0 * pi )
	1158	ENDIF
	1159
	1160	!
	1161	!-- Check time step and cfl_factor
	1162	IF ( dt /= -1.0 ) THEN
	1163	IF ( dt <= 0.0 .AND. dt /= -1.0 ) THEN
[215]	1164	WRITE( message_string, * ) 'dt = ', dt , ' <= 0.0'
[226]	1165	CALL message( 'check_parameters', 'PA0044', 1, 2, 0, 6, 0 )
[1]	1166	ENDIF
	1167	dt_3d = dt
	1168	dt_fixed = .TRUE.
	1169	ENDIF
	1170
	1171	IF ( cfl_factor <= 0.0 .OR. cfl_factor > 1.0 ) THEN
	1172	IF ( cfl_factor == -1.0 ) THEN
	1173	IF ( momentum_advec == 'ups-scheme' .OR. &
	1174	scalar_advec == 'ups-scheme' ) THEN
	1175	cfl_factor = 0.8
	1176	ELSE
	1177	IF ( timestep_scheme == 'runge-kutta-2' ) THEN
	1178	cfl_factor = 0.8
	1179	ELSEIF ( timestep_scheme == 'runge-kutta-3' ) THEN
	1180	cfl_factor = 0.9
	1181	ELSE
	1182	cfl_factor = 0.1
	1183	ENDIF
	1184	ENDIF
	1185	ELSE
[215]	1186	WRITE( message_string, * ) 'cfl_factor = ', cfl_factor, &
	1187	' out of range & 0.0 < cfl_factor <= 1.0 is required'
[226]	1188	CALL message( 'check_parameters', 'PA0045', 1, 2, 0, 6, 0 )
[1]	1189	ENDIF
	1190	ENDIF
	1191
	1192	!
	1193	!-- Store simulated time at begin
	1194	simulated_time_at_begin = simulated_time
	1195
	1196	!
[291]	1197	!-- Store reference time for coupled runs and change the coupling flag,
	1198	!-- if ...
	1199	IF ( simulated_time == 0.0 ) THEN
	1200	IF ( coupling_start_time == 0.0 ) THEN
	1201	time_since_reference_point = 0.0
	1202	ELSEIF ( time_since_reference_point < 0.0 ) THEN
	1203	run_coupled = .FALSE.
	1204	ENDIF
	1205	ENDIF
	1206
	1207	!
[1]	1208	!-- Set wind speed in the Galilei-transformed system
	1209	IF ( galilei_transformation ) THEN
	1210	IF ( use_ug_for_galilei_tr .AND. &
	1211	ug_vertical_gradient_level(1) == 0.0 .AND. &
	1212	vg_vertical_gradient_level(1) == 0.0 ) THEN
	1213	u_gtrans = ug_surface
	1214	v_gtrans = vg_surface
	1215	ELSEIF ( use_ug_for_galilei_tr .AND. &
	1216	ug_vertical_gradient_level(1) /= 0.0 ) THEN
[215]	1217	message_string = 'baroclinicity (ug) not allowed simultaneously' // &
	1218	' with galilei transformation'
[226]	1219	CALL message( 'check_parameters', 'PA0046', 1, 2, 0, 6, 0 )
[1]	1220	ELSEIF ( use_ug_for_galilei_tr .AND. &
	1221	vg_vertical_gradient_level(1) /= 0.0 ) THEN
[215]	1222	message_string = 'baroclinicity (vg) not allowed simultaneously' // &
	1223	' with galilei transformation'
[226]	1224	CALL message( 'check_parameters', 'PA0047', 1, 2, 0, 6, 0 )
[1]	1225	ELSE
[215]	1226	message_string = 'variable translation speed used for galilei-' // &
	1227	'transformation, which may cause & instabilities in stably ' // &
	1228	'stratified regions'
[226]	1229	CALL message( 'check_parameters', 'PA0048', 0, 1, 0, 6, 0 )
[1]	1230	ENDIF
	1231	ENDIF
	1232
	1233	!
	1234	!-- In case of using a prandtl-layer, calculated (or prescribed) surface
	1235	!-- fluxes have to be used in the diffusion-terms
	1236	IF ( prandtl_layer ) use_surface_fluxes = .TRUE.
	1237
	1238	!
	1239	!-- Check boundary conditions and set internal variables:
	1240	!-- Lateral boundary conditions
[73]	1241	IF ( bc_lr /= 'cyclic' .AND. bc_lr /= 'dirichlet/radiation' .AND. &
	1242	bc_lr /= 'radiation/dirichlet' ) THEN
[215]	1243	message_string = 'unknown boundary condition: bc_lr = "' // &
	1244	TRIM( bc_lr ) // '"'
[226]	1245	CALL message( 'check_parameters', 'PA0049', 1, 2, 0, 6, 0 )
[1]	1246	ENDIF
[73]	1247	IF ( bc_ns /= 'cyclic' .AND. bc_ns /= 'dirichlet/radiation' .AND. &
	1248	bc_ns /= 'radiation/dirichlet' ) THEN
[215]	1249	message_string = 'unknown boundary condition: bc_ns = "' // &
	1250	TRIM( bc_ns ) // '"'
[226]	1251	CALL message( 'check_parameters', 'PA0050', 1, 2, 0, 6, 0 )
[1]	1252	ENDIF
	1253
	1254	!
[366]	1255	!-- Internal variables used for speed optimization in if clauses
[707]	1256	IF ( bc_lr /= 'cyclic' ) bc_lr_cyc = .FALSE.
	1257	IF ( bc_lr == 'dirichlet/radiation' ) bc_lr_dirrad = .TRUE.
	1258	IF ( bc_lr == 'radiation/dirichlet' ) bc_lr_raddir = .TRUE.
	1259	IF ( bc_ns /= 'cyclic' ) bc_ns_cyc = .FALSE.
	1260	IF ( bc_ns == 'dirichlet/radiation' ) bc_ns_dirrad = .TRUE.
	1261	IF ( bc_ns == 'radiation/dirichlet' ) bc_ns_raddir = .TRUE.
[366]	1262
	1263	!
[1]	1264	!-- Non-cyclic lateral boundaries require the multigrid method and Piascek-
[667]	1265	!-- Willimas or Wicker - Skamarock advection scheme. Several schemes
	1266	!-- and tools do not work with non-cyclic boundary conditions.
[1]	1267	IF ( bc_lr /= 'cyclic' .OR. bc_ns /= 'cyclic' ) THEN
	1268	IF ( psolver /= 'multigrid' ) THEN
[215]	1269	message_string = 'non-cyclic lateral boundaries do not allow ' // &
	1270	'psolver = "' // TRIM( psolver ) // '"'
[226]	1271	CALL message( 'check_parameters', 'PA0051', 1, 2, 0, 6, 0 )
[1]	1272	ENDIF
[667]	1273	IF ( momentum_advec /= 'pw-scheme' .AND. &
	1274	momentum_advec /= 'ws-scheme') THEN
[215]	1275	message_string = 'non-cyclic lateral boundaries do not allow ' // &
	1276	'momentum_advec = "' // TRIM( momentum_advec ) // '"'
[226]	1277	CALL message( 'check_parameters', 'PA0052', 1, 2, 0, 6, 0 )
[1]	1278	ENDIF
[667]	1279	IF ( scalar_advec /= 'pw-scheme' .AND. &
	1280	scalar_advec /= 'ws-scheme' ) THEN
[215]	1281	message_string = 'non-cyclic lateral boundaries do not allow ' // &
	1282	'scalar_advec = "' // TRIM( scalar_advec ) // '"'
[226]	1283	CALL message( 'check_parameters', 'PA0053', 1, 2, 0, 6, 0 )
[1]	1284	ENDIF
[667]	1285	IF ( (scalar_advec == 'ws-scheme' .OR. momentum_advec == 'ws-scheme' ) &
	1286	.AND. loop_optimization == 'vector' ) THEN
	1287	message_string = 'non-cyclic lateral boundaries do not allow ' // &
	1288	'loop_optimization = vector and ' // &
	1289	'scalar_advec = "' // TRIM( scalar_advec ) // '"'
	1290	! The error message number still needs modification.
	1291	CALL message( 'check_parameters', 'PA0342', 1, 2, 0, 6, 0 )
	1292	END IF
[1]	1293	IF ( galilei_transformation ) THEN
[215]	1294	message_string = 'non-cyclic lateral boundaries do not allow ' // &
	1295	'galilei_transformation = .T.'
[226]	1296	CALL message( 'check_parameters', 'PA0054', 1, 2, 0, 6, 0 )
[1]	1297	ENDIF
	1298	ENDIF
	1299
	1300	!
	1301	!-- Bottom boundary condition for the turbulent Kinetic energy
	1302	IF ( bc_e_b == 'neumann' ) THEN
	1303	ibc_e_b = 1
	1304	IF ( adjust_mixing_length .AND. prandtl_layer ) THEN
[215]	1305	message_string = 'adjust_mixing_length = TRUE and bc_e_b = "neumann"'
[226]	1306	CALL message( 'check_parameters', 'PA0055', 0, 1, 0, 6, 0 )
[1]	1307	ENDIF
	1308	ELSEIF ( bc_e_b == '(u)*2+neumann' ) THEN
	1309	ibc_e_b = 2
	1310	IF ( .NOT. adjust_mixing_length .AND. prandtl_layer ) THEN
[215]	1311	message_string = 'adjust_mixing_length = FALSE and bc_e_b = "' // &
	1312	TRIM( bc_e_b ) // '"'
[226]	1313	CALL message( 'check_parameters', 'PA0056', 0, 1, 0, 6, 0 )
[1]	1314	ENDIF
	1315	IF ( .NOT. prandtl_layer ) THEN
	1316	bc_e_b = 'neumann'
	1317	ibc_e_b = 1
[215]	1318	message_string = 'boundary condition bc_e_b changed to "' // &
	1319	TRIM( bc_e_b ) // '"'
[226]	1320	CALL message( 'check_parameters', 'PA0057', 0, 1, 0, 6, 0 )
[1]	1321	ENDIF
	1322	ELSE
[215]	1323	message_string = 'unknown boundary condition: bc_e_b = "' // &
	1324	TRIM( bc_e_b ) // '"'
[226]	1325	CALL message( 'check_parameters', 'PA0058', 1, 2, 0, 6, 0 )
[1]	1326	ENDIF
	1327
	1328	!
	1329	!-- Boundary conditions for perturbation pressure
	1330	IF ( bc_p_b == 'dirichlet' ) THEN
	1331	ibc_p_b = 0
	1332	ELSEIF ( bc_p_b == 'neumann' ) THEN
	1333	ibc_p_b = 1
	1334	ELSEIF ( bc_p_b == 'neumann+inhomo' ) THEN
	1335	ibc_p_b = 2
	1336	ELSE
[215]	1337	message_string = 'unknown boundary condition: bc_p_b = "' // &
	1338	TRIM( bc_p_b ) // '"'
[226]	1339	CALL message( 'check_parameters', 'PA0059', 1, 2, 0, 6, 0 )
[1]	1340	ENDIF
	1341	IF ( ibc_p_b == 2 .AND. .NOT. prandtl_layer ) THEN
[215]	1342	message_string = 'boundary condition: bc_p_b = "' // TRIM( bc_p_b ) // &
	1343	'" not allowed with prandtl_layer = .FALSE.'
[226]	1344	CALL message( 'check_parameters', 'PA0060', 1, 2, 0, 6, 0 )
[1]	1345	ENDIF
	1346	IF ( bc_p_t == 'dirichlet' ) THEN
	1347	ibc_p_t = 0
	1348	ELSEIF ( bc_p_t == 'neumann' ) THEN
	1349	ibc_p_t = 1
	1350	ELSE
[215]	1351	message_string = 'unknown boundary condition: bc_p_t = "' // &
	1352	TRIM( bc_p_t ) // '"'
[226]	1353	CALL message( 'check_parameters', 'PA0061', 1, 2, 0, 6, 0 )
[1]	1354	ENDIF
	1355
	1356	!
	1357	!-- Boundary conditions for potential temperature
[102]	1358	IF ( coupling_mode == 'atmosphere_to_ocean' ) THEN
	1359	ibc_pt_b = 2
[1]	1360	ELSE
[102]	1361	IF ( bc_pt_b == 'dirichlet' ) THEN
	1362	ibc_pt_b = 0
	1363	ELSEIF ( bc_pt_b == 'neumann' ) THEN
	1364	ibc_pt_b = 1
	1365	ELSE
[215]	1366	message_string = 'unknown boundary condition: bc_pt_b = "' // &
	1367	TRIM( bc_pt_b ) // '"'
[226]	1368	CALL message( 'check_parameters', 'PA0062', 1, 2, 0, 6, 0 )
[1]	1369	ENDIF
	1370	ENDIF
[102]	1371
[1]	1372	IF ( bc_pt_t == 'dirichlet' ) THEN
	1373	ibc_pt_t = 0
	1374	ELSEIF ( bc_pt_t == 'neumann' ) THEN
	1375	ibc_pt_t = 1
[19]	1376	ELSEIF ( bc_pt_t == 'initial_gradient' ) THEN
	1377	ibc_pt_t = 2
[1]	1378	ELSE
[215]	1379	message_string = 'unknown boundary condition: bc_pt_t = "' // &
	1380	TRIM( bc_pt_t ) // '"'
[226]	1381	CALL message( 'check_parameters', 'PA0063', 1, 2, 0, 6, 0 )
[1]	1382	ENDIF
	1383
[20]	1384	IF ( surface_heatflux == 9999999.9 ) constant_heatflux = .FALSE.
	1385	IF ( top_heatflux == 9999999.9 ) constant_top_heatflux = .FALSE.
[103]	1386	IF ( top_momentumflux_u /= 9999999.9 .AND. &
	1387	top_momentumflux_v /= 9999999.9 ) THEN
	1388	constant_top_momentumflux = .TRUE.
	1389	ELSEIF ( .NOT. ( top_momentumflux_u == 9999999.9 .AND. &
[215]	1390	top_momentumflux_v == 9999999.9 ) ) THEN
	1391	message_string = 'both, top_momentumflux_u AND top_momentumflux_v ' // &
	1392	'must be set'
[226]	1393	CALL message( 'check_parameters', 'PA0064', 1, 2, 0, 6, 0 )
[103]	1394	ENDIF
[1]	1395
	1396	!
	1397	!-- A given surface temperature implies Dirichlet boundary condition for
	1398	!-- temperature. In this case specification of a constant heat flux is
	1399	!-- forbidden.
	1400	IF ( ibc_pt_b == 0 .AND. constant_heatflux .AND. &
	1401	surface_heatflux /= 0.0 ) THEN
[215]	1402	message_string = 'boundary_condition: bc_pt_b = "' // TRIM( bc_pt_b ) //&
	1403	'& is not allowed with constant_heatflux = .TRUE.'
[226]	1404	CALL message( 'check_parameters', 'PA0065', 1, 2, 0, 6, 0 )
[1]	1405	ENDIF
	1406	IF ( constant_heatflux .AND. pt_surface_initial_change /= 0.0 ) THEN
[215]	1407	WRITE ( message_string, * ) 'constant_heatflux = .TRUE. is not allo', &
	1408	'wed with pt_surface_initial_change (/=0) = ', &
	1409	pt_surface_initial_change
[226]	1410	CALL message( 'check_parameters', 'PA0066', 1, 2, 0, 6, 0 )
[1]	1411	ENDIF
	1412
	1413	!
[19]	1414	!-- A given temperature at the top implies Dirichlet boundary condition for
	1415	!-- temperature. In this case specification of a constant heat flux is
	1416	!-- forbidden.
	1417	IF ( ibc_pt_t == 0 .AND. constant_top_heatflux .AND. &
	1418	top_heatflux /= 0.0 ) THEN
[215]	1419	message_string = 'boundary_condition: bc_pt_t = "' // TRIM( bc_pt_t ) //&
	1420	'" is not allowed with constant_top_heatflux = .TRUE.'
[226]	1421	CALL message( 'check_parameters', 'PA0067', 1, 2, 0, 6, 0 )
[19]	1422	ENDIF
	1423
	1424	!
[95]	1425	!-- Boundary conditions for salinity
	1426	IF ( ocean ) THEN
	1427	IF ( bc_sa_t == 'dirichlet' ) THEN
	1428	ibc_sa_t = 0
	1429	ELSEIF ( bc_sa_t == 'neumann' ) THEN
	1430	ibc_sa_t = 1
	1431	ELSE
[215]	1432	message_string = 'unknown boundary condition: bc_sa_t = "' // &
	1433	TRIM( bc_sa_t ) // '"'
[226]	1434	CALL message( 'check_parameters', 'PA0068', 1, 2, 0, 6, 0 )
[95]	1435	ENDIF
	1436
	1437	IF ( top_salinityflux == 9999999.9 ) constant_top_salinityflux = .FALSE.
[97]	1438	IF ( ibc_sa_t == 1 .AND. top_salinityflux == 9999999.9 ) THEN
[215]	1439	message_string = 'boundary condition: bc_sa_t = "' // &
	1440	TRIM( bc_sa_t ) // '" requires to set ' // &
	1441	'top_salinityflux'
[226]	1442	CALL message( 'check_parameters', 'PA0069', 1, 2, 0, 6, 0 )
[97]	1443	ENDIF
[95]	1444
	1445	!
	1446	!-- A fixed salinity at the top implies Dirichlet boundary condition for
	1447	!-- salinity. In this case specification of a constant salinity flux is
	1448	!-- forbidden.
	1449	IF ( ibc_sa_t == 0 .AND. constant_top_salinityflux .AND. &
	1450	top_salinityflux /= 0.0 ) THEN
[215]	1451	message_string = 'boundary condition: bc_sa_t = "' // &
	1452	TRIM( bc_sa_t ) // '" is not allowed with ' // &
	1453	'constant_top_salinityflux = .TRUE.'
[226]	1454	CALL message( 'check_parameters', 'PA0070', 1, 2, 0, 6, 0 )
[95]	1455	ENDIF
	1456
	1457	ENDIF
	1458
	1459	!
[75]	1460	!-- In case of humidity or passive scalar, set boundary conditions for total
[1]	1461	!-- water content / scalar
[75]	1462	IF ( humidity .OR. passive_scalar ) THEN
	1463	IF ( humidity ) THEN
[1]	1464	sq = 'q'
	1465	ELSE
	1466	sq = 's'
	1467	ENDIF
	1468	IF ( bc_q_b == 'dirichlet' ) THEN
	1469	ibc_q_b = 0
	1470	ELSEIF ( bc_q_b == 'neumann' ) THEN
	1471	ibc_q_b = 1
	1472	ELSE
[215]	1473	message_string = 'unknown boundary condition: bc_' // TRIM( sq ) // &
	1474	'_b ="' // TRIM( bc_q_b ) // '"'
[226]	1475	CALL message( 'check_parameters', 'PA0071', 1, 2, 0, 6, 0 )
[1]	1476	ENDIF
	1477	IF ( bc_q_t == 'dirichlet' ) THEN
	1478	ibc_q_t = 0
	1479	ELSEIF ( bc_q_t == 'neumann' ) THEN
	1480	ibc_q_t = 1
	1481	ELSE
[215]	1482	message_string = 'unknown boundary condition: bc_' // TRIM( sq ) // &
	1483	'_t ="' // TRIM( bc_q_t ) // '"'
[226]	1484	CALL message( 'check_parameters', 'PA0072', 1, 2, 0, 6, 0 )
[1]	1485	ENDIF
	1486
[600]	1487	IF ( surface_waterflux == 9999999.9 ) constant_waterflux = .FALSE.
[1]	1488
	1489	!
	1490	!-- A given surface humidity implies Dirichlet boundary condition for
[75]	1491	!-- humidity. In this case specification of a constant water flux is
[1]	1492	!-- forbidden.
	1493	IF ( ibc_q_b == 0 .AND. constant_waterflux ) THEN
[215]	1494	message_string = 'boundary condition: bc_' // TRIM( sq ) // '_b ' // &
	1495	'= "' // TRIM( bc_q_b ) // '" is not allowed wi' // &
	1496	'th prescribed surface flux'
[226]	1497	CALL message( 'check_parameters', 'PA0073', 1, 2, 0, 6, 0 )
[1]	1498	ENDIF
	1499	IF ( constant_waterflux .AND. q_surface_initial_change /= 0.0 ) THEN
[215]	1500	WRITE( message_string, * ) 'a prescribed surface flux is not allo', &
	1501	'wed with ', sq, '_surface_initial_change (/=0) = ', &
	1502	q_surface_initial_change
[226]	1503	CALL message( 'check_parameters', 'PA0074', 1, 2, 0, 6, 0 )
[1]	1504	ENDIF
	1505
	1506	ENDIF
	1507
	1508	!
	1509	!-- Boundary conditions for horizontal components of wind speed
	1510	IF ( bc_uv_b == 'dirichlet' ) THEN
	1511	ibc_uv_b = 0
	1512	ELSEIF ( bc_uv_b == 'neumann' ) THEN
	1513	ibc_uv_b = 1
	1514	IF ( prandtl_layer ) THEN
[215]	1515	message_string = 'boundary condition: bc_uv_b = "' // &
	1516	TRIM( bc_uv_b ) // '" is not allowed with prandtl_layer = .TRUE.'
[226]	1517	CALL message( 'check_parameters', 'PA0075', 1, 2, 0, 6, 0 )
[1]	1518	ENDIF
	1519	ELSE
[215]	1520	message_string = 'unknown boundary condition: bc_uv_b = "' // &
	1521	TRIM( bc_uv_b ) // '"'
[226]	1522	CALL message( 'check_parameters', 'PA0076', 1, 2, 0, 6, 0 )
[1]	1523	ENDIF
[667]	1524	!
	1525	!-- In case of coupled simulations u and v at the ground in atmosphere will be
	1526	!-- assigned with the u and v values of the ocean surface
	1527	IF ( coupling_mode == 'atmosphere_to_ocean' ) THEN
	1528	ibc_uv_b = 2
	1529	ENDIF
[215]	1530
[108]	1531	IF ( coupling_mode == 'ocean_to_atmosphere' ) THEN
	1532	bc_uv_t = 'neumann'
[1]	1533	ibc_uv_t = 1
	1534	ELSE
[132]	1535	IF ( bc_uv_t == 'dirichlet' .OR. bc_uv_t == 'dirichlet_0' ) THEN
[108]	1536	ibc_uv_t = 0
	1537	ELSEIF ( bc_uv_t == 'neumann' ) THEN
	1538	ibc_uv_t = 1
	1539	ELSE
[215]	1540	message_string = 'unknown boundary condition: bc_uv_t = "' // &
	1541	TRIM( bc_uv_t ) // '"'
[226]	1542	CALL message( 'check_parameters', 'PA0077', 1, 2, 0, 6, 0 )
[1]	1543	ENDIF
	1544	ENDIF
	1545
	1546	!
	1547	!-- Compute and check, respectively, the Rayleigh Damping parameter
	1548	IF ( rayleigh_damping_factor == -1.0 ) THEN
	1549	IF ( momentum_advec == 'ups-scheme' ) THEN
	1550	rayleigh_damping_factor = 0.01
	1551	ELSE
	1552	rayleigh_damping_factor = 0.0
	1553	ENDIF
	1554	ELSE
	1555	IF ( rayleigh_damping_factor < 0.0 .OR. rayleigh_damping_factor > 1.0 ) &
	1556	THEN
[215]	1557	WRITE( message_string, * ) 'rayleigh_damping_factor = ', &
	1558	rayleigh_damping_factor, ' out of range [0.0,1.0]'
[226]	1559	CALL message( 'check_parameters', 'PA0078', 1, 2, 0, 6, 0 )
[1]	1560	ENDIF
	1561	ENDIF
	1562
	1563	IF ( rayleigh_damping_height == -1.0 ) THEN
[108]	1564	IF ( .NOT. ocean ) THEN
	1565	rayleigh_damping_height = 0.66666666666 * zu(nzt)
	1566	ELSE
	1567	rayleigh_damping_height = 0.66666666666 * zu(nzb)
	1568	ENDIF
[1]	1569	ELSE
[108]	1570	IF ( .NOT. ocean ) THEN
	1571	IF ( rayleigh_damping_height < 0.0 .OR. &
	1572	rayleigh_damping_height > zu(nzt) ) THEN
[215]	1573	WRITE( message_string, * ) 'rayleigh_damping_height = ', &
	1574	rayleigh_damping_height, ' out of range [0.0,', zu(nzt), ']'
[226]	1575	CALL message( 'check_parameters', 'PA0079', 1, 2, 0, 6, 0 )
[1]	1576	ENDIF
[108]	1577	ELSE
	1578	IF ( rayleigh_damping_height > 0.0 .OR. &
	1579	rayleigh_damping_height < zu(nzb) ) THEN
[215]	1580	WRITE( message_string, * ) 'rayleigh_damping_height = ', &
	1581	rayleigh_damping_height, ' out of range [0.0,', zu(nzb), ']'
[226]	1582	CALL message( 'check_parameters', 'PA0079', 1, 2, 0, 6, 0 )
[108]	1583	ENDIF
[1]	1584	ENDIF
	1585	ENDIF
	1586
	1587	!
	1588	!-- Check limiters for Upstream-Spline scheme
	1589	IF ( overshoot_limit_u < 0.0 .OR. overshoot_limit_v < 0.0 .OR. &
	1590	overshoot_limit_w < 0.0 .OR. overshoot_limit_pt < 0.0 .OR. &
	1591	overshoot_limit_e < 0.0 ) THEN
[215]	1592	message_string = 'overshoot_limit_... < 0.0 is not allowed'
[226]	1593	CALL message( 'check_parameters', 'PA0080', 1, 2, 0, 6, 0 )
[1]	1594	ENDIF
	1595	IF ( ups_limit_u < 0.0 .OR. ups_limit_v < 0.0 .OR. ups_limit_w < 0.0 .OR. &
	1596	ups_limit_pt < 0.0 .OR. ups_limit_e < 0.0 ) THEN
[215]	1597	message_string = 'ups_limit_... < 0.0 is not allowed'
[226]	1598	CALL message( 'check_parameters', 'PA0081', 1, 2, 0, 6, 0 )
[1]	1599	ENDIF
	1600
	1601	!
	1602	!-- Check number of chosen statistic regions. More than 10 regions are not
	1603	!-- allowed, because so far no more than 10 corresponding output files can
	1604	!-- be opened (cf. check_open)
	1605	IF ( statistic_regions > 9 .OR. statistic_regions < 0 ) THEN
[215]	1606	WRITE ( message_string, * ) 'number of statistic_regions = ', &
	1607	statistic_regions+1, ' but only 10 regions are allowed'
[226]	1608	CALL message( 'check_parameters', 'PA0082', 1, 2, 0, 6, 0 )
[1]	1609	ENDIF
	1610	IF ( normalizing_region > statistic_regions .OR. &
	1611	normalizing_region < 0) THEN
[215]	1612	WRITE ( message_string, * ) 'normalizing_region = ', &
	1613	normalizing_region, ' must be >= 0 and <= ',statistic_regions, &
	1614	' (value of statistic_regions)'
[226]	1615	CALL message( 'check_parameters', 'PA0083', 1, 2, 0, 6, 0 )
[1]	1616	ENDIF
	1617
	1618	!
[116]	1619	!-- Check the interval for sorting particles.
	1620	!-- Using particles as cloud droplets requires sorting after each timestep.
	1621	IF ( dt_sort_particles /= 0.0 .AND. cloud_droplets ) THEN
	1622	dt_sort_particles = 0.0
[215]	1623	message_string = 'dt_sort_particles is reset to 0.0 because of cloud' //&
	1624	'_droplets = .TRUE.'
[226]	1625	CALL message( 'check_parameters', 'PA0084', 0, 1, 0, 6, 0 )
[116]	1626	ENDIF
	1627
	1628	!
[1]	1629	!-- Set the default intervals for data output, if necessary
	1630	!-- NOTE: dt_dosp has already been set in package_parin
	1631	IF ( dt_data_output /= 9999999.9 ) THEN
	1632	IF ( dt_dopr == 9999999.9 ) dt_dopr = dt_data_output
	1633	IF ( dt_dopts == 9999999.9 ) dt_dopts = dt_data_output
	1634	IF ( dt_do2d_xy == 9999999.9 ) dt_do2d_xy = dt_data_output
	1635	IF ( dt_do2d_xz == 9999999.9 ) dt_do2d_xz = dt_data_output
	1636	IF ( dt_do2d_yz == 9999999.9 ) dt_do2d_yz = dt_data_output
	1637	IF ( dt_do3d == 9999999.9 ) dt_do3d = dt_data_output
	1638	IF ( dt_data_output_av == 9999999.9 ) dt_data_output_av = dt_data_output
[564]	1639	DO mid = 1, max_masks
[410]	1640	IF ( dt_domask(mid) == 9999999.9 ) dt_domask(mid) = dt_data_output
	1641	ENDDO
[1]	1642	ENDIF
	1643
	1644	!
	1645	!-- Set the default skip time intervals for data output, if necessary
	1646	IF ( skip_time_dopr == 9999999.9 ) &
	1647	skip_time_dopr = skip_time_data_output
	1648	IF ( skip_time_dosp == 9999999.9 ) &
	1649	skip_time_dosp = skip_time_data_output
	1650	IF ( skip_time_do2d_xy == 9999999.9 ) &
	1651	skip_time_do2d_xy = skip_time_data_output
	1652	IF ( skip_time_do2d_xz == 9999999.9 ) &
	1653	skip_time_do2d_xz = skip_time_data_output
	1654	IF ( skip_time_do2d_yz == 9999999.9 ) &
	1655	skip_time_do2d_yz = skip_time_data_output
	1656	IF ( skip_time_do3d == 9999999.9 ) &
	1657	skip_time_do3d = skip_time_data_output
	1658	IF ( skip_time_data_output_av == 9999999.9 ) &
	1659	skip_time_data_output_av = skip_time_data_output
[564]	1660	DO mid = 1, max_masks
[410]	1661	IF ( skip_time_domask(mid) == 9999999.9 ) &
	1662	skip_time_domask(mid) = skip_time_data_output
	1663	ENDDO
[1]	1664
	1665	!
	1666	!-- Check the average intervals (first for 3d-data, then for profiles and
	1667	!-- spectra)
	1668	IF ( averaging_interval > dt_data_output_av ) THEN
[215]	1669	WRITE( message_string, * ) 'averaging_interval = ', &
	1670	averaging_interval, ' must be <= dt_data_output = ', dt_data_output
[226]	1671	CALL message( 'check_parameters', 'PA0085', 1, 2, 0, 6, 0 )
[1]	1672	ENDIF
	1673
	1674	IF ( averaging_interval_pr == 9999999.9 ) THEN
	1675	averaging_interval_pr = averaging_interval
	1676	ENDIF
	1677
	1678	IF ( averaging_interval_pr > dt_dopr ) THEN
[215]	1679	WRITE( message_string, * ) 'averaging_interval_pr = ', &
	1680	averaging_interval_pr, ' must be <= dt_dopr = ', dt_dopr
[226]	1681	CALL message( 'check_parameters', 'PA0086', 1, 2, 0, 6, 0 )
[1]	1682	ENDIF
	1683
	1684	IF ( averaging_interval_sp == 9999999.9 ) THEN
	1685	averaging_interval_sp = averaging_interval
	1686	ENDIF
	1687
	1688	IF ( averaging_interval_sp > dt_dosp ) THEN
[215]	1689	WRITE( message_string, * ) 'averaging_interval_sp = ', &
	1690	averaging_interval_sp, ' must be <= dt_dosp = ', dt_dosp
[226]	1691	CALL message( 'check_parameters', 'PA0087', 1, 2, 0, 6, 0 )
[1]	1692	ENDIF
	1693
	1694	!
	1695	!-- Set the default interval for profiles entering the temporal average
	1696	IF ( dt_averaging_input_pr == 9999999.9 ) THEN
	1697	dt_averaging_input_pr = dt_averaging_input
	1698	ENDIF
	1699
	1700	!
	1701	!-- Set the default interval for the output of timeseries to a reasonable
	1702	!-- value (tries to minimize the number of calls of flow_statistics)
	1703	IF ( dt_dots == 9999999.9 ) THEN
	1704	IF ( averaging_interval_pr == 0.0 ) THEN
	1705	dt_dots = MIN( dt_run_control, dt_dopr )
	1706	ELSE
	1707	dt_dots = MIN( dt_run_control, dt_averaging_input_pr )
	1708	ENDIF
	1709	ENDIF
	1710
	1711	!
	1712	!-- Check the sample rate for averaging (first for 3d-data, then for profiles)
	1713	IF ( dt_averaging_input > averaging_interval ) THEN
[215]	1714	WRITE( message_string, * ) 'dt_averaging_input = ', &
	1715	dt_averaging_input, ' must be <= averaging_interval = ', &
	1716	averaging_interval
[226]	1717	CALL message( 'check_parameters', 'PA0088', 1, 2, 0, 6, 0 )
[1]	1718	ENDIF
	1719
	1720	IF ( dt_averaging_input_pr > averaging_interval_pr ) THEN
[215]	1721	WRITE( message_string, * ) 'dt_averaging_input_pr = ', &
	1722	dt_averaging_input_pr, ' must be <= averaging_interval_pr = ', &
	1723	averaging_interval_pr
[226]	1724	CALL message( 'check_parameters', 'PA0089', 1, 2, 0, 6, 0 )
[1]	1725	ENDIF
	1726
	1727	!
[72]	1728	!-- Set the default value for the integration interval of precipitation amount
	1729	IF ( precipitation ) THEN
	1730	IF ( precipitation_amount_interval == 9999999.9 ) THEN
	1731	precipitation_amount_interval = dt_do2d_xy
	1732	ELSE
	1733	IF ( precipitation_amount_interval > dt_do2d_xy ) THEN
[215]	1734	WRITE( message_string, * ) 'precipitation_amount_interval = ', &
	1735	precipitation_amount_interval, ' must not be larger than ', &
	1736	'dt_do2d_xy = ', dt_do2d_xy
[226]	1737	CALL message( 'check_parameters', 'PA0090', 1, 2, 0, 6, 0 )
[72]	1738	ENDIF
	1739	ENDIF
	1740	ENDIF
	1741
	1742	!
[1]	1743	!-- Determine the number of output profiles and check whether they are
	1744	!-- permissible
	1745	DO WHILE ( data_output_pr(dopr_n+1) /= ' ' )
	1746
	1747	dopr_n = dopr_n + 1
	1748	i = dopr_n
	1749
	1750	!
	1751	!-- Determine internal profile number (for hom, homs)
	1752	!-- and store height levels
	1753	SELECT CASE ( TRIM( data_output_pr(i) ) )
	1754
	1755	CASE ( 'u', '#u' )
	1756	dopr_index(i) = 1
[87]	1757	dopr_unit(i) = 'm/s'
[1]	1758	hom(:,2,1,:) = SPREAD( zu, 2, statistic_regions+1 )
	1759	IF ( data_output_pr(i)(1:1) == '#' ) THEN
	1760	dopr_initial_index(i) = 5
	1761	hom(:,2,5,:) = SPREAD( zu, 2, statistic_regions+1 )
	1762	data_output_pr(i) = data_output_pr(i)(2:)
	1763	ENDIF
	1764
	1765	CASE ( 'v', '#v' )
	1766	dopr_index(i) = 2
[87]	1767	dopr_unit(i) = 'm/s'
	1768	hom(:,2,2,:) = SPREAD( zu, 2, statistic_regions+1 )
[1]	1769	IF ( data_output_pr(i)(1:1) == '#' ) THEN
	1770	dopr_initial_index(i) = 6
	1771	hom(:,2,6,:) = SPREAD( zu, 2, statistic_regions+1 )
	1772	data_output_pr(i) = data_output_pr(i)(2:)
	1773	ENDIF
	1774
	1775	CASE ( 'w' )
	1776	dopr_index(i) = 3
[87]	1777	dopr_unit(i) = 'm/s'
	1778	hom(:,2,3,:) = SPREAD( zw, 2, statistic_regions+1 )
[1]	1779
	1780	CASE ( 'pt', '#pt' )
	1781	IF ( .NOT. cloud_physics ) THEN
	1782	dopr_index(i) = 4
[87]	1783	dopr_unit(i) = 'K'
[1]	1784	hom(:,2,4,:) = SPREAD( zu, 2, statistic_regions+1 )
	1785	IF ( data_output_pr(i)(1:1) == '#' ) THEN
	1786	dopr_initial_index(i) = 7
	1787	hom(:,2,7,:) = SPREAD( zu, 2, statistic_regions+1 )
[87]	1788	hom(nzb,2,7,:) = 0.0 ! because zu(nzb) is negative
[1]	1789	data_output_pr(i) = data_output_pr(i)(2:)
	1790	ENDIF
	1791	ELSE
	1792	dopr_index(i) = 43
[87]	1793	dopr_unit(i) = 'K'
[1]	1794	hom(:,2,43,:) = SPREAD( zu, 2, statistic_regions+1 )
	1795	IF ( data_output_pr(i)(1:1) == '#' ) THEN
	1796	dopr_initial_index(i) = 28
	1797	hom(:,2,28,:) = SPREAD( zu, 2, statistic_regions+1 )
[87]	1798	hom(nzb,2,28,:) = 0.0 ! because zu(nzb) is negative
[1]	1799	data_output_pr(i) = data_output_pr(i)(2:)
	1800	ENDIF
	1801	ENDIF
	1802
	1803	CASE ( 'e' )
	1804	dopr_index(i) = 8
[87]	1805	dopr_unit(i) = 'm2/s2'
[1]	1806	hom(:,2,8,:) = SPREAD( zu, 2, statistic_regions+1 )
	1807	hom(nzb,2,8,:) = 0.0
	1808
	1809	CASE ( 'km', '#km' )
	1810	dopr_index(i) = 9
[87]	1811	dopr_unit(i) = 'm2/s'
[1]	1812	hom(:,2,9,:) = SPREAD( zu, 2, statistic_regions+1 )
	1813	hom(nzb,2,9,:) = 0.0
	1814	IF ( data_output_pr(i)(1:1) == '#' ) THEN
	1815	dopr_initial_index(i) = 23
	1816	hom(:,2,23,:) = hom(:,2,9,:)
	1817	data_output_pr(i) = data_output_pr(i)(2:)
	1818	ENDIF
	1819
	1820	CASE ( 'kh', '#kh' )
	1821	dopr_index(i) = 10
[87]	1822	dopr_unit(i) = 'm2/s'
[1]	1823	hom(:,2,10,:) = SPREAD( zu, 2, statistic_regions+1 )
	1824	hom(nzb,2,10,:) = 0.0
	1825	IF ( data_output_pr(i)(1:1) == '#' ) THEN
	1826	dopr_initial_index(i) = 24
	1827	hom(:,2,24,:) = hom(:,2,10,:)
	1828	data_output_pr(i) = data_output_pr(i)(2:)
	1829	ENDIF
	1830
	1831	CASE ( 'l', '#l' )
	1832	dopr_index(i) = 11
[87]	1833	dopr_unit(i) = 'm'
[1]	1834	hom(:,2,11,:) = SPREAD( zu, 2, statistic_regions+1 )
	1835	hom(nzb,2,11,:) = 0.0
	1836	IF ( data_output_pr(i)(1:1) == '#' ) THEN
	1837	dopr_initial_index(i) = 25
	1838	hom(:,2,25,:) = hom(:,2,11,:)
	1839	data_output_pr(i) = data_output_pr(i)(2:)
	1840	ENDIF
	1841
	1842	CASE ( 'w"u"' )
	1843	dopr_index(i) = 12
[87]	1844	dopr_unit(i) = 'm2/s2'
[1]	1845	hom(:,2,12,:) = SPREAD( zw, 2, statistic_regions+1 )
	1846	IF ( prandtl_layer ) hom(nzb,2,12,:) = zu(1)
	1847
	1848	CASE ( 'wu' )
	1849	dopr_index(i) = 13
[87]	1850	dopr_unit(i) = 'm2/s2'
[1]	1851	hom(:,2,13,:) = SPREAD( zw, 2, statistic_regions+1 )
	1852
	1853	CASE ( 'w"v"' )
	1854	dopr_index(i) = 14
[87]	1855	dopr_unit(i) = 'm2/s2'
[1]	1856	hom(:,2,14,:) = SPREAD( zw, 2, statistic_regions+1 )
	1857	IF ( prandtl_layer ) hom(nzb,2,14,:) = zu(1)
	1858
	1859	CASE ( 'wv' )
	1860	dopr_index(i) = 15
[87]	1861	dopr_unit(i) = 'm2/s2'
[1]	1862	hom(:,2,15,:) = SPREAD( zw, 2, statistic_regions+1 )
	1863
	1864	CASE ( 'w"pt"' )
	1865	dopr_index(i) = 16
[87]	1866	dopr_unit(i) = 'K m/s'
[1]	1867	hom(:,2,16,:) = SPREAD( zw, 2, statistic_regions+1 )
	1868
	1869	CASE ( 'wpt' )
	1870	dopr_index(i) = 17
[87]	1871	dopr_unit(i) = 'K m/s'
[1]	1872	hom(:,2,17,:) = SPREAD( zw, 2, statistic_regions+1 )
	1873
	1874	CASE ( 'wpt' )
	1875	dopr_index(i) = 18
[87]	1876	dopr_unit(i) = 'K m/s'
[1]	1877	hom(:,2,18,:) = SPREAD( zw, 2, statistic_regions+1 )
	1878
	1879	CASE ( 'wu' )
	1880	dopr_index(i) = 19
[87]	1881	dopr_unit(i) = 'm2/s2'
[1]	1882	hom(:,2,19,:) = SPREAD( zw, 2, statistic_regions+1 )
	1883	IF ( prandtl_layer ) hom(nzb,2,19,:) = zu(1)
	1884
	1885	CASE ( 'wv' )
	1886	dopr_index(i) = 20
[87]	1887	dopr_unit(i) = 'm2/s2'
[1]	1888	hom(:,2,20,:) = SPREAD( zw, 2, statistic_regions+1 )
	1889	IF ( prandtl_layer ) hom(nzb,2,20,:) = zu(1)
	1890
	1891	CASE ( 'wptBC' )
	1892	dopr_index(i) = 21
[87]	1893	dopr_unit(i) = 'K m/s'
[1]	1894	hom(:,2,21,:) = SPREAD( zw, 2, statistic_regions+1 )
	1895
	1896	CASE ( 'wptBC' )
	1897	dopr_index(i) = 22
[87]	1898	dopr_unit(i) = 'K m/s'
[1]	1899	hom(:,2,22,:) = SPREAD( zw, 2, statistic_regions+1 )
	1900
[96]	1901	CASE ( 'sa', '#sa' )
	1902	IF ( .NOT. ocean ) THEN
[215]	1903	message_string = 'data_output_pr = ' // &
	1904	TRIM( data_output_pr(i) ) // ' is not imp' // &
	1905	'lemented for ocean = .FALSE.'
[226]	1906	CALL message( 'check_parameters', 'PA0091', 1, 2, 0, 6, 0 )
[96]	1907	ELSE
	1908	dopr_index(i) = 23
	1909	dopr_unit(i) = 'psu'
	1910	hom(:,2,23,:) = SPREAD( zu, 2, statistic_regions+1 )
	1911	IF ( data_output_pr(i)(1:1) == '#' ) THEN
	1912	dopr_initial_index(i) = 26
	1913	hom(:,2,26,:) = SPREAD( zu, 2, statistic_regions+1 )
	1914	hom(nzb,2,26,:) = 0.0 ! weil zu(nzb) negativ ist
	1915	data_output_pr(i) = data_output_pr(i)(2:)
	1916	ENDIF
	1917	ENDIF
	1918
[1]	1919	CASE ( 'u*2' )
	1920	dopr_index(i) = 30
[87]	1921	dopr_unit(i) = 'm2/s2'
[1]	1922	hom(:,2,30,:) = SPREAD( zu, 2, statistic_regions+1 )
	1923
	1924	CASE ( 'v*2' )
	1925	dopr_index(i) = 31
[87]	1926	dopr_unit(i) = 'm2/s2'
[1]	1927	hom(:,2,31,:) = SPREAD( zu, 2, statistic_regions+1 )
	1928
	1929	CASE ( 'w*2' )
	1930	dopr_index(i) = 32
[87]	1931	dopr_unit(i) = 'm2/s2'
[1]	1932	hom(:,2,32,:) = SPREAD( zw, 2, statistic_regions+1 )
	1933
	1934	CASE ( 'pt*2' )
	1935	dopr_index(i) = 33
[87]	1936	dopr_unit(i) = 'K2'
[1]	1937	hom(:,2,33,:) = SPREAD( zu, 2, statistic_regions+1 )
	1938
	1939	CASE ( 'e*' )
	1940	dopr_index(i) = 34
[87]	1941	dopr_unit(i) = 'm2/s2'
[1]	1942	hom(:,2,34,:) = SPREAD( zu, 2, statistic_regions+1 )
	1943
	1944	CASE ( 'w2pt' )
	1945	dopr_index(i) = 35
[87]	1946	dopr_unit(i) = 'K m2/s2'
[1]	1947	hom(:,2,35,:) = SPREAD( zw, 2, statistic_regions+1 )
	1948
	1949	CASE ( 'wpt2' )
	1950	dopr_index(i) = 36
[87]	1951	dopr_unit(i) = 'K2 m/s'
[1]	1952	hom(:,2,36,:) = SPREAD( zw, 2, statistic_regions+1 )
	1953
	1954	CASE ( 'we' )
	1955	dopr_index(i) = 37
[87]	1956	dopr_unit(i) = 'm3/s3'
[1]	1957	hom(:,2,37,:) = SPREAD( zw, 2, statistic_regions+1 )
	1958
	1959	CASE ( 'w*3' )
	1960	dopr_index(i) = 38
[87]	1961	dopr_unit(i) = 'm3/s3'
[1]	1962	hom(:,2,38,:) = SPREAD( zw, 2, statistic_regions+1 )
	1963
	1964	CASE ( 'Sw' )
	1965	dopr_index(i) = 39
[89]	1966	dopr_unit(i) = 'none'
[1]	1967	hom(:,2,39,:) = SPREAD( zw, 2, statistic_regions+1 )
	1968
[232]	1969	CASE ( 'p' )
	1970	dopr_index(i) = 40
	1971	dopr_unit(i) = 'Pa'
	1972	hom(:,2,40,:) = SPREAD( zu, 2, statistic_regions+1 )
	1973
[1]	1974	CASE ( 'q', '#q' )
[108]	1975	IF ( .NOT. humidity ) THEN
[215]	1976	message_string = 'data_output_pr = ' // &
	1977	TRIM( data_output_pr(i) ) // ' is not imp' // &
	1978	'lemented for humidity = .FALSE.'
[226]	1979	CALL message( 'check_parameters', 'PA0092', 1, 2, 0, 6, 0 )
[1]	1980	ELSE
	1981	dopr_index(i) = 41
[87]	1982	dopr_unit(i) = 'kg/kg'
	1983	hom(:,2,41,:) = SPREAD( zu, 2, statistic_regions+1 )
[1]	1984	IF ( data_output_pr(i)(1:1) == '#' ) THEN
	1985	dopr_initial_index(i) = 26
	1986	hom(:,2,26,:) = SPREAD( zu, 2, statistic_regions+1 )
	1987	hom(nzb,2,26,:) = 0.0 ! weil zu(nzb) negativ ist
	1988	data_output_pr(i) = data_output_pr(i)(2:)
	1989	ENDIF
	1990	ENDIF
	1991
	1992	CASE ( 's', '#s' )
	1993	IF ( .NOT. passive_scalar ) THEN
[215]	1994	message_string = 'data_output_pr = ' // &
	1995	TRIM( data_output_pr(i) ) // ' is not imp' // &
	1996	'lemented for passive_scalar = .FALSE.'
[226]	1997	CALL message( 'check_parameters', 'PA0093', 1, 2, 0, 6, 0 )
[1]	1998	ELSE
	1999	dopr_index(i) = 41
[87]	2000	dopr_unit(i) = 'kg/m3'
	2001	hom(:,2,41,:) = SPREAD( zu, 2, statistic_regions+1 )
[1]	2002	IF ( data_output_pr(i)(1:1) == '#' ) THEN
	2003	dopr_initial_index(i) = 26
	2004	hom(:,2,26,:) = SPREAD( zu, 2, statistic_regions+1 )
	2005	hom(nzb,2,26,:) = 0.0 ! weil zu(nzb) negativ ist
	2006	data_output_pr(i) = data_output_pr(i)(2:)
	2007	ENDIF
	2008	ENDIF
	2009
	2010	CASE ( 'qv', '#qv' )
	2011	IF ( .NOT. cloud_physics ) THEN
	2012	dopr_index(i) = 41
[87]	2013	dopr_unit(i) = 'kg/kg'
	2014	hom(:,2,41,:) = SPREAD( zu, 2, statistic_regions+1 )
[1]	2015	IF ( data_output_pr(i)(1:1) == '#' ) THEN
	2016	dopr_initial_index(i) = 26
	2017	hom(:,2,26,:) = SPREAD( zu, 2, statistic_regions+1 )
	2018	hom(nzb,2,26,:) = 0.0 ! weil zu(nzb) negativ ist
	2019	data_output_pr(i) = data_output_pr(i)(2:)
	2020	ENDIF
	2021	ELSE
	2022	dopr_index(i) = 42
[87]	2023	dopr_unit(i) = 'kg/kg'
	2024	hom(:,2,42,:) = SPREAD( zu, 2, statistic_regions+1 )
[1]	2025	IF ( data_output_pr(i)(1:1) == '#' ) THEN
	2026	dopr_initial_index(i) = 27
	2027	hom(:,2,27,:) = SPREAD( zu, 2, statistic_regions+1 )
	2028	hom(nzb,2,27,:) = 0.0 ! weil zu(nzb) negativ ist
	2029	data_output_pr(i) = data_output_pr(i)(2:)
	2030	ENDIF
	2031	ENDIF
	2032
	2033	CASE ( 'lpt', '#lpt' )
	2034	IF ( .NOT. cloud_physics ) THEN
[215]	2035	message_string = 'data_output_pr = ' // &
	2036	TRIM( data_output_pr(i) ) // ' is not imp' // &
	2037	'lemented for cloud_physics = .FALSE.'
[226]	2038	CALL message( 'check_parameters', 'PA0094', 1, 2, 0, 6, 0 )
[1]	2039	ELSE
	2040	dopr_index(i) = 4
[87]	2041	dopr_unit(i) = 'K'
[1]	2042	hom(:,2,4,:) = SPREAD( zu, 2, statistic_regions+1 )
	2043	IF ( data_output_pr(i)(1:1) == '#' ) THEN
	2044	dopr_initial_index(i) = 7
	2045	hom(:,2,7,:) = SPREAD( zu, 2, statistic_regions+1 )
	2046	hom(nzb,2,7,:) = 0.0 ! weil zu(nzb) negativ ist
	2047	data_output_pr(i) = data_output_pr(i)(2:)
	2048	ENDIF
	2049	ENDIF
	2050
	2051	CASE ( 'vpt', '#vpt' )
	2052	dopr_index(i) = 44
[87]	2053	dopr_unit(i) = 'K'
	2054	hom(:,2,44,:) = SPREAD( zu, 2, statistic_regions+1 )
[1]	2055	IF ( data_output_pr(i)(1:1) == '#' ) THEN
	2056	dopr_initial_index(i) = 29
	2057	hom(:,2,29,:) = SPREAD( zu, 2, statistic_regions+1 )
	2058	hom(nzb,2,29,:) = 0.0 ! weil zu(nzb) negativ ist
	2059	data_output_pr(i) = data_output_pr(i)(2:)
	2060	ENDIF
	2061
	2062	CASE ( 'w"vpt"' )
	2063	dopr_index(i) = 45
[87]	2064	dopr_unit(i) = 'K m/s'
[1]	2065	hom(:,2,45,:) = SPREAD( zw, 2, statistic_regions+1 )
	2066
	2067	CASE ( 'wvpt' )
	2068	dopr_index(i) = 46
[87]	2069	dopr_unit(i) = 'K m/s'
[1]	2070	hom(:,2,46,:) = SPREAD( zw, 2, statistic_regions+1 )
	2071
	2072	CASE ( 'wvpt' )
	2073	dopr_index(i) = 47
[87]	2074	dopr_unit(i) = 'K m/s'
[1]	2075	hom(:,2,47,:) = SPREAD( zw, 2, statistic_regions+1 )
	2076
	2077	CASE ( 'w"q"' )
[108]	2078	IF ( .NOT. humidity ) THEN
[215]	2079	message_string = 'data_output_pr = ' // &
	2080	TRIM( data_output_pr(i) ) // ' is not imp' // &
	2081	'lemented for humidity = .FALSE.'
[226]	2082	CALL message( 'check_parameters', 'PA0092', 1, 2, 0, 6, 0 )
[1]	2083	ELSE
	2084	dopr_index(i) = 48
[87]	2085	dopr_unit(i) = 'kg/kg m/s'
[1]	2086	hom(:,2,48,:) = SPREAD( zw, 2, statistic_regions+1 )
	2087	ENDIF
	2088
	2089	CASE ( 'wq' )
[108]	2090	IF ( .NOT. humidity ) THEN
[215]	2091	message_string = 'data_output_pr = ' // &
	2092	TRIM( data_output_pr(i) ) // ' is not imp' // &
	2093	'lemented for humidity = .FALSE.'
[226]	2094	CALL message( 'check_parameters', 'PA0092', 1, 2, 0, 6, 0 )
[1]	2095	ELSE
	2096	dopr_index(i) = 49
[87]	2097	dopr_unit(i) = 'kg/kg m/s'
[1]	2098	hom(:,2,49,:) = SPREAD( zw, 2, statistic_regions+1 )
	2099	ENDIF
	2100
	2101	CASE ( 'wq' )
[108]	2102	IF ( .NOT. humidity ) THEN
[215]	2103	message_string = 'data_output_pr = ' // &
	2104	TRIM( data_output_pr(i) ) // ' is not imp' // &
	2105	'lemented for humidity = .FALSE.'
[226]	2106	CALL message( 'check_parameters', 'PA0092', 1, 2, 0, 6, 0 )
[1]	2107	ELSE
	2108	dopr_index(i) = 50
[87]	2109	dopr_unit(i) = 'kg/kg m/s'
[1]	2110	hom(:,2,50,:) = SPREAD( zw, 2, statistic_regions+1 )
	2111	ENDIF
	2112
	2113	CASE ( 'w"s"' )
	2114	IF ( .NOT. passive_scalar ) THEN
[215]	2115	message_string = 'data_output_pr = ' // &
	2116	TRIM( data_output_pr(i) ) // ' is not imp' // &
	2117	'lemented for passive_scalar = .FALSE.'
[226]	2118	CALL message( 'check_parameters', 'PA0093', 1, 2, 0, 6, 0 )
[1]	2119	ELSE
	2120	dopr_index(i) = 48
[87]	2121	dopr_unit(i) = 'kg/m3 m/s'
[1]	2122	hom(:,2,48,:) = SPREAD( zw, 2, statistic_regions+1 )
	2123	ENDIF
	2124
	2125	CASE ( 'ws' )
	2126	IF ( .NOT. passive_scalar ) THEN
[215]	2127	message_string = 'data_output_pr = ' // &
	2128	TRIM( data_output_pr(i) ) // ' is not imp' // &
	2129	'lemented for passive_scalar = .FALSE.'
[226]	2130	CALL message( 'check_parameters', 'PA0093', 1, 2, 0, 6, 0 )
[1]	2131	ELSE
	2132	dopr_index(i) = 49
[87]	2133	dopr_unit(i) = 'kg/m3 m/s'
[1]	2134	hom(:,2,49,:) = SPREAD( zw, 2, statistic_regions+1 )
	2135	ENDIF
	2136
	2137	CASE ( 'ws' )
	2138	IF ( .NOT. passive_scalar ) THEN
[215]	2139	message_string = 'data_output_pr = ' // &
	2140	TRIM( data_output_pr(i) ) // ' is not imp' // &
	2141	'lemented for passive_scalar = .FALSE.'
[226]	2142	CALL message( 'check_parameters', 'PA0093', 1, 2, 0, 6, 0 )
[1]	2143	ELSE
	2144	dopr_index(i) = 50
[87]	2145	dopr_unit(i) = 'kg/m3 m/s'
[1]	2146	hom(:,2,50,:) = SPREAD( zw, 2, statistic_regions+1 )
	2147	ENDIF
	2148
	2149	CASE ( 'w"qv"' )
[75]	2150	IF ( humidity .AND. .NOT. cloud_physics ) &
[1]	2151	THEN
	2152	dopr_index(i) = 48
[87]	2153	dopr_unit(i) = 'kg/kg m/s'
[1]	2154	hom(:,2,48,:) = SPREAD( zw, 2, statistic_regions+1 )
[75]	2155	ELSEIF( humidity .AND. cloud_physics ) THEN
[1]	2156	dopr_index(i) = 51
[87]	2157	dopr_unit(i) = 'kg/kg m/s'
[1]	2158	hom(:,2,51,:) = SPREAD( zw, 2, statistic_regions+1 )
	2159	ELSE
[215]	2160	message_string = 'data_output_pr = ' // &
	2161	TRIM( data_output_pr(i) ) // ' is not imp' // &
	2162	'lemented for cloud_physics = .FALSE. an&' // &
	2163	'd humidity = .FALSE.'
[226]	2164	CALL message( 'check_parameters', 'PA0095', 1, 2, 0, 6, 0 )
[1]	2165	ENDIF
	2166
	2167	CASE ( 'wqv' )
[75]	2168	IF ( humidity .AND. .NOT. cloud_physics ) &
[1]	2169	THEN
	2170	dopr_index(i) = 49
[87]	2171	dopr_unit(i) = 'kg/kg m/s'
[1]	2172	hom(:,2,49,:) = SPREAD( zw, 2, statistic_regions+1 )
[75]	2173	ELSEIF( humidity .AND. cloud_physics ) THEN
[1]	2174	dopr_index(i) = 52
[87]	2175	dopr_unit(i) = 'kg/kg m/s'
[1]	2176	hom(:,2,52,:) = SPREAD( zw, 2, statistic_regions+1 )
	2177	ELSE
[215]	2178	message_string = 'data_output_pr = ' // &
	2179	TRIM( data_output_pr(i) ) // ' is not imp' // &
	2180	'lemented for cloud_physics = .FALSE. an&' // &
	2181	'd humidity = .FALSE.'
[226]	2182	CALL message( 'check_parameters', 'PA0095', 1, 2, 0, 6, 0 )
[1]	2183	ENDIF
	2184
	2185	CASE ( 'wqv' )
[75]	2186	IF ( humidity .AND. .NOT. cloud_physics ) &
[1]	2187	THEN
	2188	dopr_index(i) = 50
[87]	2189	dopr_unit(i) = 'kg/kg m/s'
[1]	2190	hom(:,2,50,:) = SPREAD( zw, 2, statistic_regions+1 )
[75]	2191	ELSEIF( humidity .AND. cloud_physics ) THEN
[1]	2192	dopr_index(i) = 53
[87]	2193	dopr_unit(i) = 'kg/kg m/s'
[1]	2194	hom(:,2,53,:) = SPREAD( zw, 2, statistic_regions+1 )
	2195	ELSE
[215]	2196	message_string = 'data_output_pr = ' // &
	2197	TRIM( data_output_pr(i) ) // ' is not imp' // &
	2198	'lemented for cloud_physics = .FALSE. an&' // &
	2199	'd humidity = .FALSE.'
[226]	2200	CALL message( 'check_parameters', 'PA0095', 1, 2, 0, 6, 0 )
[1]	2201	ENDIF
	2202
	2203	CASE ( 'ql' )
	2204	IF ( .NOT. cloud_physics .AND. .NOT. cloud_droplets ) THEN
[215]	2205	message_string = 'data_output_pr = ' // &
	2206	TRIM( data_output_pr(i) ) // ' is not imp' // &
	2207	'lemented for cloud_physics = .FALSE. or' // &
	2208	'&cloud_droplets = .FALSE.'
[226]	2209	CALL message( 'check_parameters', 'PA0096', 1, 2, 0, 6, 0 )
[1]	2210	ELSE
	2211	dopr_index(i) = 54
[87]	2212	dopr_unit(i) = 'kg/kg'
[1]	2213	hom(:,2,54,:) = SPREAD( zu, 2, statistic_regions+1 )
	2214	ENDIF
	2215
[524]	2216	CASE ( 'wuu*:dz' )
[1]	2217	dopr_index(i) = 55
[87]	2218	dopr_unit(i) = 'm2/s3'
[1]	2219	hom(:,2,55,:) = SPREAD( zu, 2, statistic_regions+1 )
	2220
[524]	2221	CASE ( 'wp:dz' )
[1]	2222	dopr_index(i) = 56
[87]	2223	dopr_unit(i) = 'm2/s3'
[106]	2224	hom(:,2,56,:) = SPREAD( zw, 2, statistic_regions+1 )
[1]	2225
[524]	2226	CASE ( 'w"e:dz' )
[1]	2227	dopr_index(i) = 57
[87]	2228	dopr_unit(i) = 'm2/s3'
[1]	2229	hom(:,2,57,:) = SPREAD( zu, 2, statistic_regions+1 )
	2230
[667]	2231
[1]	2232	CASE ( 'u"pt"' )
	2233	dopr_index(i) = 58
[87]	2234	dopr_unit(i) = 'K m/s'
[1]	2235	hom(:,2,58,:) = SPREAD( zu, 2, statistic_regions+1 )
	2236
	2237	CASE ( 'upt' )
	2238	dopr_index(i) = 59
[87]	2239	dopr_unit(i) = 'K m/s'
[1]	2240	hom(:,2,59,:) = SPREAD( zu, 2, statistic_regions+1 )
	2241
	2242	CASE ( 'upt_t' )
	2243	dopr_index(i) = 60
[87]	2244	dopr_unit(i) = 'K m/s'
[1]	2245	hom(:,2,60,:) = SPREAD( zu, 2, statistic_regions+1 )
	2246
	2247	CASE ( 'v"pt"' )
	2248	dopr_index(i) = 61
[87]	2249	dopr_unit(i) = 'K m/s'
[1]	2250	hom(:,2,61,:) = SPREAD( zu, 2, statistic_regions+1 )
	2251
	2252	CASE ( 'vpt' )
	2253	dopr_index(i) = 62
[87]	2254	dopr_unit(i) = 'K m/s'
[1]	2255	hom(:,2,62,:) = SPREAD( zu, 2, statistic_regions+1 )
	2256
	2257	CASE ( 'vpt_t' )
	2258	dopr_index(i) = 63
[87]	2259	dopr_unit(i) = 'K m/s'
[1]	2260	hom(:,2,63,:) = SPREAD( zu, 2, statistic_regions+1 )
	2261
[96]	2262	CASE ( 'rho' )
[388]	2263	IF ( .NOT. ocean ) THEN
	2264	message_string = 'data_output_pr = ' // &
	2265	TRIM( data_output_pr(i) ) // ' is not imp' // &
	2266	'lemented for ocean = .FALSE.'
	2267	CALL message( 'check_parameters', 'PA0091', 1, 2, 0, 6, 0 )
	2268	ELSE
	2269	dopr_index(i) = 64
	2270	dopr_unit(i) = 'kg/m3'
	2271	hom(:,2,64,:) = SPREAD( zu, 2, statistic_regions+1 )
	2272	ENDIF
[1]	2273
[96]	2274	CASE ( 'w"sa"' )
	2275	IF ( .NOT. ocean ) THEN
[215]	2276	message_string = 'data_output_pr = ' // &
	2277	TRIM( data_output_pr(i) ) // ' is not imp' // &
	2278	'lemented for ocean = .FALSE.'
[226]	2279	CALL message( 'check_parameters', 'PA0091', 1, 2, 0, 6, 0 )
[96]	2280	ELSE
	2281	dopr_index(i) = 65
	2282	dopr_unit(i) = 'psu m/s'
	2283	hom(:,2,65,:) = SPREAD( zw, 2, statistic_regions+1 )
	2284	ENDIF
	2285
	2286	CASE ( 'wsa' )
	2287	IF ( .NOT. ocean ) THEN
[215]	2288	message_string = 'data_output_pr = ' // &
	2289	TRIM( data_output_pr(i) ) // ' is not imp' // &
	2290	'lemented for ocean = .FALSE.'
[226]	2291	CALL message( 'check_parameters', 'PA0091', 1, 2, 0, 6, 0 )
[96]	2292	ELSE
	2293	dopr_index(i) = 66
	2294	dopr_unit(i) = 'psu m/s'
	2295	hom(:,2,66,:) = SPREAD( zw, 2, statistic_regions+1 )
	2296	ENDIF
	2297
	2298	CASE ( 'wsa' )
	2299	IF ( .NOT. ocean ) THEN
[215]	2300	message_string = 'data_output_pr = ' // &
	2301	TRIM( data_output_pr(i) ) // ' is not imp' // &
	2302	'lemented for ocean = .FALSE.'
[226]	2303	CALL message( 'check_parameters', 'PA0091', 1, 2, 0, 6, 0 )
[96]	2304	ELSE
	2305	dopr_index(i) = 67
	2306	dopr_unit(i) = 'psu m/s'
	2307	hom(:,2,67,:) = SPREAD( zw, 2, statistic_regions+1 )
	2308	ENDIF
	2309
[106]	2310	CASE ( 'wp' )
	2311	dopr_index(i) = 68
	2312	dopr_unit(i) = 'm3/s3'
	2313	hom(:,2,68,:) = SPREAD( zu, 2, statistic_regions+1 )
[96]	2314
[106]	2315	CASE ( 'w"e' )
	2316	dopr_index(i) = 69
	2317	dopr_unit(i) = 'm3/s3'
	2318	hom(:,2,69,:) = SPREAD( zu, 2, statistic_regions+1 )
	2319
[197]	2320	CASE ( 'q*2' )
	2321	IF ( .NOT. humidity ) THEN
[215]	2322	message_string = 'data_output_pr = ' // &
	2323	TRIM( data_output_pr(i) ) // ' is not imp' // &
	2324	'lemented for humidity = .FALSE.'
[226]	2325	CALL message( 'check_parameters', 'PA0092', 1, 2, 0, 6, 0 )
[197]	2326	ELSE
	2327	dopr_index(i) = 70
	2328	dopr_unit(i) = 'kg2/kg2'
	2329	hom(:,2,70,:) = SPREAD( zu, 2, statistic_regions+1 )
	2330	ENDIF
[106]	2331
[388]	2332	CASE ( 'prho' )
	2333	IF ( .NOT. ocean ) THEN
	2334	message_string = 'data_output_pr = ' // &
	2335	TRIM( data_output_pr(i) ) // ' is not imp' // &
	2336	'lemented for ocean = .FALSE.'
	2337	CALL message( 'check_parameters', 'PA0091', 1, 2, 0, 6, 0 )
	2338	ELSE
	2339	dopr_index(i) = 71
	2340	dopr_unit(i) = 'kg/m3'
	2341	hom(:,2,71,:) = SPREAD( zu, 2, statistic_regions+1 )
	2342	ENDIF
	2343
	2344	CASE ( 'hyp' )
	2345	dopr_index(i) = 72
[531]	2346	dopr_unit(i) = 'dbar'
[388]	2347	hom(:,2,72,:) = SPREAD( zu, 2, statistic_regions+1 )
	2348
[1]	2349	CASE DEFAULT
[87]	2350
	2351	CALL user_check_data_output_pr( data_output_pr(i), i, unit )
	2352
	2353	IF ( unit == 'illegal' ) THEN
[215]	2354	IF ( data_output_pr_user(1) /= ' ' ) THEN
	2355	message_string = 'illegal value for data_output_pr or ' // &
	2356	'data_output_pr_user = "' // &
	2357	TRIM( data_output_pr(i) ) // '"'
[226]	2358	CALL message( 'check_parameters', 'PA0097', 1, 2, 0, 6, 0 )
[215]	2359	ELSE
	2360	message_string = 'illegal value for data_output_pr = "' // &
	2361	TRIM( data_output_pr(i) ) // '"'
[226]	2362	CALL message( 'check_parameters', 'PA0098', 1, 2, 0, 6, 0 )
[87]	2363	ENDIF
[1]	2364	ENDIF
	2365
	2366	END SELECT
[667]	2367
[1]	2368	!
	2369	!-- Check to which of the predefined coordinate systems the profile belongs
	2370	DO k = 1, crmax
	2371	IF ( INDEX( cross_profiles(k), ' '//TRIM( data_output_pr(i) )//' ' ) &
	2372	/=0 ) &
	2373	THEN
	2374	dopr_crossindex(i) = k
	2375	EXIT
	2376	ENDIF
	2377	ENDDO
	2378	!
	2379	!-- Generate the text for the labels of the PROFIL output file. "-characters
	2380	!-- must be substituted, otherwise PROFIL would interpret them as TeX
	2381	!-- control characters
	2382	dopr_label(i) = data_output_pr(i)
	2383	position = INDEX( dopr_label(i) , '"' )
	2384	DO WHILE ( position /= 0 )
	2385	dopr_label(i)(position:position) = ''''
	2386	position = INDEX( dopr_label(i) , '"' )
	2387	ENDDO
	2388
	2389	ENDDO
	2390
	2391	!
	2392	!-- y-value range of the coordinate system (PROFIL).
	2393	!-- x-value range determined in plot_1d.
[94]	2394	IF ( .NOT. ocean ) THEN
	2395	cross_uymin = 0.0
	2396	IF ( z_max_do1d == -1.0 ) THEN
	2397	cross_uymax = zu(nzt+1)
	2398	ELSEIF ( z_max_do1d < zu(nzb+1) .OR. z_max_do1d > zu(nzt+1) ) THEN
[215]	2399	WRITE( message_string, * ) 'z_max_do1d = ', z_max_do1d, ' must ', &
	2400	'be >= ', zu(nzb+1), ' or <= ', zu(nzt+1)
[226]	2401	CALL message( 'check_parameters', 'PA0099', 1, 2, 0, 6, 0 )
[94]	2402	ELSE
	2403	cross_uymax = z_max_do1d
	2404	ENDIF
[1]	2405	ENDIF
	2406
	2407	!
	2408	!-- Check whether the chosen normalizing factor for the coordinate systems is
	2409	!-- permissible
	2410	DO i = 1, crmax
	2411	SELECT CASE ( TRIM( cross_normalized_x(i) ) ) ! TRIM required on IBM
	2412
	2413	CASE ( '', 'wpt0', 'ws2', 'tsw2', 'ws3', 'ws2tsw', 'wstsw2' )
	2414	j = 0
	2415
	2416	CASE DEFAULT
[215]	2417	message_string = 'unknown normalization method cross_normali' // &
	2418	'zed_x = "' // TRIM( cross_normalized_x(i) ) // &
	2419	'"'
[226]	2420	CALL message( 'check_parameters', 'PA0100', 1, 2, 0, 6, 0 )
[1]	2421
	2422	END SELECT
	2423	SELECT CASE ( TRIM( cross_normalized_y(i) ) ) ! TRIM required on IBM
	2424
	2425	CASE ( '', 'z_i' )
	2426	j = 0
	2427
	2428	CASE DEFAULT
[215]	2429	message_string = 'unknown normalization method cross_normali' // &
	2430	'zed_y = "' // TRIM( cross_normalized_y(i) ) // &
	2431	'"'
[226]	2432	CALL message( 'check_parameters', 'PA0101', 1, 2, 0, 6, 0 )
[1]	2433
	2434	END SELECT
	2435	ENDDO
	2436	!
	2437	!-- Check normalized y-value range of the coordinate system (PROFIL)
	2438	IF ( z_max_do1d_normalized /= -1.0 .AND. z_max_do1d_normalized <= 0.0 ) &
	2439	THEN
[215]	2440	WRITE( message_string, * ) 'z_max_do1d_normalized = ', &
	2441	z_max_do1d_normalized, ' must be >= 0.0'
[226]	2442	CALL message( 'check_parameters', 'PA0101', 1, 2, 0, 6, 0 )
[1]	2443	ENDIF
	2444
	2445
	2446	!
	2447	!-- Append user-defined data output variables to the standard data output
	2448	IF ( data_output_user(1) /= ' ' ) THEN
	2449	i = 1
	2450	DO WHILE ( data_output(i) /= ' ' .AND. i <= 100 )
	2451	i = i + 1
	2452	ENDDO
	2453	j = 1
	2454	DO WHILE ( data_output_user(j) /= ' ' .AND. j <= 100 )
	2455	IF ( i > 100 ) THEN
[215]	2456	message_string = 'number of output quantitities given by data' // &
	2457	'_output and data_output_user exceeds the limit of 100'
[226]	2458	CALL message( 'check_parameters', 'PA0102', 1, 2, 0, 6, 0 )
[1]	2459	ENDIF
	2460	data_output(i) = data_output_user(j)
	2461	i = i + 1
	2462	j = j + 1
	2463	ENDDO
	2464	ENDIF
	2465
	2466	!
	2467	!-- Check and set steering parameters for 2d/3d data output and averaging
	2468	i = 1
	2469	DO WHILE ( data_output(i) /= ' ' .AND. i <= 100 )
	2470	!
	2471	!-- Check for data averaging
	2472	ilen = LEN_TRIM( data_output(i) )
	2473	j = 0 ! no data averaging
	2474	IF ( ilen > 3 ) THEN
	2475	IF ( data_output(i)(ilen-2:ilen) == '_av' ) THEN
	2476	j = 1 ! data averaging
	2477	data_output(i) = data_output(i)(1:ilen-3)
	2478	ENDIF
	2479	ENDIF
	2480	!
	2481	!-- Check for cross section or volume data
	2482	ilen = LEN_TRIM( data_output(i) )
	2483	k = 0 ! 3d data
	2484	var = data_output(i)(1:ilen)
	2485	IF ( ilen > 3 ) THEN
	2486	IF ( data_output(i)(ilen-2:ilen) == '_xy' .OR. &
	2487	data_output(i)(ilen-2:ilen) == '_xz' .OR. &
	2488	data_output(i)(ilen-2:ilen) == '_yz' ) THEN
	2489	k = 1 ! 2d data
	2490	var = data_output(i)(1:ilen-3)
	2491	ENDIF
	2492	ENDIF
	2493	!
	2494	!-- Check for allowed value and set units
	2495	SELECT CASE ( TRIM( var ) )
	2496
	2497	CASE ( 'e' )
	2498	IF ( constant_diffusion ) THEN
[215]	2499	message_string = 'output of "' // TRIM( var ) // '" requi' // &
	2500	'res constant_diffusion = .FALSE.'
[226]	2501	CALL message( 'check_parameters', 'PA0103', 1, 2, 0, 6, 0 )
[1]	2502	ENDIF
	2503	unit = 'm2/s2'
	2504
	2505	CASE ( 'pc', 'pr' )
	2506	IF ( .NOT. particle_advection ) THEN
[215]	2507	message_string = 'output of "' // TRIM( var ) // '" requir' // &
	2508	'es a "particles_par"-NAMELIST in the parameter file (PARIN)'
[226]	2509	CALL message( 'check_parameters', 'PA0104', 1, 2, 0, 6, 0 )
[1]	2510	ENDIF
	2511	IF ( TRIM( var ) == 'pc' ) unit = 'number'
	2512	IF ( TRIM( var ) == 'pr' ) unit = 'm'
	2513
	2514	CASE ( 'q', 'vpt' )
[75]	2515	IF ( .NOT. humidity ) THEN
[215]	2516	message_string = 'output of "' // TRIM( var ) // '" requi' // &
	2517	'res humidity = .TRUE.'
[226]	2518	CALL message( 'check_parameters', 'PA0105', 1, 2, 0, 6, 0 )
[1]	2519	ENDIF
	2520	IF ( TRIM( var ) == 'q' ) unit = 'kg/kg'
	2521	IF ( TRIM( var ) == 'vpt' ) unit = 'K'
	2522
	2523	CASE ( 'ql' )
	2524	IF ( .NOT. ( cloud_physics .OR. cloud_droplets ) ) THEN
[215]	2525	message_string = 'output of "' // TRIM( var ) // '" requi' // &
	2526	'res cloud_physics = .TRUE. or cloud_droplets = .TRUE.'
[226]	2527	CALL message( 'check_parameters', 'PA0106', 1, 2, 0, 6, 0 )
[1]	2528	ENDIF
	2529	unit = 'kg/kg'
	2530
	2531	CASE ( 'ql_c', 'ql_v', 'ql_vp' )
	2532	IF ( .NOT. cloud_droplets ) THEN
[215]	2533	message_string = 'output of "' // TRIM( var ) // '" requi' // &
	2534	'res cloud_droplets = .TRUE.'
[226]	2535	CALL message( 'check_parameters', 'PA0107', 1, 2, 0, 6, 0 )
[1]	2536	ENDIF
	2537	IF ( TRIM( var ) == 'ql_c' ) unit = 'kg/kg'
	2538	IF ( TRIM( var ) == 'ql_v' ) unit = 'm3'
	2539	IF ( TRIM( var ) == 'ql_vp' ) unit = 'none'
	2540
	2541	CASE ( 'qv' )
	2542	IF ( .NOT. cloud_physics ) THEN
[215]	2543	message_string = 'output of "' // TRIM( var ) // '" requi' // &
	2544	'res cloud_physics = .TRUE.'
[226]	2545	CALL message( 'check_parameters', 'PA0108', 1, 2, 0, 6, 0 )
[1]	2546	ENDIF
	2547	unit = 'kg/kg'
	2548
[96]	2549	CASE ( 'rho' )
	2550	IF ( .NOT. ocean ) THEN
[215]	2551	message_string = 'output of "' // TRIM( var ) // '" requi' // &
	2552	'res ocean = .TRUE.'
[226]	2553	CALL message( 'check_parameters', 'PA0109', 1, 2, 0, 6, 0 )
[96]	2554	ENDIF
	2555	unit = 'kg/m3'
	2556
[1]	2557	CASE ( 's' )
	2558	IF ( .NOT. passive_scalar ) THEN
[215]	2559	message_string = 'output of "' // TRIM( var ) // '" requi' // &
	2560	'res passive_scalar = .TRUE.'
[226]	2561	CALL message( 'check_parameters', 'PA0110', 1, 2, 0, 6, 0 )
[1]	2562	ENDIF
	2563	unit = 'conc'
	2564
[96]	2565	CASE ( 'sa' )
	2566	IF ( .NOT. ocean ) THEN
[215]	2567	message_string = 'output of "' // TRIM( var ) // '" requi' // &
	2568	'res ocean = .TRUE.'
[226]	2569	CALL message( 'check_parameters', 'PA0109', 1, 2, 0, 6, 0 )
[96]	2570	ENDIF
	2571	unit = 'psu'
	2572
[354]	2573	CASE ( 'u', 't', 'lwp', 'pra', 'prr', 'qsws', 'shf', 'z0' )
[1]	2574	IF ( k == 0 .OR. data_output(i)(ilen-2:ilen) /= '_xy' ) THEN
[215]	2575	message_string = 'illegal value for data_output: "' // &
	2576	TRIM( var ) // '" & only 2d-horizontal ' // &
	2577	'cross sections are allowed for this value'
[226]	2578	CALL message( 'check_parameters', 'PA0111', 1, 2, 0, 6, 0 )
[1]	2579	ENDIF
	2580	IF ( TRIM( var ) == 'lwp*' .AND. .NOT. cloud_physics ) THEN
[215]	2581	message_string = 'output of "' // TRIM( var ) // '" requi' // &
	2582	'res cloud_physics = .TRUE.'
[226]	2583	CALL message( 'check_parameters', 'PA0108', 1, 2, 0, 6, 0 )
[1]	2584	ENDIF
[72]	2585	IF ( TRIM( var ) == 'pra*' .AND. .NOT. precipitation ) THEN
[215]	2586	message_string = 'output of "' // TRIM( var ) // '" requi' // &
	2587	'res precipitation = .TRUE.'
[226]	2588	CALL message( 'check_parameters', 'PA0112', 1, 2, 0, 6, 0 )
[72]	2589	ENDIF
	2590	IF ( TRIM( var ) == 'pra*' .AND. j == 1 ) THEN
[215]	2591	message_string = 'temporal averaging of precipitation ' // &
	2592	'amount "' // TRIM( var ) // '" is not possible'
[226]	2593	CALL message( 'check_parameters', 'PA0113', 1, 2, 0, 6, 0 )
[72]	2594	ENDIF
	2595	IF ( TRIM( var ) == 'prr*' .AND. .NOT. precipitation ) THEN
[215]	2596	message_string = 'output of "' // TRIM( var ) // '" requi' // &
	2597	'res precipitation = .TRUE.'
[226]	2598	CALL message( 'check_parameters', 'PA0112', 1, 2, 0, 6, 0 )
[72]	2599	ENDIF
[354]	2600	IF ( TRIM( var ) == 'qsws*' .AND. .NOT. humidity ) THEN
	2601	message_string = 'output of "' // TRIM( var ) // '" requi' // &
	2602	'res humidity = .TRUE.'
	2603	CALL message( 'check_parameters', 'PA0322', 1, 2, 0, 6, 0 )
	2604	ENDIF
[72]	2605
[354]	2606	IF ( TRIM( var ) == 'lwp' ) unit = 'kg/kgm'
	2607	IF ( TRIM( var ) == 'pra*' ) unit = 'mm'
	2608	IF ( TRIM( var ) == 'prr*' ) unit = 'mm/s'
	2609	IF ( TRIM( var ) == 'qsws*' ) unit = 'kgm/kgs'
	2610	IF ( TRIM( var ) == 'shf' ) unit = 'Km/s'
	2611	IF ( TRIM( var ) == 't*' ) unit = 'K'
	2612	IF ( TRIM( var ) == 'u*' ) unit = 'm/s'
	2613	IF ( TRIM( var ) == 'z0*' ) unit = 'm'
[72]	2614
[1]	2615
	2616	CASE ( 'p', 'pt', 'u', 'v', 'w' )
	2617	IF ( TRIM( var ) == 'p' ) unit = 'Pa'
	2618	IF ( TRIM( var ) == 'pt' ) unit = 'K'
	2619	IF ( TRIM( var ) == 'u' ) unit = 'm/s'
	2620	IF ( TRIM( var ) == 'v' ) unit = 'm/s'
	2621	IF ( TRIM( var ) == 'w' ) unit = 'm/s'
	2622	CONTINUE
	2623
	2624	CASE DEFAULT
	2625	CALL user_check_data_output( var, unit )
	2626
	2627	IF ( unit == 'illegal' ) THEN
[215]	2628	IF ( data_output_user(1) /= ' ' ) THEN
	2629	message_string = 'illegal value for data_output or ' // &
	2630	'data_output_user = "' // TRIM( data_output(i) ) // '"'
[226]	2631	CALL message( 'check_parameters', 'PA0114', 1, 2, 0, 6, 0 )
[215]	2632	ELSE
	2633	message_string = 'illegal value for data_output =' // &
	2634	TRIM( data_output(i) ) // '"'
[226]	2635	CALL message( 'check_parameters', 'PA0115', 1, 2, 0, 6, 0 )
[1]	2636	ENDIF
	2637	ENDIF
	2638
	2639	END SELECT
	2640	!
	2641	!-- Set the internal steering parameters appropriately
	2642	IF ( k == 0 ) THEN
	2643	do3d_no(j) = do3d_no(j) + 1
	2644	do3d(j,do3d_no(j)) = data_output(i)
	2645	do3d_unit(j,do3d_no(j)) = unit
	2646	ELSE
	2647	do2d_no(j) = do2d_no(j) + 1
	2648	do2d(j,do2d_no(j)) = data_output(i)
	2649	do2d_unit(j,do2d_no(j)) = unit
	2650	IF ( data_output(i)(ilen-2:ilen) == '_xy' ) THEN
	2651	data_output_xy(j) = .TRUE.
	2652	ENDIF
	2653	IF ( data_output(i)(ilen-2:ilen) == '_xz' ) THEN
	2654	data_output_xz(j) = .TRUE.
	2655	ENDIF
	2656	IF ( data_output(i)(ilen-2:ilen) == '_yz' ) THEN
	2657	data_output_yz(j) = .TRUE.
	2658	ENDIF
	2659	ENDIF
	2660
	2661	IF ( j == 1 ) THEN
	2662	!
	2663	!-- Check, if variable is already subject to averaging
	2664	found = .FALSE.
	2665	DO k = 1, doav_n
	2666	IF ( TRIM( doav(k) ) == TRIM( var ) ) found = .TRUE.
	2667	ENDDO
	2668
	2669	IF ( .NOT. found ) THEN
	2670	doav_n = doav_n + 1
	2671	doav(doav_n) = var
	2672	ENDIF
	2673	ENDIF
	2674
	2675	i = i + 1
	2676	ENDDO
	2677
	2678	!
[376]	2679	!-- Averaged 2d or 3d output requires that an averaging interval has been set
	2680	IF ( doav_n > 0 .AND. averaging_interval == 0.0 ) THEN
	2681	WRITE( message_string, * ) 'output of averaged quantity "', &
	2682	TRIM( doav(1) ), '_av" requires to set a ', &
	2683	'non-zero & averaging interval'
	2684	CALL message( 'check_parameters', 'PA0323', 1, 2, 0, 6, 0 )
	2685	ENDIF
	2686
	2687	!
[308]	2688	!-- Check sectional planes and store them in one shared array
	2689	IF ( ANY( section_xy > nz + 1 ) ) THEN
	2690	WRITE( message_string, * ) 'section_xy must be <= nz + 1 = ', nz + 1
	2691	CALL message( 'check_parameters', 'PA0319', 1, 2, 0, 6, 0 )
	2692	ENDIF
	2693	IF ( ANY( section_xz > ny + 1 ) ) THEN
	2694	WRITE( message_string, * ) 'section_xz must be <= ny + 1 = ', ny + 1
	2695	CALL message( 'check_parameters', 'PA0320', 1, 2, 0, 6, 0 )
	2696	ENDIF
	2697	IF ( ANY( section_yz > nx + 1 ) ) THEN
	2698	WRITE( message_string, * ) 'section_yz must be <= nx + 1 = ', nx + 1
	2699	CALL message( 'check_parameters', 'PA0321', 1, 2, 0, 6, 0 )
	2700	ENDIF
[1]	2701	section(:,1) = section_xy
	2702	section(:,2) = section_xz
	2703	section(:,3) = section_yz
	2704
	2705	!
	2706	!-- Upper plot limit (grid point value) for 1D profiles
	2707	IF ( z_max_do1d == -1.0 ) THEN
[667]	2708
[1]	2709	nz_do1d = nzt+1
[667]	2710
[1]	2711	ELSE
	2712	DO k = nzb+1, nzt+1
	2713	nz_do1d = k
	2714	IF ( zw(k) > z_max_do1d ) EXIT
	2715	ENDDO
	2716	ENDIF
	2717
	2718	!
	2719	!-- Upper plot limit for 2D vertical sections
	2720	IF ( z_max_do2d == -1.0 ) z_max_do2d = zu(nzt)
	2721	IF ( z_max_do2d < zu(nzb+1) .OR. z_max_do2d > zu(nzt) ) THEN
[215]	2722	WRITE( message_string, * ) 'z_max_do2d = ', z_max_do2d, &
	2723	' must be >= ', zu(nzb+1), '(zu(nzb+1)) and <= ', zu(nzt), &
	2724	' (zu(nzt))'
[226]	2725	CALL message( 'check_parameters', 'PA0116', 1, 2, 0, 6, 0 )
[1]	2726	ENDIF
	2727
	2728	!
	2729	!-- Upper plot limit for 3D arrays
	2730	IF ( nz_do3d == -9999 ) nz_do3d = nzt + 1
	2731
	2732	!
	2733	!-- Determine and check accuracy for compressed 3D plot output
	2734	IF ( do3d_compress ) THEN
	2735	!
	2736	!-- Compression only permissible on T3E machines
	2737	IF ( host(1:3) /= 't3e' ) THEN
[215]	2738	message_string = 'do3d_compress = .TRUE. not allowed on host "' // &
	2739	TRIM( host ) // '"'
[226]	2740	CALL message( 'check_parameters', 'PA0117', 1, 2, 0, 6, 0 )
[1]	2741	ENDIF
	2742
	2743	i = 1
	2744	DO WHILE ( do3d_comp_prec(i) /= ' ' )
	2745
	2746	ilen = LEN_TRIM( do3d_comp_prec(i) )
	2747	IF ( LLT( do3d_comp_prec(i)(ilen:ilen), '0' ) .OR. &
	2748	LGT( do3d_comp_prec(i)(ilen:ilen), '9' ) ) THEN
[215]	2749	WRITE( message_string, * ) 'illegal precision: do3d_comp_prec', &
	2750	'(', i, ') = "', TRIM(do3d_comp_prec(i)),'"'
[226]	2751	CALL message( 'check_parameters', 'PA0118', 1, 2, 0, 6, 0 )
[1]	2752	ENDIF
	2753
	2754	prec = IACHAR( do3d_comp_prec(i)(ilen:ilen) ) - IACHAR( '0' )
	2755	var = do3d_comp_prec(i)(1:ilen-1)
	2756
	2757	SELECT CASE ( var )
	2758
	2759	CASE ( 'u' )
	2760	j = 1
	2761	CASE ( 'v' )
	2762	j = 2
	2763	CASE ( 'w' )
	2764	j = 3
	2765	CASE ( 'p' )
	2766	j = 4
	2767	CASE ( 'pt' )
	2768	j = 5
	2769
	2770	CASE DEFAULT
[215]	2771	WRITE( message_string, * ) 'unknown variable "', &
	2772	TRIM( do3d_comp_prec(i) ), '" given for do3d_comp_prec(', &
	2773	i, ')'
[226]	2774	CALL message( 'check_parameters', 'PA0119', 1, 2, 0, 6, 0 )
[1]	2775
	2776	END SELECT
	2777
	2778	plot_3d_precision(j)%precision = prec
	2779	i = i + 1
	2780
	2781	ENDDO
	2782	ENDIF
	2783
	2784	!
	2785	!-- Check the data output format(s)
	2786	IF ( data_output_format(1) == ' ' ) THEN
	2787	!
	2788	!-- Default value
	2789	netcdf_output = .TRUE.
	2790	ELSE
	2791	i = 1
	2792	DO WHILE ( data_output_format(i) /= ' ' )
	2793
	2794	SELECT CASE ( data_output_format(i) )
	2795
	2796	CASE ( 'netcdf' )
	2797	netcdf_output = .TRUE.
	2798	CASE ( 'iso2d' )
	2799	iso2d_output = .TRUE.
	2800	CASE ( 'profil' )
	2801	profil_output = .TRUE.
	2802	CASE ( 'avs' )
	2803	avs_output = .TRUE.
	2804
	2805	CASE DEFAULT
[215]	2806	message_string = 'unknown value for data_output_format "' // &
	2807	TRIM( data_output_format(i) ) // '"'
[226]	2808	CALL message( 'check_parameters', 'PA0120', 1, 2, 0, 6, 0 )
[1]	2809
	2810	END SELECT
	2811
	2812	i = i + 1
	2813	IF ( i > 10 ) EXIT
	2814
	2815	ENDDO
	2816
	2817	ENDIF
	2818
	2819	!
[410]	2820	!-- Check mask conditions
[553]	2821	DO mid = 1, max_masks
[567]	2822	IF ( data_output_masks(mid,1) /= ' ' .OR. &
	2823	data_output_masks_user(mid,1) /= ' ' ) THEN
[553]	2824	masks = masks + 1
	2825	ENDIF
	2826	ENDDO
	2827
[410]	2828	IF ( masks < 0 .OR. masks > max_masks ) THEN
	2829	WRITE( message_string, * ) 'illegal value: masks must be >= 0 and ', &
	2830	'<= ', max_masks, ' (=max_masks)'
[564]	2831	CALL message( 'check_parameters', 'PA0325', 1, 2, 0, 6, 0 )
[410]	2832	ENDIF
	2833	IF ( masks > 0 ) THEN
	2834	mask_scale(1) = mask_scale_x
	2835	mask_scale(2) = mask_scale_y
	2836	mask_scale(3) = mask_scale_z
	2837	IF ( ANY( mask_scale <= 0.0 ) ) THEN
	2838	WRITE( message_string, * ) &
	2839	'illegal value: mask_scale_x, mask_scale_y and mask_scale_z', &
	2840	'must be > 0.0'
[564]	2841	CALL message( 'check_parameters', 'PA0326', 1, 2, 0, 6, 0 )
[410]	2842	ENDIF
	2843	!
	2844	!-- Generate masks for masked data output
	2845	CALL init_masks
	2846	ENDIF
	2847
	2848	!
[493]	2849	!-- Check the NetCDF data format
	2850	IF ( netcdf_data_format > 2 ) THEN
	2851	#if defined( __netcdf4 )
	2852	CONTINUE
	2853	#else
	2854	message_string = 'NetCDF: NetCDF4 format requested but no ' // &
	2855	'cpp-directive __netcdf4 given & switch ' // &
	2856	'back to 64-bit offset format'
	2857	CALL message( 'check_parameters', 'PA0171', 0, 1, 0, 6, 0 )
	2858	netcdf_data_format = 2
	2859	#endif
	2860	ENDIF
	2861
	2862	!
[667]	2863
[1]	2864	!-- Check netcdf precison
	2865	ldum = .FALSE.
	2866	CALL define_netcdf_header( 'ch', ldum, 0 )
	2867
	2868	!
	2869	!-- Check, whether a constant diffusion coefficient shall be used
	2870	IF ( km_constant /= -1.0 ) THEN
	2871	IF ( km_constant < 0.0 ) THEN
[215]	2872	WRITE( message_string, * ) 'km_constant = ', km_constant, ' < 0.0'
[226]	2873	CALL message( 'check_parameters', 'PA0121', 1, 2, 0, 6, 0 )
[1]	2874	ELSE
	2875	IF ( prandtl_number < 0.0 ) THEN
[215]	2876	WRITE( message_string, * ) 'prandtl_number = ', prandtl_number, &
	2877	' < 0.0'
[226]	2878	CALL message( 'check_parameters', 'PA0122', 1, 2, 0, 6, 0 )
[1]	2879	ENDIF
	2880	constant_diffusion = .TRUE.
	2881
	2882	IF ( prandtl_layer ) THEN
[215]	2883	message_string = 'prandtl_layer is not allowed with fixed ' // &
	2884	'value of km'
[226]	2885	CALL message( 'check_parameters', 'PA0123', 1, 2, 0, 6, 0 )
[1]	2886	ENDIF
	2887	ENDIF
	2888	ENDIF
	2889
	2890	!
	2891	!-- In case of non-cyclic lateral boundaries, set the default maximum value
	2892	!-- for the horizontal diffusivity used within the outflow damping layer,
	2893	!-- and check/set the width of the damping layer
	2894	IF ( bc_lr /= 'cyclic' ) THEN
	2895	IF ( km_damp_max == -1.0 ) THEN
	2896	km_damp_max = 0.5 * dx
	2897	ENDIF
	2898	IF ( outflow_damping_width == -1.0 ) THEN
	2899	outflow_damping_width = MIN( 20, nx/2 )
	2900	ENDIF
	2901	IF ( outflow_damping_width <= 0 .OR. outflow_damping_width > nx ) THEN
[215]	2902	message_string = 'outflow_damping width out of range'
[226]	2903	CALL message( 'check_parameters', 'PA0124', 1, 2, 0, 6, 0 )
[1]	2904	ENDIF
	2905	ENDIF
	2906
	2907	IF ( bc_ns /= 'cyclic' ) THEN
	2908	IF ( km_damp_max == -1.0 ) THEN
	2909	km_damp_max = 0.5 * dy
	2910	ENDIF
	2911	IF ( outflow_damping_width == -1.0 ) THEN
	2912	outflow_damping_width = MIN( 20, ny/2 )
	2913	ENDIF
	2914	IF ( outflow_damping_width <= 0 .OR. outflow_damping_width > ny ) THEN
[215]	2915	message_string = 'outflow_damping width out of range'
[226]	2916	CALL message( 'check_parameters', 'PA0124', 1, 2, 0, 6, 0 )
[1]	2917	ENDIF
	2918	ENDIF
	2919
	2920	!
	2921	!-- Check value range for rif
	2922	IF ( rif_min >= rif_max ) THEN
[215]	2923	WRITE( message_string, * ) 'rif_min = ', rif_min, ' must be less ', &
	2924	'than rif_max = ', rif_max
[226]	2925	CALL message( 'check_parameters', 'PA0125', 1, 2, 0, 6, 0 )
[1]	2926	ENDIF
	2927
	2928	!
	2929	!-- Determine upper and lower hight level indices for random perturbations
[97]	2930	IF ( disturbance_level_b == -9999999.9 ) THEN
	2931	IF ( ocean ) THEN
	2932	disturbance_level_b = zu((nzt*2)/3)
	2933	disturbance_level_ind_b = ( nzt * 2 ) / 3
	2934	ELSE
	2935	disturbance_level_b = zu(nzb+3)
	2936	disturbance_level_ind_b = nzb + 3
	2937	ENDIF
[1]	2938	ELSEIF ( disturbance_level_b < zu(3) ) THEN
[215]	2939	WRITE( message_string, * ) 'disturbance_level_b = ', &
	2940	disturbance_level_b, ' must be >= ', zu(3), '(zu(3))'
[226]	2941	CALL message( 'check_parameters', 'PA0126', 1, 2, 0, 6, 0 )
[1]	2942	ELSEIF ( disturbance_level_b > zu(nzt-2) ) THEN
[215]	2943	WRITE( message_string, * ) 'disturbance_level_b = ', &
	2944	disturbance_level_b, ' must be <= ', zu(nzt-2), '(zu(nzt-2))'
[226]	2945	CALL message( 'check_parameters', 'PA0127', 1, 2, 0, 6, 0 )
[1]	2946	ELSE
	2947	DO k = 3, nzt-2
	2948	IF ( disturbance_level_b <= zu(k) ) THEN
	2949	disturbance_level_ind_b = k
	2950	EXIT
	2951	ENDIF
	2952	ENDDO
	2953	ENDIF
	2954
[97]	2955	IF ( disturbance_level_t == -9999999.9 ) THEN
	2956	IF ( ocean ) THEN
	2957	disturbance_level_t = zu(nzt-3)
	2958	disturbance_level_ind_t = nzt - 3
	2959	ELSE
	2960	disturbance_level_t = zu(nzt/3)
	2961	disturbance_level_ind_t = nzt / 3
	2962	ENDIF
[1]	2963	ELSEIF ( disturbance_level_t > zu(nzt-2) ) THEN
[215]	2964	WRITE( message_string, * ) 'disturbance_level_t = ', &
	2965	disturbance_level_t, ' must be <= ', zu(nzt-2), '(zu(nzt-2))'
[226]	2966	CALL message( 'check_parameters', 'PA0128', 1, 2, 0, 6, 0 )
[1]	2967	ELSEIF ( disturbance_level_t < disturbance_level_b ) THEN
[215]	2968	WRITE( message_string, * ) 'disturbance_level_t = ', &
	2969	disturbance_level_t, ' must be >= disturbance_level_b = ', &
	2970	disturbance_level_b
[226]	2971	CALL message( 'check_parameters', 'PA0129', 1, 2, 0, 6, 0 )
[1]	2972	ELSE
	2973	DO k = 3, nzt-2
	2974	IF ( disturbance_level_t <= zu(k) ) THEN
	2975	disturbance_level_ind_t = k
	2976	EXIT
	2977	ENDIF
	2978	ENDDO
	2979	ENDIF
	2980
	2981	!
	2982	!-- Check again whether the levels determined this way are ok.
	2983	!-- Error may occur at automatic determination and too few grid points in
	2984	!-- z-direction.
	2985	IF ( disturbance_level_ind_t < disturbance_level_ind_b ) THEN
[215]	2986	WRITE( message_string, * ) 'disturbance_level_ind_t = ', &
	2987	disturbance_level_ind_t, ' must be >= disturbance_level_b = ', &
	2988	disturbance_level_b
[226]	2989	CALL message( 'check_parameters', 'PA0130', 1, 2, 0, 6, 0 )
[1]	2990	ENDIF
	2991
	2992	!
	2993	!-- Determine the horizontal index range for random perturbations.
	2994	!-- In case of non-cyclic horizontal boundaries, no perturbations are imposed
	2995	!-- near the inflow and the perturbation area is further limited to ...(1)
	2996	!-- after the initial phase of the flow.
	2997	dist_nxl = 0; dist_nxr = nx
	2998	dist_nys = 0; dist_nyn = ny
	2999	IF ( bc_lr /= 'cyclic' ) THEN
	3000	IF ( inflow_disturbance_begin == -1 ) THEN
	3001	inflow_disturbance_begin = MIN( 10, nx/2 )
	3002	ENDIF
	3003	IF ( inflow_disturbance_begin < 0 .OR. inflow_disturbance_begin > nx )&
	3004	THEN
[215]	3005	message_string = 'inflow_disturbance_begin out of range'
[226]	3006	CALL message( 'check_parameters', 'PA0131', 1, 2, 0, 6, 0 )
[1]	3007	ENDIF
	3008	IF ( inflow_disturbance_end == -1 ) THEN
	3009	inflow_disturbance_end = MIN( 100, 3*nx/4 )
	3010	ENDIF
	3011	IF ( inflow_disturbance_end < 0 .OR. inflow_disturbance_end > nx ) &
	3012	THEN
[215]	3013	message_string = 'inflow_disturbance_end out of range'
[226]	3014	CALL message( 'check_parameters', 'PA0132', 1, 2, 0, 6, 0 )
[1]	3015	ENDIF
	3016	ELSEIF ( bc_ns /= 'cyclic' ) THEN
	3017	IF ( inflow_disturbance_begin == -1 ) THEN
	3018	inflow_disturbance_begin = MIN( 10, ny/2 )
	3019	ENDIF
	3020	IF ( inflow_disturbance_begin < 0 .OR. inflow_disturbance_begin > ny )&
	3021	THEN
[215]	3022	message_string = 'inflow_disturbance_begin out of range'
[226]	3023	CALL message( 'check_parameters', 'PA0131', 1, 2, 0, 6, 0 )
[1]	3024	ENDIF
	3025	IF ( inflow_disturbance_end == -1 ) THEN
	3026	inflow_disturbance_end = MIN( 100, 3*ny/4 )
	3027	ENDIF
	3028	IF ( inflow_disturbance_end < 0 .OR. inflow_disturbance_end > ny ) &
	3029	THEN
[215]	3030	message_string = 'inflow_disturbance_end out of range'
[226]	3031	CALL message( 'check_parameters', 'PA0132', 1, 2, 0, 6, 0 )
[1]	3032	ENDIF
	3033	ENDIF
	3034
[73]	3035	IF ( bc_lr == 'radiation/dirichlet' ) THEN
[1]	3036	dist_nxr = nx - inflow_disturbance_begin
	3037	dist_nxl(1) = nx - inflow_disturbance_end
[73]	3038	ELSEIF ( bc_lr == 'dirichlet/radiation' ) THEN
[1]	3039	dist_nxl = inflow_disturbance_begin
	3040	dist_nxr(1) = inflow_disturbance_end
[73]	3041	ENDIF
	3042	IF ( bc_ns == 'dirichlet/radiation' ) THEN
[1]	3043	dist_nyn = ny - inflow_disturbance_begin
	3044	dist_nys(1) = ny - inflow_disturbance_end
[73]	3045	ELSEIF ( bc_ns == 'radiation/dirichlet' ) THEN
[1]	3046	dist_nys = inflow_disturbance_begin
	3047	dist_nyn(1) = inflow_disturbance_end
	3048	ENDIF
	3049
	3050	!
[151]	3051	!-- A turbulent inflow requires Dirichlet conditions at the respective inflow
	3052	!-- boundary (so far, a turbulent inflow is realized from the left side only)
	3053	IF ( turbulent_inflow .AND. bc_lr /= 'dirichlet/radiation' ) THEN
[215]	3054	message_string = 'turbulent_inflow = .T. requires a Dirichlet ' // &
	3055	'condition at the inflow boundary'
[226]	3056	CALL message( 'check_parameters', 'PA0133', 1, 2, 0, 6, 0 )
[151]	3057	ENDIF
	3058
	3059	!
	3060	!-- In case of turbulent inflow calculate the index of the recycling plane
	3061	IF ( turbulent_inflow ) THEN
	3062	IF ( recycling_width == 9999999.9 ) THEN
	3063	!
	3064	!-- Set the default value for the width of the recycling domain
	3065	recycling_width = 0.1 * nx * dx
	3066	ELSE
	3067	IF ( recycling_width < dx .OR. recycling_width > nx * dx ) THEN
[215]	3068	WRITE( message_string, * ) 'illegal value for recycling_width:', &
	3069	' ', recycling_width
[226]	3070	CALL message( 'check_parameters', 'PA0134', 1, 2, 0, 6, 0 )
[151]	3071	ENDIF
	3072	ENDIF
	3073	!
	3074	!-- Calculate the index
	3075	recycling_plane = recycling_width / dx
	3076	ENDIF
	3077
	3078	!
[1]	3079	!-- Check random generator
	3080	IF ( random_generator /= 'system-specific' .AND. &
	3081	random_generator /= 'numerical-recipes' ) THEN
[215]	3082	message_string = 'unknown random generator: random_generator = "' // &
	3083	TRIM( random_generator ) // '"'
[226]	3084	CALL message( 'check_parameters', 'PA0135', 1, 2, 0, 6, 0 )
[1]	3085	ENDIF
	3086
	3087	!
	3088	!-- Determine damping level index for 1D model
	3089	IF ( INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 ) THEN
	3090	IF ( damp_level_1d == -1.0 ) THEN
	3091	damp_level_1d = zu(nzt+1)
	3092	damp_level_ind_1d = nzt + 1
	3093	ELSEIF ( damp_level_1d < 0.0 .OR. damp_level_1d > zu(nzt+1) ) THEN
[215]	3094	WRITE( message_string, * ) 'damp_level_1d = ', damp_level_1d, &
	3095	' must be > 0.0 and < ', zu(nzt+1), '(zu(nzt+1))'
[226]	3096	CALL message( 'check_parameters', 'PA0136', 1, 2, 0, 6, 0 )
[1]	3097	ELSE
	3098	DO k = 1, nzt+1
	3099	IF ( damp_level_1d <= zu(k) ) THEN
	3100	damp_level_ind_1d = k
	3101	EXIT
	3102	ENDIF
	3103	ENDDO
	3104	ENDIF
	3105	ENDIF
[215]	3106
[1]	3107	!
	3108	!-- Check some other 1d-model parameters
	3109	IF ( TRIM( mixing_length_1d ) /= 'as_in_3d_model' .AND. &
	3110	TRIM( mixing_length_1d ) /= 'blackadar' ) THEN
[215]	3111	message_string = 'mixing_length_1d = "' // TRIM( mixing_length_1d ) // &
	3112	'" is unknown'
[226]	3113	CALL message( 'check_parameters', 'PA0137', 1, 2, 0, 6, 0 )
[1]	3114	ENDIF
	3115	IF ( TRIM( dissipation_1d ) /= 'as_in_3d_model' .AND. &
	3116	TRIM( dissipation_1d ) /= 'detering' ) THEN
[215]	3117	message_string = 'dissipation_1d = "' // TRIM( dissipation_1d ) // &
	3118	'" is unknown'
[226]	3119	CALL message( 'check_parameters', 'PA0138', 1, 2, 0, 6, 0 )
[1]	3120	ENDIF
	3121
	3122	!
	3123	!-- Set time for the next user defined restart (time_restart is the
	3124	!-- internal parameter for steering restart events)
	3125	IF ( restart_time /= 9999999.9 ) THEN
[291]	3126	IF ( restart_time > time_since_reference_point ) THEN
	3127	time_restart = restart_time
	3128	ENDIF
[1]	3129	ELSE
	3130	!
	3131	!-- In case of a restart run, set internal parameter to default (no restart)
	3132	!-- if the NAMELIST-parameter restart_time is omitted
	3133	time_restart = 9999999.9
	3134	ENDIF
	3135
	3136	!
	3137	!-- Set default value of the time needed to terminate a model run
	3138	IF ( termination_time_needed == -1.0 ) THEN
	3139	IF ( host(1:3) == 'ibm' ) THEN
	3140	termination_time_needed = 300.0
	3141	ELSE
	3142	termination_time_needed = 35.0
	3143	ENDIF
	3144	ENDIF
	3145
	3146	!
	3147	!-- Check the time needed to terminate a model run
	3148	IF ( host(1:3) == 't3e' ) THEN
	3149	!
	3150	!-- Time needed must be at least 30 seconds on all CRAY machines, because
	3151	!-- MPP_TREMAIN gives the remaining CPU time only in steps of 30 seconds
	3152	IF ( termination_time_needed <= 30.0 ) THEN
[215]	3153	WRITE( message_string, * ) 'termination_time_needed = ', &
	3154	termination_time_needed, ' must be > 30.0 on host "', &
	3155	TRIM( host ), '"'
[226]	3156	CALL message( 'check_parameters', 'PA0139', 1, 2, 0, 6, 0 )
[1]	3157	ENDIF
	3158	ELSEIF ( host(1:3) == 'ibm' ) THEN
	3159	!
	3160	!-- On IBM-regatta machines the time should be at least 300 seconds,
	3161	!-- because the job time consumed before executing palm (for compiling,
	3162	!-- copying of files, etc.) has to be regarded
	3163	IF ( termination_time_needed < 300.0 ) THEN
[215]	3164	WRITE( message_string, * ) 'termination_time_needed = ', &
	3165	termination_time_needed, ' should be >= 300.0 on host "', &
	3166	TRIM( host ), '"'
[226]	3167	CALL message( 'check_parameters', 'PA0140', 1, 2, 0, 6, 0 )
[1]	3168	ENDIF
	3169	ENDIF
	3170
[217]	3171	!
[240]	3172	!-- Check pressure gradient conditions
	3173	IF ( dp_external .AND. conserve_volume_flow ) THEN
[388]	3174	WRITE( message_string, * ) 'Both dp_external and conserve_volume_flo', &
	3175	'w are .TRUE. but one of them must be .FALSE.'
[240]	3176	CALL message( 'check_parameters', 'PA0150', 1, 2, 0, 6, 0 )
	3177	ENDIF
	3178	IF ( dp_external ) THEN
	3179	IF ( dp_level_b < zu(nzb) .OR. dp_level_b > zu(nzt) ) THEN
	3180	WRITE( message_string, * ) 'dp_level_b = ', dp_level_b, ' is out ', &
	3181	' of range'
	3182	CALL message( 'check_parameters', 'PA0151', 1, 2, 0, 6, 0 )
	3183	ENDIF
	3184	IF ( .NOT. ANY( dpdxy /= 0.0 ) ) THEN
[388]	3185	WRITE( message_string, * ) 'dp_external is .TRUE. but dpdxy is ze', &
	3186	'ro, i.e. the external pressure gradient & will not be applied'
[240]	3187	CALL message( 'check_parameters', 'PA0152', 0, 1, 0, 6, 0 )
	3188	ENDIF
	3189	ENDIF
	3190	IF ( ANY( dpdxy /= 0.0 ) .AND. .NOT. dp_external ) THEN
	3191	WRITE( message_string, * ) 'dpdxy is nonzero but dp_external is ', &
	3192	'.FALSE., i.e. the external pressure gradient & will not be applied'
	3193	CALL message( 'check_parameters', 'PA0153', 0, 1, 0, 6, 0 )
	3194	ENDIF
[241]	3195	IF ( conserve_volume_flow ) THEN
	3196	IF ( TRIM( conserve_volume_flow_mode ) == 'default' ) THEN
[667]	3197
	3198	conserve_volume_flow_mode = 'initial_profiles'
	3199
[241]	3200	ELSEIF ( TRIM( conserve_volume_flow_mode ) /= 'initial_profiles' .AND. &
	3201	TRIM( conserve_volume_flow_mode ) /= 'inflow_profile' .AND. &
	3202	TRIM( conserve_volume_flow_mode ) /= 'bulk_velocity' ) THEN
	3203	WRITE( message_string, * ) 'unknown conserve_volume_flow_mode: ', &
	3204	conserve_volume_flow_mode
	3205	CALL message( 'check_parameters', 'PA0154', 1, 2, 0, 6, 0 )
	3206	ENDIF
[667]	3207	IF ( (bc_lr /= 'cyclic' .OR. bc_ns /= 'cyclic') .AND. &
	3208	TRIM( conserve_volume_flow_mode ) == 'bulk_velocity' ) THEN
	3209	WRITE( message_string, * ) 'non-cyclic boundary conditions ', &
	3210	'require conserve_volume_flow_mode = ''initial_profiles'''
[241]	3211	CALL message( 'check_parameters', 'PA0155', 1, 2, 0, 6, 0 )
	3212	ENDIF
	3213	IF ( bc_lr == 'cyclic' .AND. bc_ns == 'cyclic' .AND. &
	3214	TRIM( conserve_volume_flow_mode ) == 'inflow_profile' ) THEN
	3215	WRITE( message_string, * ) 'cyclic boundary conditions ', &
[667]	3216	'require conserve_volume_flow_mode = ''initial_profiles''', &
[241]	3217	' or ''bulk_velocity'''
	3218	CALL message( 'check_parameters', 'PA0156', 1, 2, 0, 6, 0 )
	3219	ENDIF
	3220	ENDIF
	3221	IF ( ( u_bulk /= 0.0 .OR. v_bulk /= 0.0 ) .AND. &
	3222	( .NOT. conserve_volume_flow .OR. &
	3223	TRIM( conserve_volume_flow_mode ) /= 'bulk_velocity' ) ) THEN
	3224	WRITE( message_string, * ) 'nonzero bulk velocity requires ', &
[667]	3225	'conserve_volume_flow = .T. and ', &
[241]	3226	'conserve_volume_flow_mode = ''bulk_velocity'''
	3227	CALL message( 'check_parameters', 'PA0157', 1, 2, 0, 6, 0 )
	3228	ENDIF
[240]	3229
	3230	!
[264]	3231	!-- Check particle attributes
	3232	IF ( particle_color /= 'none' ) THEN
	3233	IF ( particle_color /= 'absuv' .AND. particle_color /= 'pt*' .AND. &
	3234	particle_color /= 'z' ) THEN
	3235	message_string = 'illegal value for parameter particle_color: ' // &
	3236	TRIM( particle_color)
	3237	CALL message( 'check_parameters', 'PA0313', 1, 2, 0, 6, 0 )
	3238	ELSE
	3239	IF ( color_interval(2) <= color_interval(1) ) THEN
	3240	message_string = 'color_interval(2) <= color_interval(1)'
	3241	CALL message( 'check_parameters', 'PA0315', 1, 2, 0, 6, 0 )
	3242	ENDIF
	3243	ENDIF
	3244	ENDIF
	3245
	3246	IF ( particle_dvrpsize /= 'none' ) THEN
	3247	IF ( particle_dvrpsize /= 'absw' ) THEN
	3248	message_string = 'illegal value for parameter particle_dvrpsize:' // &
	3249	' ' // TRIM( particle_color)
	3250	CALL message( 'check_parameters', 'PA0314', 1, 2, 0, 6, 0 )
	3251	ELSE
	3252	IF ( dvrpsize_interval(2) <= dvrpsize_interval(1) ) THEN
	3253	message_string = 'dvrpsize_interval(2) <= dvrpsize_interval(1)'
	3254	CALL message( 'check_parameters', 'PA0316', 1, 2, 0, 6, 0 )
	3255	ENDIF
	3256	ENDIF
	3257	ENDIF
	3258
	3259	!
[217]	3260	!-- Check &userpar parameters
	3261	CALL user_check_parameters
[1]	3262
[217]	3263
[667]	3264
[1]	3265	END SUBROUTINE check_parameters

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

Context Navigation

source: palm/trunk/SOURCE/check_parameters.f90 @ 707

Download in other formats: