Home

Context Navigation

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

Last change on this file since 72 was 72, checked in by raasch, 18 years ago
preliminary changes for precipitation output
Property svn:keywords set to `Id`
File size: 86.3 KB

Line
1	SUBROUTINE check_parameters
2
3	!------------------------------------------------------------------------------!
4	! Actual revisions:
5	! -----------------
6	! "by_user" allowed as initializing action, -data_output_ts,
7	! leapfrog with non-flat topography not allowed any more, loop_optimization
8	! and pt_reference are checked,
9	! output of precipitation amount/rate and roughnes length + check
10	! possible negative humidities are avoided in initial profile
11	!
12	! Former revisions:
13	! -----------------
14	! $Id: check_parameters.f90 72 2007-03-19 08:20:46Z raasch $
15	!
16	! 20 2007-02-26 00:12:32Z raasch
17	! Temperature and humidity gradients at top are now calculated for nzt+1,
18	! top_heatflux and respective boundary condition bc_pt_t is checked
19	!
20	! RCS Log replace by Id keyword, revision history cleaned up
21	!
22	! Revision 1.61 2006/08/04 14:20:25 raasch
23	! do2d_unit and do3d_unit now defined as 2d-arrays, check of
24	! use_upstream_for_tke, default value for dt_dopts,
25	! generation of file header moved from routines palm and header to here
26	!
27	! Revision 1.1 1997/08/26 06:29:23 raasch
28	! Initial revision
29	!
30	!
31	! Description:
32	! ------------
33	! Check control parameters and deduce further quantities.
34	!------------------------------------------------------------------------------!
35
36	USE arrays_3d
37	USE constants
38	USE control_parameters
39	USE grid_variables
40	USE indices
41	USE model_1d
42	USE netcdf_control
43	USE particle_attributes
44	USE pegrid
45	USE profil_parameter
46	USE statistics
47	USE transpose_indices
48
49	IMPLICIT NONE
50
51	CHARACTER (LEN=1) :: sq
52	CHARACTER (LEN=6) :: var
53	CHARACTER (LEN=7) :: unit
54	CHARACTER (LEN=8) :: date
55	CHARACTER (LEN=10) :: time
56	CHARACTER (LEN=100) :: action
57
58	INTEGER :: i, ilen, intervals, iter, j, k, nnxh, nnyh, position, prec
59	LOGICAL :: found, ldum
60	REAL :: gradient, maxn, maxp
61
62
63	!
64	!-- Warning, if host is not set
65	IF ( host(1:1) == ' ' ) THEN
66	IF ( myid == 0 ) THEN
67	PRINT*, '+++ WARNING: check_parameters:'
68	PRINT*, ' "host" is not set. Please check that environment', &
69	' variable "localhost"'
70	PRINT*, ' is set before running PALM'
71	ENDIF
72	ENDIF
73
74	!
75	!-- Generate the file header which is used as a header for most of PALM's
76	!-- output files
77	CALL DATE_AND_TIME( date, time )
78	run_date = date(7:8)//'-'//date(5:6)//'-'//date(3:4)
79	run_time = time(1:2)//':'//time(3:4)//':'//time(5:6)
80
81	WRITE ( run_description_header, '(A,2X,A,A,A,I2.2,2X,A,A,2X,A,1X,A)' ) &
82	TRIM( version ),'run: ', TRIM( run_identifier ), '.', &
83	runnr, 'host: ', TRIM( host ), run_date, run_time
84
85	!
86	!-- Check the general loop optimization method
87	IF ( loop_optimization == 'default' ) THEN
88	IF ( host(1:3) == 'nec' ) THEN
89	loop_optimization = 'vector'
90	ELSE
91	loop_optimization = 'cache'
92	ENDIF
93	ENDIF
94	IF ( loop_optimization /= 'noopt' .AND. loop_optimization /= 'cache' &
95	.AND. loop_optimization /= 'vector' ) THEN
96	IF ( myid == 0 ) THEN
97	PRINT*, '+++ check_parameters:'
98	PRINT*, ' illegal value given for loop_optimization: ', &
99	TRIM( loop_optimization )
100	ENDIF
101	CALL local_stop
102	ENDIF
103
104	!
105	!-- Check topography setting (check for illegal parameter combinations)
106	IF ( topography /= 'flat' ) THEN
107	action = ' '
108	IF ( scalar_advec /= 'pw-scheme' ) THEN
109	WRITE( action, '(A,A)' ) 'scalar_advec = ', scalar_advec
110	ENDIF
111	IF ( momentum_advec /= 'pw-scheme' ) THEN
112	WRITE( action, '(A,A)' ) 'momentum_advec = ', momentum_advec
113	ENDIF
114	IF ( timestep_scheme(1:8) == 'leapfrog' ) THEN
115	WRITE( action, '(A,A)' ) 'timestep_scheme = ', timestep_scheme
116	ENDIF
117	IF ( psolver == 'multigrid' .OR. psolver == 'sor' ) THEN
118	WRITE( action, '(A,A)' ) 'psolver = ', psolver
119	ENDIF
120	IF ( sloping_surface ) THEN
121	WRITE( action, '(A)' ) 'sloping surface = .TRUE.'
122	ENDIF
123	IF ( galilei_transformation ) THEN
124	WRITE( action, '(A)' ) 'galilei_transformation = .TRUE.'
125	ENDIF
126	IF ( cloud_physics ) THEN
127	WRITE( action, '(A)' ) 'cloud_physics = .TRUE.'
128	ENDIF
129	IF ( cloud_droplets ) THEN
130	WRITE( action, '(A)' ) 'cloud_droplets = .TRUE.'
131	ENDIF
132	IF ( moisture ) THEN
133	WRITE( action, '(A)' ) 'moisture = .TRUE.'
134	ENDIF
135	IF ( .NOT. prandtl_layer ) THEN
136	WRITE( action, '(A)' ) 'prandtl_layer = .FALSE.'
137	ENDIF
138	IF ( action /= ' ' ) THEN
139	IF ( myid == 0 ) THEN
140	PRINT*, '+++ check_parameters:'
141	PRINT*, ' a non-flat topography does not allow ', TRIM( action )
142	ENDIF
143	CALL local_stop
144	ENDIF
145	ENDIF
146	!
147	!-- Check whether there are any illegal values
148	!-- Pressure solver:
149	IF ( psolver /= 'poisfft' .AND. psolver /= 'poisfft_hybrid' .AND. &
150	psolver /= 'sor' .AND. psolver /= 'multigrid' ) THEN
151	IF ( myid == 0 ) THEN
152	PRINT*, '+++ check_parameters:'
153	PRINT*, ' unknown solver for perturbation pressure: psolver=', &
154	psolver
155	ENDIF
156	CALL local_stop
157	ENDIF
158
159	#if defined( __parallel )
160	IF ( psolver == 'poisfft_hybrid' .AND. pdims(2) /= 1 ) THEN
161	IF ( myid == 0 ) THEN
162	PRINT*, '+++ check_parameters:'
163	PRINT*, ' psolver="', TRIM( psolver ), '" only works for a ', &
164	'1d domain-decomposition along x'
165	PRINT*, ' please do not set npey/=1 in the parameter file'
166	ENDIF
167	CALL local_stop
168	ENDIF
169	IF ( ( psolver == 'poisfft_hybrid' .OR. psolver == 'multigrid' ) .AND. &
170	( nxra > nxr .OR. nyna > nyn .OR. nza > nz ) ) THEN
171	IF ( myid == 0 ) THEN
172	PRINT*, '+++ check_parameters:'
173	PRINT*, ' psolver="', TRIM( psolver ), '" does not work for ', &
174	'subdomains with unequal size'
175	PRINT*, ' please set grid_matching = ''strict'' in the parameter',&
176	' file'
177	ENDIF
178	CALL local_stop
179	ENDIF
180	#else
181	IF ( psolver == 'poisfft_hybrid' ) THEN
182	IF ( myid == 0 ) THEN
183	PRINT*, '+++ check_parameters:'
184	PRINT*, ' psolver="', TRIM( psolver ), '" only works for a ', &
185	'parallel environment'
186	ENDIF
187	CALL local_stop
188	ENDIF
189	#endif
190
191	IF ( psolver == 'multigrid' ) THEN
192	IF ( cycle_mg == 'w' ) THEN
193	gamma_mg = 2
194	ELSEIF ( cycle_mg == 'v' ) THEN
195	gamma_mg = 1
196	ELSE
197	IF ( myid == 0 ) THEN
198	PRINT*, '+++ check_parameters:'
199	PRINT*, ' unknown multigrid cycle: cycle_mg=', cycle_mg
200	ENDIF
201	CALL local_stop
202	ENDIF
203	ENDIF
204
205	IF ( fft_method /= 'singleton-algorithm' .AND. &
206	fft_method /= 'temperton-algorithm' .AND. &
207	fft_method /= 'system-specific' ) THEN
208	IF ( myid == 0 ) THEN
209	PRINT*, '+++ check_parameters:'
210	PRINT*, ' unknown fft-algorithm: fft_method=', fft_method
211	ENDIF
212	CALL local_stop
213	ENDIF
214
215	!
216	!-- Advection schemes:
217	IF ( momentum_advec /= 'pw-scheme' .AND. momentum_advec /= 'ups-scheme' ) &
218	THEN
219	IF ( myid == 0 ) PRINT*, '+++ check_parameters: unknown advection ', &
220	'scheme: momentum_advec=', momentum_advec
221	CALL local_stop
222	ENDIF
223	IF ( ( momentum_advec == 'ups-scheme' .OR. scalar_advec == 'ups-scheme' )&
224	.AND. timestep_scheme /= 'euler' ) THEN
225	IF ( myid == 0 ) PRINT*, '+++ check_parameters: momentum_advec=', &
226	momentum_advec, ' is not allowed with ', &
227	'timestep_scheme=', timestep_scheme
228	CALL local_stop
229	ENDIF
230
231	IF ( scalar_advec /= 'pw-scheme' .AND. scalar_advec /= 'bc-scheme' .AND.&
232	scalar_advec /= 'ups-scheme' ) THEN
233	IF ( myid == 0 ) PRINT*, '+++ check_parameters: unknown advection ', &
234	'scheme: scalar_advec=', scalar_advec
235	CALL local_stop
236	ENDIF
237
238	IF ( use_sgs_for_particles .AND. .NOT. use_upstream_for_tke ) THEN
239	use_upstream_for_tke = .TRUE.
240	IF ( myid == 0 ) THEN
241	PRINT*, '+++ WARNING: check_parameters: use_upstream_for_tke set ', &
242	'.TRUE. because use_sgs_for_particles = .TRUE.'
243	ENDIF
244	ENDIF
245
246	IF ( use_upstream_for_tke .AND. timestep_scheme(1:8) == 'leapfrog' ) THEN
247	IF ( myid == 0 ) THEN
248	PRINT*, '+++ check_parameters: use_upstream_for_tke = .TRUE. ', &
249	'not allowed with timestep_scheme=', timestep_scheme
250	ENDIF
251	CALL local_stop
252	ENDIF
253
254	!
255	!-- Timestep schemes:
256	SELECT CASE ( TRIM( timestep_scheme ) )
257
258	CASE ( 'euler' )
259	intermediate_timestep_count_max = 1
260	asselin_filter_factor = 0.0
261
262	CASE ( 'leapfrog', 'leapfrog+euler' )
263	intermediate_timestep_count_max = 1
264
265	CASE ( 'runge-kutta-2' )
266	intermediate_timestep_count_max = 2
267	asselin_filter_factor = 0.0
268
269	CASE ( 'runge-kutta-3' )
270	intermediate_timestep_count_max = 3
271	asselin_filter_factor = 0.0
272
273	CASE DEFAULT
274	IF ( myid == 0 ) PRINT*, '+++ check_parameters: unknown timestep ',&
275	'scheme: timestep_scheme=', timestep_scheme
276	CALL local_stop
277
278	END SELECT
279
280	IF ( scalar_advec == 'ups-scheme' .AND. timestep_scheme(1:5) == 'runge' )&
281	THEN
282	IF ( myid == 0 ) THEN
283	PRINT*, '+++ check_parameters: scalar advection scheme "', &
284	TRIM( scalar_advec ), '"'
285	PRINT*, ' does not work with timestep_scheme "', &
286	TRIM( timestep_scheme ), '"'
287	ENDIF
288	CALL local_stop
289	ENDIF
290
291	IF ( momentum_advec /= 'pw-scheme' .AND. timestep_scheme(1:5) == 'runge' ) &
292	THEN
293	IF ( myid == 0 ) THEN
294	PRINT*, '+++ check_parameters: momentum advection scheme "', &
295	TRIM( momentum_advec ), '"'
296	PRINT*, ' does not work with timestep_scheme "', &
297	TRIM( timestep_scheme ), '"'
298	ENDIF
299	CALL local_stop
300	ENDIF
301
302
303	IF ( initializing_actions == ' ' ) THEN
304	IF ( myid == 0 ) THEN
305	PRINT*, '+++ check parameters:'
306	PRINT*, ' no value found for initializing_actions'
307	ENDIF
308	CALL local_stop
309	ENDIF
310
311	IF ( TRIM( initializing_actions ) /= 'read_restart_data' ) THEN
312	!
313	!-- No model continuation run; several initialising actions are possible
314	action = initializing_actions
315	DO WHILE ( TRIM( action ) /= '' )
316	position = INDEX( action, ' ' )
317	SELECT CASE ( action(1:position-1) )
318
319	CASE ( 'set_constant_profiles', 'set_1d-model_profiles', &
320	'by_user', 'initialize_vortex', 'initialize_ptanom' )
321	action = action(position+1:)
322
323	CASE DEFAULT
324	IF ( myid == 0 ) PRINT*, '+++ check_parameters: initializi', &
325	'ng_action unkown or not allowed: action = "', &
326	TRIM(action), '"'
327	CALL local_stop
328
329	END SELECT
330	ENDDO
331	ENDIF
332	IF ( INDEX( initializing_actions, 'set_constant_profiles' ) /= 0 .AND. &
333	INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 ) THEN
334	IF ( myid == 0 ) PRINT*, '+++ check_parameters: initializing_actions', &
335	'"set_constant_profiles" and "set_1d-model_profiles" are not', &
336	' allowed simultaneously'
337	CALL local_stop
338	ENDIF
339	IF ( INDEX( initializing_actions, 'set_constant_profiles' ) /= 0 .AND. &
340	INDEX( initializing_actions, 'by_user' ) /= 0 ) THEN
341	IF ( myid == 0 ) PRINT*, '+++ check_parameters: initializing_actions', &
342	'"set_constant_profiles" and "by_user" are not', &
343	' allowed simultaneously'
344	CALL local_stop
345	ENDIF
346	IF ( INDEX( initializing_actions, 'by_user' ) /= 0 .AND. &
347	INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 ) THEN
348	IF ( myid == 0 ) PRINT*, '+++ check_parameters: initializing_actions', &
349	'"by_user" and "set_1d-model_profiles" are not', &
350	' allowed simultaneously'
351	CALL local_stop
352	ENDIF
353
354	IF ( cloud_physics .AND. .NOT. moisture ) THEN
355	IF ( myid == 0 ) PRINT*, '+++ check_parameters: cloud_physics =', &
356	cloud_physics, ' is not allowed with ', &
357	'moisture =', moisture
358	CALL local_stop
359	ENDIF
360
361	IF ( precipitation .AND. .NOT. cloud_physics ) THEN
362	IF ( myid == 0 ) PRINT*, '+++ check_parameters: precipitation =', &
363	precipitation, ' is not allowed with ', &
364	'cloud_physics =', cloud_physics
365	CALL local_stop
366	ENDIF
367
368	IF ( moisture .AND. sloping_surface ) THEN
369	IF ( myid == 0 ) PRINT*, '+++ check_parameters: moisture = TRUE', &
370	'and hang = TRUE are not', &
371	' allowed simultaneously'
372	CALL local_stop
373	ENDIF
374
375	IF ( moisture .AND. scalar_advec == 'ups-scheme' ) THEN
376	IF ( myid == 0 ) PRINT*, '+++ check_parameters: UPS-scheme', &
377	'is not implemented for moisture'
378	CALL local_stop
379	ENDIF
380
381	IF ( passive_scalar .AND. moisture ) THEN
382	IF ( myid == 0 ) PRINT*, '+++ check_parameters: moisture = TRUE and', &
383	'passive_scalar = TRUE is not allowed ', &
384	'simultaneously'
385	CALL local_stop
386	ENDIF
387
388	IF ( passive_scalar .AND. scalar_advec == 'ups-scheme' ) THEN
389	IF ( myid == 0 ) PRINT*, '+++ check_parameters: UPS-scheme', &
390	'is not implemented for passive_scalar'
391	CALL local_stop
392	ENDIF
393
394	IF ( grid_matching /= 'strict' .AND. grid_matching /= 'match' ) THEN
395	IF ( myid == 0 ) PRINT*, '+++ check_parameters: illegal value "', &
396	TRIM( grid_matching ), &
397	'" found for parameter grid_matching'
398	CALL local_stop
399	ENDIF
400
401	!
402	!-- In case of no model continuation run, check initialising parameters and
403	!-- deduce further quantities
404	IF ( TRIM( initializing_actions ) /= 'read_restart_data' ) THEN
405
406	!
407	!-- Initial profiles for 1D and 3D model, respectively
408	u_init = ug_surface
409	v_init = vg_surface
410	pt_init = pt_surface
411	IF ( moisture ) q_init = q_surface
412	IF ( passive_scalar ) q_init = s_surface
413
414	!
415	!--
416	!-- If required, compute initial profile of the geostrophic wind
417	!-- (component ug)
418	i = 1
419	gradient = 0.0
420	ug_vertical_gradient_level_ind(1) = 0
421	ug(0) = ug_surface
422	DO k = 1, nzt+1
423	IF ( ug_vertical_gradient_level(i) < zu(k) .AND. &
424	ug_vertical_gradient_level(i) >= 0.0 ) THEN
425	gradient = ug_vertical_gradient(i) / 100.0
426	ug_vertical_gradient_level_ind(i) = k - 1
427	i = i + 1
428	IF ( i > 10 ) THEN
429	IF ( myid == 0 ) THEN
430	PRINT*, '+++ check_parameters: upper bound 10 of array', &
431	' "ug_vertical_gradient_level_ind" exceeded'
432	ENDIF
433	CALL local_stop
434	ENDIF
435	ENDIF
436	IF ( gradient /= 0.0 ) THEN
437	IF ( k /= 1 ) THEN
438	ug(k) = ug(k-1) + dzu(k) * gradient
439	ELSE
440	ug(k) = ug_surface + 0.5 * dzu(k) * gradient
441	ENDIF
442	ELSE
443	ug(k) = ug(k-1)
444	ENDIF
445	ENDDO
446
447	u_init = ug
448
449	!
450	!-- In case of no given gradients for ug, choose a vanishing gradient
451	IF ( ug_vertical_gradient_level(1) == -1.0 ) THEN
452	ug_vertical_gradient_level(1) = 0.0
453	ENDIF
454
455	!
456	!--
457	!-- If required, compute initial profile of the geostrophic wind
458	!-- (component vg)
459	i = 1
460	gradient = 0.0
461	vg_vertical_gradient_level_ind(1) = 0
462	vg(0) = vg_surface
463	DO k = 1, nzt+1
464	IF ( vg_vertical_gradient_level(i) < zu(k) .AND. &
465	vg_vertical_gradient_level(i) >= 0.0 ) THEN
466	gradient = vg_vertical_gradient(i) / 100.0
467	vg_vertical_gradient_level_ind(i) = k - 1
468	i = i + 1
469	IF ( i > 10 ) THEN
470	IF ( myid == 0 ) THEN
471	PRINT*, '+++ check_parameters: upper bound 10 of array', &
472	' "vg_vertical_gradient_level_ind" exceeded'
473	ENDIF
474	CALL local_stop
475	ENDIF
476	ENDIF
477	IF ( gradient /= 0.0 ) THEN
478	IF ( k /= 1 ) THEN
479	vg(k) = vg(k-1) + dzu(k) * gradient
480	ELSE
481	vg(k) = vg_surface + 0.5 * dzu(k) * gradient
482	ENDIF
483	ELSE
484	vg(k) = vg(k-1)
485	ENDIF
486	ENDDO
487
488	v_init = vg
489
490	!
491	!-- In case of no given gradients for vg, choose a vanishing gradient
492	IF ( vg_vertical_gradient_level(1) == -1.0 ) THEN
493	vg_vertical_gradient_level(1) = 0.0
494	ENDIF
495
496	!
497	!-- If required, compute initial temperature profile using the given
498	!-- temperature gradient
499	i = 1
500	gradient = 0.0
501	pt_vertical_gradient_level_ind(1) = 0
502	DO k = 1, nzt+1
503	IF ( pt_vertical_gradient_level(i) < zu(k) .AND. &
504	pt_vertical_gradient_level(i) >= 0.0 ) THEN
505	gradient = pt_vertical_gradient(i) / 100.0
506	pt_vertical_gradient_level_ind(i) = k - 1
507	i = i + 1
508	IF ( i > 10 ) THEN
509	IF ( myid == 0 ) THEN
510	PRINT*, '+++ check_parameters: upper bound 10 of array', &
511	' "pt_vertical_gradient_level_ind" exceeded'
512	ENDIF
513	CALL local_stop
514	ENDIF
515	ENDIF
516	IF ( gradient /= 0.0 ) THEN
517	IF ( k /= 1 ) THEN
518	pt_init(k) = pt_init(k-1) + dzu(k) * gradient
519	ELSE
520	pt_init(k) = pt_init(k-1) + 0.5 * dzu(k) * gradient
521	ENDIF
522	ELSE
523	pt_init(k) = pt_init(k-1)
524	ENDIF
525	ENDDO
526
527	!
528	!-- In case of no given temperature gradients, choose gradient of neutral
529	!-- stratification
530	IF ( pt_vertical_gradient_level(1) == -1.0 ) THEN
531	pt_vertical_gradient_level(1) = 0.0
532	ENDIF
533
534	!
535	!-- Store temperature gradient at the top boundary for possile Neumann
536	!-- boundary condition
537	bc_pt_t_val = ( pt_init(nzt+1) - pt_init(nzt) ) / dzu(nzt+1)
538
539	!
540	!-- If required, compute initial humidity or scalar profile using the given
541	!-- humidity/scalar gradient. In case of scalar transport, initially store
542	!-- values of the scalar parameters on humidity parameters
543	IF ( passive_scalar ) THEN
544	bc_q_b = bc_s_b
545	bc_q_t = bc_s_t
546	q_surface = s_surface
547	q_surface_initial_change = s_surface_initial_change
548	q_vertical_gradient = s_vertical_gradient
549	q_vertical_gradient_level = s_vertical_gradient_level
550	surface_waterflux = surface_scalarflux
551	ENDIF
552
553	IF ( moisture .OR. passive_scalar ) THEN
554
555	i = 1
556	gradient = 0.0
557	q_vertical_gradient_level_ind(1) = 0
558	DO k = 1, nzt+1
559	IF ( q_vertical_gradient_level(i) < zu(k) .AND. &
560	q_vertical_gradient_level(i) >= 0.0 ) THEN
561	gradient = q_vertical_gradient(i) / 100.0
562	q_vertical_gradient_level_ind(i) = k - 1
563	i = i + 1
564	IF ( i > 10 ) THEN
565	IF ( myid == 0 ) THEN
566	PRINT*, '+++ check_parameters: upper bound 10 of arr', &
567	'ay "q_vertical_gradient_level_ind" exceeded'
568	ENDIF
569	CALL local_stop
570	ENDIF
571	ENDIF
572	IF ( gradient /= 0.0 ) THEN
573	IF ( k /= 1 ) THEN
574	q_init(k) = q_init(k-1) + dzu(k) * gradient
575	ELSE
576	q_init(k) = q_init(k-1) + 0.5 * dzu(k) * gradient
577	ENDIF
578	ELSE
579	q_init(k) = q_init(k-1)
580	ENDIF
581	!
582	!-- Avoid negative humidities
583	IF ( q_init(k) < 0.0 ) THEN
584	q_init(k) = 0.0
585	ENDIF
586	ENDDO
587
588	!
589	!-- In case of no given humidity gradients, choose zero gradient
590	!-- conditions
591	IF ( q_vertical_gradient_level(1) == -1.0 ) THEN
592	q_vertical_gradient_level(1) = 0.0
593	ENDIF
594
595	!
596	!-- Store humidity gradient at the top boundary for possile Neumann
597	!-- boundary condition
598	bc_q_t_val = ( q_init(nzt+1) - q_init(nzt) ) / dzu(nzt+1)
599
600	ENDIF
601
602	ENDIF
603
604	!
605	!-- Compute Coriolis parameter
606	f = 2.0 * omega * SIN( phi / 180.0 * pi )
607	fs = 2.0 * omega * COS( phi / 180.0 * pi )
608
609	!
610	!-- Reference temperature to be used in buoyancy terms
611	IF ( pt_reference /= 9999999.9 ) use_pt_reference = .TRUE.
612
613	!
614	!-- In case of a given slope, compute the relevant quantities
615	IF ( alpha_surface /= 0.0 ) THEN
616	IF ( ABS( alpha_surface ) > 90.0 ) THEN
617	IF ( myid == 0 ) PRINT*, '+++ check_parameters: ABS( alpha_surface',&
618	'=', alpha_surface, ' ) must be < 90.0'
619	CALL local_stop
620	ENDIF
621	sloping_surface = .TRUE.
622	cos_alpha_surface = COS( alpha_surface / 180.0 * pi )
623	sin_alpha_surface = SIN( alpha_surface / 180.0 * pi )
624	ENDIF
625
626	!
627	!-- Check time step and cfl_factor
628	IF ( dt /= -1.0 ) THEN
629	IF ( dt <= 0.0 .AND. dt /= -1.0 ) THEN
630	IF ( myid == 0 ) PRINT*, '+++ check_parameters: dt=', dt, ' <= 0.0'
631	CALL local_stop
632	ENDIF
633	dt_3d = dt
634	dt_fixed = .TRUE.
635	ENDIF
636
637	IF ( cfl_factor <= 0.0 .OR. cfl_factor > 1.0 ) THEN
638	IF ( cfl_factor == -1.0 ) THEN
639	IF ( momentum_advec == 'ups-scheme' .OR. &
640	scalar_advec == 'ups-scheme' ) THEN
641	cfl_factor = 0.8
642	ELSE
643	IF ( timestep_scheme == 'runge-kutta-2' ) THEN
644	cfl_factor = 0.8
645	ELSEIF ( timestep_scheme == 'runge-kutta-3' ) THEN
646	cfl_factor = 0.9
647	ELSE
648	cfl_factor = 0.1
649	ENDIF
650	ENDIF
651	ELSE
652	IF ( myid == 0 ) THEN
653	PRINT*, '+++ check_parameters: cfl_factor=', cfl_factor, &
654	' out of range'
655	PRINT*, '+++ 0.0 < cfl_factor <= 1.0 is required'
656	ENDIF
657	CALL local_stop
658	ENDIF
659	ENDIF
660
661	!
662	!-- Store simulated time at begin
663	simulated_time_at_begin = simulated_time
664
665	!
666	!-- Set wind speed in the Galilei-transformed system
667	IF ( galilei_transformation ) THEN
668	IF ( use_ug_for_galilei_tr .AND. &
669	ug_vertical_gradient_level(1) == 0.0 .AND. &
670	vg_vertical_gradient_level(1) == 0.0 ) THEN
671	u_gtrans = ug_surface
672	v_gtrans = vg_surface
673	ELSEIF ( use_ug_for_galilei_tr .AND. &
674	ug_vertical_gradient_level(1) /= 0.0 ) THEN
675	IF ( myid == 0 ) THEN
676	PRINT*, '+++ check_parameters:'
677	PRINT*, ' baroclinicity (ug) not allowed'
678	PRINT*, ' simultaneously with galilei transformation'
679	ENDIF
680	CALL local_stop
681	ELSEIF ( use_ug_for_galilei_tr .AND. &
682	vg_vertical_gradient_level(1) /= 0.0 ) THEN
683	IF ( myid == 0 ) THEN
684	PRINT*, '+++ check_parameters:'
685	PRINT*, ' baroclinicity (vg) not allowed'
686	PRINT*, ' simultaneously with galilei transformation'
687	ENDIF
688	CALL local_stop
689	ELSE
690	IF ( myid == 0 ) THEN
691	PRINT*, '+++ WARNING: check_parameters:'
692	PRINT*, ' variable translation speed used for galilei-tran' // &
693	'sformation, which'
694	PRINT*, ' may cause instabilities in stably stratified regions'
695	ENDIF
696	ENDIF
697	ENDIF
698
699	!
700	!-- In case of using a prandtl-layer, calculated (or prescribed) surface
701	!-- fluxes have to be used in the diffusion-terms
702	IF ( prandtl_layer ) use_surface_fluxes = .TRUE.
703
704	!
705	!-- Check boundary conditions and set internal variables:
706	!-- Lateral boundary conditions
707	IF ( bc_lr /= 'cyclic' .AND. bc_lr /= 'dirichlet/neumann' .AND. &
708	bc_lr /= 'neumann/dirichlet' ) THEN
709	IF ( myid == 0 ) THEN
710	PRINT*, '+++ check_parameters:'
711	PRINT*, ' unknown boundary condition: bc_lr = ', bc_lr
712	ENDIF
713	CALL local_stop
714	ENDIF
715	IF ( bc_ns /= 'cyclic' .AND. bc_ns /= 'dirichlet/neumann' .AND. &
716	bc_ns /= 'neumann/dirichlet' ) THEN
717	IF ( myid == 0 ) THEN
718	PRINT*, '+++ check_parameters:'
719	PRINT*, ' unknown boundary condition: bc_ns = ', bc_ns
720	ENDIF
721	CALL local_stop
722	ENDIF
723
724	!
725	!-- Non-cyclic lateral boundaries require the multigrid method and Piascek-
726	!-- Willimas advection scheme. Several schemes and tools do not work with
727	!-- non-cyclic boundary conditions.
728	IF ( bc_lr /= 'cyclic' .OR. bc_ns /= 'cyclic' ) THEN
729	IF ( psolver /= 'multigrid' ) THEN
730	IF ( myid == 0 ) THEN
731	PRINT*, '+++ check_parameters:'
732	PRINT*, ' non-cyclic lateral boundaries do not allow', &
733	' psolver = ', psolver
734	ENDIF
735	CALL local_stop
736	ENDIF
737	IF ( momentum_advec /= 'pw-scheme' ) THEN
738	IF ( myid == 0 ) THEN
739	PRINT*, '+++ check_parameters:'
740	PRINT*, ' non-cyclic lateral boundaries do not allow', &
741	' momentum_advec = ', momentum_advec
742	ENDIF
743	CALL local_stop
744	ENDIF
745	IF ( scalar_advec /= 'pw-scheme' ) THEN
746	IF ( myid == 0 ) THEN
747	PRINT*, '+++ check_parameters:'
748	PRINT*, ' non-cyclic lateral boundaries do not allow', &
749	' scalar_advec = ', scalar_advec
750	ENDIF
751	CALL local_stop
752	ENDIF
753	IF ( galilei_transformation ) THEN
754	IF ( myid == 0 ) THEN
755	PRINT*, '+++ check_parameters:'
756	PRINT*, ' non-cyclic lateral boundaries do not allow', &
757	' galilei_transformation = .T.'
758	ENDIF
759	CALL local_stop
760	ENDIF
761	IF ( conserve_volume_flow ) THEN
762	IF ( myid == 0 ) THEN
763	PRINT*, '+++ check_parameters:'
764	PRINT*, ' non-cyclic lateral boundaries do not allow', &
765	' conserve_volume_flow = .T.'
766	ENDIF
767	CALL local_stop
768	ENDIF
769	ENDIF
770
771	!
772	!-- Bottom boundary condition for the turbulent Kinetic energy
773	IF ( bc_e_b == 'neumann' ) THEN
774	ibc_e_b = 1
775	IF ( adjust_mixing_length .AND. prandtl_layer ) THEN
776	IF ( myid == 0 ) THEN
777	PRINT*, '+++ WARNING: check_parameters:'
778	PRINT*, ' adjust_mixing_length = TRUE and bc_e_b = ', bc_e_b
779	ENDIF
780	ENDIF
781	ELSEIF ( bc_e_b == '(u)*2+neumann' ) THEN
782	ibc_e_b = 2
783	IF ( .NOT. adjust_mixing_length .AND. prandtl_layer ) THEN
784	IF ( myid == 0 ) THEN
785	PRINT*, '+++ WARNING: check_parameters:'
786	PRINT*, ' adjust_mixing_length = FALSE and bc_e_b = ', bc_e_b
787	ENDIF
788	ENDIF
789	IF ( .NOT. prandtl_layer ) THEN
790	bc_e_b = 'neumann'
791	ibc_e_b = 1
792	IF ( myid == 0 ) THEN
793	PRINT*, '+++ WARNING: check_parameters:'
794	PRINT*, ' boundary condition bc_e_b changed to "', bc_e_b, '"'
795	ENDIF
796	ENDIF
797	ELSE
798	IF ( myid == 0 ) THEN
799	PRINT*, '+++ check_parameters:'
800	PRINT*, ' unknown boundary condition: bc_e_b = ', bc_e_b
801	ENDIF
802	CALL local_stop
803	ENDIF
804
805	!
806	!-- Boundary conditions for perturbation pressure
807	IF ( bc_p_b == 'dirichlet' ) THEN
808	ibc_p_b = 0
809	ELSEIF ( bc_p_b == 'neumann' ) THEN
810	ibc_p_b = 1
811	ELSEIF ( bc_p_b == 'neumann+inhomo' ) THEN
812	ibc_p_b = 2
813	ELSE
814	IF ( myid == 0 ) THEN
815	PRINT*, '+++ check_parameters:'
816	PRINT*, ' unknown boundary condition: bc_p_b = ', bc_p_b
817	ENDIF
818	CALL local_stop
819	ENDIF
820	IF ( ibc_p_b == 2 .AND. .NOT. prandtl_layer ) THEN
821	IF ( myid == 0 ) THEN
822	PRINT*, '+++ check_parameters:'
823	PRINT*, ' boundary condition: bc_p_b = ', TRIM( bc_p_b ), &
824	' not allowed with'
825	PRINT*, ' prandtl_layer = .FALSE.'
826	ENDIF
827	CALL local_stop
828	ENDIF
829	IF ( bc_p_t == 'dirichlet' ) THEN
830	ibc_p_t = 0
831	ELSEIF ( bc_p_t == 'neumann' ) THEN
832	ibc_p_t = 1
833	ELSE
834	IF ( myid == 0 ) THEN
835	PRINT*, '+++ check_parameters:'
836	PRINT*, ' unknown boundary condition: bc_p_t = ', bc_p_t
837	ENDIF
838	CALL local_stop
839	ENDIF
840
841	!
842	!-- Boundary conditions for potential temperature
843	IF ( bc_pt_b == 'dirichlet' ) THEN
844	ibc_pt_b = 0
845	ELSEIF ( bc_pt_b == 'neumann' ) THEN
846	ibc_pt_b = 1
847	ELSE
848	IF ( myid == 0 ) THEN
849	PRINT*, '+++ check_parameters:'
850	PRINT*, ' unknown boundary condition: bc_pt_b = ', bc_pt_b
851	ENDIF
852	CALL local_stop
853	ENDIF
854	IF ( bc_pt_t == 'dirichlet' ) THEN
855	ibc_pt_t = 0
856	ELSEIF ( bc_pt_t == 'neumann' ) THEN
857	ibc_pt_t = 1
858	ELSEIF ( bc_pt_t == 'initial_gradient' ) THEN
859	ibc_pt_t = 2
860	ELSE
861	IF ( myid == 0 ) THEN
862	PRINT*, '+++ check_parameters:'
863	PRINT*, ' unknown boundary condition: bc_pt_t = ', bc_pt_t
864	ENDIF
865	CALL local_stop
866	ENDIF
867
868	IF ( surface_heatflux == 9999999.9 ) constant_heatflux = .FALSE.
869	IF ( top_heatflux == 9999999.9 ) constant_top_heatflux = .FALSE.
870
871	!
872	!-- A given surface temperature implies Dirichlet boundary condition for
873	!-- temperature. In this case specification of a constant heat flux is
874	!-- forbidden.
875	IF ( ibc_pt_b == 0 .AND. constant_heatflux .AND. &
876	surface_heatflux /= 0.0 ) THEN
877	IF ( myid == 0 ) THEN
878	PRINT*, '+++ check_parameters:'
879	PRINT*, ' boundary_condition: bc_pt_b = ', bc_pt_b
880	PRINT*, ' is not allowed with constant_heatflux = .TRUE.'
881	ENDIF
882	CALL local_stop
883	ENDIF
884	IF ( constant_heatflux .AND. pt_surface_initial_change /= 0.0 ) THEN
885	IF ( myid == 0 ) THEN
886	PRINT*, '+++ check_parameters: constant_heatflux = .TRUE. is not'
887	PRINT*, ' allowed with pt_surface_initial_change (/=0) = ', &
888	pt_surface_initial_change
889	ENDIF
890	CALL local_stop
891	ENDIF
892
893	!
894	!-- A given temperature at the top implies Dirichlet boundary condition for
895	!-- temperature. In this case specification of a constant heat flux is
896	!-- forbidden.
897	IF ( ibc_pt_t == 0 .AND. constant_top_heatflux .AND. &
898	top_heatflux /= 0.0 ) THEN
899	IF ( myid == 0 ) THEN
900	PRINT*, '+++ check_parameters:'
901	PRINT*, ' boundary_condition: bc_pt_t = ', bc_pt_t
902	PRINT*, ' is not allowed with constant_top_heatflux = .TRUE.'
903	ENDIF
904	CALL local_stop
905	ENDIF
906
907	!
908	!-- In case of moisture or passive scalar, set boundary conditions for total
909	!-- water content / scalar
910	IF ( moisture .OR. passive_scalar ) THEN
911	IF ( moisture ) THEN
912	sq = 'q'
913	ELSE
914	sq = 's'
915	ENDIF
916	IF ( bc_q_b == 'dirichlet' ) THEN
917	ibc_q_b = 0
918	ELSEIF ( bc_q_b == 'neumann' ) THEN
919	ibc_q_b = 1
920	ELSE
921	IF ( myid == 0 ) THEN
922	PRINT*, '+++ check_parameters:'
923	PRINT*, ' unknown boundary condition: bc_', sq, '_b = ', bc_q_b
924	ENDIF
925	CALL local_stop
926	ENDIF
927	IF ( bc_q_t == 'dirichlet' ) THEN
928	ibc_q_t = 0
929	ELSEIF ( bc_q_t == 'neumann' ) THEN
930	ibc_q_t = 1
931	ELSE
932	IF ( myid == 0 ) THEN
933	PRINT*, '+++ check_parameters:'
934	PRINT*, ' unknown boundary condition: bc_', sq, '_t = ', bc_q_t
935	ENDIF
936	CALL local_stop
937	ENDIF
938
939	IF ( surface_waterflux == 0.0 ) constant_waterflux = .FALSE.
940
941	!
942	!-- A given surface humidity implies Dirichlet boundary condition for
943	!-- moisture. In this case specification of a constant water flux is
944	!-- forbidden.
945	IF ( ibc_q_b == 0 .AND. constant_waterflux ) THEN
946	IF ( myid == 0 ) THEN
947	PRINT*, '+++ check_parameters:'
948	PRINT*, ' boundary_condition: bc_', sq, '_b = ', bc_q_b
949	PRINT*, ' is not allowed with prescribed surface flux'
950	ENDIF
951	CALL local_stop
952	ENDIF
953	IF ( constant_waterflux .AND. q_surface_initial_change /= 0.0 ) THEN
954	IF ( myid == 0 ) THEN
955	PRINT*, '+++ check_parameters: a prescribed surface flux is not'
956	PRINT*, ' allowed with ', sq, '_surface_initial_change (/=0)', &
957	' = ', q_surface_initial_change
958	ENDIF
959	CALL local_stop
960	ENDIF
961
962	ENDIF
963
964	!
965	!-- Boundary conditions for horizontal components of wind speed
966	IF ( bc_uv_b == 'dirichlet' ) THEN
967	ibc_uv_b = 0
968	ELSEIF ( bc_uv_b == 'neumann' ) THEN
969	ibc_uv_b = 1
970	IF ( prandtl_layer ) THEN
971	IF ( myid == 0 ) THEN
972	PRINT*, '+++ check_parameters:'
973	PRINT*, ' boundary condition: bc_uv_b = ', TRIM( bc_uv_b ), &
974	' is not allowed with'
975	PRINT*, ' prandtl_layer = .TRUE.'
976	ENDIF
977	CALL local_stop
978	ENDIF
979	ELSE
980	IF ( myid == 0 ) THEN
981	PRINT*, '+++ check_parameters:'
982	PRINT*, ' unknown boundary condition: bc_uv_b = ', bc_uv_b
983	ENDIF
984	CALL local_stop
985	ENDIF
986	IF ( bc_uv_t == 'dirichlet' ) THEN
987	ibc_uv_t = 0
988	ELSEIF ( bc_uv_t == 'neumann' ) THEN
989	ibc_uv_t = 1
990	ELSE
991	IF ( myid == 0 ) THEN
992	PRINT*, '+++ check_parameters:'
993	PRINT*, ' unknown boundary condition: bc_uv_t = ', bc_uv_t
994	ENDIF
995	CALL local_stop
996	ENDIF
997
998	!
999	!-- Compute and check, respectively, the Rayleigh Damping parameter
1000	IF ( rayleigh_damping_factor == -1.0 ) THEN
1001	IF ( momentum_advec == 'ups-scheme' ) THEN
1002	rayleigh_damping_factor = 0.01
1003	ELSE
1004	rayleigh_damping_factor = 0.0
1005	ENDIF
1006	ELSE
1007	IF ( rayleigh_damping_factor < 0.0 .OR. rayleigh_damping_factor > 1.0 ) &
1008	THEN
1009	IF ( myid == 0 ) THEN
1010	PRINT*, '+++ check_parameters:'
1011	PRINT*, ' rayleigh_damping_factor = ', rayleigh_damping_factor,&
1012	' out of range [0.0,1.0]'
1013	ENDIF
1014	CALL local_stop
1015	ENDIF
1016	ENDIF
1017
1018	IF ( rayleigh_damping_height == -1.0 ) THEN
1019	rayleigh_damping_height = 0.66666666666 * zu(nzt)
1020	ELSE
1021	IF ( rayleigh_damping_height < 0.0 .OR. &
1022	rayleigh_damping_height > zu(nzt) ) THEN
1023	IF ( myid == 0 ) THEN
1024	PRINT*, '+++ check_parameters:'
1025	PRINT*, ' rayleigh_damping_height = ', rayleigh_damping_height,&
1026	' out of range [0.0,', zu(nzt), ']'
1027	ENDIF
1028	CALL local_stop
1029	ENDIF
1030	ENDIF
1031
1032	!
1033	!-- Check limiters for Upstream-Spline scheme
1034	IF ( overshoot_limit_u < 0.0 .OR. overshoot_limit_v < 0.0 .OR. &
1035	overshoot_limit_w < 0.0 .OR. overshoot_limit_pt < 0.0 .OR. &
1036	overshoot_limit_e < 0.0 ) THEN
1037	IF ( myid == 0 ) THEN
1038	PRINT*, '+++ check_parameters:'
1039	PRINT*, ' overshoot_limit_... < 0.0 is not allowed'
1040	ENDIF
1041	CALL local_stop
1042	ENDIF
1043	IF ( ups_limit_u < 0.0 .OR. ups_limit_v < 0.0 .OR. ups_limit_w < 0.0 .OR. &
1044	ups_limit_pt < 0.0 .OR. ups_limit_e < 0.0 ) THEN
1045	IF ( myid == 0 ) THEN
1046	PRINT*, '+++ check_parameters:'
1047	PRINT*, ' ups_limit_... < 0.0 is not allowed'
1048	ENDIF
1049	CALL local_stop
1050	ENDIF
1051
1052	!
1053	!-- Check number of chosen statistic regions. More than 10 regions are not
1054	!-- allowed, because so far no more than 10 corresponding output files can
1055	!-- be opened (cf. check_open)
1056	IF ( statistic_regions > 9 .OR. statistic_regions < 0 ) THEN
1057	IF ( myid == 0 ) THEN
1058	PRINT*, '+++ check_parameters: Number of statistic_regions = ', &
1059	statistic_regions+1
1060	PRINT*, ' Only 10 regions are allowed'
1061	ENDIF
1062	CALL local_stop
1063	ENDIF
1064	IF ( normalizing_region > statistic_regions .OR. &
1065	normalizing_region < 0) THEN
1066	IF ( myid == 0 ) THEN
1067	PRINT*, '+++ check_parameters: normalizing_region = ', &
1068	normalizing_region, ' is unknown'
1069	PRINT*, ' Must be <= ', statistic_regions
1070	ENDIF
1071	CALL local_stop
1072	ENDIF
1073
1074	!
1075	!-- Set the default intervals for data output, if necessary
1076	!-- NOTE: dt_dosp has already been set in package_parin
1077	IF ( dt_data_output /= 9999999.9 ) THEN
1078	IF ( dt_dopr == 9999999.9 ) dt_dopr = dt_data_output
1079	IF ( dt_dopts == 9999999.9 ) dt_dopts = dt_data_output
1080	IF ( dt_do2d_xy == 9999999.9 ) dt_do2d_xy = dt_data_output
1081	IF ( dt_do2d_xz == 9999999.9 ) dt_do2d_xz = dt_data_output
1082	IF ( dt_do2d_yz == 9999999.9 ) dt_do2d_yz = dt_data_output
1083	IF ( dt_do3d == 9999999.9 ) dt_do3d = dt_data_output
1084	IF ( dt_data_output_av == 9999999.9 ) dt_data_output_av = dt_data_output
1085	ENDIF
1086
1087	!
1088	!-- Set the default skip time intervals for data output, if necessary
1089	IF ( skip_time_dopr == 9999999.9 ) &
1090	skip_time_dopr = skip_time_data_output
1091	IF ( skip_time_dosp == 9999999.9 ) &
1092	skip_time_dosp = skip_time_data_output
1093	IF ( skip_time_do2d_xy == 9999999.9 ) &
1094	skip_time_do2d_xy = skip_time_data_output
1095	IF ( skip_time_do2d_xz == 9999999.9 ) &
1096	skip_time_do2d_xz = skip_time_data_output
1097	IF ( skip_time_do2d_yz == 9999999.9 ) &
1098	skip_time_do2d_yz = skip_time_data_output
1099	IF ( skip_time_do3d == 9999999.9 ) &
1100	skip_time_do3d = skip_time_data_output
1101	IF ( skip_time_data_output_av == 9999999.9 ) &
1102	skip_time_data_output_av = skip_time_data_output
1103
1104	!
1105	!-- Check the average intervals (first for 3d-data, then for profiles and
1106	!-- spectra)
1107	IF ( averaging_interval > dt_data_output_av ) THEN
1108	IF ( myid == 0 ) THEN
1109	PRINT*, '+++ check_parameters: average_interval=', &
1110	averaging_interval, ' must be <= dt_data_output=', &
1111	dt_data_output
1112	ENDIF
1113	CALL local_stop
1114	ENDIF
1115
1116	IF ( averaging_interval_pr == 9999999.9 ) THEN
1117	averaging_interval_pr = averaging_interval
1118	ENDIF
1119
1120	IF ( averaging_interval_pr > dt_dopr ) THEN
1121	IF ( myid == 0 ) THEN
1122	PRINT*, '+++ check_parameters: averaging_interval_pr=', &
1123	averaging_interval_pr, ' must be <= dt_dopr=', dt_dopr
1124	ENDIF
1125	CALL local_stop
1126	ENDIF
1127
1128	IF ( averaging_interval_sp == 9999999.9 ) THEN
1129	averaging_interval_sp = averaging_interval
1130	ENDIF
1131
1132	IF ( averaging_interval_sp > dt_dosp ) THEN
1133	IF ( myid == 0 ) THEN
1134	PRINT*, '+++ check_parameters: averaging_interval_sp=', &
1135	averaging_interval_sp, ' must be <= dt_dosp=', dt_dosp
1136	ENDIF
1137	CALL local_stop
1138	ENDIF
1139
1140	!
1141	!-- Set the default interval for profiles entering the temporal average
1142	IF ( dt_averaging_input_pr == 9999999.9 ) THEN
1143	dt_averaging_input_pr = dt_averaging_input
1144	ENDIF
1145
1146	!
1147	!-- Set the default interval for the output of timeseries to a reasonable
1148	!-- value (tries to minimize the number of calls of flow_statistics)
1149	IF ( dt_dots == 9999999.9 ) THEN
1150	IF ( averaging_interval_pr == 0.0 ) THEN
1151	dt_dots = MIN( dt_run_control, dt_dopr )
1152	ELSE
1153	dt_dots = MIN( dt_run_control, dt_averaging_input_pr )
1154	ENDIF
1155	ENDIF
1156
1157	!
1158	!-- Check the sample rate for averaging (first for 3d-data, then for profiles)
1159	IF ( dt_averaging_input > averaging_interval ) THEN
1160	IF ( myid == 0 ) THEN
1161	PRINT*, '+++ check_parameters: dt_averaging_input=', &
1162	dt_averaging_input, ' must be <= averaging_interval=', &
1163	averaging_interval
1164	ENDIF
1165	CALL local_stop
1166	ENDIF
1167
1168	IF ( dt_averaging_input_pr > averaging_interval_pr ) THEN
1169	IF ( myid == 0 ) THEN
1170	PRINT*, '+++ check_parameters: dt_averaging_input_pr=', &
1171	dt_averaging_input_pr, &
1172	' must be <= averaging_interval_pr=', &
1173	averaging_interval_pr
1174	ENDIF
1175	CALL local_stop
1176	ENDIF
1177
1178	!
1179	!-- Set the default value for the integration interval of precipitation amount
1180	IF ( precipitation ) THEN
1181	IF ( precipitation_amount_interval == 9999999.9 ) THEN
1182	precipitation_amount_interval = dt_do2d_xy
1183	ELSE
1184	IF ( precipitation_amount_interval > dt_do2d_xy ) THEN
1185	IF ( myid == 0 ) PRINT*, '+++ check_parameters: ', &
1186	'precipitation_amount_interval =', &
1187	precipitation_amount_interval, &
1188	' must not be larger than dt_do2d_xy', &
1189	' = ', dt_do2d_xy
1190	CALL local_stop
1191	ENDIF
1192	ENDIF
1193	ENDIF
1194
1195	!
1196	!-- Determine the number of output profiles and check whether they are
1197	!-- permissible
1198	DO WHILE ( data_output_pr(dopr_n+1) /= ' ' )
1199
1200	dopr_n = dopr_n + 1
1201	i = dopr_n
1202
1203	!
1204	!-- Determine internal profile number (for hom, homs)
1205	!-- and store height levels
1206	SELECT CASE ( TRIM( data_output_pr(i) ) )
1207
1208	CASE ( 'u', '#u' )
1209	dopr_index(i) = 1
1210	hom(:,2,1,:) = SPREAD( zu, 2, statistic_regions+1 )
1211	IF ( data_output_pr(i)(1:1) == '#' ) THEN
1212	dopr_initial_index(i) = 5
1213	hom(:,2,5,:) = SPREAD( zu, 2, statistic_regions+1 )
1214	data_output_pr(i) = data_output_pr(i)(2:)
1215	ENDIF
1216
1217	CASE ( 'v', '#v' )
1218	dopr_index(i) = 2
1219	hom(:,2,2,:) = SPREAD( zu, 2, statistic_regions+1 )
1220	IF ( data_output_pr(i)(1:1) == '#' ) THEN
1221	dopr_initial_index(i) = 6
1222	hom(:,2,6,:) = SPREAD( zu, 2, statistic_regions+1 )
1223	data_output_pr(i) = data_output_pr(i)(2:)
1224	ENDIF
1225
1226	CASE ( 'w' )
1227	dopr_index(i) = 3
1228	hom(:,2,3,:) = SPREAD( zw, 2, statistic_regions+1 )
1229
1230	CASE ( 'pt', '#pt' )
1231	IF ( .NOT. cloud_physics ) THEN
1232	dopr_index(i) = 4
1233	hom(:,2,4,:) = SPREAD( zu, 2, statistic_regions+1 )
1234	IF ( data_output_pr(i)(1:1) == '#' ) THEN
1235	dopr_initial_index(i) = 7
1236	hom(:,2,7,:) = SPREAD( zu, 2, statistic_regions+1 )
1237	hom(nzb,2,7,:) = 0.0 ! weil zu(nzb) negativ ist
1238	data_output_pr(i) = data_output_pr(i)(2:)
1239	ENDIF
1240	ELSE
1241	dopr_index(i) = 43
1242	hom(:,2,43,:) = SPREAD( zu, 2, statistic_regions+1 )
1243	IF ( data_output_pr(i)(1:1) == '#' ) THEN
1244	dopr_initial_index(i) = 28
1245	hom(:,2,28,:) = SPREAD( zu, 2, statistic_regions+1 )
1246	hom(nzb,2,28,:) = 0.0 ! weil zu(nzb) negativ ist
1247	data_output_pr(i) = data_output_pr(i)(2:)
1248	ENDIF
1249	ENDIF
1250
1251	CASE ( 'e' )
1252	dopr_index(i) = 8
1253	hom(:,2,8,:) = SPREAD( zu, 2, statistic_regions+1 )
1254	hom(nzb,2,8,:) = 0.0
1255
1256	CASE ( 'km', '#km' )
1257	dopr_index(i) = 9
1258	hom(:,2,9,:) = SPREAD( zu, 2, statistic_regions+1 )
1259	hom(nzb,2,9,:) = 0.0
1260	IF ( data_output_pr(i)(1:1) == '#' ) THEN
1261	dopr_initial_index(i) = 23
1262	hom(:,2,23,:) = hom(:,2,9,:)
1263	data_output_pr(i) = data_output_pr(i)(2:)
1264	ENDIF
1265
1266	CASE ( 'kh', '#kh' )
1267	dopr_index(i) = 10
1268	hom(:,2,10,:) = SPREAD( zu, 2, statistic_regions+1 )
1269	hom(nzb,2,10,:) = 0.0
1270	IF ( data_output_pr(i)(1:1) == '#' ) THEN
1271	dopr_initial_index(i) = 24
1272	hom(:,2,24,:) = hom(:,2,10,:)
1273	data_output_pr(i) = data_output_pr(i)(2:)
1274	ENDIF
1275
1276	CASE ( 'l', '#l' )
1277	dopr_index(i) = 11
1278	hom(:,2,11,:) = SPREAD( zu, 2, statistic_regions+1 )
1279	hom(nzb,2,11,:) = 0.0
1280	IF ( data_output_pr(i)(1:1) == '#' ) THEN
1281	dopr_initial_index(i) = 25
1282	hom(:,2,25,:) = hom(:,2,11,:)
1283	data_output_pr(i) = data_output_pr(i)(2:)
1284	ENDIF
1285
1286	CASE ( 'w"u"' )
1287	dopr_index(i) = 12
1288	hom(:,2,12,:) = SPREAD( zw, 2, statistic_regions+1 )
1289	IF ( prandtl_layer ) hom(nzb,2,12,:) = zu(1)
1290
1291	CASE ( 'wu' )
1292	dopr_index(i) = 13
1293	hom(:,2,13,:) = SPREAD( zw, 2, statistic_regions+1 )
1294
1295	CASE ( 'w"v"' )
1296	dopr_index(i) = 14
1297	hom(:,2,14,:) = SPREAD( zw, 2, statistic_regions+1 )
1298	IF ( prandtl_layer ) hom(nzb,2,14,:) = zu(1)
1299
1300	CASE ( 'wv' )
1301	dopr_index(i) = 15
1302	hom(:,2,15,:) = SPREAD( zw, 2, statistic_regions+1 )
1303
1304	CASE ( 'w"pt"' )
1305	dopr_index(i) = 16
1306	hom(:,2,16,:) = SPREAD( zw, 2, statistic_regions+1 )
1307
1308	CASE ( 'wpt' )
1309	dopr_index(i) = 17
1310	hom(:,2,17,:) = SPREAD( zw, 2, statistic_regions+1 )
1311
1312	CASE ( 'wpt' )
1313	dopr_index(i) = 18
1314	hom(:,2,18,:) = SPREAD( zw, 2, statistic_regions+1 )
1315
1316	CASE ( 'wu' )
1317	dopr_index(i) = 19
1318	hom(:,2,19,:) = SPREAD( zw, 2, statistic_regions+1 )
1319	IF ( prandtl_layer ) hom(nzb,2,19,:) = zu(1)
1320
1321	CASE ( 'wv' )
1322	dopr_index(i) = 20
1323	hom(:,2,20,:) = SPREAD( zw, 2, statistic_regions+1 )
1324	IF ( prandtl_layer ) hom(nzb,2,20,:) = zu(1)
1325
1326	CASE ( 'wptBC' )
1327	dopr_index(i) = 21
1328	hom(:,2,21,:) = SPREAD( zw, 2, statistic_regions+1 )
1329
1330	CASE ( 'wptBC' )
1331	dopr_index(i) = 22
1332	hom(:,2,22,:) = SPREAD( zw, 2, statistic_regions+1 )
1333
1334	CASE ( 'u*2' )
1335	dopr_index(i) = 30
1336	hom(:,2,30,:) = SPREAD( zu, 2, statistic_regions+1 )
1337
1338	CASE ( 'v*2' )
1339	dopr_index(i) = 31
1340	hom(:,2,31,:) = SPREAD( zu, 2, statistic_regions+1 )
1341
1342	CASE ( 'w*2' )
1343	dopr_index(i) = 32
1344	hom(:,2,32,:) = SPREAD( zw, 2, statistic_regions+1 )
1345
1346	CASE ( 'pt*2' )
1347	dopr_index(i) = 33
1348	hom(:,2,33,:) = SPREAD( zu, 2, statistic_regions+1 )
1349
1350	CASE ( 'e*' )
1351	dopr_index(i) = 34
1352	hom(:,2,34,:) = SPREAD( zu, 2, statistic_regions+1 )
1353
1354	CASE ( 'w2pt' )
1355	dopr_index(i) = 35
1356	hom(:,2,35,:) = SPREAD( zw, 2, statistic_regions+1 )
1357
1358	CASE ( 'wpt2' )
1359	dopr_index(i) = 36
1360	hom(:,2,36,:) = SPREAD( zw, 2, statistic_regions+1 )
1361
1362	CASE ( 'we' )
1363	dopr_index(i) = 37
1364	hom(:,2,37,:) = SPREAD( zw, 2, statistic_regions+1 )
1365
1366	CASE ( 'w*3' )
1367	dopr_index(i) = 38
1368	hom(:,2,38,:) = SPREAD( zw, 2, statistic_regions+1 )
1369
1370	CASE ( 'Sw' )
1371	dopr_index(i) = 39
1372	hom(:,2,39,:) = SPREAD( zw, 2, statistic_regions+1 )
1373
1374	CASE ( 'q', '#q' )
1375	IF ( .NOT. cloud_physics ) THEN
1376	IF ( myid == 0 ) THEN
1377	PRINT*, '+++ check_parameters: data_output_pr = ', &
1378	data_output_pr(i), &
1379	' is not implemented for cloud_physics = FALSE'
1380	ENDIF
1381	CALL local_stop
1382	ELSE
1383	dopr_index(i) = 41
1384	hom(:,2,41,:) = SPREAD( zu, 2, statistic_regions+1 )
1385	IF ( data_output_pr(i)(1:1) == '#' ) THEN
1386	dopr_initial_index(i) = 26
1387	hom(:,2,26,:) = SPREAD( zu, 2, statistic_regions+1 )
1388	hom(nzb,2,26,:) = 0.0 ! weil zu(nzb) negativ ist
1389	data_output_pr(i) = data_output_pr(i)(2:)
1390	ENDIF
1391	ENDIF
1392
1393	CASE ( 's', '#s' )
1394	IF ( .NOT. passive_scalar ) THEN
1395	IF ( myid == 0 ) THEN
1396	PRINT*, '+++ check_parameters: data_output_pr = ', &
1397	data_output_pr(i), &
1398	' is not implemented for passive_scalar = FALSE'
1399	ENDIF
1400	CALL local_stop
1401	ELSE
1402	dopr_index(i) = 41
1403	hom(:,2,41,:) = SPREAD( zu, 2, statistic_regions+1 )
1404	IF ( data_output_pr(i)(1:1) == '#' ) THEN
1405	dopr_initial_index(i) = 26
1406	hom(:,2,26,:) = SPREAD( zu, 2, statistic_regions+1 )
1407	hom(nzb,2,26,:) = 0.0 ! weil zu(nzb) negativ ist
1408	data_output_pr(i) = data_output_pr(i)(2:)
1409	ENDIF
1410	ENDIF
1411
1412	CASE ( 'qv', '#qv' )
1413	IF ( .NOT. cloud_physics ) THEN
1414	dopr_index(i) = 41
1415	hom(:,2,41,:) = SPREAD( zu, 2, statistic_regions+1 )
1416	IF ( data_output_pr(i)(1:1) == '#' ) THEN
1417	dopr_initial_index(i) = 26
1418	hom(:,2,26,:) = SPREAD( zu, 2, statistic_regions+1 )
1419	hom(nzb,2,26,:) = 0.0 ! weil zu(nzb) negativ ist
1420	data_output_pr(i) = data_output_pr(i)(2:)
1421	ENDIF
1422	ELSE
1423	dopr_index(i) = 42
1424	hom(:,2,42,:) = SPREAD( zu, 2, statistic_regions+1 )
1425	IF ( data_output_pr(i)(1:1) == '#' ) THEN
1426	dopr_initial_index(i) = 27
1427	hom(:,2,27,:) = SPREAD( zu, 2, statistic_regions+1 )
1428	hom(nzb,2,27,:) = 0.0 ! weil zu(nzb) negativ ist
1429	data_output_pr(i) = data_output_pr(i)(2:)
1430	ENDIF
1431	ENDIF
1432
1433	CASE ( 'lpt', '#lpt' )
1434	IF ( .NOT. cloud_physics ) THEN
1435	IF ( myid == 0 ) THEN
1436	PRINT*, '+++ check_parameters: data_output_pr = ', &
1437	data_output_pr(i), &
1438	' is not implemented for cloud_physics = FALSE'
1439	ENDIF
1440	CALL local_stop
1441	ELSE
1442	dopr_index(i) = 4
1443	hom(:,2,4,:) = SPREAD( zu, 2, statistic_regions+1 )
1444	IF ( data_output_pr(i)(1:1) == '#' ) THEN
1445	dopr_initial_index(i) = 7
1446	hom(:,2,7,:) = SPREAD( zu, 2, statistic_regions+1 )
1447	hom(nzb,2,7,:) = 0.0 ! weil zu(nzb) negativ ist
1448	data_output_pr(i) = data_output_pr(i)(2:)
1449	ENDIF
1450	ENDIF
1451
1452	CASE ( 'vpt', '#vpt' )
1453	dopr_index(i) = 44
1454	hom(:,2,44,:) = SPREAD( zu, 2, statistic_regions+1 )
1455	IF ( data_output_pr(i)(1:1) == '#' ) THEN
1456	dopr_initial_index(i) = 29
1457	hom(:,2,29,:) = SPREAD( zu, 2, statistic_regions+1 )
1458	hom(nzb,2,29,:) = 0.0 ! weil zu(nzb) negativ ist
1459	data_output_pr(i) = data_output_pr(i)(2:)
1460	ENDIF
1461
1462	CASE ( 'w"vpt"' )
1463	dopr_index(i) = 45
1464	hom(:,2,45,:) = SPREAD( zw, 2, statistic_regions+1 )
1465
1466	CASE ( 'wvpt' )
1467	dopr_index(i) = 46
1468	hom(:,2,46,:) = SPREAD( zw, 2, statistic_regions+1 )
1469
1470	CASE ( 'wvpt' )
1471	dopr_index(i) = 47
1472	hom(:,2,47,:) = SPREAD( zw, 2, statistic_regions+1 )
1473
1474	CASE ( 'w"q"' )
1475	IF ( .NOT. cloud_physics ) THEN
1476	IF ( myid == 0 ) THEN
1477	PRINT*, '+++ check_parameters: data_output_pr = ', &
1478	data_output_pr(i), &
1479	' is not implemented for cloud_physics = FALSE'
1480	ENDIF
1481	CALL local_stop
1482	ELSE
1483	dopr_index(i) = 48
1484	hom(:,2,48,:) = SPREAD( zw, 2, statistic_regions+1 )
1485	ENDIF
1486
1487	CASE ( 'wq' )
1488	IF ( .NOT. cloud_physics ) THEN
1489	IF ( myid == 0 ) THEN
1490	PRINT*, '+++ check_parameters: data_output_pr = ', &
1491	data_output_pr(i), &
1492	' is not implemented for cloud_physics = FALSE'
1493	ENDIF
1494	CALL local_stop
1495	ELSE
1496	dopr_index(i) = 49
1497	hom(:,2,49,:) = SPREAD( zw, 2, statistic_regions+1 )
1498	ENDIF
1499
1500	CASE ( 'wq' )
1501	IF ( .NOT. cloud_physics ) THEN
1502	IF ( myid == 0 ) THEN
1503	PRINT*, '+++ check_parameters: data_output_pr = ', &
1504	data_output_pr(i), &
1505	' is not implemented for cloud_physics = FALSE'
1506	ENDIF
1507	CALL local_stop
1508	ELSE
1509	dopr_index(i) = 50
1510	hom(:,2,50,:) = SPREAD( zw, 2, statistic_regions+1 )
1511	ENDIF
1512
1513	CASE ( 'w"s"' )
1514	IF ( .NOT. passive_scalar ) THEN
1515	IF ( myid == 0 ) THEN
1516	PRINT*, '+++ check_parameters: data_output_pr = ', &
1517	data_output_pr(i), &
1518	' is not implemented for passive_scalar = FALSE'
1519	ENDIF
1520	CALL local_stop
1521	ELSE
1522	dopr_index(i) = 48
1523	hom(:,2,48,:) = SPREAD( zw, 2, statistic_regions+1 )
1524	ENDIF
1525
1526	CASE ( 'ws' )
1527	IF ( .NOT. passive_scalar ) THEN
1528	IF ( myid == 0 ) THEN
1529	PRINT*, '+++ check_parameters: data_output_pr = ', &
1530	data_output_pr(i), &
1531	' is not implemented for passive_scalar = FALSE'
1532	ENDIF
1533	CALL local_stop
1534	ELSE
1535	dopr_index(i) = 49
1536	hom(:,2,49,:) = SPREAD( zw, 2, statistic_regions+1 )
1537	ENDIF
1538
1539	CASE ( 'ws' )
1540	IF ( .NOT. passive_scalar ) THEN
1541	IF ( myid == 0 ) THEN
1542	PRINT*, '+++ check_parameters: data_output_pr = ', &
1543	data_output_pr(i), &
1544	' is not implemented for passive_scalar = FALSE'
1545	ENDIF
1546	CALL local_stop
1547	ELSE
1548	dopr_index(i) = 50
1549	hom(:,2,50,:) = SPREAD( zw, 2, statistic_regions+1 )
1550	ENDIF
1551
1552	CASE ( 'w"qv"' )
1553	IF ( moisture .AND. .NOT. cloud_physics ) &
1554	THEN
1555	dopr_index(i) = 48
1556	hom(:,2,48,:) = SPREAD( zw, 2, statistic_regions+1 )
1557	ELSEIF( moisture .AND. cloud_physics ) THEN
1558	dopr_index(i) = 51
1559	hom(:,2,51,:) = SPREAD( zw, 2, statistic_regions+1 )
1560	ELSE
1561	IF ( myid == 0 ) THEN
1562	PRINT*, '+++ check_parameters: data_output_pr = ', &
1563	data_output_pr(i), &
1564	' is not implemented for cloud_physics = FALSE', &
1565	' and moisture = FALSE'
1566	ENDIF
1567	CALL local_stop
1568	ENDIF
1569
1570	CASE ( 'wqv' )
1571	IF ( moisture .AND. .NOT. cloud_physics ) &
1572	THEN
1573	dopr_index(i) = 49
1574	hom(:,2,49,:) = SPREAD( zw, 2, statistic_regions+1 )
1575	ELSEIF( moisture .AND. cloud_physics ) THEN
1576	dopr_index(i) = 52
1577	hom(:,2,52,:) = SPREAD( zw, 2, statistic_regions+1 )
1578	ELSE
1579	IF ( myid == 0 ) THEN
1580	PRINT*, '+++ check_parameters: data_output_pr = ', &
1581	data_output_pr(i), &
1582	' is not implemented for cloud_physics = FALSE', &
1583	' and moisture = FALSE'
1584	ENDIF
1585	CALL local_stop
1586	ENDIF
1587
1588	CASE ( 'wqv' )
1589	IF ( moisture .AND. .NOT. cloud_physics ) &
1590	THEN
1591	dopr_index(i) = 50
1592	hom(:,2,50,:) = SPREAD( zw, 2, statistic_regions+1 )
1593	ELSEIF( moisture .AND. cloud_physics ) THEN
1594	dopr_index(i) = 53
1595	hom(:,2,53,:) = SPREAD( zw, 2, statistic_regions+1 )
1596	ELSE
1597	IF ( myid == 0 ) THEN
1598	PRINT*, '+++ check_parameters: data_output_pr = ', &
1599	data_output_pr(i), &
1600	' is not implemented for cloud_physics = FALSE', &
1601	' and moisture = FALSE'
1602	ENDIF
1603	CALL local_stop
1604	ENDIF
1605
1606	CASE ( 'ql' )
1607	IF ( .NOT. cloud_physics .AND. .NOT. cloud_droplets ) THEN
1608	IF ( myid == 0 ) THEN
1609	PRINT*, '+++ check_parameters: data_output_pr = ', &
1610	data_output_pr(i), &
1611	' is not implemented for cloud_physics = FALSE'
1612	ENDIF
1613	CALL local_stop
1614	ELSE
1615	dopr_index(i) = 54
1616	hom(:,2,54,:) = SPREAD( zu, 2, statistic_regions+1 )
1617	ENDIF
1618
1619	CASE ( 'wuu*/dz' )
1620	dopr_index(i) = 55
1621	hom(:,2,55,:) = SPREAD( zu, 2, statistic_regions+1 )
1622
1623	CASE ( 'wp/dz' )
1624	dopr_index(i) = 56
1625	hom(:,2,56,:) = SPREAD( zu, 2, statistic_regions+1 )
1626
1627	CASE ( 'w"e/dz' )
1628	dopr_index(i) = 57
1629	hom(:,2,57,:) = SPREAD( zu, 2, statistic_regions+1 )
1630
1631	CASE ( 'u"pt"' )
1632	dopr_index(i) = 58
1633	hom(:,2,58,:) = SPREAD( zu, 2, statistic_regions+1 )
1634
1635	CASE ( 'upt' )
1636	dopr_index(i) = 59
1637	hom(:,2,59,:) = SPREAD( zu, 2, statistic_regions+1 )
1638
1639	CASE ( 'upt_t' )
1640	dopr_index(i) = 60
1641	hom(:,2,60,:) = SPREAD( zu, 2, statistic_regions+1 )
1642
1643	CASE ( 'v"pt"' )
1644	dopr_index(i) = 61
1645	hom(:,2,61,:) = SPREAD( zu, 2, statistic_regions+1 )
1646
1647	CASE ( 'vpt' )
1648	dopr_index(i) = 62
1649	hom(:,2,62,:) = SPREAD( zu, 2, statistic_regions+1 )
1650
1651	CASE ( 'vpt_t' )
1652	dopr_index(i) = 63
1653	hom(:,2,63,:) = SPREAD( zu, 2, statistic_regions+1 )
1654
1655
1656	CASE DEFAULT
1657	IF ( myid == 0 ) THEN
1658	PRINT*, '+++ check_parameters: unknown output profile: ', &
1659	'data_output_pr = ', data_output_pr(i)
1660	ENDIF
1661	CALL local_stop
1662
1663	END SELECT
1664	!
1665	!-- Check to which of the predefined coordinate systems the profile belongs
1666	DO k = 1, crmax
1667	IF ( INDEX( cross_profiles(k), ' '//TRIM( data_output_pr(i) )//' ' ) &
1668	/=0 ) &
1669	THEN
1670	dopr_crossindex(i) = k
1671	EXIT
1672	ENDIF
1673	ENDDO
1674	!
1675	!-- Generate the text for the labels of the PROFIL output file. "-characters
1676	!-- must be substituted, otherwise PROFIL would interpret them as TeX
1677	!-- control characters
1678	dopr_label(i) = data_output_pr(i)
1679	position = INDEX( dopr_label(i) , '"' )
1680	DO WHILE ( position /= 0 )
1681	dopr_label(i)(position:position) = ''''
1682	position = INDEX( dopr_label(i) , '"' )
1683	ENDDO
1684
1685	ENDDO
1686
1687	!
1688	!-- y-value range of the coordinate system (PROFIL).
1689	!-- x-value range determined in plot_1d.
1690	cross_uymin = 0.0
1691	IF ( z_max_do1d == -1.0 ) THEN
1692	cross_uymax = zu(nzt+1)
1693	ELSEIF ( z_max_do1d < zu(nzb+1) .OR. z_max_do1d > zu(nzt+1) ) THEN
1694	IF ( myid == 0 ) PRINT*, '+++ check_parameters: z_max_do1d=', &
1695	z_max_do1d,' must be >= ', zu(nzb+1), &
1696	' or <= ', zu(nzt+1)
1697	CALL local_stop
1698	ELSE
1699	cross_uymax = z_max_do1d
1700	ENDIF
1701
1702	!
1703	!-- Check whether the chosen normalizing factor for the coordinate systems is
1704	!-- permissible
1705	DO i = 1, crmax
1706	SELECT CASE ( TRIM( cross_normalized_x(i) ) ) ! TRIM required on IBM
1707
1708	CASE ( '', 'wpt0', 'ws2', 'tsw2', 'ws3', 'ws2tsw', 'wstsw2' )
1709	j = 0
1710
1711	CASE DEFAULT
1712	IF ( myid == 0 ) THEN
1713	PRINT*, '+++ check_parameters: unknown normalize method'
1714	PRINT*, ' cross_normalized_x="',cross_normalized_x(i),'"'
1715	ENDIF
1716	CALL local_stop
1717
1718	END SELECT
1719	SELECT CASE ( TRIM( cross_normalized_y(i) ) ) ! TRIM required on IBM
1720
1721	CASE ( '', 'z_i' )
1722	j = 0
1723
1724	CASE DEFAULT
1725	IF ( myid == 0 ) THEN
1726	PRINT*, '+++ check_parameters: unknown normalize method'
1727	PRINT*, ' cross_normalized_y="',cross_normalized_y(i),'"'
1728	ENDIF
1729	CALL local_stop
1730
1731	END SELECT
1732	ENDDO
1733	!
1734	!-- Check normalized y-value range of the coordinate system (PROFIL)
1735	IF ( z_max_do1d_normalized /= -1.0 .AND. z_max_do1d_normalized <= 0.0 ) &
1736	THEN
1737	IF ( myid == 0 ) PRINT*,'+++ check_parameters: z_max_do1d_normalize', &
1738	'd=', z_max_do1d_normalized, ' must be >= 0.0'
1739	CALL local_stop
1740	ENDIF
1741
1742
1743	!
1744	!-- Append user-defined data output variables to the standard data output
1745	IF ( data_output_user(1) /= ' ' ) THEN
1746	i = 1
1747	DO WHILE ( data_output(i) /= ' ' .AND. i <= 100 )
1748	i = i + 1
1749	ENDDO
1750	j = 1
1751	DO WHILE ( data_output_user(j) /= ' ' .AND. j <= 100 )
1752	IF ( i > 100 ) THEN
1753	IF ( myid == 0 ) THEN
1754	PRINT*, '+++ check_parameters: number of output quantitities', &
1755	' given by data_output and data_output_user'
1756	PRINT*, ' exceeds the limit of 100'
1757	ENDIF
1758	CALL local_stop
1759	ENDIF
1760	data_output(i) = data_output_user(j)
1761	i = i + 1
1762	j = j + 1
1763	ENDDO
1764	ENDIF
1765
1766	!
1767	!-- Check and set steering parameters for 2d/3d data output and averaging
1768	i = 1
1769	DO WHILE ( data_output(i) /= ' ' .AND. i <= 100 )
1770	!
1771	!-- Check for data averaging
1772	ilen = LEN_TRIM( data_output(i) )
1773	j = 0 ! no data averaging
1774	IF ( ilen > 3 ) THEN
1775	IF ( data_output(i)(ilen-2:ilen) == '_av' ) THEN
1776	j = 1 ! data averaging
1777	data_output(i) = data_output(i)(1:ilen-3)
1778	ENDIF
1779	ENDIF
1780	!
1781	!-- Check for cross section or volume data
1782	ilen = LEN_TRIM( data_output(i) )
1783	k = 0 ! 3d data
1784	var = data_output(i)(1:ilen)
1785	IF ( ilen > 3 ) THEN
1786	IF ( data_output(i)(ilen-2:ilen) == '_xy' .OR. &
1787	data_output(i)(ilen-2:ilen) == '_xz' .OR. &
1788	data_output(i)(ilen-2:ilen) == '_yz' ) THEN
1789	k = 1 ! 2d data
1790	var = data_output(i)(1:ilen-3)
1791	ENDIF
1792	ENDIF
1793	!
1794	!-- Check for allowed value and set units
1795	SELECT CASE ( TRIM( var ) )
1796
1797	CASE ( 'e' )
1798	IF ( constant_diffusion ) THEN
1799	IF ( myid == 0 ) THEN
1800	PRINT*, '+++ check_parameters: output of "', TRIM( var ), &
1801	'" requires constant_diffusion = .FALSE.'
1802	ENDIF
1803	CALL local_stop
1804	ENDIF
1805	unit = 'm2/s2'
1806
1807	CASE ( 'pc', 'pr' )
1808	IF ( .NOT. particle_advection ) THEN
1809	IF ( myid == 0 ) THEN
1810	PRINT*, '+++ check_parameters: output of "', TRIM( var ), &
1811	'" requires particle package'
1812	PRINT*, ' (mrun-option "-p particles")'
1813	ENDIF
1814	CALL local_stop
1815	ENDIF
1816	IF ( TRIM( var ) == 'pc' ) unit = 'number'
1817	IF ( TRIM( var ) == 'pr' ) unit = 'm'
1818
1819	CASE ( 'q', 'vpt' )
1820	IF ( .NOT. moisture ) THEN
1821	IF ( myid == 0 ) THEN
1822	PRINT*, '+++ check_parameters: output of "', TRIM( var ), &
1823	'" requires moisture = .TRUE.'
1824	ENDIF
1825	CALL local_stop
1826	ENDIF
1827	IF ( TRIM( var ) == 'q' ) unit = 'kg/kg'
1828	IF ( TRIM( var ) == 'vpt' ) unit = 'K'
1829
1830	CASE ( 'ql' )
1831	IF ( .NOT. ( cloud_physics .OR. cloud_droplets ) ) THEN
1832	IF ( myid == 0 ) THEN
1833	PRINT*, '+++ check_parameters: output of "', TRIM( var ), &
1834	'" requires cloud_physics = .TRUE.'
1835	PRINT*, ' or cloud_droplets = .TRUE.'
1836	ENDIF
1837	CALL local_stop
1838	ENDIF
1839	unit = 'kg/kg'
1840
1841	CASE ( 'ql_c', 'ql_v', 'ql_vp' )
1842	IF ( .NOT. cloud_droplets ) THEN
1843	IF ( myid == 0 ) THEN
1844	PRINT*, '+++ check_parameters: output of "', TRIM( var ), &
1845	'" requires cloud_droplets = .TRUE.'
1846	ENDIF
1847	CALL local_stop
1848	ENDIF
1849	IF ( TRIM( var ) == 'ql_c' ) unit = 'kg/kg'
1850	IF ( TRIM( var ) == 'ql_v' ) unit = 'm3'
1851	IF ( TRIM( var ) == 'ql_vp' ) unit = 'none'
1852
1853	CASE ( 'qv' )
1854	IF ( .NOT. cloud_physics ) THEN
1855	IF ( myid == 0 ) THEN
1856	PRINT*, '+++ check_parameters: output of "', TRIM( var ), &
1857	'" requires cloud_physics = .TRUE.'
1858	ENDIF
1859	CALL local_stop
1860	ENDIF
1861	unit = 'kg/kg'
1862
1863	CASE ( 's' )
1864	IF ( .NOT. passive_scalar ) THEN
1865	IF ( myid == 0 ) THEN
1866	PRINT*, '+++ check_parameters: output of "', TRIM( var ), &
1867	'" requires passive_scalar = .TRUE.'
1868	ENDIF
1869	CALL local_stop
1870	ENDIF
1871	unit = 'conc'
1872
1873	CASE ( 'u', 't', 'lwp', 'pra', 'prr', 'z0' )
1874	IF ( k == 0 .OR. data_output(i)(ilen-2:ilen) /= '_xy' ) THEN
1875	IF ( myid == 0 ) THEN
1876	PRINT*, '+++ check_parameters: illegal value for data_',&
1877	'output: "', TRIM( var ), '" is only allowed'
1878	PRINT*, ' for horizontal cross section'
1879	ENDIF
1880	CALL local_stop
1881	ENDIF
1882	IF ( TRIM( var ) == 'lwp*' .AND. .NOT. cloud_physics ) THEN
1883	IF ( myid == 0 ) THEN
1884	PRINT*, '+++ check_parameters: output of "', TRIM( var ), &
1885	'" requires cloud_physics = .TRUE.'
1886	ENDIF
1887	CALL local_stop
1888	ENDIF
1889	IF ( TRIM( var ) == 'pra*' .AND. .NOT. precipitation ) THEN
1890	IF ( myid == 0 ) THEN
1891	PRINT*, '+++ check_parameters: output of "', TRIM( var ), &
1892	'" requires precipitation = .TRUE.'
1893	ENDIF
1894	CALL local_stop
1895	ENDIF
1896	IF ( TRIM( var ) == 'pra*' .AND. j == 1 ) THEN
1897	IF ( myid == 0 ) THEN
1898	PRINT*, '+++ check_parameters: temporal averaging of ', &
1899	' precipitation amount "', TRIM( var ), &
1900	'" not possible'
1901	ENDIF
1902	CALL local_stop
1903	ENDIF
1904	IF ( TRIM( var ) == 'prr*' .AND. .NOT. precipitation ) THEN
1905	IF ( myid == 0 ) THEN
1906	PRINT*, '+++ check_parameters: output of "', TRIM( var ), &
1907	'" requires precipitation = .TRUE.'
1908	ENDIF
1909	CALL local_stop
1910	ENDIF
1911
1912
1913	IF ( TRIM( var ) == 'u*' ) unit = 'm/s'
1914	IF ( TRIM( var ) == 't*' ) unit = 'K'
1915	IF ( TRIM( var ) == 'lwp' ) unit = 'kg/kgm'
1916	IF ( TRIM( var ) == 'pra*' ) unit = 'mm'
1917	IF ( TRIM( var ) == 'prr*' ) unit = 'mm/s'
1918	IF ( TRIM( var ) == 'z0*' ) unit = 'm'
1919
1920	CASE ( 'p', 'pt', 'u', 'v', 'w' )
1921	IF ( TRIM( var ) == 'p' ) unit = 'Pa'
1922	IF ( TRIM( var ) == 'pt' ) unit = 'K'
1923	IF ( TRIM( var ) == 'u' ) unit = 'm/s'
1924	IF ( TRIM( var ) == 'v' ) unit = 'm/s'
1925	IF ( TRIM( var ) == 'w' ) unit = 'm/s'
1926	CONTINUE
1927
1928	CASE DEFAULT
1929	CALL user_check_data_output( var, unit )
1930
1931	IF ( unit == 'illegal' ) THEN
1932	IF ( myid == 0 ) THEN
1933	IF ( data_output_user(1) /= ' ' ) THEN
1934	PRINT*, '+++ check_parameters: illegal value for data_',&
1935	'output or data_output_user: "', &
1936	TRIM( data_output(i) ), '"'
1937	ELSE
1938	PRINT*, '+++ check_parameters: illegal value for data_',&
1939	'output: "', TRIM( data_output(i) ), '"'
1940	ENDIF
1941	ENDIF
1942	CALL local_stop
1943	ENDIF
1944
1945	END SELECT
1946	!
1947	!-- Set the internal steering parameters appropriately
1948	IF ( k == 0 ) THEN
1949	do3d_no(j) = do3d_no(j) + 1
1950	do3d(j,do3d_no(j)) = data_output(i)
1951	do3d_unit(j,do3d_no(j)) = unit
1952	ELSE
1953	do2d_no(j) = do2d_no(j) + 1
1954	do2d(j,do2d_no(j)) = data_output(i)
1955	do2d_unit(j,do2d_no(j)) = unit
1956	IF ( data_output(i)(ilen-2:ilen) == '_xy' ) THEN
1957	data_output_xy(j) = .TRUE.
1958	ENDIF
1959	IF ( data_output(i)(ilen-2:ilen) == '_xz' ) THEN
1960	data_output_xz(j) = .TRUE.
1961	ENDIF
1962	IF ( data_output(i)(ilen-2:ilen) == '_yz' ) THEN
1963	data_output_yz(j) = .TRUE.
1964	ENDIF
1965	ENDIF
1966
1967	IF ( j == 1 ) THEN
1968	!
1969	!-- Check, if variable is already subject to averaging
1970	found = .FALSE.
1971	DO k = 1, doav_n
1972	IF ( TRIM( doav(k) ) == TRIM( var ) ) found = .TRUE.
1973	ENDDO
1974
1975	IF ( .NOT. found ) THEN
1976	doav_n = doav_n + 1
1977	doav(doav_n) = var
1978	ENDIF
1979	ENDIF
1980
1981	i = i + 1
1982	ENDDO
1983
1984	!
1985	!-- Store sectional planes in one shared array
1986	section(:,1) = section_xy
1987	section(:,2) = section_xz
1988	section(:,3) = section_yz
1989
1990	!
1991	!-- Upper plot limit (grid point value) for 1D profiles
1992	IF ( z_max_do1d == -1.0 ) THEN
1993	nz_do1d = nzt+1
1994	ELSE
1995	DO k = nzb+1, nzt+1
1996	nz_do1d = k
1997	IF ( zw(k) > z_max_do1d ) EXIT
1998	ENDDO
1999	ENDIF
2000
2001	!
2002	!-- Upper plot limit for 2D vertical sections
2003	IF ( z_max_do2d == -1.0 ) z_max_do2d = zu(nzt)
2004	IF ( z_max_do2d < zu(nzb+1) .OR. z_max_do2d > zu(nzt) ) THEN
2005	IF ( myid == 0 ) PRINT*, '+++ check_parameters: z_max_do2d=', &
2006	z_max_do2d, ' must be >= ', zu(nzb+1), &
2007	'(zu(nzb+1)) and <= ', zu(nzt), ' (zu(nzt))'
2008	CALL local_stop
2009	ENDIF
2010
2011	!
2012	!-- Upper plot limit for 3D arrays
2013	IF ( nz_do3d == -9999 ) nz_do3d = nzt + 1
2014
2015	!
2016	!-- Determine and check accuracy for compressed 3D plot output
2017	IF ( do3d_compress ) THEN
2018	!
2019	!-- Compression only permissible on T3E machines
2020	IF ( host(1:3) /= 't3e' ) THEN
2021	IF ( myid == 0 ) THEN
2022	PRINT*, '+++ check_parameters: do3d_compress = .TRUE. not allow', &
2023	'ed on host "', TRIM( host ), '"'
2024	ENDIF
2025	CALL local_stop
2026	ENDIF
2027
2028	i = 1
2029	DO WHILE ( do3d_comp_prec(i) /= ' ' )
2030
2031	ilen = LEN_TRIM( do3d_comp_prec(i) )
2032	IF ( LLT( do3d_comp_prec(i)(ilen:ilen), '0' ) .OR. &
2033	LGT( do3d_comp_prec(i)(ilen:ilen), '9' ) ) THEN
2034	IF ( myid == 0 ) THEN
2035	PRINT*, '+++ check_parameters: illegal precision: ', &
2036	'do3d_comp_prec(', i, ')="', TRIM(do3d_comp_prec(i)),'"'
2037	ENDIF
2038	CALL local_stop
2039	ENDIF
2040
2041	prec = IACHAR( do3d_comp_prec(i)(ilen:ilen) ) - IACHAR( '0' )
2042	var = do3d_comp_prec(i)(1:ilen-1)
2043
2044	SELECT CASE ( var )
2045
2046	CASE ( 'u' )
2047	j = 1
2048	CASE ( 'v' )
2049	j = 2
2050	CASE ( 'w' )
2051	j = 3
2052	CASE ( 'p' )
2053	j = 4
2054	CASE ( 'pt' )
2055	j = 5
2056
2057	CASE DEFAULT
2058	IF ( myid == 0 ) THEN
2059	PRINT*, '+++ check_parameters: unknown variable in ', &
2060	'assignment'
2061	PRINT*, ' do3d_comp_prec(', i, ')="', &
2062	TRIM( do3d_comp_prec(i) ),'"'
2063	ENDIF
2064	CALL local_stop
2065
2066	END SELECT
2067
2068	plot_3d_precision(j)%precision = prec
2069	i = i + 1
2070
2071	ENDDO
2072	ENDIF
2073
2074	!
2075	!-- Check the data output format(s)
2076	IF ( data_output_format(1) == ' ' ) THEN
2077	!
2078	!-- Default value
2079	netcdf_output = .TRUE.
2080	ELSE
2081	i = 1
2082	DO WHILE ( data_output_format(i) /= ' ' )
2083
2084	SELECT CASE ( data_output_format(i) )
2085
2086	CASE ( 'netcdf' )
2087	netcdf_output = .TRUE.
2088	CASE ( 'iso2d' )
2089	iso2d_output = .TRUE.
2090	CASE ( 'profil' )
2091	profil_output = .TRUE.
2092	CASE ( 'avs' )
2093	avs_output = .TRUE.
2094
2095	CASE DEFAULT
2096	IF ( myid == 0 ) THEN
2097	PRINT*, '+++ check_parameters:'
2098	PRINT*, ' unknown value for data_output_format "', &
2099	TRIM( data_output_format(i) ),'"'
2100	ENDIF
2101	CALL local_stop
2102
2103	END SELECT
2104
2105	i = i + 1
2106	IF ( i > 10 ) EXIT
2107
2108	ENDDO
2109
2110	ENDIF
2111
2112	!
2113	!-- Check netcdf precison
2114	ldum = .FALSE.
2115	CALL define_netcdf_header( 'ch', ldum, 0 )
2116
2117	!
2118	!-- Check, whether a constant diffusion coefficient shall be used
2119	IF ( km_constant /= -1.0 ) THEN
2120	IF ( km_constant < 0.0 ) THEN
2121	IF ( myid == 0 ) PRINT*, '+++ check_parameters: km_constant=', &
2122	km_constant, ' < 0.0'
2123	CALL local_stop
2124	ELSE
2125	IF ( prandtl_number < 0.0 ) THEN
2126	IF ( myid == 0 ) PRINT*,'+++ check_parameters: prandtl_number=',&
2127	prandtl_number, ' < 0.0'
2128	CALL local_stop
2129	ENDIF
2130	constant_diffusion = .TRUE.
2131
2132	IF ( prandtl_layer ) THEN
2133	IF ( myid == 0 ) PRINT*, '+++ check_parameters: prandtl_layer ',&
2134	'is not allowed with fixed value of km'
2135	CALL local_stop
2136	ENDIF
2137	ENDIF
2138	ENDIF
2139
2140	!
2141	!-- In case of non-cyclic lateral boundaries, set the default maximum value
2142	!-- for the horizontal diffusivity used within the outflow damping layer,
2143	!-- and check/set the width of the damping layer
2144	IF ( bc_lr /= 'cyclic' ) THEN
2145	IF ( km_damp_max == -1.0 ) THEN
2146	km_damp_max = 0.5 * dx
2147	ENDIF
2148	IF ( outflow_damping_width == -1.0 ) THEN
2149	outflow_damping_width = MIN( 20, nx/2 )
2150	ENDIF
2151	IF ( outflow_damping_width <= 0 .OR. outflow_damping_width > nx ) THEN
2152	IF ( myid == 0 ) PRINT*, '+++ check_parameters: outflow_damping w',&
2153	'idth out of range'
2154	CALL local_stop
2155	ENDIF
2156	ENDIF
2157
2158	IF ( bc_ns /= 'cyclic' ) THEN
2159	IF ( km_damp_max == -1.0 ) THEN
2160	km_damp_max = 0.5 * dy
2161	ENDIF
2162	IF ( outflow_damping_width == -1.0 ) THEN
2163	outflow_damping_width = MIN( 20, ny/2 )
2164	ENDIF
2165	IF ( outflow_damping_width <= 0 .OR. outflow_damping_width > ny ) THEN
2166	IF ( myid == 0 ) PRINT*, '+++ check_parameters: outflow_damping w',&
2167	'idth out of range'
2168	CALL local_stop
2169	ENDIF
2170	ENDIF
2171
2172	!
2173	!-- Check value range for rif
2174	IF ( rif_min >= rif_max ) THEN
2175	IF ( myid == 0 ) PRINT*, '+++ check_parameters: rif_min=', rif_min, &
2176	' must be less than rif_max=', rif_max
2177	CALL local_stop
2178	ENDIF
2179
2180	!
2181	!-- Determine upper and lower hight level indices for random perturbations
2182	IF ( disturbance_level_b == -1.0 ) THEN
2183	disturbance_level_b = zu(nzb+3)
2184	disturbance_level_ind_b = nzb + 3
2185	ELSEIF ( disturbance_level_b < zu(3) ) THEN
2186	IF ( myid == 0 ) PRINT*, '+++ check_parameters: disturbance_level_b=',&
2187	disturbance_level_b, ' must be >= ',zu(3), &
2188	'(zu(3))'
2189	CALL local_stop
2190	ELSEIF ( disturbance_level_b > zu(nzt-2) ) THEN
2191	IF ( myid == 0 ) PRINT*, '+++ check_parameters: disturbance_level_b=',&
2192	disturbance_level_b, ' must be <= ',zu(nzt-2),&
2193	'(zu(nzt-2))'
2194	CALL local_stop
2195	ELSE
2196	DO k = 3, nzt-2
2197	IF ( disturbance_level_b <= zu(k) ) THEN
2198	disturbance_level_ind_b = k
2199	EXIT
2200	ENDIF
2201	ENDDO
2202	ENDIF
2203
2204	IF ( disturbance_level_t == -1.0 ) THEN
2205	disturbance_level_t = zu(nzt/3)
2206	disturbance_level_ind_t = nzt / 3
2207	ELSEIF ( disturbance_level_t > zu(nzt-2) ) THEN
2208	IF ( myid == 0 ) PRINT*, '+++ check_parameters: disturbance_level_t=',&
2209	disturbance_level_t, ' must be <= ',zu(nzt-2),&
2210	'(zu(nzt-2))'
2211	CALL local_stop
2212	ELSEIF ( disturbance_level_t < disturbance_level_b ) THEN
2213	IF ( myid == 0 ) PRINT*, '+++ check_parameters: disturbance_level_t=',&
2214	disturbance_level_t, ' must be >= ', &
2215	'disturbance_level_b=', disturbance_level_b
2216	CALL local_stop
2217	ELSE
2218	DO k = 3, nzt-2
2219	IF ( disturbance_level_t <= zu(k) ) THEN
2220	disturbance_level_ind_t = k
2221	EXIT
2222	ENDIF
2223	ENDDO
2224	ENDIF
2225
2226	!
2227	!-- Check again whether the levels determined this way are ok.
2228	!-- Error may occur at automatic determination and too few grid points in
2229	!-- z-direction.
2230	IF ( disturbance_level_ind_t < disturbance_level_ind_b ) THEN
2231	IF ( myid == 0 ) PRINT*, '+++ check_parameters: ', &
2232	'disturbance_level_ind_t=', &
2233	disturbance_level_ind_t, ' must be >= ', &
2234	'disturbance_level_ind_b=', &
2235	disturbance_level_ind_b
2236	CALL local_stop
2237	ENDIF
2238
2239	!
2240	!-- Determine the horizontal index range for random perturbations.
2241	!-- In case of non-cyclic horizontal boundaries, no perturbations are imposed
2242	!-- near the inflow and the perturbation area is further limited to ...(1)
2243	!-- after the initial phase of the flow.
2244	dist_nxl = 0; dist_nxr = nx
2245	dist_nys = 0; dist_nyn = ny
2246	IF ( bc_lr /= 'cyclic' ) THEN
2247	IF ( inflow_disturbance_begin == -1 ) THEN
2248	inflow_disturbance_begin = MIN( 10, nx/2 )
2249	ENDIF
2250	IF ( inflow_disturbance_begin < 0 .OR. inflow_disturbance_begin > nx )&
2251	THEN
2252	IF ( myid == 0 ) PRINT*, '+++ check_parameters: inflow_disturbance',&
2253	'_begin out of range'
2254	CALL local_stop
2255	ENDIF
2256	IF ( inflow_disturbance_end == -1 ) THEN
2257	inflow_disturbance_end = MIN( 100, 3*nx/4 )
2258	ENDIF
2259	IF ( inflow_disturbance_end < 0 .OR. inflow_disturbance_end > nx ) &
2260	THEN
2261	IF ( myid == 0 ) PRINT*, '+++ check_parameters: inflow_disturbance',&
2262	'_end out of range'
2263	CALL local_stop
2264	ENDIF
2265	ELSEIF ( bc_ns /= 'cyclic' ) THEN
2266	IF ( inflow_disturbance_begin == -1 ) THEN
2267	inflow_disturbance_begin = MIN( 10, ny/2 )
2268	ENDIF
2269	IF ( inflow_disturbance_begin < 0 .OR. inflow_disturbance_begin > ny )&
2270	THEN
2271	IF ( myid == 0 ) PRINT*, '+++ check_parameters: inflow_disturbance',&
2272	'_begin out of range'
2273	CALL local_stop
2274	ENDIF
2275	IF ( inflow_disturbance_end == -1 ) THEN
2276	inflow_disturbance_end = MIN( 100, 3*ny/4 )
2277	ENDIF
2278	IF ( inflow_disturbance_end < 0 .OR. inflow_disturbance_end > ny ) &
2279	THEN
2280	IF ( myid == 0 ) PRINT*, '+++ check_parameters: inflow_disturbance',&
2281	'_end out of range'
2282	CALL local_stop
2283	ENDIF
2284	ENDIF
2285
2286	IF ( outflow_l ) THEN
2287	dist_nxr = nx - inflow_disturbance_begin
2288	dist_nxl(1) = nx - inflow_disturbance_end
2289	ELSEIF ( outflow_r ) THEN
2290	dist_nxl = inflow_disturbance_begin
2291	dist_nxr(1) = inflow_disturbance_end
2292	ELSEIF ( outflow_s ) THEN
2293	dist_nyn = ny - inflow_disturbance_begin
2294	dist_nys(1) = ny - inflow_disturbance_end
2295	ELSEIF ( outflow_n ) THEN
2296	dist_nys = inflow_disturbance_begin
2297	dist_nyn(1) = inflow_disturbance_end
2298	ENDIF
2299
2300	!
2301	!-- Check random generator
2302	IF ( random_generator /= 'system-specific' .AND. &
2303	random_generator /= 'numerical-recipes' ) THEN
2304	IF ( myid == 0 ) THEN
2305	PRINT*, '+++ check_parameters:'
2306	PRINT*, ' unknown random generator: random_generator=', &
2307	random_generator
2308	ENDIF
2309	CALL local_stop
2310	ENDIF
2311
2312	!
2313	!-- Determine damping level index for 1D model
2314	IF ( INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 ) THEN
2315	IF ( damp_level_1d == -1.0 ) THEN
2316	damp_level_1d = zu(nzt+1)
2317	damp_level_ind_1d = nzt + 1
2318	ELSEIF ( damp_level_1d < 0.0 .OR. damp_level_1d > zu(nzt+1) ) THEN
2319	IF ( myid == 0 ) PRINT*, '+++ check_parameters: damp_level_1d=', &
2320	damp_level_1d, ' must be > 0.0 and < ', &
2321	'zu(nzt+1)', zu(nzt+1)
2322	CALL local_stop
2323	ELSE
2324	DO k = 1, nzt+1
2325	IF ( damp_level_1d <= zu(k) ) THEN
2326	damp_level_ind_1d = k
2327	EXIT
2328	ENDIF
2329	ENDDO
2330	ENDIF
2331	ENDIF
2332	!
2333	!-- Check some other 1d-model parameters
2334	IF ( TRIM( mixing_length_1d ) /= 'as_in_3d_model' .AND. &
2335	TRIM( mixing_length_1d ) /= 'blackadar' ) THEN
2336	IF ( myid == 0 ) PRINT*, '+++ check_parameters: mixing_length_1d = "', &
2337	TRIM( mixing_length_1d ), '" is unknown'
2338	CALL local_stop
2339	ENDIF
2340	IF ( TRIM( dissipation_1d ) /= 'as_in_3d_model' .AND. &
2341	TRIM( dissipation_1d ) /= 'detering' ) THEN
2342	IF ( myid == 0 ) PRINT*, '+++ check_parameters: dissipation_1d = "', &
2343	TRIM( dissipation_1d ), '" is unknown'
2344	CALL local_stop
2345	ENDIF
2346
2347	!
2348	!-- Set time for the next user defined restart (time_restart is the
2349	!-- internal parameter for steering restart events)
2350	IF ( restart_time /= 9999999.9 ) THEN
2351	IF ( restart_time > simulated_time ) time_restart = restart_time
2352	ELSE
2353	!
2354	!-- In case of a restart run, set internal parameter to default (no restart)
2355	!-- if the NAMELIST-parameter restart_time is omitted
2356	time_restart = 9999999.9
2357	ENDIF
2358
2359	!
2360	!-- Set default value of the time needed to terminate a model run
2361	IF ( termination_time_needed == -1.0 ) THEN
2362	IF ( host(1:3) == 'ibm' ) THEN
2363	termination_time_needed = 300.0
2364	ELSE
2365	termination_time_needed = 35.0
2366	ENDIF
2367	ENDIF
2368
2369	!
2370	!-- Check the time needed to terminate a model run
2371	IF ( host(1:3) == 't3e' ) THEN
2372	!
2373	!-- Time needed must be at least 30 seconds on all CRAY machines, because
2374	!-- MPP_TREMAIN gives the remaining CPU time only in steps of 30 seconds
2375	IF ( termination_time_needed <= 30.0 ) THEN
2376	IF ( myid == 0 ) THEN
2377	PRINT*, '+++ check_parameters: termination_time_needed', &
2378	termination_time_needed
2379	PRINT*, ' must be > 30.0 on host "', host, &
2380	'"'
2381	ENDIF
2382	CALL local_stop
2383	ENDIF
2384	ELSEIF ( host(1:3) == 'ibm' ) THEN
2385	!
2386	!-- On IBM-regatta machines the time should be at least 300 seconds,
2387	!-- because the job time consumed before executing palm (for compiling,
2388	!-- copying of files, etc.) has to be regarded
2389	IF ( termination_time_needed < 300.0 ) THEN
2390	IF ( myid == 0 ) THEN
2391	PRINT*, '+++ WARNING: check_parameters: termination_time_', &
2392	'needed', termination_time_needed
2393	PRINT*, ' should be >= 300.0', &
2394	' on host "', host, '"'
2395	ENDIF
2396	ENDIF
2397	ENDIF
2398
2399
2400	END SUBROUTINE check_parameters

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |