Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |