Home

Context Navigation

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

Last change on this file since 2128 was 2119, checked in by raasch, 8 years ago
last commit documented
Property svn:keywords set to `Id`
File size: 99.9 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |