Home

Context Navigation

source: palm/trunk/SOURCE/modules.f90 @ 2107

Last change on this file since 2107 was 2107, checked in by kanani, 7 years ago
Preparation for doxygen comments
Property svn:keywords set to `Id`
File size: 100.1 KB

Line
1	!> @file modules.f90
2	!------------------------------------------------------------------------------!
3	! This file is part of PALM.
4	!
5	! PALM is free software: you can redistribute it and/or modify it under the
6	! terms of the GNU General Public License as published by the Free Software
7	! Foundation, either version 3 of the License, or (at your option) any later
8	! version.
9	!
10	! PALM is distributed in the hope that it will be useful, but WITHOUT ANY
11	! WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
12	! A PARTICULAR PURPOSE. See the GNU General Public License for more details.
13	!
14	! You should have received a copy of the GNU General Public License along with
15	! PALM. If not, see <http://www.gnu.org/licenses/>.
16	!
17	! Copyright 1997-2017 Leibniz Universitaet Hannover
18	!------------------------------------------------------------------------------!
19	!
20	! Current revisions:
21	! ------------------
22	! Preparation for doxygen comments (Giersch)
23	!
24	! Former revisions:
25	! -----------------
26	! $Id: modules.f90 2107 2017-01-09 12:21:49Z kanani $
27	!
28	! 2050 2016-11-08 15:00:55Z gronemeier
29	! Implement turbulent outflow condition
30	!
31	! 2037 2016-10-26 11:15:40Z knoop
32	! Anelastic approximation implemented
33	!
34	! 2031 2016-10-21 15:11:58Z knoop
35	! renamed variable rho to rho_ocean and rho_av to rho_ocean_av
36	!
37	! 2011 2016-09-19 17:29:57Z kanani
38	! +urban_surface, +lsf_exception, +varnamelength
39	!
40	! 2007 2016-08-24 15:47:17Z kanani
41	! Increased DIMENSION of data_output, data_output_user, do2d, do3d
42	!
43	! 2000 2016-08-20 18:09:15Z knoop
44	! Forced header and separation lines into 80 columns
45	!
46	! 1992 2016-08-12 15:14:59Z suehring
47	! +constant_top_scalarflux, top_scalarflux
48	! default of bc_s_t adjusted
49	!
50	! 1968 2016-07-18 12:01:49Z suehring
51	! Changed dimension for MPI-datatypes type_x_int and type_y_int
52	!
53	! 1960 2016-07-12 16:34:24Z suehring
54	! Separate humidity and passive scalar
55	! +bc_s_t_val, diss_s_s, diss_l_s, flux_s_s, flux_l_s, s, sp, s1, s2, s3, ssws_av,
56	! s_init, s_surf, sums_wsss_ws_l, ss, ssws, sswst, ts_m, wall_sflux
57	! +constant_scalarflux, ibc_s_b, ibc_s_t, s_vertical_gradient_level_ind
58	!
59	! Unused variables removed
60	! -gamma_x, gamma_y, gamma_z, var_x, var_y, var_z
61	!
62	! Change initial values (in order to unify gradient calculation):
63	! pt_vertical_gradient_level, sa_vertical_gradient_level
64	!
65	! 1957 2016-07-07 10:43:48Z suehring
66	! +fl_max, num_leg, num_var_fl, num_var_fl_user, var_fl_max, virtual_flight
67	!
68	! 1918 2016-05-27 14:35:57Z raasch
69	! default timestep switched from -1.0 to +1.0 in order to avoid wrong sign of
70	! initially calculated divergence
71	!
72	! 1906 2016-05-24 14:38:08Z suehring
73	! default value of mg_switch_to_pe0_level changed to -1
74	!
75	! 1849 2016-04-08 11:33:18Z hoffmann
76	! bfactor, mass_of_solute, molecular_weight_of_solute, molecular_weight_of_water,
77	! vanthoff moved to mod_particle_attributes.
78	! dt_micro and several cloud_parameters moved to microphysics_mod.
79	! 1d-microphysics profiles moved to microphysics_mod.
80	!
81	! 1845 2016-04-08 08:29:13Z raasch
82	! -nzb_2d
83	!
84	! 1833 2016-04-07 14:23:03Z raasch
85	! spectra parameter moved to spectra module
86	!
87	! 1831 2016-04-07 13:15:51Z hoffmann
88	! curvature_solution_effects removed
89	! turbulence renamed collision_turbulence, drizzle renamed
90	! cloud_water_sedimentation
91	!
92	! 1822 2016-04-07 07:49:42Z hoffmann
93	! icloud_scheme removed. microphysics_sat_adjust, microphysics_kessler,
94	! microphysics_seifert added.
95	!
96	! 1815 2016-04-06 13:49:59Z raasch
97	! cpp-directive for decalpha removed
98	!
99	! 1808 2016-04-05 19:44:00Z raasch
100	! MPI module used by default on all machines
101	!
102	! 1804 2016-04-05 16:30:18Z maronga
103	! Removed code for parameter file check (__check)
104	!
105	! 1788 2016-03-10 11:01:04Z maronga
106	! Added roughness length for moisture (z0q)
107	!
108	! 1786 2016-03-08 05:49:27Z raasch
109	! module spectrum moved to new separate module
110	!
111	! 1783 2016-03-06 18:36:17Z raasch
112	! netcdf variables moved to the netcdf-interface module
113	!
114	! 1779 2016-03-03 08:01:28Z raasch
115	! coupling_char extended to LEN=3
116	!
117	! 1764 2016-02-28 12:45:19Z raasch
118	! some reformatting
119	!
120	! 1762 2016-02-25 12:31:13Z hellstea
121	! +nest_* variables, size of volume_flow arrays increased by one element
122	!
123	! 1738 2015-12-18 13:56:05Z raasch
124	! +mean_surface_level_height
125	!
126	! 1695 2015-10-27 10:03:11Z maronga
127	! Removed rif (forgotten in last revision)
128	!
129	! 1693 2015-10-27 08:35:45Z maronga
130	! Renamed zp -> z_mo
131	!
132	! 1691 2015-10-26 16:17:44Z maronga
133	! Renamed Obukhov length. Added ol, removed rif. Increased number of profiles
134	! (pr_palm). Changed default values for dissipation_1d to 'detering' and
135	! (mixing_length_1d to 'blackadar'. Added most_method. rif_min and rif_max
136	! renamed to zeta_min and zeta_max and new values assigned.
137	!
138	! 1682 2015-10-07 23:56:08Z knoop
139	! Code annotations made doxygen readable
140	!
141	! 1677 2015-10-02 13:25:23Z boeske
142	! +ngp_yz_int, type_xz_int, type_yz_int
143	!
144	! 1666 2015-09-23 07:31:10Z raasch
145	! +user_interface_current_revision, user_interface_required_revision in
146	! control_parameters
147	!
148	! 1639 2015-08-31 14:46:48Z knoop
149	! Bugfix: string 'unknown' extended to match LEN=13
150	!
151	! 1575 2015-03-27 09:56:27Z raasch
152	! +ngp_xz
153	!
154	! 1560 2015-03-06 10:48:54Z keck
155	! +recycling_yshift
156	!
157	! 1557 2015-03-05 16:43:04Z suehring
158	! +monotonic_adjustment
159	!
160	! 1551 2015-03-03 14:18:16Z maronga
161	! Increased pr_palm to 120. Increased length of dots_unit and dots_label to 13
162	! digits. Increased length of domask, do2d, and do3d to 20 digits.
163	!
164	! 1496 2014-12-02 17:25:50Z maronga
165	! Renamed "radiation" -> "cloud_top_radiation"
166	!
167	! 1484 2014-10-21 10:53:05Z kanani
168	! Changes due to new module structure of the plant canopy model:
169	! canopy-model related parameters/variables moved to module
170	! plant_canopy_model_mod
171	!
172	! 1468 2014-09-24 14:06:57Z maronga
173	! Adapted for use on up to 6-digit processor cores.
174	! Increased identifier string length for user-defined quantities to 20.
175	!
176	! 1450 2014-08-21 07:31:51Z heinze
177	! ntnudge from 100 to 1000 increased to allow longer simulations
178	!
179	! 1431 2014-07-15 14:47:17Z suehring
180	! +var_d
181	!
182	! 1429 2014-07-15 12:53:45Z knoop
183	! +ensemble_member_nr to prepare the random_generator for ensemble runs
184	!
185	! 1382 2014-04-30 12:15:41Z boeske
186	! Renamed variables which store large scale forcing tendencies
187	! pt_lsa -> td_lsa_lpt, pt_subs -> td_sub_lpt,
188	! q_lsa -> td_lsa_q, q_subs -> td_sub_q
189	!
190	! 1365 2014-04-22 15:03:56Z boeske
191	! Usage of large scale forcing enabled:
192	! increase pr_palm from 90 to 100 to allow for more standard profiles
193	! + ngp_sums_ls, pt_lsa, pt_subs, q_lsa, q_subs, tmp_tnudge, sums_ls_l,
194	! use_subsidence_tendencies
195	!
196	! 1361 2014-04-16 15:17:48Z hoffmann
197	! tend_* removed
198	! call_microphysics_at_all_substeps added
199	! default of drizzle set to true
200	!
201	! 1359 2014-04-11 17:15:14Z hoffmann
202	! particle_attributes moved to mod_particle_attributes.f90
203	!
204	! 1353 2014-04-08 15:21:23Z heinze
205	! REAL constants provided with KIND-attribute
206	!
207	! 1327 2014-03-21 11:00:16Z raasch
208	! REAL constants defined as wp-kind
209	! -avs_output, data_output_format, do3d_compress, iso2d_output, netcdf_output
210	!
211	! 1320 2014-03-20 08:40:49Z raasch
212	! ONLY-attribute added to USE-statements,
213	! kind-parameters added to all INTEGER and REAL declaration statements,
214	! kinds are defined in new module kinds,
215	! old module precision_kind is removed,
216	! revision history before 2012 removed,
217	! comment fields (!:) to be used for variable explanations added to
218	! all variable declaration statements
219	!
220	! 1318 2014-03-17 13:35:16Z raasch
221	! module cpulog moved to new separate module-file
222	! interface for cpu_log removed
223	!
224	! 1314 2014-03-14 18:25:17Z suehring
225	! + log_z_z0, number_of_sublayers, z0_av_global
226	! 1308 2014-03-13 14:58:42Z fricke
227	! +ntdim_2d_xy, ntdim_2d_xz, ntdim_2d_yz, ntdim_3d
228	!
229	! 1257 2013-11-08 15:18:40Z raasch
230	! set default values for grid indices of maximum velocity components
231	! u\|v\|w_max_ijk
232	!
233	! 1241 2013-10-30 11:36:58Z heinze
234	! Usage of nudging enabled
235	! +nudging, ntnudge, ptnudge, qnudge, tnudge, unudge, vnudge, wnudge
236	! increase pr_palm from 80 to 90 to allow for more standard profiles
237	!
238	! Enable prescribed time depenend surface fluxes and geostrophic wind read in
239	! from external file LSF_DATA
240	! +large_scale_forcing, lsf_surf, lsf_vert, nlsf, time_surf, shf_surf, qsws_surf,
241	! pt_surf, q_surf, p_surf, time_vert, ug_vert, vg_vert, wsubs_vert
242	!
243	! 1221 2013-09-10 08:59:13Z raasch
244	! wall_flags_0 changed to 32bit int, +wall_flags_00,
245	! +rflags_s_inner, rflags_invers
246	!
247	! 1216 2013-08-26 09:31:42Z raasch
248	! +transpose_compute_overlap,
249	! several variables are now defined in the serial (non-parallel) case also
250	!
251	! 1212 2013-08-15 08:46:27Z raasch
252	! +tri
253	!
254	! 1179 2013-06-14 05:57:58Z raasch
255	! +reference_state, ref_state, use_initial_profile_as_reference, vpt_reference,
256	! use_reference renamed use_single_reference_value
257	!
258	! 1159 2013-05-21 11:58:22Z fricke
259	! -bc_lr_dirneu, bc_lr_neudir, bc_ns_dirneu, bc_ns_neudir
260	! +use_cmax
261	!
262	! 1128 2013-04-12 06:19:32Z raasch
263	! +background_communication, i_left, i_right, j_north, j_south, req, req_count,
264	! send_receive, sendrecv_in_background, wait_stat
265	!
266	! 1115 2013-03-26 18:16:16Z hoffmann
267	! unused variables removed
268	!
269	! 1113 2013-03-10 02:48:14Z raasch
270	! +on_device
271	!
272	! 1111 2013-03-08 23:54:10Z raasch
273	! +tric, nr_timesteps_this_run
274	!
275	! 1106 2013-03-04 05:31:38Z raasch
276	! array_kind renamed precision_kind, pdims defined in serial code
277	! bugfix: default value assigned to coupling_start_time
278	!
279	! 1095 2013-02-03 02:21:01Z raasch
280	! FORTRAN error in r1092 removed
281	!
282	! 1092 2013-02-02 11:24:22Z raasch
283	! character length in some derived type changed for better alignment
284	!
285	! 1065 2012-11-22 17:42:36Z hoffmann
286	! + c_sedimentation, limiter_sedimentation, turbulence, a_1, a_2, a_3, b_1, b_2,
287	! + b_3, c_1, c_2, c_3, beta_cc
288	!
289	! bottom boundary condition of qr, nr changed from Dirichlet to Neumann
290	!
291	! 1053 2012-11-13 17:11:03Z hoffmann
292	! necessary expansions according to the two new prognostic equations (nr, qr)
293	! of the two-moment cloud physics scheme:
294	! +_init, flux_l_, diss_l_, flux_s_, diss_s_, sws, swst, tend_, , _p
295	! +t_m, _1, _2, _3, _av, bc__b, bc__t, ibc__b, ibc__t, bc__t_val,
296	! +_vertical_gradient, _surface_initial_change, *_vertical_gradient_level,
297	! +_vertical_gradient_level_ind, _surface, constant_waterflux_*,
298	! +cloud_scheme, icloud_scheme, surface_waterflux_, sums_wss_ws_l, wall_*flux
299	!
300	! constants for the two-moment scheme:
301	! +a_vent, a_term, b_vent, b_term, c_evap, c_term, cof, eps_sb, k_cc, k_cr, k_rr,
302	! +k_br, kappa_rr, kin_vis_air, mu_constant_value, nc, pirho_l, dpirho_l, rho_1,
303	! +schmidt, schmidt_p_1d3, stp, x0, xmin, xmax, dt_precipitation, w_precipitation
304	!
305	! steering parameters for the two_moment scheme:
306	! +mu_constant, ventilation_effect
307	!
308	! 1036 2012-10-22 13:43:42Z raasch
309	! code put under GPL (PALM 3.9)
310	!
311	! 1031 2012-10-19 14:35:30Z raasch
312	! +output_format_netcdf
313	!
314	! 1015 2012-09-27 09:23:24Z raasch
315	! +acc_rank, num_acc_per_node,
316	! -adjust_mixing_length
317	!
318	! 1010 2012-09-20 07:59:54Z raasch
319	! pointer free version can be generated with cpp switch __nopointer
320	!
321	! 1003 2012-09-14 14:35:53Z raasch
322	! -grid_matching, nxa, nya, etc., nnx_pe, nny_pe, spl_*
323	!
324	! 1001 2012-09-13 14:08:46Z raasch
325	! -asselin_filter_factor, cut_spline_overshoot, dt_changed, last_dt_change,
326	! last_dt_change_1d, long_filter_factor, overshoot_limit_, ups_limit_
327	! several pointer/target arrays converted to normal ones
328	!
329	! 996 2012-09-07 10:41:47Z raasch
330	! -use_prior_plot1d_parameters
331	!
332	! 978 2012-08-09 08:28:32Z fricke
333	! +c_u_m, c_u_m_l, c_v_m, c_v_m_l, c_w_m, c_w_m_l,
334	! +bc_lr_dirneu, bc_lr_neudir, bc_ns_dirneu, bc_ns_neudir
335	! -km_damp_x, km_damp_y, km_damp_max, outflow_damping_width
336	! +z0h, z0h_av, z0h_factor, z0h1d
337	! +ptdf_x, ptdf_y, pt_damping_width, pt_damping_factor
338	!
339	! 964 2012-07-26 09:14:24Z raasch
340	! -cross_linecolors, cross_linestyles, cross_normalized_x, cross_normx_factor,
341	! cross_normalized_y, cross_normy_factor, cross_pnc_local,
342	! cross_profile_numbers, cross_profile_number_counter, cross_uxmax,
343	! cross_uxmax_computed, cross_uxmax_normalized,
344	! cross_uxmax_normalized_computed, cross_uxmin, cross_uxmin_computed,
345	! cross_uxmin_normalized, cross_uxmin_normalized_computed, cross_uymax,
346	! cross_uymin, cross_xtext, dopr_crossindex, dopr_label, linecolors, linestyles,
347	! nz_do1d, profil_output, z_max_do1d, z_max_do1d_normalized
348	!
349	! 951 2012-07-19 14:22:52Z hoffmann
350	! changing profile_columns and profile_rows
351	!
352	! 940 2012-07-09 14:31:00Z raasch
353	! +neutral
354	!
355	! 927 2012-06-06 19:15:04Z raasch
356	! +masking_method
357	!
358	! 880 2012-04-13 06:28:59Z raasch
359	! gathered_size, subdomain_size moved to control_parameters
360	!
361	! 866 2012-03-28 06:44:41Z raasch
362	! interface for global_min_max changed
363	!
364	! 861 2012-03-26 14:18:34Z suehring
365	! +wall_flags_0
366	! -boundary_flags
367	! +nzb_max
368	! +adv_sca_1, +adv_mom_1
369	!
370	! 849 2012-03-15 10:35:09Z raasch
371	! +deleted_particles, deleted_tails, tr.._count_sum, tr.._count_recv_sum in
372	! particle_attributes,
373	! +de_dx, de_dy, de_dz in arrays_3d,
374	! first_call_advec_particles renamed first_call_lpm
375	!
376	! 828 2012-02-21 12:00:36Z raasch
377	! +dissipation_classes, radius_classes, use_kernel_tables,
378	! particle feature color renamed class
379	!
380	! 825 2012-02-19 03:03:44Z raasch
381	! +bfactor, curvature_solution_effects, eps_ros, molecular_weight_of_water,
382	! vanthoff, -b_cond in cloud_parameters,
383	! dopts_num increased to 29, particle attributes speed_x\|y\|z_sgs renamed
384	! rvar1\|2\|3
385	! wang_collision_kernel and turbulence_effects_on_collision replaced by
386	! collision_kernel, hall_kernel, palm_kernel, wang_kernel
387	!
388	! 809 2012-01-30 13:32:58Z marongas
389	! Bugfix: replaced .AND. and .NOT. with && and ! in the preprocessor directives
390	!
391	! 807 2012-01-25 11:53:51Z maronga
392	! New cpp directive "__check" implemented which is used by check_namelist_files.
393	! New parameter check_restart has been defined which is needed by
394	! check_namelist_files only.
395	!
396	! 805 2012-01-17 15:53:28Z franke
397	! Bugfix collective_wait must be out of parallel branch for runs in serial mode
398	!
399	! 801 2012-01-10 17:30:36Z suehring
400	! Dimesion of sums_wsus_ws_l, ! sums_wsvs_ws_l, sums_us2_ws_l, sums_vs2_ws_l,
401	! sums_ws2_ws_l, sums_wspts_ws_l, sums_wsqs_ws_l, sums_wssas_ws_l increased.
402	! for thread-safe summation in advec_ws.
403	!
404	! RCS Log replace by Id keyword, revision history cleaned up
405	!
406	! Revision 1.95 2007/02/11 13:18:30 raasch
407	! version 3.1b (last under RCS control)
408	!
409	! Revision 1.1 1997/07/24 11:21:26 raasch
410	! Initial revision
411	!
412	!
413	!------------------------------------------------------------------------------!
414	! Description:
415	! ------------
416	!> Definition of all variables
417	!>
418	!> @todo Add description for each variable
419	!------------------------------------------------------------------------------!
420
421
422	!------------------------------------------------------------------------------!
423	! Description:
424	! ------------
425	!> Definition of variables for special advection schemes.
426	!------------------------------------------------------------------------------!
427	MODULE advection
428
429	USE kinds
430
431	REAL(wp), DIMENSION(:), ALLOCATABLE :: aex !<
432	REAL(wp), DIMENSION(:), ALLOCATABLE :: bex !<
433	REAL(wp), DIMENSION(:), ALLOCATABLE :: dex !<
434	REAL(wp), DIMENSION(:), ALLOCATABLE :: eex !<
435
436	SAVE
437
438	END MODULE advection
439
440
441	!------------------------------------------------------------------------------!
442	! Description:
443	! ------------
444	!> Definition of all arrays defined on the computational grid.
445	!------------------------------------------------------------------------------!
446	MODULE arrays_3d
447
448	USE kinds
449
450	REAL(wp), DIMENSION(:), ALLOCATABLE :: c_u_m !<
451	REAL(wp), DIMENSION(:), ALLOCATABLE :: c_u_m_l !<
452	REAL(wp), DIMENSION(:), ALLOCATABLE :: c_v_m !<
453	REAL(wp), DIMENSION(:), ALLOCATABLE :: c_v_m_l !<
454	REAL(wp), DIMENSION(:), ALLOCATABLE :: c_w_m !<
455	REAL(wp), DIMENSION(:), ALLOCATABLE :: c_w_m_l !<
456	REAL(wp), DIMENSION(:), ALLOCATABLE :: ddzu !<
457	REAL(wp), DIMENSION(:), ALLOCATABLE :: ddzu_pres !<
458	REAL(wp), DIMENSION(:), ALLOCATABLE :: dd2zu !<
459	REAL(wp), DIMENSION(:), ALLOCATABLE :: dzu !<
460	REAL(wp), DIMENSION(:), ALLOCATABLE :: ddzw !<
461	REAL(wp), DIMENSION(:), ALLOCATABLE :: dzw !<
462	REAL(wp), DIMENSION(:), ALLOCATABLE :: hyp !<
463	REAL(wp), DIMENSION(:), ALLOCATABLE :: inflow_damping_factor !<
464	REAL(wp), DIMENSION(:), ALLOCATABLE :: l_grid !<
465	REAL(wp), DIMENSION(:), ALLOCATABLE :: ptdf_x !<
466	REAL(wp), DIMENSION(:), ALLOCATABLE :: ptdf_y !<
467	REAL(wp), DIMENSION(:), ALLOCATABLE :: p_surf !<
468	REAL(wp), DIMENSION(:), ALLOCATABLE :: pt_surf !<
469	REAL(wp), DIMENSION(:), ALLOCATABLE :: pt_init !<
470	REAL(wp), DIMENSION(:), ALLOCATABLE :: qsws_surf !<
471	REAL(wp), DIMENSION(:), ALLOCATABLE :: q_init !<
472	REAL(wp), DIMENSION(:), ALLOCATABLE :: q_surf !<
473	REAL(wp), DIMENSION(:), ALLOCATABLE :: rdf !<
474	REAL(wp), DIMENSION(:), ALLOCATABLE :: rdf_sc !<
475	REAL(wp), DIMENSION(:), ALLOCATABLE :: ref_state !<
476	REAL(wp), DIMENSION(:), ALLOCATABLE :: s_init !<
477	REAL(wp), DIMENSION(:), ALLOCATABLE :: s_surf !<
478	REAL(wp), DIMENSION(:), ALLOCATABLE :: sa_init !<
479	REAL(wp), DIMENSION(:), ALLOCATABLE :: shf_surf !<
480	REAL(wp), DIMENSION(:), ALLOCATABLE :: timenudge !<
481	REAL(wp), DIMENSION(:), ALLOCATABLE :: time_surf !<
482	REAL(wp), DIMENSION(:), ALLOCATABLE :: time_vert !<
483	REAL(wp), DIMENSION(:), ALLOCATABLE :: tmp_tnudge !<
484	REAL(wp), DIMENSION(:), ALLOCATABLE :: ug !<
485	REAL(wp), DIMENSION(:), ALLOCATABLE :: u_init !<
486	REAL(wp), DIMENSION(:), ALLOCATABLE :: u_nzb_p1_for_vfc !<
487	REAL(wp), DIMENSION(:), ALLOCATABLE :: vg !<
488	REAL(wp), DIMENSION(:), ALLOCATABLE :: v_init !<
489	REAL(wp), DIMENSION(:), ALLOCATABLE :: v_nzb_p1_for_vfc !<
490	REAL(wp), DIMENSION(:), ALLOCATABLE :: w_subs !<
491	REAL(wp), DIMENSION(:), ALLOCATABLE :: zu !<
492	REAL(wp), DIMENSION(:), ALLOCATABLE :: zw !<
493
494	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: c_u !<
495	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: c_v !<
496	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: c_w !<
497	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: diss_s_e !<
498	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: diss_s_nr !<
499	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: diss_s_pt !<
500	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: diss_s_q !<
501	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: diss_s_qr !<
502	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: diss_s_s !<
503	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: diss_s_sa !<
504	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: diss_s_u !<
505	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: diss_s_v !<
506	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: diss_s_w !<
507	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: dzu_mg !<
508	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: dzw_mg !<
509	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: flux_s_e !<
510	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: flux_s_nr !<
511	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: flux_s_pt !<
512	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: flux_s_q !<
513	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: flux_s_qr !<
514	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: flux_s_s !<
515	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: flux_s_sa !<
516	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: flux_s_u !<
517	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: flux_s_v !<
518	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: flux_s_w !<
519	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: f1_mg !<
520	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: f2_mg !<
521	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: f3_mg !<
522	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: mean_inflow_profiles !<
523	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: nrs !<
524	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: nrsws !<
525	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: nrswst !<
526	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: ol !< Obukhov length
527	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: precipitation_amount !<
528	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: precipitation_rate !<
529	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: ptnudge !<
530	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: pt_slope_ref !<
531	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: qnudge !<
532	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: qs !<
533	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: qsws !<
534	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: qswst !<
535	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: qswst_remote !<
536	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: qrs !<
537	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: qrsws !<
538	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: qrswst !<
539	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: saswsb !<
540	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: saswst !<
541	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: shf !<
542	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: ss !<
543	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: ssws !<
544	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: sswst !<
545	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: tnudge !<
546	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: td_lsa_lpt !<
547	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: td_lsa_q !<
548	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: td_sub_lpt !<
549	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: td_sub_q !<
550	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: total_2d_a !<
551	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: total_2d_o !<
552	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: ts !<
553	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: tswst !<
554	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: ug_vert !<
555	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: unudge !<
556	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: us !<
557	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: usws !<
558	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: uswst !<
559	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: vnudge !<
560	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: vg_vert !<
561	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: vsws !<
562	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: vswst !<
563	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: wnudge !<
564	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: wsubs_vert !<
565	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: z0 !< roughness length for momentum
566	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: z0h !< roughness length for heat
567	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: z0q !< roughness length for moisture
568
569	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: d !<
570	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: de_dx !<
571	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: de_dy !<
572	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: de_dz !<
573	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: diss !<
574	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: diss_l_e !<
575	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: diss_l_nr !<
576	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: diss_l_pt !<
577	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: diss_l_q !<
578	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: diss_l_qr !<
579	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: diss_l_s !<
580	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: diss_l_sa !<
581	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: diss_l_u !<
582	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: diss_l_v !<
583	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: diss_l_w !<
584	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: flux_l_e !<
585	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: flux_l_nr !<
586	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: flux_l_pt !<
587	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: flux_l_q !<
588	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: flux_l_qr !<
589	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: flux_l_s !<
590	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: flux_l_sa !<
591	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: flux_l_u !<
592	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: flux_l_v !<
593	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: flux_l_w !<
594	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: kh !<
595	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: km !<
596	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: l_wall !<
597	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: prr !<
598	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: p_loc !<
599	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: tend !<
600	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: tric !<
601	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: u_m_l !<
602	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: u_m_n !<
603	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: u_m_r !<
604	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: u_m_s !<
605	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: v_m_l !<
606	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: v_m_n !<
607	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: v_m_r !<
608	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: v_m_s !<
609	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: w_m_l !<
610	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: w_m_n !<
611	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: w_m_r !<
612	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: w_m_s !<
613
614	#if defined( __nopointer )
615	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: e !<
616	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: e_p !<
617	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: nr !<
618	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: nr_p !<
619	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: p !<
620	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: prho !<
621	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: pt !<
622	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: pt_p !<
623	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: q !<
624	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: q_p !<
625	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: qc !<
626	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: ql !<
627	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: ql_c !<
628	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: ql_v !<
629	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: ql_vp !<
630	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: qr !<
631	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: qr_p !<
632	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: rho_ocean !<
633	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: s !<
634	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: s_p !<
635	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: sa !<
636	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: sa_p !<
637	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: te_m !<
638	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: tnr_m !<
639	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: tpt_m !<
640	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: tq_m !<
641	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: tqr_m !<
642	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: ts_m !<
643	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: tsa_m !<
644	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: tu_m !<
645	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: tv_m !<
646	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: tw_m !<
647	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: u !<
648	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: u_p !<
649	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: v !<
650	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: v_p !<
651	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: vpt !<
652	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: w !<
653	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: w_p !<
654	#else
655	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: e_1 !<
656	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: e_2 !<
657	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: e_3 !<
658	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: p !<
659	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: prho_1 !<
660	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: nr_1 !<
661	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: nr_2 !<
662	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: nr_3 !<
663	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: pt_1 !<
664	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: pt_2 !<
665	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: pt_3 !<
666	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: q_1 !<
667	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: q_2 !<
668	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: q_3 !<
669	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: qc_1 !<
670	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: ql_v !<
671	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: ql_vp !<
672	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: ql_1 !<
673	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: ql_2 !<
674	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: qr_1 !<
675	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: qr_2 !<
676	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: qr_3 !<
677	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: rho_1 !<
678	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: s_1 !<
679	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: s_2 !<
680	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: s_3 !<
681	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: sa_1 !<
682	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: sa_2 !<
683	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: sa_3 !<
684	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: u_1 !<
685	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: u_2 !<
686	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: u_3 !<
687	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: v_1 !<
688	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: v_2 !<
689	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: v_3 !<
690	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: vpt_1 !<
691	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: w_1 !<
692	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: w_2 !<
693	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: w_3 !<
694
695	REAL(wp), DIMENSION(:,:,:), POINTER :: e !<
696	REAL(wp), DIMENSION(:,:,:), POINTER :: e_p !<
697	REAL(wp), DIMENSION(:,:,:), POINTER :: nr !<
698	REAL(wp), DIMENSION(:,:,:), POINTER :: nr_p !<
699	REAL(wp), DIMENSION(:,:,:), POINTER :: prho !<
700	REAL(wp), DIMENSION(:,:,:), POINTER :: pt !<
701	REAL(wp), DIMENSION(:,:,:), POINTER :: pt_p !<
702	REAL(wp), DIMENSION(:,:,:), POINTER :: q !<
703	REAL(wp), DIMENSION(:,:,:), POINTER :: q_p !<
704	REAL(wp), DIMENSION(:,:,:), POINTER :: qc !<
705	REAL(wp), DIMENSION(:,:,:), POINTER :: ql !<
706	REAL(wp), DIMENSION(:,:,:), POINTER :: ql_c !<
707	REAL(wp), DIMENSION(:,:,:), POINTER :: qr !<
708	REAL(wp), DIMENSION(:,:,:), POINTER :: qr_p !<
709	REAL(wp), DIMENSION(:,:,:), POINTER :: rho_ocean !<
710	REAL(wp), DIMENSION(:,:,:), POINTER :: s !<
711	REAL(wp), DIMENSION(:,:,:), POINTER :: s_p !<
712	REAL(wp), DIMENSION(:,:,:), POINTER :: sa !<
713	REAL(wp), DIMENSION(:,:,:), POINTER :: sa_p !<
714	REAL(wp), DIMENSION(:,:,:), POINTER :: te_m !<
715	REAL(wp), DIMENSION(:,:,:), POINTER :: tnr_m !<
716	REAL(wp), DIMENSION(:,:,:), POINTER :: tpt_m !<
717	REAL(wp), DIMENSION(:,:,:), POINTER :: tq_m !<
718	REAL(wp), DIMENSION(:,:,:), POINTER :: tqr_m !<
719	REAL(wp), DIMENSION(:,:,:), POINTER :: ts_m !<
720	REAL(wp), DIMENSION(:,:,:), POINTER :: tsa_m !<
721	REAL(wp), DIMENSION(:,:,:), POINTER :: tu_m !<
722	REAL(wp), DIMENSION(:,:,:), POINTER :: tv_m !<
723	REAL(wp), DIMENSION(:,:,:), POINTER :: tw_m !<
724	REAL(wp), DIMENSION(:,:,:), POINTER :: u !<
725	REAL(wp), DIMENSION(:,:,:), POINTER :: u_p !<
726	REAL(wp), DIMENSION(:,:,:), POINTER :: v !<
727	REAL(wp), DIMENSION(:,:,:), POINTER :: v_p !<
728	REAL(wp), DIMENSION(:,:,:), POINTER :: vpt !<
729	REAL(wp), DIMENSION(:,:,:), POINTER :: w !<
730	REAL(wp), DIMENSION(:,:,:), POINTER :: w_p !<
731	#endif
732
733	REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: rif_wall !<
734	REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: tri !<
735
736	REAL(wp), DIMENSION(:), ALLOCATABLE :: rho_air !< air density profile on the uv grid
737	REAL(wp), DIMENSION(:), ALLOCATABLE :: rho_air_zw !< air density profile on the w grid
738	REAL(wp), DIMENSION(:), ALLOCATABLE :: drho_air !< inverse air density profile on the uv grid
739	REAL(wp), DIMENSION(:), ALLOCATABLE :: drho_air_zw !< inverse air density profile on the w grid
740
741	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: rho_air_mg !< air density profiles on the uv grid for multigrid
742	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: rho_air_zw_mg !< air density profiles on the w grid for multigrid
743
744	REAL(wp), DIMENSION(:), ALLOCATABLE :: heatflux_input_conversion !< conversion factor array for heatflux input
745	REAL(wp), DIMENSION(:), ALLOCATABLE :: waterflux_input_conversion !< conversion factor array for waterflux input
746	REAL(wp), DIMENSION(:), ALLOCATABLE :: momentumflux_input_conversion !< conversion factor array for momentumflux input
747	REAL(wp), DIMENSION(:), ALLOCATABLE :: heatflux_output_conversion !< conversion factor array for heatflux output
748	REAL(wp), DIMENSION(:), ALLOCATABLE :: waterflux_output_conversion !< conversion factor array for waterflux output
749	REAL(wp), DIMENSION(:), ALLOCATABLE :: momentumflux_output_conversion !< conversion factor array for momentumflux output
750
751	SAVE
752
753	END MODULE arrays_3d
754
755
756	!------------------------------------------------------------------------------!
757	! Description:
758	! ------------
759	!> Definition of variables needed for time-averaging of 2d/3d data.
760	!------------------------------------------------------------------------------!
761	MODULE averaging
762
763	USE kinds
764
765	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: lwp_av !< Avg. liquid water path
766	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: precipitation_rate_av !< Avg. of precipitation rate
767	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: ol_av !< Avg. of Obukhov length
768	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: qsws_av !< Avg. of surface moisture flux
769	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: ssws_av !< Avg. of surface scalar flux
770	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: shf_av !< Avg. of surface heat flux
771	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: ts_av !< Avg. of characteristic temperature scale
772	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: us_av !< Avg. of friction velocity
773	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: z0_av !< Avg. of roughness length for momentum
774	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: z0h_av !< Avg. of roughness length for heat
775	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: z0q_av !< Avg. of roughness length for moisture
776
777	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: e_av !<
778	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: lpt_av !<
779	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: nr_av !<
780	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: p_av !<
781	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: pc_av !<
782	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: pr_av !<
783	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: prr_av !<
784	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: pt_av !<
785	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: q_av !<
786	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: qc_av !<
787	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: ql_av !<
788	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: ql_c_av !<
789	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: ql_v_av !<
790	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: ql_vp_av !<
791	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: qr_av !<
792	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: qv_av !<
793	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: rho_ocean_av !<
794	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: s_av !<
795	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: sa_av !<
796	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: u_av !<
797	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: v_av !<
798	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: vpt_av !<
799	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: w_av !<
800
801	END MODULE averaging
802
803
804	!------------------------------------------------------------------------------!
805	! Description:
806	! ------------
807	!> Definition of variables and constants for cloud physics.
808	!------------------------------------------------------------------------------!
809	MODULE cloud_parameters
810
811	USE kinds
812
813	REAL(wp) :: cp = 1005.0_wp !< heat capacity of dry air (J kg-1 K-1)
814	REAL(wp) :: l_v = 2.5E+06_wp !< latent heat of vaporization (J kg-1)
815	REAL(wp) :: l_d_cp !< l_v / cp
816	REAL(wp) :: l_d_r !< l_v / r_d
817	REAL(wp) :: l_d_rv !< l_v / r_v
818	REAL(wp) :: rho_l = 1.0E3_wp !< density of water (kg m-3)
819	REAL(wp) :: r_d = 287.0_wp !< sp. gas const. dry air (J kg-1 K-1)
820	REAL(wp) :: r_v = 461.51_wp !< sp. gas const. water vapor (J kg-1 K-1)
821
822
823	REAL(wp), DIMENSION(:), ALLOCATABLE :: hyrho !<
824	REAL(wp), DIMENSION(:), ALLOCATABLE :: pt_d_t !<
825	REAL(wp), DIMENSION(:), ALLOCATABLE :: t_d_pt !<
826
827	SAVE
828
829	END MODULE cloud_parameters
830
831
832	!------------------------------------------------------------------------------!
833	! Description:
834	! ------------
835	!> Definition of general constants.
836	!------------------------------------------------------------------------------!
837	MODULE constants
838
839	USE kinds
840
841	REAL(wp) :: pi = 3.141592654_wp !<
842	REAL(wp) :: adv_mom_1 !<
843	REAL(wp) :: adv_mom_3 !<
844	REAL(wp) :: adv_mom_5 !<
845	REAL(wp) :: adv_sca_1 !<
846	REAL(wp) :: adv_sca_3 !<
847	REAL(wp) :: adv_sca_5 !<
848
849	SAVE
850
851	END MODULE constants
852
853
854	!------------------------------------------------------------------------------!
855	! Description:
856	! ------------
857	!> Definition of parameters for program control
858	!------------------------------------------------------------------------------!
859	MODULE control_parameters
860
861	USE kinds
862
863	TYPE plot_precision
864	CHARACTER (LEN=8) :: variable !<
865	INTEGER(iwp) :: precision !<
866	END TYPE plot_precision
867
868	TYPE(plot_precision), DIMENSION(100) :: plot_3d_precision = & !<
869	(/ plot_precision( 'u', 2 ), plot_precision( 'v', 2 ), &
870	plot_precision( 'w', 2 ), plot_precision( 'p', 5 ), &
871	plot_precision( 'pt', 2 ), &
872	( plot_precision( ' ', 1 ), i9 = 1,95 ) /)
873
874	TYPE file_status
875	LOGICAL :: opened !<
876	LOGICAL :: opened_before !<
877	END TYPE file_status
878
879	INTEGER, PARAMETER :: mask_xyz_dimension = 100 !<
880	INTEGER, PARAMETER :: max_masks = 50 !<
881	INTEGER(iwp), PARAMETER :: varnamelength = 30 !< length of output variable names
882
883	TYPE(file_status), DIMENSION(200+2*max_masks) :: & !<
884	openfile = file_status(.FALSE.,.FALSE.)
885
886	CHARACTER (LEN=1) :: cycle_mg = 'w' !<
887	CHARACTER (LEN=1) :: timestep_reason = ' ' !<
888	CHARACTER (LEN=3) :: coupling_char = '' !<
889	CHARACTER (LEN=5) :: write_binary = 'false' !<
890	CHARACTER (LEN=8) :: most_method = 'lookup' !< NAMELIST parameter defining method to be used to calculate Okukhov length,
891	CHARACTER (LEN=8) :: run_date !<
892	CHARACTER (LEN=8) :: run_time !<
893	CHARACTER (LEN=9) :: simulated_time_chr !<
894	CHARACTER (LEN=11) :: topography_grid_convention = ' ' !<
895	CHARACTER (LEN=12) :: version = ' ' !<
896	CHARACTER (LEN=12) :: revision = ' ' !<
897	CHARACTER (LEN=12) :: user_interface_current_revision = ' ' !<
898	CHARACTER (LEN=12) :: user_interface_required_revision = ' ' !<
899	CHARACTER (LEN=16) :: conserve_volume_flow_mode = 'default' !<
900	CHARACTER (LEN=16) :: loop_optimization = 'default' !<
901	CHARACTER (LEN=16) :: momentum_advec = 'ws-scheme' !<
902	CHARACTER (LEN=16) :: psolver = 'poisfft' !<
903	CHARACTER (LEN=16) :: scalar_advec = 'ws-scheme' !<
904	CHARACTER (LEN=20) :: approximation = 'boussinesq' !<
905	CHARACTER (LEN=40) :: flux_input_mode = 'approximation-specific' !<
906	CHARACTER (LEN=40) :: flux_output_mode = 'approximation-specific' !<
907	CHARACTER (LEN=20) :: bc_e_b = 'neumann' !<
908	CHARACTER (LEN=20) :: bc_lr = 'cyclic' !<
909	CHARACTER (LEN=20) :: bc_ns = 'cyclic' !<
910	CHARACTER (LEN=20) :: bc_p_b = 'neumann' !<
911	CHARACTER (LEN=20) :: bc_p_t = 'dirichlet' !<
912	CHARACTER (LEN=20) :: bc_pt_b = 'dirichlet' !<
913	CHARACTER (LEN=20) :: bc_pt_t = 'initial_gradient' !<
914	CHARACTER (LEN=20) :: bc_q_b = 'dirichlet' !<
915	CHARACTER (LEN=20) :: bc_q_t = 'neumann' !<
916	CHARACTER (LEN=20) :: bc_s_b = 'dirichlet' !<
917	CHARACTER (LEN=20) :: bc_s_t = 'initial_gradient' !<
918	CHARACTER (LEN=20) :: bc_sa_t = 'neumann' !<
919	CHARACTER (LEN=20) :: bc_uv_b = 'dirichlet' !<
920	CHARACTER (LEN=20) :: bc_uv_t = 'dirichlet' !<
921	CHARACTER (LEN=20) :: cloud_scheme = 'saturation_adjust' !<
922	CHARACTER (LEN=20) :: coupling_mode = 'uncoupled' !<
923	CHARACTER (LEN=20) :: coupling_mode_remote = 'uncoupled' !<
924	CHARACTER (LEN=20) :: dissipation_1d = 'detering' !<
925	CHARACTER (LEN=20) :: fft_method = 'system-specific' !<
926	CHARACTER (LEN=20) :: mixing_length_1d = 'blackadar' !<
927	CHARACTER (LEN=20) :: random_generator = 'numerical-recipes' !<
928	CHARACTER (LEN=20) :: reference_state = 'initial_profile' !<
929	CHARACTER (LEN=20) :: return_addres !<
930	CHARACTER (LEN=20) :: return_username !<
931	CHARACTER (LEN=20) :: timestep_scheme = 'runge-kutta-3' !<
932	CHARACTER (LEN=40) :: avs_data_file !<
933	CHARACTER (LEN=40) :: topography = 'flat' !<
934	CHARACTER (LEN=64) :: host = ' ' !<
935	CHARACTER (LEN=80) :: log_message !<
936	CHARACTER (LEN=80) :: run_identifier !<
937	CHARACTER (LEN=100) :: initializing_actions = ' ' !<
938	CHARACTER (LEN=110) :: run_description_header !<
939	CHARACTER (LEN=1000) :: message_string = ' ' !<
940
941	CHARACTER (LEN=varnamelength), DIMENSION(500) :: data_output = ' ' !<
942	CHARACTER (LEN=varnamelength), DIMENSION(500) :: data_output_user = ' ' !<
943	CHARACTER (LEN=varnamelength), DIMENSION(500) :: doav = ' ' !<
944
945	CHARACTER (LEN=varnamelength), DIMENSION(max_masks,100) :: data_output_masks = ' ' !<
946	CHARACTER (LEN=varnamelength), DIMENSION(max_masks,100) :: data_output_masks_user = ' ' !<
947
948	CHARACTER (LEN=varnamelength), DIMENSION(300) :: data_output_pr = ' ' !<
949
950	CHARACTER (LEN=varnamelength), DIMENSION(200) :: data_output_pr_user = ' ' !<
951
952	CHARACTER (LEN=varnamelength), DIMENSION(max_masks,0:1,100) :: domask = ' ' !<
953
954	CHARACTER (LEN=varnamelength), DIMENSION(0:1,500) :: do2d = ' ' !<
955	CHARACTER (LEN=varnamelength), DIMENSION(0:1,500) :: do3d = ' ' !<
956
957	INTEGER(iwp), PARAMETER :: fl_max = 100 !<
958	INTEGER(iwp), PARAMETER :: var_fl_max = 20 !<
959
960	INTEGER(iwp) :: abort_mode = 1 !<
961	INTEGER(iwp) :: average_count_pr = 0 !<
962	INTEGER(iwp) :: average_count_3d = 0 !<
963	INTEGER(iwp) :: current_timestep_number = 0 !<
964	INTEGER(iwp) :: coupling_topology = 0 !<
965	INTEGER(iwp) :: dist_range = 0 !<
966	INTEGER(iwp) :: disturbance_level_ind_b !<
967	INTEGER(iwp) :: disturbance_level_ind_t !<
968	INTEGER(iwp) :: doav_n = 0 !<
969	INTEGER(iwp) :: dopr_n = 0 !<
970	INTEGER(iwp) :: dopr_time_count = 0 !<
971	INTEGER(iwp) :: dopts_time_count = 0 !<
972	INTEGER(iwp) :: dots_time_count = 0 !<
973	INTEGER(iwp) :: do2d_xy_n = 0 !<
974	INTEGER(iwp) :: do2d_xz_n = 0 !<
975	INTEGER(iwp) :: do2d_yz_n = 0 !<
976	INTEGER(iwp) :: do3d_avs_n = 0 !<
977	INTEGER(iwp) :: dp_level_ind_b = 0 !<
978	INTEGER(iwp) :: dvrp_filecount = 0 !<
979	INTEGER(iwp) :: dz_stretch_level_index !<
980	INTEGER(iwp) :: ensemble_member_nr = 0 !<
981	INTEGER(iwp) :: gamma_mg !<
982	INTEGER(iwp) :: gathered_size !<
983	INTEGER(iwp) :: grid_level !<
984	INTEGER(iwp) :: ibc_e_b !<
985	INTEGER(iwp) :: ibc_p_b !<
986	INTEGER(iwp) :: ibc_p_t !<
987	INTEGER(iwp) :: ibc_pt_b !<
988	INTEGER(iwp) :: ibc_pt_t !<
989	INTEGER(iwp) :: ibc_q_b !<
990	INTEGER(iwp) :: ibc_q_t !<
991	INTEGER(iwp) :: ibc_s_b !<
992	INTEGER(iwp) :: ibc_s_t !<
993	INTEGER(iwp) :: ibc_sa_t !<
994	INTEGER(iwp) :: ibc_uv_b !<
995	INTEGER(iwp) :: ibc_uv_t !<
996	INTEGER(iwp) :: inflow_disturbance_begin = -1 !<
997	INTEGER(iwp) :: inflow_disturbance_end = -1 !<
998	INTEGER(iwp) :: intermediate_timestep_count !<
999	INTEGER(iwp) :: intermediate_timestep_count_max !<
1000	INTEGER(iwp) :: io_group = 0 !<
1001	INTEGER(iwp) :: io_blocks = 1 !<
1002	INTEGER(iwp) :: iran = -1234567 !<
1003	INTEGER(iwp) :: masks = 0 !<
1004	INTEGER(iwp) :: maximum_grid_level !<
1005	INTEGER(iwp) :: maximum_parallel_io_streams = -1 !<
1006	INTEGER(iwp) :: max_pr_user = 0 !<
1007	INTEGER(iwp) :: mgcycles = 0 !<
1008	INTEGER(iwp) :: mg_cycles = -1 !<
1009	INTEGER(iwp) :: mg_switch_to_pe0_level = -1 !<
1010	INTEGER(iwp) :: mid !<
1011	INTEGER(iwp) :: nlsf = 1000 !<
1012	INTEGER(iwp) :: ntnudge = 1000 !<
1013	INTEGER(iwp) :: ngsrb = 2 !<
1014	INTEGER(iwp) :: nr_timesteps_this_run = 0 !<
1015	INTEGER(iwp) :: nsor = 20 !<
1016	INTEGER(iwp) :: nsor_ini = 100 !<
1017	INTEGER(iwp) :: n_sor !<
1018	INTEGER(iwp) :: normalizing_region = 0 !<
1019	INTEGER(iwp) :: num_leg=0 !<
1020	INTEGER(iwp) :: num_var_fl !<
1021	INTEGER(iwp) :: num_var_fl_user=0 !<
1022	INTEGER(iwp) :: nz_do3d = -9999 !<
1023	INTEGER(iwp) :: prt_time_count = 0 !<
1024	INTEGER(iwp) :: recycling_plane !<
1025	INTEGER(iwp) :: runnr = 0 !<
1026	INTEGER(iwp) :: skip_do_avs = 0 !<
1027	INTEGER(iwp) :: subdomain_size !<
1028	INTEGER(iwp) :: terminate_coupled = 0 !<
1029	INTEGER(iwp) :: terminate_coupled_remote = 0 !<
1030	INTEGER(iwp) :: timestep_count = 0 !<
1031
1032	INTEGER(iwp) :: dist_nxl(0:1) !<
1033	INTEGER(iwp) :: dist_nxr(0:1) !<
1034	INTEGER(iwp) :: dist_nyn(0:1) !<
1035	INTEGER(iwp) :: dist_nys(0:1) !<
1036	INTEGER(iwp) :: do2d_no(0:1) = 0 !<
1037	INTEGER(iwp) :: do2d_xy_time_count(0:1) !<
1038	INTEGER(iwp) :: do2d_xz_time_count(0:1) !<
1039	INTEGER(iwp) :: do2d_yz_time_count(0:1) !<
1040	INTEGER(iwp) :: do3d_no(0:1) = 0 !<
1041	INTEGER(iwp) :: do3d_time_count(0:1) !<
1042	INTEGER(iwp) :: domask_no(max_masks,0:1) = 0 !<
1043	INTEGER(iwp) :: domask_time_count(max_masks,0:1) !<
1044	INTEGER(iwp) :: mask_size(max_masks,3) = -1 !<
1045	INTEGER(iwp) :: mask_size_l(max_masks,3) = -1 !<
1046	INTEGER(iwp) :: mask_start_l(max_masks,3) = -1 !<
1047	INTEGER(iwp) :: pt_vertical_gradient_level_ind(10) = -9999 !<
1048	INTEGER(iwp) :: q_vertical_gradient_level_ind(10) = -9999 !<
1049	INTEGER(iwp) :: s_vertical_gradient_level_ind(10) = -9999 !<
1050	INTEGER(iwp) :: sa_vertical_gradient_level_ind(10) = -9999 !<
1051	INTEGER(iwp) :: section(100,3) !<
1052	INTEGER(iwp) :: section_xy(100) = -9999 !<
1053	INTEGER(iwp) :: section_xz(100) = -9999 !<
1054	INTEGER(iwp) :: section_yz(100) = -9999 !<
1055	INTEGER(iwp) :: ug_vertical_gradient_level_ind(10) = -9999 !<
1056	INTEGER(iwp) :: vg_vertical_gradient_level_ind(10) = -9999 !<
1057	INTEGER(iwp) :: subs_vertical_gradient_level_i(10) = -9999 !<
1058
1059
1060	INTEGER(iwp), DIMENSION(0:1) :: ntdim_2d_xy !<
1061	INTEGER(iwp), DIMENSION(0:1) :: ntdim_2d_xz !<
1062	INTEGER(iwp), DIMENSION(0:1) :: ntdim_2d_yz !<
1063	INTEGER(iwp), DIMENSION(0:1) :: ntdim_3d !<
1064
1065	INTEGER(iwp), DIMENSION(:), ALLOCATABLE :: grid_level_count !<
1066
1067	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: mask_i !<
1068	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: mask_j !<
1069	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: mask_k !<
1070	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: mask_i_global !<
1071	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: mask_j_global !<
1072	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: mask_k_global !<
1073
1074	LOGICAL :: bc_lr_cyc =.TRUE. !<
1075	LOGICAL :: bc_lr_dirrad = .FALSE. !<
1076	LOGICAL :: bc_lr_raddir = .FALSE. !<
1077	LOGICAL :: bc_ns_cyc = .TRUE. !<
1078	LOGICAL :: bc_ns_dirrad = .FALSE. !<
1079	LOGICAL :: bc_ns_raddir = .FALSE. !<
1080	LOGICAL :: call_microphysics_at_all_substeps = .FALSE. !<
1081	LOGICAL :: call_psolver_at_all_substeps = .TRUE. !<
1082	LOGICAL :: cloud_droplets = .FALSE. !<
1083	LOGICAL :: cloud_physics = .FALSE. !<
1084	LOGICAL :: cloud_top_radiation = .FALSE. !<
1085	LOGICAL :: conserve_volume_flow = .FALSE. !<
1086	LOGICAL :: constant_diffusion = .FALSE. !<
1087	LOGICAL :: constant_flux_layer = .TRUE. !<
1088	LOGICAL :: constant_heatflux = .TRUE. !<
1089	LOGICAL :: constant_top_heatflux = .TRUE. !<
1090	LOGICAL :: constant_top_momentumflux = .FALSE. !<
1091	LOGICAL :: constant_top_salinityflux = .TRUE. !<
1092	LOGICAL :: constant_top_scalarflux = .TRUE. !<
1093	LOGICAL :: constant_scalarflux = .TRUE. !<
1094	LOGICAL :: constant_waterflux = .TRUE. !<
1095	LOGICAL :: create_disturbances = .TRUE. !<
1096	LOGICAL :: data_output_2d_on_each_pe = .TRUE. !<
1097	LOGICAL :: dissipation_control = .FALSE. !<
1098	LOGICAL :: disturbance_created = .FALSE. !<
1099	LOGICAL :: do2d_at_begin = .FALSE. !<
1100	LOGICAL :: do3d_at_begin = .FALSE. !<
1101	LOGICAL :: do_sum = .FALSE. !<
1102	LOGICAL :: dp_external = .FALSE. !<
1103	LOGICAL :: dp_smooth = .FALSE. !<
1104	LOGICAL :: dt_fixed = .FALSE. !<
1105	LOGICAL :: dt_3d_reached !<
1106	LOGICAL :: dt_3d_reached_l !<
1107	LOGICAL :: exchange_mg = .FALSE. !<
1108	LOGICAL :: first_call_lpm = .TRUE. !<
1109	LOGICAL :: force_print_header = .FALSE. !<
1110	LOGICAL :: galilei_transformation = .FALSE. !<
1111	LOGICAL :: humidity = .FALSE. !<
1112	LOGICAL :: humidity_remote = .FALSE. !<
1113	LOGICAL :: inflow_l = .FALSE. !<
1114	LOGICAL :: inflow_n = .FALSE. !<
1115	LOGICAL :: inflow_r = .FALSE. !<
1116	LOGICAL :: inflow_s = .FALSE. !<
1117	LOGICAL :: large_scale_forcing = .FALSE. !<
1118	LOGICAL :: large_scale_subsidence = .FALSE. !<
1119	LOGICAL :: lsf_exception = .FALSE. !< temporary flag for use of lsf with buildings on flat terrain
1120	LOGICAL :: lsf_surf = .TRUE. !<
1121	LOGICAL :: lsf_vert = .TRUE. !<
1122	LOGICAL :: lptnudge = .FALSE. !<
1123	LOGICAL :: lqnudge = .FALSE. !<
1124	LOGICAL :: lunudge = .FALSE. !<
1125	LOGICAL :: lvnudge = .FALSE. !<
1126	LOGICAL :: lwnudge = .FALSE. !<
1127	LOGICAL :: masking_method = .FALSE. !<
1128	LOGICAL :: microphysics_sat_adjust = .FALSE. !<
1129	LOGICAL :: microphysics_kessler = .FALSE. !<
1130	LOGICAL :: microphysics_seifert = .FALSE. !<
1131	LOGICAL :: mg_switch_to_pe0 = .FALSE. !<
1132	LOGICAL :: monotonic_adjustment = .FALSE. !<
1133	LOGICAL :: nest_bound_l = .FALSE. !< nested boundary on left side
1134	LOGICAL :: nest_bound_n = .FALSE. !< nested boundary on north side
1135	LOGICAL :: nest_bound_r = .FALSE. !< nested boundary on right side
1136	LOGICAL :: nest_bound_s = .FALSE. !< nested boundary on south side
1137	LOGICAL :: nest_domain = .FALSE. !< domain is nested into a parent domain
1138	LOGICAL :: neutral = .FALSE. !<
1139	LOGICAL :: nudging = .FALSE. !<
1140	LOGICAL :: ocean = .FALSE. !<
1141	LOGICAL :: on_device = .FALSE. !<
1142	LOGICAL :: outflow_l = .FALSE. !<
1143	LOGICAL :: outflow_n = .FALSE. !<
1144	LOGICAL :: outflow_r = .FALSE. !<
1145	LOGICAL :: outflow_s = .FALSE. !<
1146	LOGICAL :: passive_scalar = .FALSE. !<
1147	LOGICAL :: precipitation = .FALSE. !<
1148	LOGICAL :: random_heatflux = .FALSE. !<
1149	LOGICAL :: recycling_yshift = .FALSE. !<
1150	LOGICAL :: run_control_header = .FALSE. !<
1151	LOGICAL :: run_coupled = .TRUE. !<
1152	LOGICAL :: scalar_rayleigh_damping = .TRUE. !<
1153	LOGICAL :: sloping_surface = .FALSE. !<
1154	LOGICAL :: stop_dt = .FALSE. !<
1155	LOGICAL :: synchronous_exchange = .FALSE. !<
1156	LOGICAL :: terminate_run = .FALSE. !<
1157	LOGICAL :: transpose_compute_overlap = .FALSE. !<
1158	LOGICAL :: turbulent_inflow = .FALSE. !<
1159	LOGICAL :: turbulent_outflow = .FALSE. !< flag for turbulent outflow condition
1160	LOGICAL :: urban_surface = .FALSE. !< flag for urban surface model
1161	LOGICAL :: use_cmax = .TRUE. !<
1162	LOGICAL :: use_initial_profile_as_reference = .FALSE. !<
1163	LOGICAL :: use_prescribed_profile_data = .FALSE. !<
1164	LOGICAL :: use_single_reference_value = .FALSE. !<
1165	LOGICAL :: use_subsidence_tendencies = .FALSE. !<
1166	LOGICAL :: use_surface_fluxes = .FALSE. !<
1167	LOGICAL :: use_top_fluxes = .FALSE. !<
1168	LOGICAL :: use_ug_for_galilei_tr = .TRUE. !<
1169	LOGICAL :: use_upstream_for_tke = .FALSE. !<
1170	LOGICAL :: virtual_flight = .FALSE. !< flag for virtual flight model
1171	LOGICAL :: wall_adjustment = .TRUE. !<
1172	LOGICAL :: ws_scheme_sca = .FALSE. !<
1173	LOGICAL :: ws_scheme_mom = .FALSE. !<
1174
1175	LOGICAL :: data_output_xy(0:1) = .FALSE. !<
1176	LOGICAL :: data_output_xz(0:1) = .FALSE. !<
1177	LOGICAL :: data_output_yz(0:1) = .FALSE. !<
1178
1179	REAL(wp) :: advected_distance_x = 0.0_wp !<
1180	REAL(wp) :: advected_distance_y = 0.0_wp !<
1181	REAL(wp) :: alpha_surface = 0.0_wp !<
1182	REAL(wp) :: atmos_ocean_sign = 1.0_wp !<
1183	REAL(wp) :: averaging_interval = 0.0_wp !<
1184	REAL(wp) :: averaging_interval_pr = 9999999.9_wp !<
1185	REAL(wp) :: bc_pt_t_val !<
1186	REAL(wp) :: bc_q_t_val !<
1187	REAL(wp) :: bc_s_t_val !<
1188	REAL(wp) :: bottom_salinityflux = 0.0_wp !<
1189	REAL(wp) :: building_height = 50.0_wp !<
1190	REAL(wp) :: building_length_x = 50.0_wp !<
1191	REAL(wp) :: building_length_y = 50.0_wp !<
1192	REAL(wp) :: building_wall_left = 9999999.9_wp !<
1193	REAL(wp) :: building_wall_south = 9999999.9_wp !<
1194	REAL(wp) :: canyon_height = 50.0_wp !<
1195	REAL(wp) :: canyon_width_x = 9999999.9_wp !<
1196	REAL(wp) :: canyon_width_y = 9999999.9_wp !<
1197	REAL(wp) :: canyon_wall_left = 9999999.9_wp !<
1198	REAL(wp) :: canyon_wall_south = 9999999.9_wp !<
1199	REAL(wp) :: cfl_factor = -1.0_wp !<
1200	REAL(wp) :: cos_alpha_surface !<
1201	REAL(wp) :: coupling_start_time = 0.0_wp !<
1202	REAL(wp) :: disturbance_amplitude = 0.25_wp !<
1203	REAL(wp) :: disturbance_energy_limit = 0.01_wp !<
1204	REAL(wp) :: disturbance_level_b = -9999999.9_wp !<
1205	REAL(wp) :: disturbance_level_t = -9999999.9_wp !<
1206	REAL(wp) :: dp_level_b = 0.0_wp !<
1207	REAL(wp) :: dt = -1.0_wp !<
1208	REAL(wp) :: dt_averaging_input = 0.0_wp !<
1209	REAL(wp) :: dt_averaging_input_pr = 9999999.9_wp !<
1210	REAL(wp) :: dt_coupling = 9999999.9_wp !<
1211	REAL(wp) :: dt_data_output = 9999999.9_wp !<
1212	REAL(wp) :: dt_data_output_av = 9999999.9_wp !<
1213	REAL(wp) :: dt_disturb = 9999999.9_wp !<
1214	REAL(wp) :: dt_dopr = 9999999.9_wp !<
1215	REAL(wp) :: dt_dopr_listing = 9999999.9_wp !<
1216	REAL(wp) :: dt_dopts = 9999999.9_wp !<
1217	REAL(wp) :: dt_dots = 9999999.9_wp !<
1218	REAL(wp) :: dt_do2d_xy = 9999999.9_wp !<
1219	REAL(wp) :: dt_do2d_xz = 9999999.9_wp !<
1220	REAL(wp) :: dt_do2d_yz = 9999999.9_wp !<
1221	REAL(wp) :: dt_do3d = 9999999.9_wp !<
1222	REAL(wp) :: dt_dvrp = 9999999.9_wp !<
1223	REAL(wp) :: dt_max = 20.0_wp !<
1224	REAL(wp) :: dt_restart = 9999999.9_wp !<
1225	REAL(wp) :: dt_run_control = 60.0_wp !<
1226	REAL(wp) :: dt_3d = 1.0_wp !<
1227	REAL(wp) :: dz = -1.0_wp !<
1228	REAL(wp) :: dz_max = 9999999.9_wp !<
1229	REAL(wp) :: dz_stretch_factor = 1.08_wp !<
1230	REAL(wp) :: dz_stretch_level = 100000.0_wp !<
1231	REAL(wp) :: e_init = 0.0_wp !<
1232	REAL(wp) :: e_min = 0.0_wp !<
1233	REAL(wp) :: end_time = 0.0_wp !<
1234	REAL(wp) :: f = 0.0_wp !<
1235	REAL(wp) :: fs = 0.0_wp !<
1236	REAL(wp) :: g = 9.81_wp !<
1237	REAL(wp) :: inflow_damping_height = 9999999.9_wp !<
1238	REAL(wp) :: inflow_damping_width = 9999999.9_wp !<
1239	REAL(wp) :: kappa = 0.4_wp !<
1240	REAL(wp) :: km_constant = -1.0_wp !<
1241	REAL(wp) :: mask_scale_x = 1.0_wp !<
1242	REAL(wp) :: mask_scale_y = 1.0_wp !<
1243	REAL(wp) :: mask_scale_z = 1.0_wp !<
1244	REAL(wp) :: maximum_cpu_time_allowed = 0.0_wp !<
1245	REAL(wp) :: molecular_viscosity = 1.461E-5_wp !<
1246	REAL(wp) :: old_dt = 1.0E-10_wp !<
1247	REAL(wp) :: omega = 7.29212E-5_wp !<
1248	REAL(wp) :: omega_sor = 1.8_wp !<
1249	REAL(wp) :: outflow_source_plane = -9999999.9_wp !< x-position of outflow-source plane (turbulent outflow method)
1250	REAL(wp) :: particle_maximum_age = 9999999.9_wp !<
1251	REAL(wp) :: phi = 55.0_wp !<
1252	REAL(wp) :: prandtl_number = 1.0_wp !<
1253	REAL(wp) :: precipitation_amount_interval = 9999999.9_wp !<
1254	REAL(wp) :: prho_reference !<
1255	REAL(wp) :: pt_damping_factor = 0.0_wp !<
1256	REAL(wp) :: pt_damping_width = 0.0_wp !<
1257	REAL(wp) :: pt_reference = 9999999.9_wp !<
1258	REAL(wp) :: pt_slope_offset = 0.0_wp !<
1259	REAL(wp) :: pt_surface = 300.0_wp !<
1260	REAL(wp) :: pt_surface_initial_change = 0.0_wp !<
1261	REAL(wp) :: q_surface = 0.0_wp !<
1262	REAL(wp) :: q_surface_initial_change = 0.0_wp !<
1263	REAL(wp) :: rayleigh_damping_factor = -1.0_wp !<
1264	REAL(wp) :: rayleigh_damping_height = -1.0_wp !<
1265	REAL(wp) :: recycling_width = 9999999.9_wp !<
1266	REAL(wp) :: residual_limit = 1.0E-4_wp !<
1267	REAL(wp) :: restart_time = 9999999.9_wp !<
1268	REAL(wp) :: rho_reference !<
1269	REAL(wp) :: rho_surface !<
1270	REAL(wp) :: roughness_length = 0.1_wp !<
1271	REAL(wp) :: sa_surface = 35.0_wp !<
1272	REAL(wp) :: simulated_time = 0.0_wp !<
1273	REAL(wp) :: simulated_time_at_begin !<
1274	REAL(wp) :: sin_alpha_surface !<
1275	REAL(wp) :: skip_time_data_output = 0.0_wp !<
1276	REAL(wp) :: skip_time_data_output_av = 9999999.9_wp !<
1277	REAL(wp) :: skip_time_dopr = 9999999.9_wp !<
1278	REAL(wp) :: skip_time_do2d_xy = 9999999.9_wp !<
1279	REAL(wp) :: skip_time_do2d_xz = 9999999.9_wp !<
1280	REAL(wp) :: skip_time_do2d_yz = 9999999.9_wp !<
1281	REAL(wp) :: skip_time_do3d = 9999999.9_wp !<
1282	REAL(wp) :: surface_heatflux = 9999999.9_wp !<
1283	REAL(wp) :: surface_pressure = 1013.25_wp !<
1284	REAL(wp) :: surface_scalarflux = 9999999.9_wp !<
1285	REAL(wp) :: surface_waterflux = 9999999.9_wp !<
1286	REAL(wp) :: s_surface = 0.0_wp !<
1287	REAL(wp) :: s_surface_initial_change = 0.0_wp !<
1288	REAL(wp) :: termination_time_needed = -1.0_wp !<
1289	REAL(wp) :: time_coupling = 0.0_wp !<
1290	REAL(wp) :: time_disturb = 0.0_wp !<
1291	REAL(wp) :: time_dopr = 0.0_wp !<
1292	REAL(wp) :: time_dopr_av = 0.0_wp !<
1293	REAL(wp) :: time_dopr_listing = 0.0_wp !<
1294	REAL(wp) :: time_dopts = 0.0_wp !<
1295	REAL(wp) :: time_dosp = 0.0_wp !<
1296	REAL(wp) :: time_dosp_av = 0.0_wp !<
1297	REAL(wp) :: time_dots = 0.0_wp !<
1298	REAL(wp) :: time_do2d_xy = 0.0_wp !<
1299	REAL(wp) :: time_do2d_xz = 0.0_wp !<
1300	REAL(wp) :: time_do2d_yz = 0.0_wp !<
1301	REAL(wp) :: time_do3d = 0.0_wp !<
1302	REAL(wp) :: time_do_av = 0.0_wp !<
1303	REAL(wp) :: time_do_sla = 0.0_wp !<
1304	REAL(wp) :: time_dvrp = 0.0_wp !<
1305	REAL(wp) :: time_restart = 9999999.9_wp !<
1306	REAL(wp) :: time_run_control = 0.0_wp !<
1307	REAL(wp) :: time_since_reference_point !<
1308	REAL(wp) :: top_heatflux = 9999999.9_wp !<
1309	REAL(wp) :: top_momentumflux_u = 9999999.9_wp !<
1310	REAL(wp) :: top_momentumflux_v = 9999999.9_wp !<
1311	REAL(wp) :: top_salinityflux = 9999999.9_wp !<
1312	REAL(wp) :: top_scalarflux = 9999999.9_wp !<
1313	REAL(wp) :: ug_surface = 0.0_wp !<
1314	REAL(wp) :: u_bulk = 0.0_wp !<
1315	REAL(wp) :: u_gtrans = 0.0_wp !<
1316	REAL(wp) :: vg_surface = 0.0_wp !<
1317	REAL(wp) :: vpt_reference = 9999999.9_wp !<
1318	REAL(wp) :: v_bulk = 0.0_wp !<
1319	REAL(wp) :: v_gtrans = 0.0_wp !<
1320	REAL(wp) :: wall_adjustment_factor = 1.8_wp !<
1321	REAL(wp) :: zeta_max = 20.0_wp !< Maximum value of zeta = z/L
1322	REAL(wp) :: zeta_min = -9990.0_wp !< Minimum value of zeta = z/L
1323	REAL(wp) :: z_max_do2d = -1.0_wp !<
1324	REAL(wp) :: z0h_factor = 1.0_wp !<
1325
1326	REAL(wp) :: do2d_xy_last_time(0:1) = -1.0_wp !<
1327	REAL(wp) :: do2d_xz_last_time(0:1) = -1.0_wp !<
1328	REAL(wp) :: do2d_yz_last_time(0:1) = -1.0_wp !<
1329	REAL(wp) :: dpdxy(1:2) = 0.0_wp !<
1330	REAL(wp) :: dt_domask(max_masks) = 9999999.9_wp !<
1331	REAL(wp) :: mask_scale(3) !<
1332	REAL(wp) :: pt_vertical_gradient(10) = 0.0_wp !<
1333	REAL(wp) :: pt_vertical_gradient_level(10) = -1.0_wp !<
1334	REAL(wp) :: q_vertical_gradient(10) = 0.0_wp !<
1335	REAL(wp) :: q_vertical_gradient_level(10) = -1.0_wp !<
1336	REAL(wp) :: s_vertical_gradient(10) = 0.0_wp !<
1337	REAL(wp) :: s_vertical_gradient_level(10) = -1.0_wp !<
1338	REAL(wp) :: sa_vertical_gradient(10) = 0.0_wp !<
1339	REAL(wp) :: sa_vertical_gradient_level(10) = -1.0_wp !<
1340	REAL(wp) :: skip_time_domask(max_masks) = 9999999.9_wp !<
1341	REAL(wp) :: threshold(20) = 0.0_wp !<
1342	REAL(wp) :: time_domask(max_masks) = 0.0_wp !<
1343	REAL(wp) :: tsc(10) = (/ 1.0_wp, 1.0_wp, 0.0_wp, 0.0_wp, & !<
1344	0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp /)
1345	REAL(wp) :: u_profile(100) = 9999999.9_wp !<
1346	REAL(wp) :: uv_heights(100) = 9999999.9_wp !<
1347	REAL(wp) :: v_profile(100) = 9999999.9_wp !<
1348	REAL(wp) :: ug_vertical_gradient(10) = 0.0_wp !<
1349	REAL(wp) :: ug_vertical_gradient_level(10) = -9999999.9_wp !<
1350	REAL(wp) :: vg_vertical_gradient(10) = 0.0_wp !<
1351	REAL(wp) :: vg_vertical_gradient_level(10) = -9999999.9_wp !<
1352	REAL(wp) :: volume_flow(1:3) = 0.0_wp !<
1353	REAL(wp) :: volume_flow_area(1:3) = 0.0_wp !<
1354	REAL(wp) :: volume_flow_initial(1:3) = 0.0_wp !<
1355	REAL(wp) :: wall_heatflux(0:4) = 0.0_wp !<
1356	REAL(wp) :: wall_humidityflux(0:4) = 0.0_wp !<
1357	REAL(wp) :: wall_nrflux(0:4) = 0.0_wp !<
1358	REAL(wp) :: wall_qflux(0:4) = 0.0_wp !<
1359	REAL(wp) :: wall_qrflux(0:4) = 0.0_wp !<
1360	REAL(wp) :: wall_salinityflux(0:4) = 0.0_wp !<
1361	REAL(wp) :: wall_sflux(0:4) = 0.0_wp !<
1362	REAL(wp) :: wall_scalarflux(0:4) = 0.0_wp !<
1363	REAL(wp) :: subs_vertical_gradient(10) = 0.0_wp !<
1364	REAL(wp) :: subs_vertical_gradient_level(10) = -9999999.9_wp !<
1365
1366	REAL(wp), DIMENSION(:), ALLOCATABLE :: dp_smooth_factor !<
1367
1368	REAL(wp), DIMENSION(max_masks,mask_xyz_dimension) :: mask_x = -1.0_wp !<
1369	REAL(wp), DIMENSION(max_masks,mask_xyz_dimension) :: mask_y = -1.0_wp !<
1370	REAL(wp), DIMENSION(max_masks,mask_xyz_dimension) :: mask_z = -1.0_wp !<
1371
1372	REAL(wp), DIMENSION(max_masks,3) :: mask_x_loop = -1.0_wp !<
1373	REAL(wp), DIMENSION(max_masks,3) :: mask_y_loop = -1.0_wp !<
1374	REAL(wp), DIMENSION(max_masks,3) :: mask_z_loop = -1.0_wp !<
1375
1376	!
1377	!-- internal mask arrays ("mask,dimension,selection")
1378	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: mask !<
1379	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: mask_loop !<
1380
1381	SAVE
1382
1383	END MODULE control_parameters
1384
1385
1386	!------------------------------------------------------------------------------!
1387	! Description:
1388	! ------------
1389	!> Definition of variables used with dvrp-software.
1390	!------------------------------------------------------------------------------!
1391	MODULE dvrp_variables
1392
1393	USE kinds
1394
1395	CHARACTER (LEN=10) :: dvrp_output = 'rtsp' !<
1396	CHARACTER (LEN=10) :: particle_color = 'none' !<
1397	CHARACTER (LEN=10) :: particle_dvrpsize = 'none' !<
1398
1399	CHARACTER (LEN=20), DIMENSION(10) :: mode_dvrp = & !<
1400	(/ ( ' ', i9 = 1,10 ) /)
1401
1402	CHARACTER (LEN=80) :: dvrp_directory = 'default' !<
1403	CHARACTER (LEN=80) :: dvrp_file = 'default' !<
1404	CHARACTER (LEN=80) :: dvrp_host = 'origin.rvs.uni-hannover.de' !<
1405	CHARACTER (LEN=80) :: dvrp_password = '********' !<
1406	CHARACTER (LEN=80) :: dvrp_username = ' ' !<
1407
1408	INTEGER(iwp) :: cluster_size = 1 !<
1409	INTEGER(iwp) :: dvrp_colortable_entries = 4 !<
1410	INTEGER(iwp) :: dvrp_colortable_entries_prt = 22 !<
1411	INTEGER(iwp) :: islice_dvrp !<
1412	INTEGER(iwp) :: nx_dvrp !<
1413	INTEGER(iwp) :: nxl_dvrp !<
1414	INTEGER(iwp) :: nxr_dvrp !<
1415	INTEGER(iwp) :: ny_dvrp !<
1416	INTEGER(iwp) :: nyn_dvrp !<
1417	INTEGER(iwp) :: nys_dvrp !<
1418	INTEGER(iwp) :: nz_dvrp, pathlines_fadeintime = 5 !<
1419	INTEGER(iwp) :: pathlines_fadeouttime = 5 !<
1420	INTEGER(iwp) :: pathlines_linecount = 1000 !<
1421	INTEGER(iwp) :: pathlines_maxhistory = 40 !<
1422	INTEGER(iwp) :: pathlines_wavecount = 10 !<
1423	INTEGER(iwp) :: pathlines_wavetime = 50 !<
1424	INTEGER(iwp) :: vc_gradient_normals = 0 !<
1425	INTEGER(iwp) :: vc_mode = 0 !<
1426	INTEGER(iwp) :: vc_size_x = 2 !<
1427	INTEGER(iwp) :: vc_size_y = 2 !<
1428	INTEGER(iwp) :: vc_size_z = 2 !<
1429
1430	INTEGER(iwp), DIMENSION(10) :: slicer_position_dvrp !<
1431
1432	LOGICAL :: cyclic_dvrp = .FALSE. !<
1433	LOGICAL :: dvrp_overlap !<
1434	LOGICAL :: dvrp_total_overlap !<
1435	LOGICAL :: local_dvrserver_running !<
1436	LOGICAL :: lock_steering_update = .FALSE. !<
1437	LOGICAL :: use_seperate_pe_for_dvrp_output = .FALSE. !<
1438
1439	REAL(wp) :: clip_dvrp_l = 9999999.9_wp !<
1440	REAL(wp) :: clip_dvrp_n = 9999999.9_wp !<
1441	REAL(wp) :: clip_dvrp_r = 9999999.9_wp !<
1442	REAL(wp) :: clip_dvrp_s = 9999999.9_wp !<
1443	REAL(wp) :: superelevation = 1.0_wp !<
1444	REAL(wp) :: superelevation_x = 1.0_wp !<
1445	REAL(wp) :: superelevation_y = 1.0_wp !<
1446	REAL(wp) :: vc_alpha = 38.0_wp !<
1447
1448	REAL(wp), DIMENSION(2) :: color_interval = (/ 0.0_wp, 1.0_wp /) !<
1449	REAL(wp), DIMENSION(2) :: dvrpsize_interval = (/ 0.0_wp, 1.0_wp /) !<
1450
1451	REAL(wp), DIMENSION(3) :: groundplate_color = (/ 0.0_wp, 0.6_wp, 0.0_wp /) !<
1452	REAL(wp), DIMENSION(3) :: topography_color = (/ 0.8_wp, 0.7_wp, 0.6_wp /) !<
1453
1454	REAL(wp), DIMENSION(2,10) :: slicer_range_limits_dvrp !<
1455
1456	REAL(wp), DIMENSION(3,10) :: isosurface_color !<
1457
1458	REAL(sp), DIMENSION(2,100) :: interval_values_dvrp !<
1459	REAL(sp), DIMENSION(2,100) :: interval_values_dvrp_prt !<
1460	REAL(sp), DIMENSION(2,100) :: interval_h_dvrp !<
1461	REAL(sp), DIMENSION(2,100) :: interval_h_dvrp_prt !<
1462	REAL(sp), DIMENSION(2,100) :: interval_l_dvrp = 0.5_sp !<
1463	REAL(sp), DIMENSION(2,100) :: interval_l_dvrp_prt = 0.5_sp !<
1464	REAL(sp), DIMENSION(2,100) :: interval_s_dvrp = 1.0_sp !<
1465	REAL(sp), DIMENSION(2,100) :: interval_s_dvrp_prt = 1.0_sp !<
1466	REAL(sp), DIMENSION(2,100) :: interval_a_dvrp = 0.0_sp !<
1467	REAL(sp), DIMENSION(2,100) :: interval_a_dvrp_prt = 0.0_sp !<
1468
1469	DATA slicer_range_limits_dvrp / -1.0_wp, 1.0_wp, -1.0_wp, 1.0_wp, -1.0_wp, 1.0_wp, & !<
1470	-1.0_wp, 1.0_wp, -1.0_wp, 1.0_wp, -1.0_wp, 1.0_wp, &
1471	-1.0_wp, 1.0_wp, -1.0_wp, 1.0_wp, -1.0_wp, 1.0_wp, &
1472	-1.0_wp, 1.0_wp /
1473
1474	DATA isosurface_color / 0.9_wp, 0.9_wp, 0.9_wp, 0.8_wp, 0.1_wp, 0.1_wp, 0.1_wp, 0.1_wp, 0.8_wp, & !<
1475	0.1_wp, 0.8_wp, 0.1_wp, 0.6_wp, 0.1_wp, 0.1_wp, 0.1_wp, 0.1_wp, 0.6_wp, &
1476	0.1_wp, 0.6_wp, 0.1_wp, 0.4_wp, 0.1_wp, 0.1_wp, 0.1_wp, 0.1_wp, 0.4_wp, &
1477	0.1_wp, 0.4_wp, 0.1_wp /
1478
1479	DATA interval_h_dvrp / 270.0_wp, 225.0_wp, 225.0_wp, 180.0_wp, 70.0_wp, 25.0_wp, & !<
1480	25.0_wp, -25.0_wp, 192 * 0.0_wp /
1481
1482	DATA interval_h_dvrp_prt / 270.0_wp, 225.0_wp, 225.0_wp, 180.0_wp, 70.0_wp, 25.0_wp, & !<
1483	25.0_wp, -25.0_wp, 192 * 0.0_wp /
1484
1485	REAL(sp), DIMENSION(:), ALLOCATABLE :: xcoor_dvrp !<
1486	REAL(sp), DIMENSION(:), ALLOCATABLE :: ycoor_dvrp !<
1487	REAL(sp), DIMENSION(:), ALLOCATABLE :: zcoor_dvrp !<
1488
1489	TYPE steering
1490	CHARACTER (LEN=24) :: name !<
1491	REAL(sp) :: min !<
1492	REAL(sp) :: max !<
1493	INTEGER(iwp) :: imin !<
1494	INTEGER(iwp) :: imax !<
1495	END TYPE steering
1496
1497	TYPE(steering), DIMENSION(:), ALLOCATABLE :: steering_dvrp !<
1498
1499	SAVE
1500
1501	END MODULE dvrp_variables
1502
1503
1504	!------------------------------------------------------------------------------!
1505	! Description:
1506	! ------------
1507	!> Definition of grid spacings.
1508	!------------------------------------------------------------------------------!
1509	MODULE grid_variables
1510
1511	USE kinds
1512
1513	REAL(wp) :: ddx !<
1514	REAL(wp) :: ddx2 !<
1515	REAL(wp) :: dx = 1.0_wp !<
1516	REAL(wp) :: dx2 !<
1517	REAL(wp) :: ddy !<
1518	REAL(wp) :: ddy2 !<
1519	REAL(wp) :: dy = 1.0_wp !<
1520	REAL(wp) :: dy2 !<
1521
1522	REAL(wp), DIMENSION(:), ALLOCATABLE :: ddx2_mg !<
1523	REAL(wp), DIMENSION(:), ALLOCATABLE :: ddy2_mg !<
1524
1525	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: fwxm !<
1526	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: fwxp !<
1527	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: fwym !<
1528	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: fwyp !<
1529	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: fxm !<
1530	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: fxp !<
1531	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: fym !<
1532	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: fyp !<
1533	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: wall_e_x !<
1534	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: wall_e_y !<
1535	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: wall_u !<
1536	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: wall_v !<
1537	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: wall_w_x !<
1538	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: wall_w_y !<
1539	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: zu_s_inner !<
1540	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: zw_w_inner !<
1541
1542	SAVE
1543
1544	END MODULE grid_variables
1545
1546
1547	!------------------------------------------------------------------------------!
1548	! Description:
1549	! ------------
1550	!> Definition of array bounds, number of gridpoints, and wall flag arrays.
1551	!------------------------------------------------------------------------------!
1552	MODULE indices
1553
1554	USE kinds
1555
1556	INTEGER(iwp) :: i_left !<
1557	INTEGER(iwp) :: i_right !<
1558	INTEGER(iwp) :: j_north !<
1559	INTEGER(iwp) :: j_south !<
1560	INTEGER(iwp) :: nbgp = 3 !<
1561	INTEGER(iwp) :: ngp_sums !<
1562	INTEGER(iwp) :: ngp_sums_ls !<
1563	INTEGER(iwp) :: nnx !<
1564	INTEGER(iwp) :: nx = 0 !<
1565	INTEGER(iwp) :: nx_a !<
1566	INTEGER(iwp) :: nx_o !<
1567	INTEGER(iwp) :: nxl !<
1568	INTEGER(iwp) :: nxlg !<
1569	INTEGER(iwp) :: nxlu !<
1570	INTEGER(iwp) :: nxr !<
1571	INTEGER(iwp) :: nxrg !<
1572	INTEGER(iwp) :: nx_on_file !<
1573	INTEGER(iwp) :: nny !<
1574	INTEGER(iwp) :: ny = 0 !<
1575	INTEGER(iwp) :: ny_a !<
1576	INTEGER(iwp) :: ny_o !<
1577	INTEGER(iwp) :: nyn !<
1578	INTEGER(iwp) :: nyng !<
1579	INTEGER(iwp) :: nys !<
1580	INTEGER(iwp) :: nysg !<
1581	INTEGER(iwp) :: nysv !<
1582	INTEGER(iwp) :: ny_on_file !<
1583	INTEGER(iwp) :: nnz !<
1584	INTEGER(iwp) :: nz = 0 !<
1585	INTEGER(iwp) :: nzb !<
1586	INTEGER(iwp) :: nzb_diff !<
1587	INTEGER(iwp) :: nzb_max !<
1588	INTEGER(iwp) :: nzt !<
1589	INTEGER(iwp) :: nzt_diff !<
1590
1591	INTEGER(idp), DIMENSION(:), ALLOCATABLE :: ngp_3d !<
1592	INTEGER(idp), DIMENSION(:), ALLOCATABLE :: ngp_3d_inner !< ! need to have 64 bit for grids > 2E9
1593
1594	INTEGER(iwp), DIMENSION(:), ALLOCATABLE :: ngp_2dh !<
1595	INTEGER(iwp), DIMENSION(:), ALLOCATABLE :: nxl_mg !<
1596	INTEGER(iwp), DIMENSION(:), ALLOCATABLE :: nxr_mg !<
1597	INTEGER(iwp), DIMENSION(:), ALLOCATABLE :: nyn_mg !<
1598	INTEGER(iwp), DIMENSION(:), ALLOCATABLE :: nys_mg !<
1599	INTEGER(iwp), DIMENSION(:), ALLOCATABLE :: nzt_mg !<
1600
1601
1602	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: ngp_2dh_outer !<
1603	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: ngp_2dh_s_inner !<
1604	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: mg_loc_ind !<
1605	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: nzb_diff_s_inner !<
1606	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: nzb_diff_s_outer !<
1607	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: nzb_diff_u !<
1608	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: nzb_diff_v !<
1609	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: nzb_inner !<
1610	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: nzb_outer !<
1611	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: nzb_s_inner !<
1612	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: nzb_s_outer !<
1613	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: nzb_u_inner !<
1614	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: nzb_u_outer !<
1615	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: nzb_v_inner !<
1616	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: nzb_v_outer !<
1617	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: nzb_w_inner !<
1618	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: nzb_w_outer !<
1619
1620	INTEGER(iwp), DIMENSION(:,:,:), POINTER :: flags !<
1621
1622	INTEGER(iwp), DIMENSION(:,:,:), ALLOCATABLE :: wall_flags_0 !<
1623	INTEGER(iwp), DIMENSION(:,:,:), ALLOCATABLE :: wall_flags_00 !<
1624
1625	INTEGER(iwp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: wall_flags_1 !<
1626	INTEGER(iwp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: wall_flags_2 !<
1627	INTEGER(iwp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: wall_flags_3 !<
1628	INTEGER(iwp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: wall_flags_4 !<
1629	INTEGER(iwp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: wall_flags_5 !<
1630	INTEGER(iwp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: wall_flags_6 !<
1631	INTEGER(iwp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: wall_flags_7 !<
1632	INTEGER(iwp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: wall_flags_8 !<
1633	INTEGER(iwp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: wall_flags_9 !<
1634	INTEGER(iwp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: wall_flags_10 !<
1635
1636	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: rflags_s_inner !<
1637	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: rflags_invers !<
1638
1639	SAVE
1640
1641	END MODULE indices
1642
1643
1644	!------------------------------------------------------------------------------!
1645	! Description:
1646	! ------------
1647	!> Interfaces for special subroutines which use optional parameters.
1648	!------------------------------------------------------------------------------!
1649	MODULE interfaces
1650
1651	INTERFACE
1652
1653	!------------------------------------------------------------------------------!
1654	! Description:
1655	! ------------
1656	!> @todo Missing subroutine description.
1657	!------------------------------------------------------------------------------!
1658	SUBROUTINE global_min_max ( i1, i2, j1, j2, k1, k2, feld, mode, offset, &
1659	wert, wert_ijk, wert1, wert1_ijk )
1660
1661	USE kinds
1662
1663	CHARACTER (LEN=*), INTENT(IN) :: mode !<
1664	INTEGER(iwp), INTENT(IN) :: i1 !<
1665	INTEGER(iwp), INTENT(IN) :: i2 !<
1666	INTEGER(iwp), INTENT(IN) :: j1 !<
1667	INTEGER(iwp), INTENT(IN) :: j2 !<
1668	INTEGER(iwp), INTENT(IN) :: k1 !<
1669	INTEGER(iwp), INTENT(IN) :: k2 !<
1670	INTEGER(iwp) :: wert_ijk(3) !<
1671	INTEGER(iwp), OPTIONAL :: wert1_ijk(3) !<
1672	REAL(wp) :: offset !<
1673	REAL(wp) :: wert !<
1674	REAL(wp), OPTIONAL :: wert1 !<
1675	REAL(wp), INTENT(IN) :: feld(i1:i2,j1:j2,k1:k2) !<
1676
1677	END SUBROUTINE global_min_max
1678
1679	END INTERFACE
1680
1681	SAVE
1682
1683	END MODULE interfaces
1684
1685
1686	!------------------------------------------------------------------------------!
1687	! Description:
1688	! ------------
1689	!> Interfaces for subroutines with pointer arguments called in
1690	!> prognostic_equations.
1691	!------------------------------------------------------------------------------!
1692	MODULE pointer_interfaces
1693
1694	INTERFACE
1695
1696	!------------------------------------------------------------------------------!
1697	! Description:
1698	! ------------
1699	!> @todo Missing subroutine description.
1700	!------------------------------------------------------------------------------!
1701	SUBROUTINE advec_s_bc( sk, sk_char )
1702
1703	USE kinds
1704
1705	CHARACTER (LEN=*), INTENT(IN) :: sk_char !<
1706	#if defined( __nopointer )
1707	REAL(wp), DIMENSION(:,:,:) :: sk !<
1708	#else
1709	REAL(wp), DIMENSION(:,:,:), POINTER :: sk !<
1710	#endif
1711	END SUBROUTINE advec_s_bc
1712
1713	END INTERFACE
1714
1715	SAVE
1716
1717	END MODULE pointer_interfaces
1718
1719
1720	!------------------------------------------------------------------------------!
1721	! Description:
1722	! ------------
1723	!> Definition of variables for the 1D-model.
1724	!------------------------------------------------------------------------------!
1725	MODULE model_1d
1726
1727	USE kinds
1728
1729	INTEGER(iwp) :: current_timestep_number_1d = 0 !<
1730	INTEGER(iwp) :: damp_level_ind_1d !<
1731
1732	LOGICAL :: run_control_header_1d = .FALSE. !<
1733	LOGICAL :: stop_dt_1d = .FALSE. !<
1734
1735	REAL(wp) :: damp_level_1d = -1.0_wp !<
1736	REAL(wp) :: dt_1d = 60.0_wp !<
1737	REAL(wp) :: dt_max_1d = 300.0_wp !<
1738	REAL(wp) :: dt_pr_1d = 9999999.9_wp !<
1739	REAL(wp) :: dt_run_control_1d = 60.0_wp !<
1740	REAL(wp) :: end_time_1d = 864000.0_wp !<
1741	REAL(wp) :: old_dt_1d = 1.0E-10_wp !<
1742	REAL(wp) :: qs1d !<
1743	REAL(wp) :: simulated_time_1d = 0.0_wp !<
1744	REAL(wp) :: time_pr_1d = 0.0_wp !<
1745	REAL(wp) :: time_run_control_1d = 0.0_wp !<
1746	REAL(wp) :: ts1d !<
1747	REAL(wp) :: us1d !<
1748	REAL(wp) :: usws1d !<
1749	REAL(wp) :: vsws1d !<
1750	REAL(wp) :: z01d !<
1751	REAL(wp) :: z0h1d !<
1752
1753
1754	REAL(wp), DIMENSION(:), ALLOCATABLE :: e1d !<
1755	REAL(wp), DIMENSION(:), ALLOCATABLE :: e1d_p !<
1756	REAL(wp), DIMENSION(:), ALLOCATABLE :: kh1d !<
1757	REAL(wp), DIMENSION(:), ALLOCATABLE :: km1d !<
1758	REAL(wp), DIMENSION(:), ALLOCATABLE :: l_black !<
1759	REAL(wp), DIMENSION(:), ALLOCATABLE :: l1d !<
1760	REAL(wp), DIMENSION(:), ALLOCATABLE :: rif1d !<
1761	REAL(wp), DIMENSION(:), ALLOCATABLE :: te_e !<
1762	REAL(wp), DIMENSION(:), ALLOCATABLE :: te_em !<
1763	REAL(wp), DIMENSION(:), ALLOCATABLE :: te_u !<
1764	REAL(wp), DIMENSION(:), ALLOCATABLE :: te_um !<
1765	REAL(wp), DIMENSION(:), ALLOCATABLE :: te_v !<
1766	REAL(wp), DIMENSION(:), ALLOCATABLE :: te_vm !<
1767	REAL(wp), DIMENSION(:), ALLOCATABLE :: u1d !<
1768	REAL(wp), DIMENSION(:), ALLOCATABLE :: u1d_p !<
1769	REAL(wp), DIMENSION(:), ALLOCATABLE :: v1d !<
1770	REAL(wp), DIMENSION(:), ALLOCATABLE :: v1d_p !<
1771
1772	SAVE
1773
1774	END MODULE model_1d
1775
1776
1777
1778
1779	!------------------------------------------------------------------------------!
1780	! Description:
1781	! ------------
1782	!> Definition of variables which define processor topology and the exchange of
1783	!> ghost point layers. This modules must be placed in all routines which contain
1784	!> MPI-calls.
1785	!------------------------------------------------------------------------------!
1786	MODULE pegrid
1787
1788	USE kinds
1789
1790	#if defined( __parallel )
1791	#if defined( __mpifh )
1792	INCLUDE "mpif.h"
1793	#else
1794	USE MPI
1795	#endif
1796	#endif
1797	CHARACTER(LEN=2) :: send_receive = 'al'
1798	CHARACTER(LEN=7) :: myid_char = ''
1799
1800	INTEGER(iwp) :: acc_rank !<
1801	INTEGER(iwp) :: comm1dx !<
1802	INTEGER(iwp) :: comm1dy !<
1803	INTEGER(iwp) :: comm2d !<
1804	INTEGER(iwp) :: comm_inter !<
1805	INTEGER(iwp) :: comm_palm !<
1806	INTEGER(iwp) :: id_inflow = 0 !<
1807	INTEGER(iwp) :: id_outflow = 0 !< myidx of procs at outflow (turbulent outflow method)
1808	INTEGER(iwp) :: id_outflow_source = 0 !< myidx of procs including ouflow source plane (turbulent outflow method)
1809	INTEGER(iwp) :: id_recycling = 0 !<
1810	INTEGER(iwp) :: ierr !<
1811	INTEGER(iwp) :: myid = 0 !<
1812	INTEGER(iwp) :: myidx = 0 !<
1813	INTEGER(iwp) :: myidy = 0 !<
1814	INTEGER(iwp) :: ndim = 2 !<
1815	INTEGER(iwp) :: ngp_a !<
1816	INTEGER(iwp) :: ngp_o !<
1817	INTEGER(iwp) :: ngp_xy !<
1818	INTEGER(iwp) :: ngp_y !<
1819	INTEGER(iwp) :: npex = -1 !<
1820	INTEGER(iwp) :: npey = -1 !<
1821	INTEGER(iwp) :: numprocs = 1 !<
1822	INTEGER(iwp) :: numprocs_previous_run = -1 !<
1823	INTEGER(iwp) :: num_acc_per_node = 0 !<
1824	INTEGER(iwp) :: pleft !<
1825	INTEGER(iwp) :: pnorth !<
1826	INTEGER(iwp) :: pright !<
1827	INTEGER(iwp) :: psouth !<
1828	INTEGER(iwp) :: req_count = 0 !<
1829	INTEGER(iwp) :: sendrecvcount_xy !<
1830	INTEGER(iwp) :: sendrecvcount_yz !<
1831	INTEGER(iwp) :: sendrecvcount_zx !<
1832	INTEGER(iwp) :: sendrecvcount_zyd !<
1833	INTEGER(iwp) :: sendrecvcount_yxd !<
1834	INTEGER(iwp) :: target_id !<
1835	INTEGER(iwp) :: tasks_per_node = -9999 !<
1836	INTEGER(iwp) :: threads_per_task = 1 !<
1837	INTEGER(iwp) :: type_x !<
1838	INTEGER(iwp) :: type_xy !<
1839	INTEGER(iwp) :: type_y !<
1840
1841	INTEGER(iwp) :: pdims(2) = 1 !<
1842	INTEGER(iwp) :: req(100) !<
1843
1844	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: hor_index_bounds !<
1845	INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: hor_index_bounds_previous_run !<
1846
1847	LOGICAL :: background_communication =.FALSE. !<
1848	LOGICAL :: collective_wait = .FALSE. !<
1849	LOGICAL :: sendrecv_in_background = .FALSE. !<
1850
1851	#if defined( __parallel )
1852	#if defined( __mpi2 )
1853	CHARACTER (LEN=MPI_MAX_PORT_NAME) :: port_name !<
1854	#endif
1855
1856	INTEGER(iwp) :: ibuf(12) !<
1857	INTEGER(iwp) :: pcoord(2) !<
1858	INTEGER(iwp) :: status(MPI_STATUS_SIZE) !<
1859
1860	INTEGER(iwp), DIMENSION(MPI_STATUS_SIZE,100) :: wait_stat !<
1861
1862	INTEGER(iwp) :: ngp_yz_int !<
1863	INTEGER(iwp) :: type_xz_int !<
1864	INTEGER(iwp) :: type_yz_int !<
1865
1866	INTEGER(iwp), DIMENSION(:), ALLOCATABLE :: ngp_xz !<
1867	INTEGER(iwp), DIMENSION(:), ALLOCATABLE :: ngp_yz !<
1868	INTEGER(iwp), DIMENSION(:), ALLOCATABLE :: type_x_int !<
1869	INTEGER(iwp), DIMENSION(:), ALLOCATABLE :: type_xz !<
1870	INTEGER(iwp), DIMENSION(:), ALLOCATABLE :: type_y_int !<
1871	INTEGER(iwp), DIMENSION(:), ALLOCATABLE :: type_yz !<
1872
1873	LOGICAL :: left_border_pe = .FALSE. !<
1874	LOGICAL :: north_border_pe = .FALSE. !<
1875	LOGICAL :: reorder = .TRUE. !<
1876	LOGICAL :: right_border_pe = .FALSE. !<
1877	LOGICAL :: south_border_pe = .FALSE. !<
1878
1879	LOGICAL, DIMENSION(2) :: cyclic = (/ .TRUE. , .TRUE. /) !<
1880	LOGICAL, DIMENSION(2) :: remain_dims !<
1881	#endif
1882
1883	SAVE
1884
1885	END MODULE pegrid
1886
1887
1888	!------------------------------------------------------------------------------!
1889	! Description:
1890	! ------------
1891	!> Definition of variables which control PROFIL-output.
1892	!------------------------------------------------------------------------------!
1893	MODULE profil_parameter
1894
1895	USE kinds
1896
1897	INTEGER(iwp), PARAMETER :: crmax = 100 !<
1898
1899	CHARACTER (LEN=20), DIMENSION(20) :: cross_ts_profiles = & !<
1900	(/ ' E E* ', ' dt ', &
1901	' u* w* ', ' th* ', &
1902	' umax vmax wmax ', ' div_old div_new ', &
1903	' z_i_wpt z_i_pt ', ' w"pt"0 w"pt" wpt ', &
1904	' pt(0) pt(zp) ', ' splux spluy spluz ', &
1905	' L ', &
1906	( ' ', i9 = 1, 9 ) /)
1907
1908	CHARACTER (LEN=100), DIMENSION(crmax) :: cross_profiles = & !<
1909	(/ ' u v ', &
1910	' pt ', &
1911	' w"pt" wpt wptBC wpt wptBC ', &
1912	' w"u" wu wu w"v" wv wv ', &
1913	' km kh ', &
1914	' l ', &
1915	( ' ', i9 = 1, 94 ) /)
1916
1917	INTEGER(iwp) :: profile_columns = 2 !<
1918	INTEGER(iwp) :: profile_rows = 3 !<
1919	INTEGER(iwp) :: profile_number = 0 !<
1920
1921	INTEGER(iwp) :: cross_ts_numbers(crmax,crmax) = 0 !<
1922	INTEGER(iwp) :: cross_ts_number_count(crmax) = 0 !<
1923	INTEGER(iwp) :: dopr_index(300) = 0 !<
1924	INTEGER(iwp) :: dopr_initial_index(300) = 0 !<
1925	INTEGER(iwp) :: dots_crossindex(100) = 0 !<
1926	INTEGER(iwp) :: dots_index(100) = 0 !<
1927
1928
1929	REAL(wp) :: cross_ts_uymax(20) = & !<
1930	(/ 999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp, &
1931	999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp, &
1932	999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp, &
1933	999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp /)
1934	REAL(wp) :: cross_ts_uymax_computed(20) = 999.999_wp !<
1935	REAL(wp) :: cross_ts_uymin(20) = & !<
1936	(/ 999.999_wp, 999.999_wp, 999.999_wp, -5.000_wp, 999.999_wp, &
1937	999.999_wp, 0.000_wp, 999.999_wp, 999.999_wp, 999.999_wp, &
1938	999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp, &
1939	999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp, 999.999_wp /)
1940	REAL(wp) :: cross_ts_uymin_computed(20) = 999.999_wp !<
1941
1942	SAVE
1943
1944	END MODULE profil_parameter
1945
1946	!------------------------------------------------------------------------------!
1947	! Description:
1948	! ------------
1949	!> Definition of statistical quantities, e.g. global sums.
1950	!------------------------------------------------------------------------------!
1951	MODULE statistics
1952
1953	USE kinds
1954
1955	CHARACTER (LEN=40) :: region(0:9) !<
1956
1957	INTEGER(iwp) :: pr_palm = 130 !<
1958	INTEGER(iwp) :: statistic_regions = 0 !<
1959
1960	INTEGER(iwp) :: u_max_ijk(3) = -1 !<
1961	INTEGER(iwp) :: v_max_ijk(3) = -1 !<
1962	INTEGER(iwp) :: w_max_ijk(3) = -1 !<
1963
1964	LOGICAL :: flow_statistics_called = .FALSE. !<
1965
1966	REAL(wp) :: u_max !<
1967	REAL(wp) :: v_max !<
1968	REAL(wp) :: w_max !<
1969
1970	REAL(wp), DIMENSION(:), ALLOCATABLE :: mean_surface_level_height !<
1971	REAL(wp), DIMENSION(:), ALLOCATABLE :: sums_divnew_l !<
1972	REAL(wp), DIMENSION(:), ALLOCATABLE :: sums_divold_l !<
1973	REAL(wp), DIMENSION(:), ALLOCATABLE :: weight_substep !<
1974	REAL(wp), DIMENSION(:), ALLOCATABLE :: weight_pres !<
1975
1976	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: sums !<
1977	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: sums_wsts_bc_l !<
1978	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: ts_value !<
1979	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: sums_wsus_ws_l !<
1980	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: sums_wsvs_ws_l !<
1981	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: sums_us2_ws_l !<
1982	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: sums_vs2_ws_l !<
1983	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: sums_ws2_ws_l !<
1984	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: sums_wsnrs_ws_l !<
1985	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: sums_wspts_ws_l !<
1986	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: sums_wsqs_ws_l !<
1987	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: sums_wsqrs_ws_l !<
1988	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: sums_wssas_ws_l !<
1989	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: sums_wsss_ws_l !<
1990	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: sums_ls_l !<
1991
1992	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: hom_sum !<
1993	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: rmask !<
1994	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: sums_l !<
1995	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: sums_l_l !<
1996	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: sums_up_fraction_l !<
1997
1998	REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: hom !<
1999
2000	SAVE
2001
2002	END MODULE statistics
2003
2004
2005
2006	!------------------------------------------------------------------------------!
2007	! Description:
2008	! ------------
2009	!> Definition of indices for transposed arrays.
2010	!------------------------------------------------------------------------------!
2011	MODULE transpose_indices
2012
2013	USE kinds
2014
2015	INTEGER(iwp) :: nxl_y !<
2016	INTEGER(iwp) :: nxl_yd !<
2017	INTEGER(iwp) :: nxl_z !<
2018	INTEGER(iwp) :: nxr_y !<
2019	INTEGER(iwp) :: nxr_yd !<
2020	INTEGER(iwp) :: nxr_z !<
2021	INTEGER(iwp) :: nyn_x !<
2022	INTEGER(iwp) :: nyn_z !<
2023	INTEGER(iwp) :: nys_x !<
2024	INTEGER(iwp) :: nys_z !<
2025	INTEGER(iwp) :: nzb_x !<
2026	INTEGER(iwp) :: nzb_y !<
2027	INTEGER(iwp) :: nzb_yd !<
2028	INTEGER(iwp) :: nzt_x !<
2029	INTEGER(iwp) :: nzt_y !<
2030	INTEGER(iwp) :: nzt_yd !<
2031
2032	SAVE
2033
2034	END MODULE transpose_indices

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |