Home

Context Navigation

nesting_offl_mod.f90 @ 4875

Last change on this file since 4875 was 4843, checked in by raasch, 4 years ago
local namelist parameter added to switch off the module although the respective module namelist appears in the namelist file, further copyright updates
Property svn:keywords set to `Id`
File size: 174.2 KB

Rev	Line
[3347]	1	!> @file nesting_offl_mod.f90
[4834]	2	!--------------------------------------------------------------------------------------------------!
[3347]	3	! This file is part of the PALM model system.
	4	!
[4834]	5	! PALM is free software: you can redistribute it and/or modify it under the terms of the GNU General
	6	! Public License as published by the Free Software Foundation, either version 3 of the License, or
	7	! (at your option) any later version.
[3347]	8	!
[4834]	9	! PALM is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the
	10	! implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General
	11	! Public License for more details.
[3347]	12	!
[4834]	13	! You should have received a copy of the GNU General Public License along with PALM. If not, see
	14	! <http://www.gnu.org/licenses/>.
[3347]	15	!
[4828]	16	! Copyright 1997-2021 Leibniz Universitaet Hannover
[4834]	17	!--------------------------------------------------------------------------------------------------!
[3347]	18	!
	19	! Current revisions:
	20	! ------------------
[4724]	21	!
	22	!
[3347]	23	! Former revisions:
	24	! -----------------
[3413]	25	! $Id: nesting_offl_mod.f90 4843 2021-01-15 15:22:11Z scharf $
[4843]	26	! local namelist parameter added to switch off the module although the respective module namelist
	27	! appears in the namelist file
	28	!
	29	! 4842 2021-01-14 10:42:28Z raasch
[4842]	30	! reading of namelist file and actions in case of namelist errors revised so that statement labels
	31	! and goto statements are not required any more
	32	!
	33	! 4834 2021-01-07 10:28:00Z raasch
[4834]	34	! file re-formatted to follow the PALM coding standard
	35	!
	36	! 4828 2021-01-05 11:21:41Z Giersch
[4795]	37	! Bugfix in obtaining the correct timestamp in case of restart runs
	38	!
[4834]	39	! 4724 2020-10-06 17:20:39Z suehring $
[4724]	40	! - Enable LOD=1 input of boundary conditions
	41	! - Minor bugfix - add missing initialization of the top boundary
[4834]	42	!
[4724]	43	! 4582 2020-06-29 09:22:11Z suehring
[4582]	44	! Remove unused variable
[4834]	45	!
[4582]	46	! 4581 2020-06-29 08:49:58Z suehring
[4834]	47	! Omit explicit pressure forcing via geostrophic wind components in case of mesoscale nesting.
	48	!
[4581]	49	! 4561 2020-06-12 07:05:56Z suehring
[4564]	50	! Adapt mass-flux correction also for the anelastic approximation
[4834]	51	!
[4564]	52	! 4561 2020-06-12 07:05:56Z suehring
[4457]	53	! use statement for exchange horiz added
[4834]	54	!
[4457]	55	! 4360 2020-01-07 11:25:50Z suehring
[4834]	56	! Bugfix, time coordinate is relative to origin_time rather than to 00:00:00 UTC.
	57	!
[4358]	58	! 4346 2019-12-18 11:55:56Z motisi
[4834]	59	! Introduction of wall_flags_total_0, which currently sets bits based on static topography
	60	! information used in wall_flags_static_0
	61	!
[4346]	62	! 4329 2019-12-10 15:46:36Z motisi
[4329]	63	! Renamed wall_flags_0 to wall_flags_static_0
[4834]	64	!
[4329]	65	! 4286 2019-10-30 16:01:14Z resler
[4286]	66	! Fix wrong checks of time from dynamic driver in nesting_offl_mod
[4834]	67	!
[4286]	68	! 4273 2019-10-24 13:40:54Z monakurppa
[4273]	69	! Add a logical switch nesting_offline_chem
[4834]	70	!
[4273]	71	! 4270 2019-10-23 10:46:20Z monakurppa
[4270]	72	! Implement offline nesting for salsa variables.
[4834]	73	!
[4270]	74	! 4231 2019-09-12 11:22:00Z suehring
[4231]	75	! Bugfix in array deallocation
[4834]	76	!
[4231]	77	! 4230 2019-09-11 13:58:14Z suehring
[4834]	78	! Update mean chemistry profiles. These are also used for rayleigh damping.
	79	!
[4230]	80	! 4227 2019-09-10 18:04:34Z gronemeier
	81	! implement new palm_date_time_mod
	82	!
[4226]	83	! - Data input moved into nesting_offl_mod
	84	! - check rephrased
	85	! - time variable is now relative to time_utc_init
	86	! - Define module specific data type for offline nesting in nesting_offl_mod
[4834]	87	!
[4226]	88	! 4182 2019-08-22 15:20:23Z scharf
[4182]	89	! Corrected "Former revisions" section
[4834]	90	!
[4182]	91	! 4169 2019-08-19 13:54:35Z suehring
[4169]	92	! Additional check added.
[4834]	93	!
[4169]	94	! 4168 2019-08-16 13:50:17Z suehring
[4168]	95	! Replace function get_topography_top_index by topo_top_ind
[4834]	96	!
[4168]	97	! 4125 2019-07-29 13:31:44Z suehring
[4834]	98	! In order to enable netcdf parallel access, allocate dummy arrays for the lateral boundary data on
	99	! cores that actually do not belong to these boundaries.
	100	!
[4125]	101	! 4079 2019-07-09 18:04:41Z suehring
[4834]	102	! - Set boundary condition for w at nzt+1 at the lateral boundaries, even though these won't enter
	103	! the numerical solution. However, due to the mass conservation these values might some up to very
	104	! large values which will occur in the run-control file
[4079]	105	! - Bugfix in offline nesting of chemical species
	106	! - Do not set Neumann conditions for TKE and passive scalar
[4834]	107	!
[4079]	108	! 4022 2019-06-12 11:52:39Z suehring
[4834]	109	! Detection of boundary-layer depth in stable boundary layer on basis of boundary data improved
[4022]	110	! Routine for boundary-layer depth calculation renamed and made public
[4834]	111	!
[4022]	112	! 3987 2019-05-22 09:52:13Z kanani
[3987]	113	! Introduce alternative switch for debug output during timestepping
[4834]	114	!
[3987]	115	! 3964 2019-05-09 09:48:32Z suehring
[4834]	116	! Ensure that veloctiy term in calculation of bulk Richardson number does not become zero
	117	!
[3964]	118	! 3937 2019-04-29 15:09:07Z suehring
[4834]	119	! Set boundary conditon on upper-left and upper-south grid point for the u- and v-component,
	120	! respectively.
	121	!
[3937]	122	! 3891 2019-04-12 17:52:01Z suehring
[4834]	123	! Bugfix, do not overwrite lateral and top boundary data in case of restart runs.
	124	!
[3891]	125	! 3885 2019-04-11 11:29:34Z kanani
[4834]	126	! Changes related to global restructuring of location messages and introduction of additional debug
	127	! messages
	128	!
	129	!
	130	! Do local data exchange for chemistry variables only when boundary data is coming from dynamic file
	131	!
[3858]	132	! 3737 2019-02-12 16:57:06Z suehring
[3737]	133	! Introduce mesoscale nesting for chemical species
[4834]	134	!
[3737]	135	! 3705 2019-01-29 19:56:39Z suehring
[3705]	136	! Formatting adjustments
[4834]	137	!
[3705]	138	! 3704 2019-01-29 19:51:41Z suehring
[4834]	139	! Check implemented for offline nesting in child domain
	140	!
	141	! Initial Revision:
[4182]	142	! - separate offline nesting from large_scale_nudging_mod
	143	! - revise offline nesting, adjust for usage of synthetic turbulence generator
[4834]	144	! - adjust Rayleigh damping depending on the time-depending atmospheric conditions
	145	!
	146	!
[3347]	147	! Description:
	148	! ------------
[4834]	149	!> Offline nesting in larger-scale models. Boundary conditions for the simulation are read from
	150	!> NetCDF file and are prescribed onto the respective arrays.
	151	!> Further, a mass-flux correction is performed to maintain the mass balance.
	152	!--------------------------------------------------------------------------------------------------!
[3347]	153	MODULE nesting_offl_mod
	154
[4834]	155	USE arrays_3d, &
	156	ONLY: diss, &
	157	drho_air_zw, &
	158	dzw, &
	159	e, &
	160	pt, &
	161	pt_init, &
	162	q, &
	163	q_init, &
	164	rdf, &
	165	rdf_sc, &
	166	rho_air, &
	167	rho_air_zw, &
	168	s, &
	169	u, &
	170	u_init, &
	171	ug, &
	172	v, &
	173	v_init, &
	174	vg, &
	175	w, &
	176	zu, &
[4226]	177	zw
	178
[4834]	179	USE basic_constants_and_equations_mod, &
	180	ONLY: g, &
	181	pi
[4273]	182
[4834]	183	USE chem_modules, &
[4273]	184	ONLY: chem_species, nesting_offline_chem
[3347]	185
[4834]	186	USE control_parameters, &
	187	ONLY: air_chemistry, &
	188	bc_dirichlet_l, &
	189	bc_dirichlet_n, &
	190	bc_dirichlet_r, &
	191	bc_dirichlet_s, &
	192	coupling_char, &
	193	constant_diffusion, &
	194	child_domain, &
	195	debug_output_timestep, &
	196	dt_3d, &
	197	dz, &
	198	end_time, &
	199	humidity, &
	200	initializing_actions, &
	201	message_string, &
	202	nesting_offline, &
	203	neutral, &
	204	passive_scalar, &
	205	rans_mode, &
	206	rans_tke_e, &
	207	rayleigh_damping_factor, &
	208	rayleigh_damping_height, &
	209	salsa, &
	210	spinup_time, &
	211	time_since_reference_point, &
[4169]	212	volume_flow
[4273]	213
[4834]	214	USE cpulog, &
	215	ONLY: cpu_log, &
	216	log_point, &
[4226]	217	log_point_s
[3347]	218
	219	USE grid_variables
	220
[4834]	221	USE indices, &
	222	ONLY: nbgp, nx, nxl, nxlg, nxlu, nxr, nxrg, ny, nys, nysv, nysg, nyn, nyng, nzb, nz, nzt, &
	223	topo_top_ind, wall_flags_total_0
[3347]	224
	225	USE kinds
	226
[4834]	227	USE netcdf_data_input_mod, &
	228	ONLY: char_fill, &
	229	char_lod, &
	230	check_existence, &
	231	close_input_file, &
	232	get_attribute, &
	233	get_dimension_length, &
	234	get_variable, &
	235	get_variable_pr, &
	236	input_pids_dynamic, &
	237	inquire_num_variables, &
	238	inquire_variable_names, &
	239	input_file_dynamic, &
	240	num_var_pids, &
	241	open_read_file, &
[4226]	242	pids_id
	243
[3347]	244	USE pegrid
	245
[4834]	246	USE salsa_mod, &
	247	ONLY: salsa_nesting_offl_bc, &
	248	salsa_nesting_offl_init, &
[4270]	249	salsa_nesting_offl_input
	250
[4226]	251	IMPLICIT NONE
	252
	253	!
	254	!-- Define data type for nesting in larger-scale models like COSMO.
	255	!-- Data type comprises u, v, w, pt, and q at lateral and top boundaries.
	256	TYPE nest_offl_type
	257
[4834]	258	CHARACTER(LEN=16) :: char_l = 'ls_forcing_left_' !< leading substring for variables at left boundary
	259	CHARACTER(LEN=17) :: char_n = 'ls_forcing_north_' !< leading substring for variables at north boundary
	260	CHARACTER(LEN=17) :: char_r = 'ls_forcing_right_' !< leading substring for variables at right boundary
	261	CHARACTER(LEN=17) :: char_s = 'ls_forcing_south_' !< leading substring for variables at south boundary
	262	CHARACTER(LEN=15) :: char_t = 'ls_forcing_top_' !< leading substring for variables at top boundary
[4226]	263
	264	CHARACTER(LEN=100), DIMENSION(:), ALLOCATABLE :: var_names !< list of variable in dynamic input file
[4834]	265	CHARACTER(LEN=100), DIMENSION(:), ALLOCATABLE :: var_names_chem_l !< names of mesoscale nested chemistry variables at left
	266	!< boundary
	267	CHARACTER(LEN=100), DIMENSION(:), ALLOCATABLE :: var_names_chem_n !< names of mesoscale nested chemistry variables at north
	268	!< boundary
	269	CHARACTER(LEN=100), DIMENSION(:), ALLOCATABLE :: var_names_chem_r !< names of mesoscale nested chemistry variables at right
	270	!< boundary
	271	CHARACTER(LEN=100), DIMENSION(:), ALLOCATABLE :: var_names_chem_s !< names of mesoscale nested chemistry variables at south
	272	!< boundary
	273	CHARACTER(LEN=100), DIMENSION(:), ALLOCATABLE :: var_names_chem_t !< names of mesoscale nested chemistry variables at top
	274	!< boundary
[4226]	275
[4724]	276	INTEGER(iwp) :: lod_east_pt = 2 !< level-of-detail of input data of potential temperature at the eastern boundary
	277	INTEGER(iwp) :: lod_east_qc = 2 !< level-of-detail of input data of cloud-water mixture fraction at the eastern boundary
	278	INTEGER(iwp) :: lod_east_qv = 2 !< level-of-detail of input data of specific humidity at the eastern boundary
	279	INTEGER(iwp) :: lod_east_u = 2 !< level-of-detail of input data of the u-component at the eastern boundary
	280	INTEGER(iwp) :: lod_east_v = 2 !< level-of-detail of input data of the v-component at the eastern boundary
	281	INTEGER(iwp) :: lod_east_w = 2 !< level-of-detail of input data of the w-component at the eastern boundary
	282	INTEGER(iwp) :: lod_north_pt = 2 !< level-of-detail of input data of potential temperature at the northern boundary
	283	INTEGER(iwp) :: lod_north_qc = 2 !< level-of-detail of input data of cloud-water mixture fraction at the northern boundary
	284	INTEGER(iwp) :: lod_north_qv = 2 !< level-of-detail of input data of specific humidity at the northern boundary
	285	INTEGER(iwp) :: lod_north_u = 2 !< level-of-detail of input data of the u-component at the northern boundary
	286	INTEGER(iwp) :: lod_north_v = 2 !< level-of-detail of input data of the v-component at the northern boundary
	287	INTEGER(iwp) :: lod_north_w = 2 !< level-of-detail of input data of the w-component at the northern boundary
	288	INTEGER(iwp) :: lod_south_pt = 2 !< level-of-detail of input data of potential temperature at the southern boundary
	289	INTEGER(iwp) :: lod_south_qc = 2 !< level-of-detail of input data of cloud-water mixture fraction at the southern boundary
	290	INTEGER(iwp) :: lod_south_qv = 2 !< level-of-detail of input data of specific humidity at the southern boundary
	291	INTEGER(iwp) :: lod_south_u = 2 !< level-of-detail of input data of the u-component at the southern boundary
	292	INTEGER(iwp) :: lod_south_v = 2 !< level-of-detail of input data of the v-component at the southern boundary
	293	INTEGER(iwp) :: lod_south_w = 2 !< level-of-detail of input data of the w-component at the southern boundary
	294	INTEGER(iwp) :: lod_top_pt = 2 !< level-of-detail of input data of potential temperature at the top boundary
	295	INTEGER(iwp) :: lod_top_qc = 2 !< level-of-detail of input data of cloud-water mixture fraction at the top boundary
	296	INTEGER(iwp) :: lod_top_qv = 2 !< level-of-detail of input data of specific humidity at the top boundary
	297	INTEGER(iwp) :: lod_top_u = 2 !< level-of-detail of input data of the u-component at the top boundary
	298	INTEGER(iwp) :: lod_top_v = 2 !< level-of-detail of input data of the v-component at the top boundary
	299	INTEGER(iwp) :: lod_top_w = 2 !< level-of-detail of input data of the w-component at the top boundary
	300	INTEGER(iwp) :: lod_west_pt = 2 !< level-of-detail of input data of potential temperature at the western boundary
	301	INTEGER(iwp) :: lod_west_qc = 2 !< level-of-detail of input data of cloud-water mixture fraction at the western boundary
	302	INTEGER(iwp) :: lod_west_qv = 2 !< level-of-detail of input data of specific humidity at the western boundary
	303	INTEGER(iwp) :: lod_west_u = 2 !< level-of-detail of input data of the u-component at the western boundary
	304	INTEGER(iwp) :: lod_west_v = 2 !< level-of-detail of input data of the v-component at the western boundary
	305	INTEGER(iwp) :: lod_west_w = 2 !< level-of-detail of input data of the w-component at the western boundary
	306	INTEGER(iwp) :: nt !< number of time levels in dynamic input file
	307	INTEGER(iwp) :: nzu !< number of vertical levels on scalar grid in dynamic input file
	308	INTEGER(iwp) :: nzw !< number of vertical levels on w grid in dynamic input file
[4795]	309	INTEGER(iwp) :: tind = 0 !< time index for reference time in mesoscale-offline nesting
	310	INTEGER(iwp) :: tind_p = 0 !< time index for following time in mesoscale-offline nesting
[4226]	311
[4724]	312	LOGICAL :: init = .FALSE. !< flag indicating that offline nesting is already initialized
[4226]	313
[4834]	314	LOGICAL, DIMENSION(:), ALLOCATABLE :: chem_from_file_l !< flags inidicating whether left boundary data for chemistry is in
	315	!< dynamic input file
	316	LOGICAL, DIMENSION(:), ALLOCATABLE :: chem_from_file_n !< flags inidicating whether north boundary data for chemistry is in
	317	!< dynamic input file
	318	LOGICAL, DIMENSION(:), ALLOCATABLE :: chem_from_file_r !< flags inidicating whether right boundary data for chemistry is in
	319	!< dynamic input file
	320	LOGICAL, DIMENSION(:), ALLOCATABLE :: chem_from_file_s !< flags inidicating whether south boundary data for chemistry is in
	321	!< dynamic input file
	322	LOGICAL, DIMENSION(:), ALLOCATABLE :: chem_from_file_t !< flags inidicating whether top boundary data for chemistry is in
	323	!< dynamic input file
[4226]	324
[4834]	325	REAL(wp), DIMENSION(:), ALLOCATABLE :: surface_pressure !< time dependent surface pressure
	326	REAL(wp), DIMENSION(:), ALLOCATABLE :: time !< time levels in dynamic input file
	327	REAL(wp), DIMENSION(:), ALLOCATABLE :: zu_atmos !< vertical levels at scalar grid in dynamic input file
	328	REAL(wp), DIMENSION(:), ALLOCATABLE :: zw_atmos !< vertical levels at w grid in dynamic input file
[4226]	329
	330	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: ug !< domain-averaged geostrophic component
	331	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: vg !< domain-averaged geostrophic component
	332
[4834]	333	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: pt_l !< potentital temperautre at left boundary
	334	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: pt_n !< potentital temperautre at north boundary
	335	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: pt_r !< potentital temperautre at right boundary
	336	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: pt_s !< potentital temperautre at south boundary
	337	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: pt_top !< potentital temperautre at top boundary
[4724]	338	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: q_l !< mixing ratio at left boundary
	339	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: q_n !< mixing ratio at north boundary
	340	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: q_r !< mixing ratio at right boundary
	341	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: q_s !< mixing ratio at south boundary
	342	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: q_top !< mixing ratio at top boundary
	343	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: u_l !< u-component at left boundary
	344	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: u_n !< u-component at north boundary
	345	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: u_r !< u-component at right boundary
	346	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: u_s !< u-component at south boundary
[4226]	347	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: u_top !< u-component at top boundary
[4724]	348	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: v_l !< v-component at left boundary
	349	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: v_n !< v-component at north boundary
	350	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: v_r !< v-component at right boundary
	351	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: v_s !< v-component at south boundary
[4226]	352	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: v_top !< v-component at top boundary
[4724]	353	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: w_l !< w-component at left boundary
	354	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: w_n !< w-component at north boundary
	355	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: w_r !< w-component at right boundary
	356	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: w_s !< w-component at south boundary
[4226]	357	REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: w_top !< w-component at top boundary
	358
[4724]	359	REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: chem_l !< chemical species at left boundary
	360	REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: chem_n !< chemical species at north boundary
	361	REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: chem_r !< chemical species at right boundary
	362	REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: chem_s !< chemical species at south boundary
	363	REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: chem_top !< chemical species at left boundary
	364
[4226]	365	END TYPE nest_offl_type
	366
[4724]	367	INTEGER(iwp) :: i_bound !< boundary grid point in x-direction for scalars, v, and w
	368	INTEGER(iwp) :: i_bound_u !< boundary grid point in x-direction for u
[4834]	369	INTEGER(iwp) :: i_end !< end index for array allocation along x-direction at norther/southern boundary
[4724]	370	INTEGER(iwp) :: i_start !< start index for array allocation along x-direction at norther/southern boundary (scalars, v, w)
	371	INTEGER(iwp) :: i_start_u !< start index for array allocation along x-direction at norther/southern boundary (u)
	372	INTEGER(iwp) :: j_bound !< boundary grid point in y-direction for scalars, u, and w
	373	INTEGER(iwp) :: j_bound_v !< boundary grid point in y-direction for v
[4834]	374	INTEGER(iwp) :: j_end !< end index for array allocation along y-direction at eastern/western boundary
[4724]	375	INTEGER(iwp) :: j_start !< start index for array allocation along y-direction at eastern/western boundary (scalars, u, w)
	376	INTEGER(iwp) :: j_start_v !< start index for array allocation along y-direction at eastern/western boundary (v)
	377	INTEGER(iwp) :: lod !< level-of-detail of lateral input data
	378
[4226]	379	REAL(wp) :: fac_dt !< interpolation factor
[4834]	380	REAL(wp) :: zi_ribulk = 0.0_wp !< boundary-layer depth according to bulk Richardson criterion, i.e. the height where Ri_bulk
	381	!< exceeds the critical bulk Richardson number of 0.2
[4226]	382
[4834]	383	TYPE(nest_offl_type) :: nest_offl !< data structure for data input at lateral and top boundaries (provided by Inifor)
	384
[3347]	385	SAVE
	386	PRIVATE
	387	!
	388	!-- Public subroutines
[4834]	389	PUBLIC nesting_offl_bc, &
	390	nesting_offl_calc_zi, &
	391	nesting_offl_check_parameters, &
	392	nesting_offl_geostrophic_wind, &
	393	nesting_offl_header, &
	394	nesting_offl_init, &
	395	nesting_offl_input, &
	396	nesting_offl_interpolation_factor, &
	397	nesting_offl_mass_conservation, &
	398	nesting_offl_parin
[3347]	399	!
	400	!-- Public variables
[4834]	401	PUBLIC zi_ribulk
[3347]	402
	403	INTERFACE nesting_offl_bc
	404	MODULE PROCEDURE nesting_offl_bc
	405	END INTERFACE nesting_offl_bc
[4834]	406
[4022]	407	INTERFACE nesting_offl_calc_zi
	408	MODULE PROCEDURE nesting_offl_calc_zi
	409	END INTERFACE nesting_offl_calc_zi
[4834]	410
[3579]	411	INTERFACE nesting_offl_check_parameters
	412	MODULE PROCEDURE nesting_offl_check_parameters
	413	END INTERFACE nesting_offl_check_parameters
[4226]	414
	415	INTERFACE nesting_offl_geostrophic_wind
	416	MODULE PROCEDURE nesting_offl_geostrophic_wind
	417	END INTERFACE nesting_offl_geostrophic_wind
[4834]	418
[3347]	419	INTERFACE nesting_offl_header
	420	MODULE PROCEDURE nesting_offl_header
	421	END INTERFACE nesting_offl_header
[4834]	422
[3347]	423	INTERFACE nesting_offl_init
	424	MODULE PROCEDURE nesting_offl_init
	425	END INTERFACE nesting_offl_init
[4226]	426
	427	INTERFACE nesting_offl_input
	428	MODULE PROCEDURE nesting_offl_input
	429	END INTERFACE nesting_offl_input
	430
	431	INTERFACE nesting_offl_interpolation_factor
	432	MODULE PROCEDURE nesting_offl_interpolation_factor
	433	END INTERFACE nesting_offl_interpolation_factor
[4834]	434
[3347]	435	INTERFACE nesting_offl_mass_conservation
	436	MODULE PROCEDURE nesting_offl_mass_conservation
	437	END INTERFACE nesting_offl_mass_conservation
[4834]	438
[3347]	439	INTERFACE nesting_offl_parin
	440	MODULE PROCEDURE nesting_offl_parin
	441	END INTERFACE nesting_offl_parin
	442
	443	CONTAINS
	444
[4834]	445	!--------------------------------------------------------------------------------------------------!
[4226]	446	! Description:
	447	! ------------
	448	!> Reads data at lateral and top boundaries derived from larger-scale model.
[4834]	449	!--------------------------------------------------------------------------------------------------!
	450	SUBROUTINE nesting_offl_input
[3347]	451
[4834]	452	INTEGER(iwp) :: n !< running index for chemistry variables
[4226]	453
	454	!
[4834]	455	!-- Initialize INIFOR forcing in first call.
	456	IF ( .NOT. nest_offl%init ) THEN
[4226]	457	#if defined ( __netcdf )
	458	!
[4834]	459	!-- Open file in read-only mode
	460	CALL open_read_file( TRIM( input_file_dynamic ) // TRIM( coupling_char ), pids_id )
[4226]	461	!
[4834]	462	!-- At first, inquire all variable names.
	463	CALL inquire_num_variables( pids_id, num_var_pids )
[4226]	464	!
[4834]	465	!-- Allocate memory to store variable names.
	466	ALLOCATE( nest_offl%var_names(1:num_var_pids) )
	467	CALL inquire_variable_names( pids_id, nest_offl%var_names )
[4226]	468	!
[4834]	469	!-- Read time dimension, allocate memory and finally read time array
	470	CALL get_dimension_length( pids_id, nest_offl%nt, 'time' )
[4226]	471
[4834]	472	IF ( check_existence( nest_offl%var_names, 'time' ) ) THEN
	473	ALLOCATE( nest_offl%time(0:nest_offl%nt-1) )
	474	CALL get_variable( pids_id, 'time', nest_offl%time )
	475	ENDIF
[4226]	476	!
[4834]	477	!-- Read vertical dimension of scalar und w grid
	478	CALL get_dimension_length( pids_id, nest_offl%nzu, 'z' )
	479	CALL get_dimension_length( pids_id, nest_offl%nzw, 'zw' )
[4226]	480
[4834]	481	IF ( check_existence( nest_offl%var_names, 'z' ) ) THEN
	482	ALLOCATE( nest_offl%zu_atmos(1:nest_offl%nzu) )
	483	CALL get_variable( pids_id, 'z', nest_offl%zu_atmos )
	484	ENDIF
	485	IF ( check_existence( nest_offl%var_names, 'zw' ) ) THEN
	486	ALLOCATE( nest_offl%zw_atmos(1:nest_offl%nzw) )
	487	CALL get_variable( pids_id, 'zw', nest_offl%zw_atmos )
	488	ENDIF
[4226]	489	!
[4834]	490	!-- Read surface pressure
	491	IF ( check_existence( nest_offl%var_names, 'surface_forcing_surface_pressure' ) ) THEN
	492	ALLOCATE( nest_offl%surface_pressure(0:nest_offl%nt-1) )
	493	CALL get_variable( pids_id, 'surface_forcing_surface_pressure', &
	494	nest_offl%surface_pressure )
	495	ENDIF
[4226]	496	!
[4834]	497	!-- Close input file
	498	CALL close_input_file( pids_id )
[4226]	499	#endif
[4834]	500	ENDIF
[4226]	501	!
[4834]	502	!-- Check if dynamic driver data input is required.
	503	IF ( nest_offl%time(nest_offl%tind_p) <= MAX( time_since_reference_point, 0.0_wp) .OR. &
	504	.NOT. nest_offl%init ) THEN
	505	CONTINUE
[4226]	506	!
[4834]	507	!-- Return otherwise
	508	ELSE
	509	RETURN
	510	ENDIF
[4226]	511	!
[4834]	512	!-- Start of CPU measurement
	513	CALL cpu_log( log_point_s(86), 'NetCDF input forcing', 'start' )
	514
[4226]	515	!
[4834]	516	!-- Obtain time index for current point in time. Note, the time coordinate in the input file is
	517	!-- always relative to the initial time in UTC, i.e. the time coordinate always starts at 0.0 even
	518	!-- if the initial UTC is e.g. 7200.0. Further, since time_since_reference_point is negativ here
	519	!-- when spinup is applied, use MAX function to obtain correct time index.
	520	nest_offl%tind = MINLOC( ABS( nest_offl%time - MAX( time_since_reference_point, 0.0_wp ) ), &
	521	DIM = 1 ) - 1
[4795]	522	!
[4834]	523	!-- Note, in case of restart runs, the time index for the boundary data may indicate a time in
	524	!-- the future. This needs to be checked and corrected.
[4795]	525	IF ( TRIM( initializing_actions ) == 'read_restart_data' .AND. &
	526	nest_offl%time(nest_offl%tind) > time_since_reference_point ) THEN
	527	nest_offl%tind = nest_offl%tind - 1
	528	ENDIF
[4834]	529	nest_offl%tind_p = nest_offl%tind + 1
[4226]	530	!
[4834]	531	!-- Open file in read-only mode
[4226]	532	#if defined ( __netcdf )
[4834]	533	CALL open_read_file( TRIM( input_file_dynamic ) // TRIM( coupling_char ), pids_id )
[4226]	534	!
	535	!-- Read geostrophic wind components
[4581]	536	! DO t = nest_offl%tind, nest_offl%tind_p
	537	! CALL get_variable_pr( pids_id, 'ls_forcing_ug', t+1, &
	538	! nest_offl%ug(t-nest_offl%tind,nzb+1:nzt) )
	539	! CALL get_variable_pr( pids_id, 'ls_forcing_vg', t+1, &
	540	! nest_offl%vg(t-nest_offl%tind,nzb+1:nzt) )
	541	! ENDDO
[4226]	542	!
[4834]	543	!-- Read data at lateral and top boundaries. Please note, at left and right domain boundary,
	544	!-- yz-layers are read for u, v, w, pt and q.
	545	!-- For the v-component, the data starts at nysv, while for the other quantities the data starts at
	546	!-- nys. This is equivalent at the north and south domain boundary for the u-component (nxlu).
	547	!-- Note, lateral data is also accessed by parallel IO, which is the reason why different arguments
	548	!-- are passed depending on the boundary control flags. Cores that do not belong to the respective
	549	!-- boundary only do a dummy read with count = 0, just in order to participate the collective
	550	!-- operation. This is because collective parallel access shows better performance than just a
	551	!-- conditional access.
	552	!-- Read data for LOD 2, i.e. time-dependent xz-, yz-, and xy-slices.
	553	IF ( lod == 2 ) THEN
	554	CALL get_variable( pids_id, 'ls_forcing_left_u', &
	555	nest_offl%u_l, & ! array to be read
	556	MERGE( nys+1, 1, bc_dirichlet_l), & ! start index y direction
	557	MERGE( nzb+1, 1, bc_dirichlet_l), & ! start index z direction
	558	MERGE( nest_offl%tind+1, 1, bc_dirichlet_l), & ! start index time dimension
	559	MERGE( nyn-nys+1, 0, bc_dirichlet_l), & ! number of elements along y
	560	MERGE( nest_offl%nzu, 0, bc_dirichlet_l), & ! number of elements alogn z
	561	MERGE( 2, 0, bc_dirichlet_l), & ! number of time steps (2 or 0)
	562	.TRUE. ) ! parallel IO when compiled accordingly
[4226]	563
[4834]	564	CALL get_variable( pids_id, 'ls_forcing_left_v', &
	565	nest_offl%v_l, &
	566	MERGE( nysv, 1, bc_dirichlet_l), &
	567	MERGE( nzb+1, 1, bc_dirichlet_l), &
	568	MERGE( nest_offl%tind+1, 1, bc_dirichlet_l), &
	569	MERGE( nyn-nysv+1, 0, bc_dirichlet_l), &
	570	MERGE( nest_offl%nzu, 0, bc_dirichlet_l), &
	571	MERGE( 2, 0, bc_dirichlet_l), &
	572	.TRUE. )
	573
	574	CALL get_variable( pids_id, 'ls_forcing_left_w', &
	575	nest_offl%w_l, &
	576	MERGE( nys+1, 1, bc_dirichlet_l), &
	577	MERGE( nzb+1, 1, bc_dirichlet_l), &
	578	MERGE( nest_offl%tind+1, 1, bc_dirichlet_l), &
	579	MERGE( nyn-nys+1, 0, bc_dirichlet_l), &
	580	MERGE( nest_offl%nzw, 0, bc_dirichlet_l), &
	581	MERGE( 2, 0, bc_dirichlet_l), &
	582	.TRUE. )
	583
	584	IF ( .NOT. neutral ) THEN
	585	CALL get_variable( pids_id, 'ls_forcing_left_pt', &
	586	nest_offl%pt_l, &
	587	MERGE( nys+1, 1, bc_dirichlet_l), &
[4724]	588	MERGE( nzb+1, 1, bc_dirichlet_l), &
	589	MERGE( nest_offl%tind+1, 1, bc_dirichlet_l), &
[4834]	590	MERGE( nyn-nys+1, 0, bc_dirichlet_l), &
[4724]	591	MERGE( nest_offl%nzu, 0, bc_dirichlet_l), &
	592	MERGE( 2, 0, bc_dirichlet_l), &
[4226]	593	.TRUE. )
[4834]	594	ENDIF
[4226]	595
[4834]	596	IF ( humidity ) THEN
	597	CALL get_variable( pids_id, 'ls_forcing_left_qv', &
	598	nest_offl%q_l, &
[4724]	599	MERGE( nys+1, 1, bc_dirichlet_l), &
	600	MERGE( nzb+1, 1, bc_dirichlet_l), &
	601	MERGE( nest_offl%tind+1, 1, bc_dirichlet_l), &
	602	MERGE( nyn-nys+1, 0, bc_dirichlet_l), &
[4834]	603	MERGE( nest_offl%nzu, 0, bc_dirichlet_l), &
[4724]	604	MERGE( 2, 0, bc_dirichlet_l), &
[4834]	605	.TRUE. )
	606	ENDIF
[4273]	607
[4834]	608	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	609	DO n = 1, UBOUND( nest_offl%var_names_chem_l, 1 )
	610	IF ( check_existence( nest_offl%var_names, nest_offl%var_names_chem_l(n) ) ) THEN
	611	CALL get_variable( pids_id, &
	612	TRIM( nest_offl%var_names_chem_l(n) ), &
	613	nest_offl%chem_l(:,:,:,n), &
	614	MERGE( nys+1, 1, bc_dirichlet_l), &
	615	MERGE( nzb+1, 1, bc_dirichlet_l), &
	616	MERGE( nest_offl%tind+1, 1, bc_dirichlet_l), &
	617	MERGE( nyn-nys+1, 0, bc_dirichlet_l), &
	618	MERGE( nest_offl%nzu, 0, bc_dirichlet_l), &
	619	MERGE( 2, 0, bc_dirichlet_l), &
	620	.TRUE. )
	621	nest_offl%chem_from_file_l(n) = .TRUE.
	622	ENDIF
	623	ENDDO
	624	ENDIF
	625	!
	626	!-- Read data for eastern boundary
	627	CALL get_variable( pids_id, 'ls_forcing_right_u', &
	628	nest_offl%u_r, &
	629	MERGE( nys+1, 1, bc_dirichlet_r), &
	630	MERGE( nzb+1, 1, bc_dirichlet_r), &
	631	MERGE( nest_offl%tind+1, 1, bc_dirichlet_r), &
	632	MERGE( nyn-nys+1, 0, bc_dirichlet_r), &
	633	MERGE( nest_offl%nzu, 0, bc_dirichlet_r), &
	634	MERGE( 2, 0, bc_dirichlet_r), &
	635	.TRUE. )
[4226]	636
[4834]	637	CALL get_variable( pids_id, 'ls_forcing_right_v', &
	638	nest_offl%v_r, &
	639	MERGE( nysv, 1, bc_dirichlet_r), &
	640	MERGE( nzb+1, 1, bc_dirichlet_r), &
	641	MERGE( nest_offl%tind+1, 1, bc_dirichlet_r), &
	642	MERGE( nyn-nysv+1, 0, bc_dirichlet_r), &
	643	MERGE( nest_offl%nzu, 0, bc_dirichlet_r), &
	644	MERGE( 2, 0, bc_dirichlet_r), &
	645	.TRUE. )
[4226]	646
[4834]	647	CALL get_variable( pids_id, 'ls_forcing_right_w', &
	648	nest_offl%w_r, &
	649	MERGE( nys+1, 1, bc_dirichlet_r), &
	650	MERGE( nzb+1, 1, bc_dirichlet_r), &
	651	MERGE( nest_offl%tind+1, 1, bc_dirichlet_r), &
	652	MERGE( nyn-nys+1, 0, bc_dirichlet_r), &
	653	MERGE( nest_offl%nzw, 0, bc_dirichlet_r), &
	654	MERGE( 2, 0, bc_dirichlet_r), &
	655	.TRUE. )
	656
	657	IF ( .NOT. neutral ) THEN
	658	CALL get_variable( pids_id, 'ls_forcing_right_pt', &
	659	nest_offl%pt_r, &
[4724]	660	MERGE( nys+1, 1, bc_dirichlet_r), &
	661	MERGE( nzb+1, 1, bc_dirichlet_r), &
	662	MERGE( nest_offl%tind+1, 1, bc_dirichlet_r), &
	663	MERGE( nyn-nys+1, 0, bc_dirichlet_r), &
	664	MERGE( nest_offl%nzu, 0, bc_dirichlet_r), &
	665	MERGE( 2, 0, bc_dirichlet_r), &
[4226]	666	.TRUE. )
[4834]	667	ENDIF
[4226]	668
[4834]	669	IF ( humidity ) THEN
	670	CALL get_variable( pids_id, 'ls_forcing_right_qv', &
	671	nest_offl%q_r, &
[4724]	672	MERGE( nys+1, 1, bc_dirichlet_r), &
	673	MERGE( nzb+1, 1, bc_dirichlet_r), &
	674	MERGE( nest_offl%tind+1, 1, bc_dirichlet_r), &
	675	MERGE( nyn-nys+1, 0, bc_dirichlet_r), &
[4834]	676	MERGE( nest_offl%nzu, 0, bc_dirichlet_r), &
[4724]	677	MERGE( 2, 0, bc_dirichlet_r), &
	678	.TRUE. )
[4834]	679	ENDIF
[4273]	680
[4834]	681	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	682	DO n = 1, UBOUND( nest_offl%var_names_chem_r, 1 )
	683	IF ( check_existence( nest_offl%var_names, nest_offl%var_names_chem_r(n) ) ) THEN
	684	CALL get_variable( pids_id, &
	685	TRIM( nest_offl%var_names_chem_r(n) ), &
	686	nest_offl%chem_r(:,:,:,n), &
	687	MERGE( nys+1, 1, bc_dirichlet_r), &
	688	MERGE( nzb+1, 1, bc_dirichlet_r), &
	689	MERGE( nest_offl%tind+1, 1, bc_dirichlet_r), &
	690	MERGE( nyn-nys+1, 0, bc_dirichlet_r), &
	691	MERGE( nest_offl%nzu, 0, bc_dirichlet_r), &
	692	MERGE( 2, 0, bc_dirichlet_r), &
	693	.TRUE. )
	694	nest_offl%chem_from_file_r(n) = .TRUE.
	695	ENDIF
	696	ENDDO
	697	ENDIF
	698	!
	699	!-- Read data for northern boundary
	700	CALL get_variable( pids_id, 'ls_forcing_north_u', &
	701	nest_offl%u_n, &
	702	MERGE( nxlu, 1, bc_dirichlet_n ), &
	703	MERGE( nzb+1, 1, bc_dirichlet_n ), &
	704	MERGE( nest_offl%tind+1, 1, bc_dirichlet_n ), &
	705	MERGE( nxr-nxlu+1, 0, bc_dirichlet_n ), &
	706	MERGE( nest_offl%nzu, 0, bc_dirichlet_n ), &
	707	MERGE( 2, 0, bc_dirichlet_n ), &
	708	.TRUE. )
[4273]	709
[4834]	710	CALL get_variable( pids_id, 'ls_forcing_north_v', &
	711	nest_offl%v_n, &
	712	MERGE( nxl+1, 1, bc_dirichlet_n ), &
	713	MERGE( nzb+1, 1, bc_dirichlet_n ), &
	714	MERGE( nest_offl%tind+1, 1, bc_dirichlet_n ), &
	715	MERGE( nxr-nxl+1, 0, bc_dirichlet_n ), &
	716	MERGE( nest_offl%nzu, 0, bc_dirichlet_n ), &
	717	MERGE( 2, 0, bc_dirichlet_n ), &
	718	.TRUE. )
[4273]	719
[4834]	720	CALL get_variable( pids_id, 'ls_forcing_north_w', &
	721	nest_offl%w_n, &
	722	MERGE( nxl+1, 1, bc_dirichlet_n ), &
	723	MERGE( nzb+1, 1, bc_dirichlet_n ), &
	724	MERGE( nest_offl%tind+1, 1, bc_dirichlet_n ), &
	725	MERGE( nxr-nxl+1, 0, bc_dirichlet_n ), &
	726	MERGE( nest_offl%nzw, 0, bc_dirichlet_n ), &
	727	MERGE( 2, 0, bc_dirichlet_n ), &
	728	.TRUE. )
[4724]	729
[4834]	730	IF ( .NOT. neutral ) THEN
	731	CALL get_variable( pids_id, 'ls_forcing_north_pt', &
	732	nest_offl%pt_n, &
[4724]	733	MERGE( nxl+1, 1, bc_dirichlet_n ), &
	734	MERGE( nzb+1, 1, bc_dirichlet_n ), &
	735	MERGE( nest_offl%tind+1, 1, bc_dirichlet_n ), &
	736	MERGE( nxr-nxl+1, 0, bc_dirichlet_n ), &
	737	MERGE( nest_offl%nzu, 0, bc_dirichlet_n ), &
	738	MERGE( 2, 0, bc_dirichlet_n ), &
[4273]	739	.TRUE. )
[4834]	740	ENDIF
	741	IF ( humidity ) THEN
	742	CALL get_variable( pids_id, 'ls_forcing_north_qv', &
	743	nest_offl%q_n, &
[4724]	744	MERGE( nxl+1, 1, bc_dirichlet_n ), &
	745	MERGE( nzb+1, 1, bc_dirichlet_n ), &
	746	MERGE( nest_offl%tind+1, 1, bc_dirichlet_n ), &
	747	MERGE( nxr-nxl+1, 0, bc_dirichlet_n ), &
[4834]	748	MERGE( nest_offl%nzu, 0, bc_dirichlet_n ), &
[4724]	749	MERGE( 2, 0, bc_dirichlet_n ), &
	750	.TRUE. )
[4834]	751	ENDIF
[4273]	752
[4834]	753	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	754	DO n = 1, UBOUND( nest_offl%var_names_chem_n, 1 )
	755	IF ( check_existence( nest_offl%var_names, nest_offl%var_names_chem_n(n) ) ) THEN
	756	CALL get_variable( pids_id, &
	757	TRIM( nest_offl%var_names_chem_n(n) ), &
	758	nest_offl%chem_n(:,:,:,n), &
	759	MERGE( nxl+1, 1, bc_dirichlet_n ), &
	760	MERGE( nzb+1, 1, bc_dirichlet_n ), &
	761	MERGE( nest_offl%tind+1, 1, bc_dirichlet_n ), &
	762	MERGE( nxr-nxl+1, 0, bc_dirichlet_n ), &
	763	MERGE( nest_offl%nzu, 0, bc_dirichlet_n ), &
	764	MERGE( 2, 0, bc_dirichlet_n ), &
	765	.TRUE. )
	766	nest_offl%chem_from_file_n(n) = .TRUE.
	767	ENDIF
	768	ENDDO
	769	ENDIF
[4724]	770	!
[4834]	771	!-- Read data for southern boundary
	772	CALL get_variable( pids_id, 'ls_forcing_south_u', &
	773	nest_offl%u_s, &
	774	MERGE( nxlu, 1, bc_dirichlet_s ), &
	775	MERGE( nzb+1, 1, bc_dirichlet_s ), &
	776	MERGE( nest_offl%tind+1, 1, bc_dirichlet_s ), &
	777	MERGE( nxr-nxlu+1, 0, bc_dirichlet_s ), &
	778	MERGE( nest_offl%nzu, 0, bc_dirichlet_s ), &
	779	MERGE( 2, 0, bc_dirichlet_s ), &
	780	.TRUE. )
[4273]	781
[4834]	782	CALL get_variable( pids_id, 'ls_forcing_south_v', &
	783	nest_offl%v_s, &
	784	MERGE( nxl+1, 1, bc_dirichlet_s ), &
	785	MERGE( nzb+1, 1, bc_dirichlet_s ), &
	786	MERGE( nest_offl%tind+1, 1, bc_dirichlet_s ), &
	787	MERGE( nxr-nxl+1, 0, bc_dirichlet_s ), &
	788	MERGE( nest_offl%nzu, 0, bc_dirichlet_s ), &
	789	MERGE( 2, 0, bc_dirichlet_s ), &
	790	.TRUE. )
	791
	792	CALL get_variable( pids_id, 'ls_forcing_south_w', &
	793	nest_offl%w_s, &
	794	MERGE( nxl+1, 1, bc_dirichlet_s ), &
	795	MERGE( nzb+1, 1, bc_dirichlet_s ), &
	796	MERGE( nest_offl%tind+1, 1, bc_dirichlet_s ), &
	797	MERGE( nxr-nxl+1, 0, bc_dirichlet_s ), &
	798	MERGE( nest_offl%nzw, 0, bc_dirichlet_s ), &
	799	MERGE( 2, 0, bc_dirichlet_s ), &
	800	.TRUE. )
	801
	802	IF ( .NOT. neutral ) THEN
	803	CALL get_variable( pids_id, 'ls_forcing_south_pt', &
	804	nest_offl%pt_s, &
[4724]	805	MERGE( nxl+1, 1, bc_dirichlet_s ), &
	806	MERGE( nzb+1, 1, bc_dirichlet_s ), &
	807	MERGE( nest_offl%tind+1, 1, bc_dirichlet_s ), &
	808	MERGE( nxr-nxl+1, 0, bc_dirichlet_s ), &
	809	MERGE( nest_offl%nzu, 0, bc_dirichlet_s ), &
	810	MERGE( 2, 0, bc_dirichlet_s ), &
[4273]	811	.TRUE. )
[4834]	812	ENDIF
	813	IF ( humidity ) THEN
	814	CALL get_variable( pids_id, 'ls_forcing_south_qv', &
	815	nest_offl%q_s, &
[4724]	816	MERGE( nxl+1, 1, bc_dirichlet_s ), &
	817	MERGE( nzb+1, 1, bc_dirichlet_s ), &
	818	MERGE( nest_offl%tind+1, 1, bc_dirichlet_s ), &
	819	MERGE( nxr-nxl+1, 0, bc_dirichlet_s ), &
[4834]	820	MERGE( nest_offl%nzu, 0, bc_dirichlet_s ), &
[4724]	821	MERGE( 2, 0, bc_dirichlet_s ), &
[4273]	822	.TRUE. )
[4834]	823	ENDIF
[4273]	824
[4834]	825	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	826	DO n = 1, UBOUND( nest_offl%var_names_chem_s, 1 )
	827	IF ( check_existence( nest_offl%var_names, nest_offl%var_names_chem_s(n) ) ) THEN
	828	CALL get_variable( pids_id, &
	829	TRIM( nest_offl%var_names_chem_s(n) ), &
	830	nest_offl%chem_s(:,:,:,n), &
	831	MERGE( nxl+1, 1, bc_dirichlet_s ), &
	832	MERGE( nzb+1, 1, bc_dirichlet_s ), &
	833	MERGE( nest_offl%tind+1, 1, bc_dirichlet_s ), &
	834	MERGE( nxr-nxl+1, 0, bc_dirichlet_s ), &
	835	MERGE( nest_offl%nzu, 0, bc_dirichlet_s ), &
	836	MERGE( 2, 0, bc_dirichlet_s ), &
	837	.TRUE. )
	838	nest_offl%chem_from_file_s(n) = .TRUE.
	839	ENDIF
	840	ENDDO
	841	ENDIF
[4226]	842	!
[4834]	843	!-- Top boundary
	844	CALL get_variable( pids_id, 'ls_forcing_top_u', &
	845	nest_offl%u_top(0:1,nys:nyn,nxlu:nxr), &
	846	nxlu, nys+1, nest_offl%tind+1, &
	847	nxr-nxlu+1, nyn-nys+1, 2, .TRUE. )
[4226]	848
[4834]	849	CALL get_variable( pids_id, 'ls_forcing_top_v', &
	850	nest_offl%v_top(0:1,nysv:nyn,nxl:nxr), &
	851	nxl+1, nysv, nest_offl%tind+1, &
	852	nxr-nxl+1, nyn-nysv+1, 2, .TRUE. )
[4273]	853
[4834]	854	CALL get_variable( pids_id, 'ls_forcing_top_w', &
	855	nest_offl%w_top(0:1,nys:nyn,nxl:nxr), &
	856	nxl+1, nys+1, nest_offl%tind+1, &
	857	nxr-nxl+1, nyn-nys+1, 2, .TRUE. )
	858
	859	IF ( .NOT. neutral ) THEN
	860	CALL get_variable( pids_id, 'ls_forcing_top_pt', &
	861	nest_offl%pt_top(0:1,nys:nyn,nxl:nxr), &
[4724]	862	nxl+1, nys+1, nest_offl%tind+1, &
[4226]	863	nxr-nxl+1, nyn-nys+1, 2, .TRUE. )
[4834]	864	ENDIF
	865	IF ( humidity ) THEN
	866	CALL get_variable( pids_id, 'ls_forcing_top_qv', &
	867	nest_offl%q_top(0:1,nys:nyn,nxl:nxr), &
	868	nxl+1, nys+1, nest_offl%tind+1, &
	869	nxr-nxl+1, nyn-nys+1, 2, .TRUE. )
	870	ENDIF
[4273]	871
[4834]	872	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	873	DO n = 1, UBOUND( nest_offl%var_names_chem_t, 1 )
	874	IF ( check_existence( nest_offl%var_names, nest_offl%var_names_chem_t(n) ) ) THEN
	875	CALL get_variable( pids_id, &
	876	TRIM( nest_offl%var_names_chem_t(n) ), &
	877	nest_offl%chem_top(0:1,nys:nyn,nxl:nxr,n), &
	878	nxl+1, nys+1, nest_offl%tind+1, &
	879	nxr-nxl+1, nyn-nys+1, 2, .TRUE. )
	880	nest_offl%chem_from_file_t(n) = .TRUE.
	881	ENDIF
	882	ENDDO
	883	ENDIF
	884	!
	885	!-- Read data for LOD 1, i.e. time-dependent profiles. In constrast to LOD 2 where the amount of IO
	886	!-- is larger, only the respective boundary processes read the data.
	887	ELSE
	888	IF ( bc_dirichlet_l ) THEN
	889	CALL get_variable( pids_id, 'ls_forcing_left_u', &
	890	nest_offl%u_l(0:1,:,1:1), & ! array to be read
	891	MERGE( nzb+1, 1, bc_dirichlet_l), & ! start index z direction
	892	MERGE( nest_offl%tind+1, 1, bc_dirichlet_l), & ! start index time dimension
	893	MERGE( nest_offl%nzu, 0, bc_dirichlet_l), & ! number of elements along z
	894	MERGE( 2, 0, bc_dirichlet_l) ) ! number of time steps (2 or 0)
	895	CALL get_variable( pids_id, 'ls_forcing_left_v', &
	896	nest_offl%v_l(0:1,:,1:1), &
	897	MERGE( nzb+1, 1, bc_dirichlet_l), &
	898	MERGE( nest_offl%tind+1, 1, bc_dirichlet_l), &
	899	MERGE( nest_offl%nzu, 0, bc_dirichlet_l), &
	900	MERGE( 2, 0, bc_dirichlet_l) )
	901	CALL get_variable( pids_id, 'ls_forcing_left_w', &
	902	nest_offl%w_l(0:1,:,1:1), &
	903	MERGE( nzb+1, 1, bc_dirichlet_l), &
	904	MERGE( nest_offl%tind+1, 1, bc_dirichlet_l), &
	905	MERGE( nest_offl%nzw, 0, bc_dirichlet_l), &
	906	MERGE( 2, 0, bc_dirichlet_l) )
[4724]	907	IF ( .NOT. neutral ) THEN
[4834]	908	CALL get_variable( pids_id, 'ls_forcing_left_pt', &
	909	nest_offl%pt_l(0:1,:,1:1), &
[4724]	910	MERGE( nzb+1, 1, bc_dirichlet_l), &
	911	MERGE( nest_offl%tind+1, 1, bc_dirichlet_l), &
	912	MERGE( nest_offl%nzu, 0, bc_dirichlet_l), &
	913	MERGE( 2, 0, bc_dirichlet_l) )
[4834]	914	ENDIF
	915	IF ( humidity ) THEN
	916	CALL get_variable( pids_id, 'ls_forcing_left_qv', &
	917	nest_offl%q_l(0:1,:,1:1), &
[4724]	918	MERGE( nzb+1, 1, bc_dirichlet_l), &
	919	MERGE( nest_offl%tind+1, 1, bc_dirichlet_l), &
[4834]	920	MERGE( nest_offl%nzu, 0, bc_dirichlet_l), &
[4724]	921	MERGE( 2, 0, bc_dirichlet_l) )
	922	ENDIF
[4834]	923	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	924	DO n = 1, UBOUND( nest_offl%var_names_chem_t, 1 )
	925	IF ( check_existence( nest_offl%var_names, nest_offl%var_names_chem_t(n) ) ) THEN
	926	CALL get_variable( pids_id, TRIM( nest_offl%var_names_chem_t(n) ), &
	927	nest_offl%chem_l(0:1,:,1:1,n), &
	928	MERGE( nzb+1, 1, bc_dirichlet_l), &
	929	MERGE( nest_offl%tind+1, 1, bc_dirichlet_l), &
	930	MERGE( nest_offl%nzu, 0, bc_dirichlet_l), &
	931	MERGE( 2, 0, bc_dirichlet_l) )
	932	nest_offl%chem_from_file_l(n) = .TRUE.
	933	ENDIF
	934	ENDDO
	935	ENDIF
	936	ENDIF
	937	IF ( bc_dirichlet_r ) THEN
	938	CALL get_variable( pids_id, 'ls_forcing_right_u', &
	939	nest_offl%u_r(0:1,:,1:1), &
	940	MERGE( nzb+1, 1, bc_dirichlet_r), &
	941	MERGE( nest_offl%tind+1, 1, bc_dirichlet_r), &
	942	MERGE( nest_offl%nzu, 0, bc_dirichlet_r), &
	943	MERGE( 2, 0, bc_dirichlet_r) )
	944	CALL get_variable( pids_id, 'ls_forcing_right_v', &
	945	nest_offl%v_r(0:1,:,1:1), &
	946	MERGE( nzb+1, 1, bc_dirichlet_r), &
	947	MERGE( nest_offl%tind+1, 1, bc_dirichlet_r), &
	948	MERGE( nest_offl%nzu, 0, bc_dirichlet_r), &
	949	MERGE( 2, 0, bc_dirichlet_r) )
	950	CALL get_variable( pids_id, 'ls_forcing_right_w', &
	951	nest_offl%w_r(0:1,:,1:1), &
	952	MERGE( nzb+1, 1, bc_dirichlet_r), &
	953	MERGE( nest_offl%tind+1, 1, bc_dirichlet_r), &
	954	MERGE( nest_offl%nzw, 0, bc_dirichlet_r), &
	955	MERGE( 2, 0, bc_dirichlet_r) )
	956	IF ( .NOT. neutral ) THEN
	957	CALL get_variable( pids_id, 'ls_forcing_right_pt', &
	958	nest_offl%pt_r(0:1,:,1:1), &
[4724]	959	MERGE( nzb+1, 1, bc_dirichlet_r), &
	960	MERGE( nest_offl%tind+1, 1, bc_dirichlet_r), &
	961	MERGE( nest_offl%nzu, 0, bc_dirichlet_r), &
	962	MERGE( 2, 0, bc_dirichlet_r) )
[4834]	963	ENDIF
	964	IF ( humidity ) THEN
	965	CALL get_variable( pids_id, 'ls_forcing_right_qv', &
	966	nest_offl%q_r(0:1,:,1:1), &
[4724]	967	MERGE( nzb+1, 1, bc_dirichlet_r), &
	968	MERGE( nest_offl%tind+1, 1, bc_dirichlet_r), &
	969	MERGE( nest_offl%nzu, 0, bc_dirichlet_r), &
	970	MERGE( 2, 0, bc_dirichlet_r) )
	971	ENDIF
[4834]	972	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	973	DO n = 1, UBOUND( nest_offl%var_names_chem_t, 1 )
	974	IF ( check_existence( nest_offl%var_names, nest_offl%var_names_chem_t(n) ) ) THEN
	975	CALL get_variable( pids_id, TRIM( nest_offl%var_names_chem_t(n) ), &
	976	nest_offl%chem_r(0:1,:,1:1,n), &
	977	MERGE( nzb+1, 1, bc_dirichlet_r), &
	978	MERGE( nest_offl%tind+1, 1, bc_dirichlet_r), &
	979	MERGE( nest_offl%nzu, 0, bc_dirichlet_r), &
	980	MERGE( 2, 0, bc_dirichlet_r) )
	981	nest_offl%chem_from_file_r(n) = .TRUE.
	982	ENDIF
	983	ENDDO
	984	ENDIF
	985	ENDIF
	986	IF ( bc_dirichlet_n ) THEN
	987	CALL get_variable( pids_id, 'ls_forcing_north_u', &
	988	nest_offl%u_n(0:1,:,1:1), &
	989	MERGE( nzb+1, 1, bc_dirichlet_n), &
	990	MERGE( nest_offl%tind+1, 1, bc_dirichlet_n), &
	991	MERGE( nest_offl%nzu, 0, bc_dirichlet_n), &
	992	MERGE( 2, 0, bc_dirichlet_n) )
	993	CALL get_variable( pids_id, 'ls_forcing_north_v', &
	994	nest_offl%v_n(0:1,:,1:1), &
	995	MERGE( nzb+1, 1, bc_dirichlet_n), &
	996	MERGE( nest_offl%tind+1, 1, bc_dirichlet_n), &
	997	MERGE( nest_offl%nzu, 0, bc_dirichlet_n), &
	998	MERGE( 2, 0, bc_dirichlet_n) )
	999	CALL get_variable( pids_id, 'ls_forcing_north_w', &
	1000	nest_offl%w_n(0:1,:,1:1), &
	1001	MERGE( nzb+1, 1, bc_dirichlet_n), &
	1002	MERGE( nest_offl%tind+1, 1, bc_dirichlet_n), &
	1003	MERGE( nest_offl%nzw, 0, bc_dirichlet_n), &
	1004	MERGE( 2, 0, bc_dirichlet_n) )
	1005	IF ( .NOT. neutral ) THEN
	1006	CALL get_variable( pids_id, 'ls_forcing_north_pt', &
	1007	nest_offl%pt_n(0:1,:,1:1), &
[4724]	1008	MERGE( nzb+1, 1, bc_dirichlet_n), &
	1009	MERGE( nest_offl%tind+1, 1, bc_dirichlet_n), &
	1010	MERGE( nest_offl%nzu, 0, bc_dirichlet_n), &
	1011	MERGE( 2, 0, bc_dirichlet_n) )
[4834]	1012	ENDIF
	1013	IF ( humidity ) THEN
	1014	CALL get_variable( pids_id, 'ls_forcing_north_qv', &
	1015	nest_offl%q_n(0:1,:,1:1), &
[4724]	1016	MERGE( nzb+1, 1, bc_dirichlet_n), &
	1017	MERGE( nest_offl%tind+1, 1, bc_dirichlet_n), &
	1018	MERGE( nest_offl%nzu, 0, bc_dirichlet_n), &
	1019	MERGE( 2, 0, bc_dirichlet_n) )
	1020	ENDIF
[4834]	1021	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	1022	DO n = 1, UBOUND( nest_offl%var_names_chem_t, 1 )
	1023	IF ( check_existence( nest_offl%var_names, nest_offl%var_names_chem_t(n) ) ) THEN
	1024	CALL get_variable( pids_id, TRIM( nest_offl%var_names_chem_t(n) ), &
	1025	nest_offl%chem_n(0:1,:,1:1,n), &
	1026	MERGE( nzb+1, 1, bc_dirichlet_n), &
	1027	MERGE( nest_offl%tind+1, 1, bc_dirichlet_n), &
	1028	MERGE( nest_offl%nzu, 0, bc_dirichlet_n), &
	1029	MERGE( 2, 0, bc_dirichlet_n) )
	1030	nest_offl%chem_from_file_n(n) = .TRUE.
	1031	ENDIF
	1032	ENDDO
	1033	ENDIF
	1034	ENDIF
	1035	IF ( bc_dirichlet_s ) THEN
	1036	CALL get_variable( pids_id, 'ls_forcing_south_u', &
	1037	nest_offl%u_s(0:1,:,1:1), &
	1038	MERGE( nzb+1, 1, bc_dirichlet_s), &
	1039	MERGE( nest_offl%tind+1, 1, bc_dirichlet_s), &
	1040	MERGE( nest_offl%nzu, 0, bc_dirichlet_s), &
	1041	MERGE( 2, 0, bc_dirichlet_s) )
	1042	CALL get_variable( pids_id, 'ls_forcing_south_v', &
	1043	nest_offl%v_s(0:1,:,1:1), &
	1044	MERGE( nzb+1, 1, bc_dirichlet_s), &
	1045	MERGE( nest_offl%tind+1, 1, bc_dirichlet_s), &
	1046	MERGE( nest_offl%nzu, 0, bc_dirichlet_s), &
	1047	MERGE( 2, 0, bc_dirichlet_s) )
	1048	CALL get_variable( pids_id, 'ls_forcing_south_w', &
	1049	nest_offl%w_s(0:1,:,1:1), &
	1050	MERGE( nzb+1, 1, bc_dirichlet_s), &
	1051	MERGE( nest_offl%tind+1, 1, bc_dirichlet_s), &
	1052	MERGE( nest_offl%nzw, 0, bc_dirichlet_s), &
	1053	MERGE( 2, 0, bc_dirichlet_s) )
	1054	IF ( .NOT. neutral ) THEN
	1055	CALL get_variable( pids_id, 'ls_forcing_south_pt', &
	1056	nest_offl%pt_s(0:1,:,1:1), &
[4724]	1057	MERGE( nzb+1, 1, bc_dirichlet_s), &
	1058	MERGE( nest_offl%tind+1, 1, bc_dirichlet_s), &
	1059	MERGE( nest_offl%nzu, 0, bc_dirichlet_s), &
	1060	MERGE( 2, 0, bc_dirichlet_s) )
[4834]	1061	ENDIF
	1062	IF ( humidity ) THEN
	1063	CALL get_variable( pids_id, 'ls_forcing_south_qv', &
	1064	nest_offl%q_s(0:1,:,1:1), &
[4724]	1065	MERGE( nzb+1, 1, bc_dirichlet_s), &
	1066	MERGE( nest_offl%tind+1, 1, bc_dirichlet_s), &
	1067	MERGE( nest_offl%nzu, 0, bc_dirichlet_s), &
	1068	MERGE( 2, 0, bc_dirichlet_s) )
	1069	ENDIF
	1070	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
[4834]	1071	DO n = 1, UBOUND( nest_offl%var_names_chem_t, 1 )
	1072	IF ( check_existence( nest_offl%var_names, nest_offl%var_names_chem_t(n) ) ) THEN
[4724]	1073	CALL get_variable( pids_id, TRIM( nest_offl%var_names_chem_t(n) ), &
[4834]	1074	nest_offl%chem_s(0:1,:,1:1,n), &
	1075	MERGE( nzb+1, 1, bc_dirichlet_s), &
	1076	MERGE( nest_offl%tind+1, 1, bc_dirichlet_s), &
	1077	MERGE( nest_offl%nzu, 0, bc_dirichlet_s), &
	1078	MERGE( 2, 0, bc_dirichlet_s) )
	1079	nest_offl%chem_from_file_s(n) = .TRUE.
[4724]	1080	ENDIF
	1081	ENDDO
	1082	ENDIF
[4226]	1083	ENDIF
[4834]	1084	!
	1085	!-- Read top boundary data, which is actually only a scalar value in the LOD 1 case.
	1086	CALL get_variable( pids_id, 'ls_forcing_top_u', &
	1087	nest_offl%u_top(0:1,1,1), & ! array to be read
	1088	nest_offl%tind+1, & ! start index in time
	1089	2 ) ! number of elements to be read
	1090	CALL get_variable( pids_id, 'ls_forcing_top_v', &
	1091	nest_offl%v_top(0:1,1,1), &
	1092	nest_offl%tind+1, &
	1093	2 )
	1094	CALL get_variable( pids_id, 'ls_forcing_top_w', &
	1095	nest_offl%w_top(0:1,1,1), &
	1096	nest_offl%tind+1, &
	1097	2 )
	1098	IF ( .NOT. neutral ) THEN
	1099	CALL get_variable( pids_id, 'ls_forcing_top_pt', &
	1100	nest_offl%pt_top(0:1,1,1), &
	1101	nest_offl%tind+1, &
	1102	2 )
	1103	ENDIF
	1104	IF ( humidity ) THEN
	1105	CALL get_variable( pids_id, 'ls_forcing_top_qv', &
	1106	nest_offl%q_top(0:1,1,1), &
	1107	nest_offl%tind+1, &
	1108	2 )
	1109	ENDIF
	1110	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	1111	DO n = 1, UBOUND( nest_offl%var_names_chem_t, 1 )
	1112	IF ( check_existence( nest_offl%var_names, nest_offl%var_names_chem_t(n) ) ) THEN
	1113	CALL get_variable( pids_id, TRIM( nest_offl%var_names_chem_t(n) ), &
	1114	nest_offl%chem_top(0:1,1,1,n), &
	1115	nest_offl%tind+1, &
	1116	2 )
	1117	nest_offl%chem_from_file_t(n) = .TRUE.
	1118	ENDIF
	1119	ENDDO
	1120	ENDIF
	1121	ENDIF
[4226]	1122
[4724]	1123
[4226]	1124	!
[4834]	1125	!-- Close input file
	1126	CALL close_input_file( pids_id )
[4226]	1127	#endif
	1128	!
[4834]	1129	!-- Set control flag to indicate that boundary data has been initially input.
	1130	nest_offl%init = .TRUE.
[4226]	1131	!
[4834]	1132	!-- Call offline nesting for salsa
	1133	IF ( salsa ) CALL salsa_nesting_offl_input
[4270]	1134	!
[4834]	1135	!-- End of CPU measurement
	1136	CALL cpu_log( log_point_s(86), 'NetCDF input forcing', 'stop' )
[4226]	1137
[4834]	1138	END SUBROUTINE nesting_offl_input
[4226]	1139
	1140
[4834]	1141	!--------------------------------------------------------------------------------------------------!
[3347]	1142	! Description:
	1143	! ------------
[4834]	1144	!> In this subroutine a constant mass within the model domain is guaranteed.
	1145	!> Larger-scale models may be based on a compressible equation system, which is not consistent with
	1146	!> PALMs incompressible equation system. In order to avoid a decrease or increase of mass during the
	1147	!> simulation, non-divergent flow through the lateral and top boundaries is compensated by the
	1148	!> vertical wind component at the top boundary.
	1149	!--------------------------------------------------------------------------------------------------!
	1150	SUBROUTINE nesting_offl_mass_conservation
[3347]	1151
[4834]	1152	INTEGER(iwp) :: i !< grid index in x-direction
	1153	INTEGER(iwp) :: j !< grid index in y-direction
	1154	INTEGER(iwp) :: k !< grid index in z-direction
[3347]	1155
[4834]	1156	REAL(wp) :: d_area_t !< inverse of the total area of the horizontal model domain
	1157	REAL(wp) :: w_correct !< vertical velocity increment required to compensate non-divergent flow through the boundaries
[3347]	1158
[4834]	1159	REAL(wp), DIMENSION(1:3) :: volume_flow_l !< local volume flow
[3987]	1160
	1161
[4834]	1162	IF ( debug_output_timestep ) CALL debug_message( 'nesting_offl_mass_conservation', 'start' )
[3347]	1163
[4834]	1164	CALL cpu_log( log_point(58), 'offline nesting', 'start' )
[3347]	1165
[4834]	1166	volume_flow = 0.0_wp
	1167	volume_flow_l = 0.0_wp
	1168
	1169	d_area_t = 1.0_wp / ( ( nx + 1 ) * dx * ( ny + 1 ) * dy )
	1170
	1171	IF ( bc_dirichlet_l ) THEN
	1172	i = nxl
	1173	DO j = nys, nyn
	1174	DO k = nzb+1, nzt
	1175	volume_flow_l(1) = volume_flow_l(1) + u(k,j,i) * dzw(k) * dy * rho_air(k) &
	1176	* MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j,i), 1 ) )
[3347]	1177	ENDDO
[4834]	1178	ENDDO
	1179	ENDIF
	1180	IF ( bc_dirichlet_r ) THEN
	1181	i = nxr+1
	1182	DO j = nys, nyn
	1183	DO k = nzb+1, nzt
	1184	volume_flow_l(1) = volume_flow_l(1) - u(k,j,i) * dzw(k) * dy * rho_air(k) &
	1185	* MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j,i), 1 ) )
[3347]	1186	ENDDO
[4834]	1187	ENDDO
	1188	ENDIF
	1189	IF ( bc_dirichlet_s ) THEN
	1190	j = nys
	1191	DO i = nxl, nxr
	1192	DO k = nzb+1, nzt
	1193	volume_flow_l(2) = volume_flow_l(2) + v(k,j,i) * dzw(k) * dx * rho_air(k) &
	1194	* MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j,i), 2 ) )
[3347]	1195	ENDDO
[4834]	1196	ENDDO
	1197	ENDIF
	1198	IF ( bc_dirichlet_n ) THEN
	1199	j = nyn+1
[3347]	1200	DO i = nxl, nxr
[4834]	1201	DO k = nzb+1, nzt
	1202	volume_flow_l(2) = volume_flow_l(2) - v(k,j,i) * dzw(k) * dx * rho_air(k) &
	1203	* MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j,i), 2 ) )
[3347]	1204	ENDDO
	1205	ENDDO
[4834]	1206	ENDIF
	1207	!
	1208	!-- Top boundary
	1209	k = nzt
	1210	DO i = nxl, nxr
	1211	DO j = nys, nyn
	1212	volume_flow_l(3) = volume_flow_l(3) - rho_air_zw(k) * w(k,j,i) * dx * dy
	1213	ENDDO
	1214	ENDDO
[3347]	1215
	1216	#if defined( __parallel )
[4834]	1217	IF ( collective_wait ) CALL MPI_BARRIER( comm2d, ierr )
	1218	CALL MPI_ALLREDUCE( volume_flow_l, volume_flow, 3, MPI_REAL, MPI_SUM, comm2d, ierr )
[3347]	1219	#else
[4834]	1220	volume_flow = volume_flow_l
[3347]	1221	#endif
	1222
[4834]	1223	w_correct = SUM( volume_flow ) * d_area_t * drho_air_zw(nzt)
[3347]	1224
[4834]	1225	DO i = nxl, nxr
	1226	DO j = nys, nyn
	1227	DO k = nzt, nzt + 1
	1228	w(k,j,i) = w(k,j,i) + w_correct &
	1229	* MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j,i), 3 ) )
[3347]	1230	ENDDO
	1231	ENDDO
[4834]	1232	ENDDO
[3347]	1233
[4834]	1234	CALL cpu_log( log_point(58), 'offline nesting', 'stop' )
[3987]	1235
[4834]	1236	IF ( debug_output_timestep ) CALL debug_message( 'nesting_offl_mass_conservation', 'end' )
[3347]	1237
[4834]	1238	END SUBROUTINE nesting_offl_mass_conservation
[3347]	1239
[4834]	1240
	1241	!--------------------------------------------------------------------------------------------------!
[3347]	1242	! Description:
	1243	! ------------
[4834]	1244	!> Set the lateral and top boundary conditions in case the PALM domain is nested offline in a
	1245	!> mesoscale model. Further, average boundary data and determine mean profiles, further used for
	1246	!> correct damping in the sponge layer.
	1247	!--------------------------------------------------------------------------------------------------!
	1248	SUBROUTINE nesting_offl_bc
[3347]	1249
[4834]	1250	USE exchange_horiz_mod, &
	1251	ONLY: exchange_horiz
[4457]	1252
[4834]	1253	INTEGER(iwp) :: i !< running index x-direction
	1254	INTEGER(iwp) :: j !< running index y-direction
	1255	INTEGER(iwp) :: k !< running index z-direction
	1256	INTEGER(iwp) :: n !< running index for chemical species
[4724]	1257
[4834]	1258	REAL(wp), DIMENSION(nzb:nzt+1) :: pt_ref !< reference profile for potential temperature
	1259	REAL(wp), DIMENSION(nzb:nzt+1) :: pt_ref_l !< reference profile for potential temperature on subdomain
	1260	REAL(wp), DIMENSION(nzb:nzt+1) :: q_ref !< reference profile for mixing ratio
	1261	REAL(wp), DIMENSION(nzb:nzt+1) :: q_ref_l !< reference profile for mixing ratio on subdomain
	1262	REAL(wp), DIMENSION(nzb:nzt+1) :: u_ref !< reference profile for u-component
	1263	REAL(wp), DIMENSION(nzb:nzt+1) :: u_ref_l !< reference profile for u-component on subdomain
	1264	REAL(wp), DIMENSION(nzb:nzt+1) :: v_ref !< reference profile for v-component
	1265	REAL(wp), DIMENSION(nzb:nzt+1) :: v_ref_l !< reference profile for v-component on subdomain
	1266	REAL(wp), DIMENSION(nzb:nzt+1) :: var_1d !< pre-interpolated profile for LOD1 mode
[3885]	1267
[4834]	1268	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: ref_chem !< reference profile for chemical species
	1269	REAL(wp), DIMENSION(:,:), ALLOCATABLE :: ref_chem_l !< reference profile for chemical species on subdomain
[4230]	1270
[4834]	1271	IF ( debug_output_timestep ) CALL debug_message( 'nesting_offl_bc', 'start' )
[3987]	1272
[4834]	1273	CALL cpu_log( log_point(58), 'offline nesting', 'start' )
[3347]	1274	!
[4834]	1275	!-- Initialize mean profiles, derived from boundary data, to zero.
	1276	pt_ref = 0.0_wp
	1277	q_ref = 0.0_wp
	1278	u_ref = 0.0_wp
	1279	v_ref = 0.0_wp
[4273]	1280
[4834]	1281	pt_ref_l = 0.0_wp
	1282	q_ref_l = 0.0_wp
	1283	u_ref_l = 0.0_wp
	1284	v_ref_l = 0.0_wp
[3347]	1285	!
[4834]	1286	!-- If required, allocate temporary arrays to compute chemistry mean profiles
	1287	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	1288	ALLOCATE( ref_chem(nzb:nzt+1,1:UBOUND( chem_species, 1 ) ) )
	1289	ALLOCATE( ref_chem_l(nzb:nzt+1,1:UBOUND( chem_species, 1 ) ) )
	1290	ref_chem = 0.0_wp
	1291	ref_chem_l = 0.0_wp
	1292	ENDIF
[4230]	1293	!
[4834]	1294	!-- Set boundary conditions of u-, v-, w-component, as well as q, and pt.
	1295	!-- Note, boundary values at the left boundary: i=-1 (v,w,pt,q) and i=0 (u), at the right boundary:
	1296	!-- i=nxr+1 (all), at the south boundary: j=-1 (u,w,pt,q) and j=0 (v), at the north boundary:
	1297	!-- j=nyn+1 (all).
	1298	!-- Please note, at the left (for u) and south (for v) boundary, values for u and v are set also at
	1299	!-- i/j=-1, since these values are used in boundary_conditions() to restore prognostic values.
	1300	!-- Further, sum up data to calculate mean profiles from boundary data, used for Rayleigh damping.
	1301	IF ( bc_dirichlet_l ) THEN
[4724]	1302	!
[4834]	1303	!-- u-component
	1304	IF ( lod == 2 ) THEN
	1305	DO j = nys, nyn
	1306	DO k = nzb+1, nzt
	1307	u(k,j,i_bound_u) = interpolate_in_time( nest_offl%u_l(0,k,j), &
	1308	nest_offl%u_l(1,k,j), &
	1309	fac_dt ) * &
	1310	MERGE( 1.0_wp, 0.0_wp, &
	1311	BTEST( wall_flags_total_0(k,j,i_bound_u), 1 ) )
[4724]	1312	ENDDO
[4834]	1313	u(:,j,i_bound_u-1) = u(:,j,i_bound_u)
	1314	u_ref_l(nzb+1:nzt) = u_ref_l(nzb+1:nzt) + u(nzb+1:nzt,j,i_bound_u)
	1315	ENDDO
	1316	ELSE
[4724]	1317	!
[4834]	1318	!-- Pre-interpolate profile before mapping onto the boundaries.
	1319	DO k = nzb+1, nzt
	1320	var_1d(k) = interpolate_in_time( nest_offl%u_l(0,k,1), &
	1321	nest_offl%u_l(1,k,1), &
	1322	fac_dt )
	1323	ENDDO
	1324	DO j = nys, nyn
	1325	u(nzb+1:nzt,j,i_bound_u) = var_1d(nzb+1:nzt) * &
[4724]	1326	MERGE( 1.0_wp, 0.0_wp, &
	1327	BTEST( wall_flags_total_0(nzb+1:nzt,j,i_bound_u), 1 ) )
[4834]	1328	u(:,j,i_bound_u-1) = u(:,j,i_bound_u)
	1329	u_ref_l(nzb+1:nzt) = u_ref_l(nzb+1:nzt) + u(nzb+1:nzt,j,i_bound_u)
	1330	ENDDO
	1331	ENDIF
	1332	!
	1333	!-- w-component
	1334	IF ( lod == 2 ) THEN
	1335	DO j = nys, nyn
	1336	DO k = nzb+1, nzt-1
	1337	w(k,j,i_bound) = interpolate_in_time( nest_offl%w_l(0,k,j), &
	1338	nest_offl%w_l(1,k,j), &
	1339	fac_dt ) * &
	1340	MERGE( 1.0_wp, 0.0_wp, &
	1341	BTEST( wall_flags_total_0(k,j,i_bound), 3 ) )
[4724]	1342	ENDDO
[4834]	1343	w(nzt,j,i_bound) = w(nzt-1,j,i_bound)
	1344	ENDDO
	1345	ELSE
	1346	DO k = nzb+1, nzt-1
	1347	var_1d(k) = interpolate_in_time( nest_offl%w_l(0,k,1), &
	1348	nest_offl%w_l(1,k,1), &
	1349	fac_dt )
	1350	ENDDO
	1351	DO j = nys, nyn
	1352	w(nzb+1:nzt-1,j,i_bound) = var_1d(nzb+1:nzt-1) * &
	1353	MERGE( 1.0_wp, 0.0_wp, &
	1354	BTEST( wall_flags_total_0(nzb+1:nzt-1,j,i_bound), 3 ) )
	1355	w(nzt,j,i_bound) = w(nzt-1,j,i_bound)
	1356	ENDDO
	1357	ENDIF
[4724]	1358	!
[4834]	1359	!-- v-component
	1360	IF ( lod == 2 ) THEN
	1361	DO j = nysv, nyn
	1362	DO k = nzb+1, nzt
	1363	v(k,j,i_bound) = interpolate_in_time( nest_offl%v_l(0,k,j), &
	1364	nest_offl%v_l(1,k,j), &
	1365	fac_dt ) * &
	1366	MERGE( 1.0_wp, 0.0_wp, &
	1367	BTEST( wall_flags_total_0(k,j,i_bound), 2 ) )
	1368	ENDDO
	1369	v_ref_l(nzb+1:nzt) = v_ref_l(nzb+1:nzt) + v(nzb+1:nzt,j,i_bound)
	1370	ENDDO
	1371	ELSE
	1372	DO k = nzb+1, nzt
	1373	var_1d(k) = interpolate_in_time( nest_offl%v_l(0,k,1), &
	1374	nest_offl%v_l(1,k,1), &
	1375	fac_dt )
	1376	ENDDO
	1377	DO j = nysv, nyn
	1378	v(nzb+1:nzt,j,i_bound) = var_1d(nzb+1:nzt) * &
	1379	MERGE( 1.0_wp, 0.0_wp, &
	1380	BTEST( wall_flags_total_0(nzb+1:nzt,j,i_bound), 2 ) )
	1381	v_ref_l(nzb+1:nzt) = v_ref_l(nzb+1:nzt) + v(nzb+1:nzt,j,i_bound)
	1382	ENDDO
	1383	ENDIF
	1384	!
	1385	!-- Potential temperature
	1386	IF ( .NOT. neutral ) THEN
[4724]	1387	IF ( lod == 2 ) THEN
	1388	DO j = nys, nyn
[4834]	1389	DO k = nzb+1, nzt
	1390	pt(k,j,i_bound) = interpolate_in_time( nest_offl%pt_l(0,k,j), &
	1391	nest_offl%pt_l(1,k,j), &
	1392	fac_dt )
[4724]	1393	ENDDO
[4834]	1394	pt_ref_l(nzb+1:nzt) = pt_ref_l(nzb+1:nzt) + pt(nzb+1:nzt,j,i_bound)
[4724]	1395	ENDDO
	1396	ELSE
[4834]	1397	DO k = nzb+1, nzt
	1398	var_1d(k) = interpolate_in_time( nest_offl%pt_l(0,k,1), &
	1399	nest_offl%pt_l(1,k,1), &
[4724]	1400	fac_dt )
[3347]	1401	ENDDO
[4724]	1402	DO j = nys, nyn
[4834]	1403	pt(nzb+1:nzt,j,i_bound) = var_1d(nzb+1:nzt)
	1404	pt_ref_l(nzb+1:nzt) = pt_ref_l(nzb+1:nzt) + pt(nzb+1:nzt,j,i_bound)
[4724]	1405	ENDDO
	1406	ENDIF
[4834]	1407	ENDIF
[4724]	1408	!
[4834]	1409	!-- Humidity
	1410	IF ( humidity ) THEN
[4724]	1411	IF ( lod == 2 ) THEN
[4834]	1412	DO j = nys, nyn
[4724]	1413	DO k = nzb+1, nzt
[4834]	1414	q(k,j,i_bound) = interpolate_in_time( nest_offl%q_l(0,k,j), &
	1415	nest_offl%q_l(1,k,j), &
	1416	fac_dt )
[4724]	1417	ENDDO
[4834]	1418	q_ref_l(nzb+1:nzt) = q_ref_l(nzb+1:nzt) + q(nzb+1:nzt,j,i_bound)
[4724]	1419	ENDDO
	1420	ELSE
[3347]	1421	DO k = nzb+1, nzt
[4834]	1422	var_1d(k) = interpolate_in_time( nest_offl%q_l(0,k,1), &
	1423	nest_offl%q_l(1,k,1), &
[4724]	1424	fac_dt )
[3347]	1425	ENDDO
[4834]	1426	DO j = nys, nyn
	1427	q(nzb+1:nzt,j,i_bound) = var_1d(nzb+1:nzt)
	1428	q_ref_l(nzb+1:nzt) = q_ref_l(nzb+1:nzt) + q(nzb+1:nzt,j,i_bound)
[4724]	1429	ENDDO
	1430	ENDIF
[4834]	1431	ENDIF
[4724]	1432	!
[4834]	1433	!-- Chemistry
	1434	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	1435	DO n = 1, UBOUND( chem_species, 1 )
	1436	IF ( nest_offl%chem_from_file_l(n) ) THEN
	1437	IF ( lod == 2 ) THEN
	1438	DO j = nys, nyn
	1439	DO k = nzb+1, nzt
	1440	chem_species(n)%conc(k,j,i_bound) = interpolate_in_time( &
[4724]	1441	nest_offl%chem_l(0,k,j,n), &
	1442	nest_offl%chem_l(1,k,j,n), &
	1443	fac_dt )
	1444	ENDDO
[4834]	1445	ref_chem_l(nzb+1:nzt,n) = ref_chem_l(nzb+1:nzt,n) &
	1446	+ chem_species(n)%conc(nzb+1:nzt,j,i_bound)
	1447	ENDDO
	1448	ELSE
	1449	DO k = nzb+1, nzt
	1450	var_1d(k) = interpolate_in_time( nest_offl%chem_l(0,k,1,n), &
	1451	nest_offl%chem_l(1,k,1,n), &
	1452	fac_dt )
	1453	ENDDO
	1454	DO j = nys, nyn
	1455	chem_species(n)%conc(nzb+1:nzt,j,i_bound) = var_1d(nzb+1:nzt)
	1456	ref_chem_l(nzb+1:nzt,n) = ref_chem_l(nzb+1:nzt,n) &
	1457	+ chem_species(n)%conc(nzb+1:nzt,j,i_bound)
	1458	ENDDO
[3737]	1459	ENDIF
[4834]	1460	ENDIF
	1461	ENDDO
[3347]	1462	ENDIF
	1463
[4834]	1464	ENDIF
	1465
	1466	IF ( bc_dirichlet_r ) THEN
[4724]	1467	!
[4834]	1468	!-- u-component
	1469	IF ( lod == 2 ) THEN
	1470	DO j = nys, nyn
[3347]	1471	DO k = nzb+1, nzt
[4834]	1472	u(k,j,i_bound_u) = interpolate_in_time( nest_offl%u_r(0,k,j), &
	1473	nest_offl%u_r(1,k,j), &
	1474	fac_dt ) * &
	1475	MERGE( 1.0_wp, 0.0_wp, &
	1476	BTEST( wall_flags_total_0(k,j,i_bound_u), 1 ) )
[3347]	1477	ENDDO
[4834]	1478	u_ref_l(nzb+1:nzt) = u_ref_l(nzb+1:nzt) + u(nzb+1:nzt,j,i_bound_u)
	1479	ENDDO
	1480	ELSE
	1481	DO k = nzb+1, nzt
	1482	var_1d(k) = interpolate_in_time( nest_offl%u_r(0,k,1), &
	1483	nest_offl%u_r(1,k,1), &
	1484	fac_dt )
	1485	ENDDO
	1486	DO j = nys, nyn
	1487	u(nzb+1:nzt,j,i_bound_u) = var_1d(nzb+1:nzt) * &
[4724]	1488	MERGE( 1.0_wp, 0.0_wp, &
	1489	BTEST( wall_flags_total_0(nzb+1:nzt,j,i_bound_u), 1 ) )
[4834]	1490	u_ref_l(nzb+1:nzt) = u_ref_l(nzb+1:nzt) + u(nzb+1:nzt,j,i_bound_u)
	1491	ENDDO
	1492	ENDIF
	1493	!
	1494	!-- w-component
	1495	IF ( lod == 2 ) THEN
	1496	DO j = nys, nyn
	1497	DO k = nzb+1, nzt-1
	1498	w(k,j,i_bound) = interpolate_in_time( nest_offl%w_r(0,k,j), &
	1499	nest_offl%w_r(1,k,j), &
	1500	fac_dt ) * &
	1501	MERGE( 1.0_wp, 0.0_wp, &
	1502	BTEST( wall_flags_total_0(k,j,i_bound), 3 ) )
[4724]	1503	ENDDO
[4834]	1504	w(nzt,j,i_bound) = w(nzt-1,j,i_bound)
	1505	ENDDO
	1506	ELSE
	1507	DO k = nzb+1, nzt-1
	1508	var_1d(k) = interpolate_in_time( nest_offl%w_r(0,k,1), &
	1509	nest_offl%w_r(1,k,1), &
	1510	fac_dt )
	1511	ENDDO
	1512	DO j = nys, nyn
	1513	w(nzb+1:nzt-1,j,i_bound) = var_1d(nzb+1:nzt-1) * &
	1514	MERGE( 1.0_wp, 0.0_wp, &
	1515	BTEST( wall_flags_total_0(nzb+1:nzt-1,j,i_bound), 3 ) )
	1516	w(nzt,j,i_bound) = w(nzt-1,j,i_bound)
	1517	ENDDO
	1518	ENDIF
[4724]	1519	!
[4834]	1520	!-- v-component
	1521	IF ( lod == 2 ) THEN
	1522	DO j = nysv, nyn
	1523	DO k = nzb+1, nzt
	1524	v(k,j,i_bound) = interpolate_in_time( nest_offl%v_r(0,k,j), &
	1525	nest_offl%v_r(1,k,j), &
	1526	fac_dt ) * &
	1527	MERGE( 1.0_wp, 0.0_wp, &
	1528	BTEST( wall_flags_total_0(k,j,i_bound), 2 ) )
	1529	ENDDO
	1530	v_ref_l(nzb+1:nzt) = v_ref_l(nzb+1:nzt) + v(nzb+1:nzt,j,i_bound)
	1531	ENDDO
	1532	ELSE
	1533	DO k = nzb+1, nzt
	1534	var_1d(k) = interpolate_in_time( nest_offl%v_r(0,k,1), &
	1535	nest_offl%v_r(1,k,1), &
	1536	fac_dt )
	1537	ENDDO
	1538	DO j = nysv, nyn
	1539	v(nzb+1:nzt,j,i_bound) = var_1d(nzb+1:nzt) * &
	1540	MERGE( 1.0_wp, 0.0_wp, &
	1541	BTEST( wall_flags_total_0(nzb+1:nzt,j,i_bound), 2 ) )
	1542	v_ref_l(nzb+1:nzt) = v_ref_l(nzb+1:nzt) + v(nzb+1:nzt,j,i_bound)
	1543	ENDDO
	1544	ENDIF
	1545	!
	1546	!-- Potential temperature
	1547	IF ( .NOT. neutral ) THEN
[4724]	1548	IF ( lod == 2 ) THEN
	1549	DO j = nys, nyn
[4834]	1550	DO k = nzb+1, nzt
	1551	pt(k,j,i_bound) = interpolate_in_time( nest_offl%pt_r(0,k,j), &
	1552	nest_offl%pt_r(1,k,j), &
	1553	fac_dt )
[4724]	1554	ENDDO
[4834]	1555	pt_ref_l(nzb+1:nzt) = pt_ref_l(nzb+1:nzt) + pt(nzb+1:nzt,j,i_bound)
[4724]	1556	ENDDO
	1557	ELSE
[4834]	1558	DO k = nzb+1, nzt
	1559	var_1d(k) = interpolate_in_time( nest_offl%pt_r(0,k,1), &
	1560	nest_offl%pt_r(1,k,1), &
[4724]	1561	fac_dt )
[3347]	1562	ENDDO
[4724]	1563	DO j = nys, nyn
[4834]	1564	pt(nzb+1:nzt,j,i_bound) = var_1d(nzb+1:nzt)
	1565	pt_ref_l(nzb+1:nzt) = pt_ref_l(nzb+1:nzt) + pt(nzb+1:nzt,j,i_bound)
[4724]	1566	ENDDO
	1567	ENDIF
[4834]	1568	ENDIF
[4724]	1569	!
[4834]	1570	!-- Humidity
	1571	IF ( humidity ) THEN
[4724]	1572	IF ( lod == 2 ) THEN
[4834]	1573	DO j = nys, nyn
[4724]	1574	DO k = nzb+1, nzt
[4834]	1575	q(k,j,i_bound) = interpolate_in_time( nest_offl%q_r(0,k,j), &
	1576	nest_offl%q_r(1,k,j), &
	1577	fac_dt )
[4724]	1578	ENDDO
[4834]	1579	q_ref_l(nzb+1:nzt) = q_ref_l(nzb+1:nzt) + q(nzb+1:nzt,j,i_bound)
[4724]	1580	ENDDO
	1581	ELSE
[3347]	1582	DO k = nzb+1, nzt
[4834]	1583	var_1d(k) = interpolate_in_time( nest_offl%q_r(0,k,1), &
	1584	nest_offl%q_r(1,k,1), &
[4724]	1585	fac_dt )
[3347]	1586	ENDDO
[4834]	1587	DO j = nys, nyn
	1588	q(nzb+1:nzt,j,i_bound) = var_1d(nzb+1:nzt)
	1589	q_ref_l(nzb+1:nzt) = q_ref_l(nzb+1:nzt) + q(nzb+1:nzt,j,i_bound)
[4724]	1590	ENDDO
	1591	ENDIF
[4834]	1592	ENDIF
[4724]	1593	!
[4834]	1594	!-- Chemistry
	1595	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	1596	DO n = 1, UBOUND( chem_species, 1 )
	1597	IF ( nest_offl%chem_from_file_r(n) ) THEN
	1598	IF ( lod == 2 ) THEN
	1599	DO j = nys, nyn
	1600	DO k = nzb+1, nzt
	1601	chem_species(n)%conc(k,j,i_bound) = interpolate_in_time( &
[4724]	1602	nest_offl%chem_r(0,k,j,n), &
	1603	nest_offl%chem_r(1,k,j,n), &
	1604	fac_dt )
	1605	ENDDO
[4834]	1606	ref_chem_l(nzb+1:nzt,n) = ref_chem_l(nzb+1:nzt,n) &
	1607	+ chem_species(n)%conc(nzb+1:nzt,j,i_bound)
	1608	ENDDO
	1609	ELSE
	1610	DO k = nzb+1, nzt
	1611	var_1d(k) = interpolate_in_time( nest_offl%chem_r(0,k,1,n), &
	1612	nest_offl%chem_r(1,k,1,n), &
	1613	fac_dt )
	1614	ENDDO
	1615	DO j = nys, nyn
	1616	chem_species(n)%conc(nzb+1:nzt,j,i_bound) = var_1d(nzb+1:nzt)
	1617	ref_chem_l(nzb+1:nzt,n) = ref_chem_l(nzb+1:nzt,n) &
	1618	+ chem_species(n)%conc(nzb+1:nzt,j,i_bound)
	1619	ENDDO
[3737]	1620	ENDIF
[4834]	1621	ENDIF
	1622	ENDDO
[3347]	1623	ENDIF
	1624
[4834]	1625	ENDIF
	1626
	1627	IF ( bc_dirichlet_n ) THEN
[4724]	1628	!
[4834]	1629	!-- v-component
	1630	IF ( lod == 2 ) THEN
	1631	DO i = nxl, nxr
[3347]	1632	DO k = nzb+1, nzt
[4834]	1633	v(k,j_bound_v,i) = interpolate_in_time( nest_offl%v_n(0,k,i), &
	1634	nest_offl%v_n(1,k,i), &
	1635	fac_dt ) * &
	1636	MERGE( 1.0_wp, 0.0_wp, &
	1637	BTEST( wall_flags_total_0(k,j_bound_v,i), 2 ) )
[3347]	1638	ENDDO
[4834]	1639	v_ref_l(nzb+1:nzt) = v_ref_l(nzb+1:nzt) + v(nzb+1:nzt,j_bound_v,i)
	1640	ENDDO
	1641	ELSE
	1642	DO k = nzb+1, nzt
	1643	var_1d(k) = interpolate_in_time( nest_offl%v_n(0,k,1), &
	1644	nest_offl%v_n(1,k,1), &
	1645	fac_dt )
	1646	ENDDO
	1647	DO i = nxl, nxr
	1648	v(nzb+1:nzt,j_bound_v,i) = var_1d(nzb+1:nzt) * &
	1649	MERGE( 1.0_wp, 0.0_wp, &
	1650	BTEST( wall_flags_total_0(nzb+1:nzt,j_bound_v,i), 2 ) )
	1651	v_ref_l(nzb+1:nzt) = v_ref_l(nzb+1:nzt) + v(nzb+1:nzt,j_bound_v,i)
	1652	ENDDO
	1653	ENDIF
[4724]	1654	!
[4834]	1655	!-- w-component
	1656	IF ( lod == 2 ) THEN
	1657	DO i = nxl, nxr
[3347]	1658	DO k = nzb+1, nzt-1
[4834]	1659	w(k,j_bound,i) = interpolate_in_time( nest_offl%w_n(0,k,i), &
	1660	nest_offl%w_n(1,k,i), &
	1661	fac_dt ) * &
	1662	MERGE( 1.0_wp, 0.0_wp, &
	1663	BTEST( wall_flags_total_0(k,j_bound,i), 3 ) )
[3347]	1664	ENDDO
[4834]	1665	w(nzt,j_bound,i) = w(nzt-1,j_bound,i)
	1666	ENDDO
	1667	ELSE
	1668	DO k = nzb+1, nzt-1
	1669	var_1d(k) = interpolate_in_time( nest_offl%w_n(0,k,1), &
	1670	nest_offl%w_n(1,k,1), &
	1671	fac_dt )
	1672	ENDDO
	1673	DO i = nxl, nxr
	1674	w(nzb+1:nzt-1,j_bound,i) = var_1d(nzb+1:nzt-1) * &
	1675	MERGE( 1.0_wp, 0.0_wp, &
	1676	BTEST( wall_flags_total_0(nzb+1:nzt-1,j_bound,i), 3 ) )
	1677	w(nzt,j_bound,i) = w(nzt-1,j_bound,i)
	1678	ENDDO
	1679	ENDIF
[4724]	1680	!
[4834]	1681	!-- u-component
	1682	IF ( lod == 2 ) THEN
	1683	DO i = nxlu, nxr
[3347]	1684	DO k = nzb+1, nzt
[4834]	1685	u(k,j_bound,i) = interpolate_in_time( nest_offl%u_n(0,k,i), &
	1686	nest_offl%u_n(1,k,i), &
	1687	fac_dt ) * &
	1688	MERGE( 1.0_wp, 0.0_wp, &
	1689	BTEST( wall_flags_total_0(k,j_bound,i), 1 ) )
[3347]	1690	ENDDO
[4834]	1691	u_ref_l(nzb+1:nzt) = u_ref_l(nzb+1:nzt) + u(nzb+1:nzt,j_bound,i)
	1692	ENDDO
	1693	ELSE
	1694	DO k = nzb+1, nzt
	1695	var_1d(k) = interpolate_in_time( nest_offl%u_n(0,k,1), &
	1696	nest_offl%u_n(1,k,1), &
	1697	fac_dt )
	1698	ENDDO
	1699	DO i = nxlu, nxr
	1700	u(nzb+1:nzt,j_bound,i) = var_1d(nzb+1:nzt) * &
	1701	MERGE( 1.0_wp, 0.0_wp, &
	1702	BTEST( wall_flags_total_0(nzb+1:nzt,j_bound,i), 1 ) )
	1703	u_ref_l(nzb+1:nzt) = u_ref_l(nzb+1:nzt) + u(nzb+1:nzt,j_bound,i)
	1704	ENDDO
[3347]	1705	ENDIF
[4724]	1706	!
[4834]	1707	!-- Potential temperature
	1708	IF ( .NOT. neutral ) THEN
[4724]	1709	IF ( lod == 2 ) THEN
	1710	DO i = nxl, nxr
	1711	DO k = nzb+1, nzt
[4834]	1712	pt(k,j_bound,i) = interpolate_in_time( nest_offl%pt_n(0,k,i), &
	1713	nest_offl%pt_n(1,k,i), &
	1714	fac_dt )
[4724]	1715	ENDDO
[4834]	1716	pt_ref_l(nzb+1:nzt) = pt_ref_l(nzb+1:nzt) + pt(nzb+1:nzt,j_bound,i)
[4724]	1717	ENDDO
	1718	ELSE
[3347]	1719	DO k = nzb+1, nzt
[4834]	1720	var_1d(k) = interpolate_in_time( nest_offl%pt_n(0,k,1), &
	1721	nest_offl%pt_n(1,k,1), &
[4724]	1722	fac_dt )
[3347]	1723	ENDDO
[4724]	1724	DO i = nxl, nxr
[4834]	1725	pt(nzb+1:nzt,j_bound,i) = var_1d(nzb+1:nzt)
	1726	pt_ref_l(nzb+1:nzt) = pt_ref_l(nzb+1:nzt) + pt(nzb+1:nzt,j_bound,i)
[4724]	1727	ENDDO
	1728	ENDIF
[4834]	1729	ENDIF
[4724]	1730	!
[4834]	1731	!-- Humidity
	1732	IF ( humidity ) THEN
[4724]	1733	IF ( lod == 2 ) THEN
	1734	DO i = nxl, nxr
	1735	DO k = nzb+1, nzt
[4834]	1736	q(k,j_bound,i) = interpolate_in_time( nest_offl%q_n(0,k,i), &
	1737	nest_offl%q_n(1,k,i), &
	1738	fac_dt )
[4724]	1739	ENDDO
[4834]	1740	q_ref_l(nzb+1:nzt) = q_ref_l(nzb+1:nzt) + q(nzb+1:nzt,j_bound,i)
[4724]	1741	ENDDO
	1742	ELSE
[3347]	1743	DO k = nzb+1, nzt
[4834]	1744	var_1d(k) = interpolate_in_time( nest_offl%q_n(0,k,1), &
	1745	nest_offl%q_n(1,k,1), &
[4724]	1746	fac_dt )
[3347]	1747	ENDDO
[4834]	1748	DO i = nxl, nxr
	1749	q(nzb+1:nzt,j_bound,i) = var_1d(nzb+1:nzt)
	1750	q_ref_l(nzb+1:nzt) = q_ref_l(nzb+1:nzt) + q(nzb+1:nzt,j_bound,i)
[4724]	1751	ENDDO
	1752	ENDIF
[4834]	1753	ENDIF
[4724]	1754	!
[4834]	1755	!-- Chemistry
	1756	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	1757	DO n = 1, UBOUND( chem_species, 1 )
	1758	IF ( nest_offl%chem_from_file_n(n) ) THEN
	1759	IF ( lod == 2 ) THEN
	1760	DO i = nxl, nxr
	1761	DO k = nzb+1, nzt
	1762	chem_species(n)%conc(k,j_bound,i) = interpolate_in_time( &
	1763	nest_offl%chem_n(0,k,i,n), &
	1764	nest_offl%chem_n(1,k,i,n), &
	1765	fac_dt )
	1766	ENDDO
	1767	ref_chem_l(nzb+1:nzt,n) = ref_chem_l(nzb+1:nzt,n) &
	1768	+ chem_species(n)%conc(nzb+1:nzt,j_bound,i)
[4724]	1769	ENDDO
[4834]	1770	ELSE
[4724]	1771	DO k = nzb+1, nzt
[4834]	1772	var_1d(k) = interpolate_in_time( nest_offl%chem_n(0,k,1,n), &
	1773	nest_offl%chem_n(1,k,1,n), &
	1774	fac_dt )
[4724]	1775	ENDDO
[4834]	1776	DO i = nxl, nxr
	1777	chem_species(n)%conc(nzb+1:nzt,j_bound,i) = var_1d(nzb+1:nzt)
	1778	ref_chem_l(nzb+1:nzt,n) = ref_chem_l(nzb+1:nzt,n) + &
	1779	chem_species(n) &
	1780	%conc(nzb+1:nzt,j_bound,i)
	1781	ENDDO
	1782	ENDIF
[4724]	1783	ENDIF
[4834]	1784	ENDDO
[3347]	1785	ENDIF
[4834]	1786	ENDIF
	1787
	1788	IF ( bc_dirichlet_s ) THEN
[3347]	1789	!
[4834]	1790	!-- v-component
[4724]	1791	IF ( lod == 2 ) THEN
[4834]	1792	DO i = nxl, nxr
	1793	DO k = nzb+1, nzt
	1794	v(k,j_bound_v,i) = interpolate_in_time( nest_offl%v_s(0,k,i), &
	1795	nest_offl%v_s(1,k,i), &
	1796	fac_dt ) * &
	1797	MERGE( 1.0_wp, 0.0_wp, &
	1798	BTEST( wall_flags_total_0(k,j_bound_v,i), 2 ) )
[4724]	1799	ENDDO
[4834]	1800	v(:,j_bound_v-1,i) = v(:,j_bound_v,i)
	1801	v_ref_l(nzb+1:nzt) = v_ref_l(nzb+1:nzt) + v(nzb+1:nzt,j_bound_v,i)
[3347]	1802	ENDDO
[4724]	1803	ELSE
[4834]	1804	DO k = nzb+1, nzt
	1805	var_1d(k) = interpolate_in_time( nest_offl%v_s(0,k,1), &
	1806	nest_offl%v_s(1,k,1), &
	1807	fac_dt )
	1808	ENDDO
	1809	DO i = nxl, nxr
	1810	v(nzb+1:nzt,j_bound_v,i) = var_1d(nzb+1:nzt) * &
[4724]	1811	MERGE( 1.0_wp, 0.0_wp, &
[4834]	1812	BTEST( wall_flags_total_0(nzb+1:nzt,j_bound_v,i), 2 ) )
	1813	v(:,j_bound_v-1,i) = v(:,j_bound_v,i)
	1814	v_ref_l(nzb+1:nzt) = v_ref_l(nzb+1:nzt) + v(nzb+1:nzt,j_bound_v,i)
	1815	ENDDO
[4724]	1816	ENDIF
[3937]	1817	!
[4834]	1818	!-- w-component
[4724]	1819	IF ( lod == 2 ) THEN
	1820	DO i = nxl, nxr
[4834]	1821	DO k = nzb+1, nzt-1
	1822	w(k,j_bound,i) = interpolate_in_time( nest_offl%w_s(0,k,i), &
	1823	nest_offl%w_s(1,k,i), &
	1824	fac_dt ) * &
	1825	MERGE( 1.0_wp, 0.0_wp, &
	1826	BTEST( wall_flags_total_0(k,j_bound,i), 3 ) )
[4724]	1827	ENDDO
[4834]	1828	w(nzt,j_bound,i) = w(nzt-1,j_bound,i)
[3347]	1829	ENDDO
[4724]	1830	ELSE
[4834]	1831	DO k = nzb+1, nzt-1
	1832	var_1d(k) = interpolate_in_time( nest_offl%w_s(0,k,1), &
	1833	nest_offl%w_s(1,k,1), &
	1834	fac_dt )
	1835	ENDDO
	1836	DO i = nxl, nxr
	1837	w(nzb+1:nzt-1,j_bound,i) = var_1d(nzb+1:nzt-1) * &
[4724]	1838	MERGE( 1.0_wp, 0.0_wp, &
[4834]	1839	BTEST( wall_flags_total_0(nzb+1:nzt-1,j_bound,i), 3 ) )
	1840	w(nzt,j_bound,i) = w(nzt-1,j_bound,i)
	1841	ENDDO
[4724]	1842	ENDIF
[3937]	1843	!
[4834]	1844	!-- u-component
[4724]	1845	IF ( lod == 2 ) THEN
[4834]	1846	DO i = nxlu, nxr
	1847	DO k = nzb+1, nzt
	1848	u(k,j_bound,i) = interpolate_in_time( nest_offl%u_s(0,k,i), &
	1849	nest_offl%u_s(1,k,i), &
	1850	fac_dt ) * &
	1851	MERGE( 1.0_wp, 0.0_wp, &
	1852	BTEST( wall_flags_total_0(k,j_bound,i), 1 ) )
[3347]	1853	ENDDO
[4834]	1854	u_ref_l(nzb+1:nzt) = u_ref_l(nzb+1:nzt) + u(nzb+1:nzt,j_bound,i)
[3347]	1855	ENDDO
[4724]	1856	ELSE
[4834]	1857	DO k = nzb+1, nzt
	1858	var_1d(k) = interpolate_in_time( nest_offl%u_s(0,k,1), &
	1859	nest_offl%u_s(1,k,1), &
	1860	fac_dt )
	1861	ENDDO
	1862	DO i = nxlu, nxr
	1863	u(nzb+1:nzt,j_bound,i) = var_1d(nzb+1:nzt) * &
[4724]	1864	MERGE( 1.0_wp, 0.0_wp, &
[4834]	1865	BTEST( wall_flags_total_0(nzb+1:nzt,j_bound,i), 1 ) )
	1866	u_ref_l(nzb+1:nzt) = u_ref_l(nzb+1:nzt) + u(nzb+1:nzt,j_bound,i)
	1867	ENDDO
[3347]	1868	ENDIF
[4724]	1869	!
[4834]	1870	!-- Potential temperature
[4724]	1871	IF ( .NOT. neutral ) THEN
	1872	IF ( lod == 2 ) THEN
	1873	DO i = nxl, nxr
[4834]	1874	DO k = nzb+1, nzt
	1875	pt(k,j_bound,i) = interpolate_in_time( nest_offl%pt_s(0,k,i), &
	1876	nest_offl%pt_s(1,k,i), &
	1877	fac_dt )
[4724]	1878	ENDDO
[4834]	1879	pt_ref_l(nzb+1:nzt) = pt_ref_l(nzb+1:nzt) + pt(nzb+1:nzt,j_bound,i)
[4724]	1880	ENDDO
	1881	ELSE
[4834]	1882	DO k = nzb+1, nzt
	1883	var_1d(k) = interpolate_in_time( nest_offl%pt_s(0,k,1), &
	1884	nest_offl%pt_s(1,k,1), &
	1885	fac_dt )
	1886	ENDDO
	1887	DO i = nxl, nxr
	1888	pt(nzb+1:nzt,j_bound,i) = var_1d(nzb+1:nzt)
	1889	pt_ref_l(nzb+1:nzt) = pt_ref_l(nzb+1:nzt) + pt(nzb+1:nzt,j_bound,i)
	1890	ENDDO
[4724]	1891	ENDIF
	1892	ENDIF
	1893	!
[4834]	1894	!-- Humidity
[3347]	1895	IF ( humidity ) THEN
[4724]	1896	IF ( lod == 2 ) THEN
	1897	DO i = nxl, nxr
[4834]	1898	DO k = nzb+1, nzt
	1899	q(k,j_bound,i) = interpolate_in_time( nest_offl%q_s(0,k,i), &
	1900	nest_offl%q_s(1,k,i), &
	1901	fac_dt )
[4724]	1902	ENDDO
[4834]	1903	q_ref_l(nzb+1:nzt) = q_ref_l(nzb+1:nzt) + q(nzb+1:nzt,j_bound,i)
[3347]	1904	ENDDO
[4724]	1905	ELSE
[4834]	1906	DO k = nzb+1, nzt
	1907	var_1d(k) = interpolate_in_time( nest_offl%q_s(0,k,1), &
	1908	nest_offl%q_s(1,k,1), &
	1909	fac_dt )
	1910	ENDDO
	1911	DO i = nxl, nxr
	1912	q(nzb+1:nzt,j_bound,i) = var_1d(nzb+1:nzt)
	1913	q_ref_l(nzb+1:nzt) = q_ref_l(nzb+1:nzt) + q(nzb+1:nzt,j_bound,i)
	1914	ENDDO
[4724]	1915	ENDIF
[3347]	1916	ENDIF
[4834]	1917	!
	1918	!-- Chemistry
[4273]	1919	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
[3737]	1920	DO n = 1, UBOUND( chem_species, 1 )
[4834]	1921	IF ( nest_offl%chem_from_file_s(n) ) THEN
[4724]	1922	IF ( lod == 2 ) THEN
	1923	DO i = nxl, nxr
[4834]	1924	DO k = nzb+1, nzt
	1925	chem_species(n)%conc(k,j_bound,i) = interpolate_in_time( &
	1926	nest_offl%chem_s(0,k,i,n), &
	1927	nest_offl%chem_s(1,k,i,n), &
	1928	fac_dt )
[4724]	1929	ENDDO
[4834]	1930	ref_chem_l(nzb+1:nzt,n) = ref_chem_l(nzb+1:nzt,n) &
	1931	+ chem_species(n)%conc(nzb+1:nzt,j_bound,i)
[3737]	1932	ENDDO
[4724]	1933	ELSE
[4834]	1934	DO k = nzb+1, nzt
	1935	var_1d(k) = interpolate_in_time( nest_offl%chem_s(0,k,1,n), &
	1936	nest_offl%chem_s(1,k,1,n), &
	1937	fac_dt )
	1938	ENDDO
	1939	DO i = nxl, nxr
	1940	chem_species(n)%conc(nzb+1:nzt,j_bound,i) = var_1d(nzb+1:nzt)
	1941	ref_chem_l(nzb+1:nzt,n) = ref_chem_l(nzb+1:nzt,n) + &
	1942	chem_species(n) &
	1943	%conc(nzb+1:nzt,j_bound,i)
	1944	ENDDO
[4724]	1945	ENDIF
[3737]	1946	ENDIF
	1947	ENDDO
	1948	ENDIF
[4834]	1949	ENDIF
[3347]	1950	!
[4834]	1951	!-- Top boundary
	1952	!-- u-component
	1953	IF ( lod == 2 ) THEN
	1954	DO i = nxlu, nxr
	1955	DO j = nys, nyn
	1956	u(nzt+1,j,i) = interpolate_in_time( nest_offl%u_top(0,j,i), &
	1957	nest_offl%u_top(1,j,i), &
	1958	fac_dt ) * &
	1959	MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(nzt+1,j,i), 1 ) )
	1960	u_ref_l(nzt+1) = u_ref_l(nzt+1) + u(nzt+1,j,i)
	1961	ENDDO
	1962	ENDDO
	1963	ELSE
	1964	var_1d(nzt+1) = interpolate_in_time( nest_offl%u_top(0,1,1), &
	1965	nest_offl%u_top(1,1,1), &
	1966	fac_dt )
	1967	u(nzt+1,nys:nyn,nxlu:nxr) = var_1d(nzt+1) * &
	1968	MERGE( 1.0_wp, 0.0_wp, &
	1969	BTEST( wall_flags_total_0(nzt+1,nys:nyn,nxlu:nxr), 1 ) )
	1970	u_ref_l(nzt+1) = u_ref_l(nzt+1) + SUM( u(nzt+1,nys:nyn,nxlu:nxr) )
	1971	ENDIF
	1972	!
	1973	!-- For left boundary set boundary condition for u-component also at top grid point.
	1974	!-- Note, this has no effect on the numeric solution, only for data output.
	1975	IF ( bc_dirichlet_l ) u(nzt+1,:,nxl) = u(nzt+1,:,nxlu)
	1976	!
	1977	!-- v-component
	1978	IF ( lod == 2 ) THEN
	1979	DO i = nxl, nxr
	1980	DO j = nysv, nyn
	1981	v(nzt+1,j,i) = interpolate_in_time( nest_offl%v_top(0,j,i), &
	1982	nest_offl%v_top(1,j,i), &
	1983	fac_dt ) * &
	1984	MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(nzt+1,j,i), 2 ) )
	1985	v_ref_l(nzt+1) = v_ref_l(nzt+1) + v(nzt+1,j,i)
	1986	ENDDO
	1987	ENDDO
	1988	ELSE
	1989	var_1d(nzt+1) = interpolate_in_time( nest_offl%v_top(0,1,1), &
	1990	nest_offl%v_top(1,1,1), &
	1991	fac_dt )
	1992	v(nzt+1,nysv:nyn,nxl:nxr) = var_1d(nzt+1) * &
	1993	MERGE( 1.0_wp, 0.0_wp, &
	1994	BTEST( wall_flags_total_0(nzt+1,nysv:nyn,nxl:nxr), 2 ) )
	1995	v_ref_l(nzt+1) = v_ref_l(nzt+1) + SUM( v(nzt+1,nysv:nyn,nxl:nxr) )
	1996	ENDIF
	1997	!
	1998	!-- For south boundary set boundary condition for v-component also at top grid point.
	1999	!-- Note, this has no effect on the numeric solution, only for data output.
	2000	IF ( bc_dirichlet_s ) v(nzt+1,nys,:) = v(nzt+1,nysv,:)
	2001	!
	2002	!-- w-component
	2003	IF ( lod == 2 ) THEN
	2004	DO i = nxl, nxr
	2005	DO j = nys, nyn
	2006	w(nzt,j,i) = interpolate_in_time( nest_offl%w_top(0,j,i), &
	2007	nest_offl%w_top(1,j,i), &
	2008	fac_dt ) * &
	2009	MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(nzt,j,i), 3 ) )
	2010	w(nzt+1,j,i) = w(nzt,j,i)
	2011	ENDDO
	2012	ENDDO
	2013	ELSE
	2014	var_1d(nzt) = interpolate_in_time( nest_offl%w_top(0,1,1), &
	2015	nest_offl%w_top(1,1,1), &
	2016	fac_dt )
	2017	w(nzt,nys:nyn,nxl:nxr) = var_1d(nzt) * &
	2018	MERGE( 1.0_wp, 0.0_wp, &
	2019	BTEST( wall_flags_total_0(nzt,nys:nyn,nxl:nxr), 3 ) )
	2020	w(nzt+1,nys:nyn,nxl:nxr) = w(nzt,nys:nyn,nxl:nxr)
	2021	ENDIF
	2022	!
	2023	!-- Potential temperture
	2024	IF ( .NOT. neutral ) THEN
	2025	IF ( lod == 2 ) THEN
	2026	DO i = nxl, nxr
	2027	DO j = nys, nyn
	2028	pt(nzt+1,j,i) = interpolate_in_time( nest_offl%pt_top(0,j,i), &
	2029	nest_offl%pt_top(1,j,i), &
	2030	fac_dt )
	2031	pt_ref_l(nzt+1) = pt_ref_l(nzt+1) + pt(nzt+1,j,i)
	2032	ENDDO
	2033	ENDDO
	2034	ELSE
	2035	var_1d(nzt+1) = interpolate_in_time( nest_offl%pt_top(0,1,1), &
	2036	nest_offl%pt_top(1,1,1), &
	2037	fac_dt )
	2038	pt(nzt+1,nys:nyn,nxl:nxr) = var_1d(nzt+1)
	2039	pt_ref_l(nzt+1) = pt_ref_l(nzt+1) + SUM( pt(nzt+1,nys:nyn,nxl:nxr) )
[3347]	2040	ENDIF
[4834]	2041	ENDIF
	2042	!
	2043	!-- humidity
	2044	IF ( humidity ) THEN
	2045	IF ( lod == 2 ) THEN
	2046	DO i = nxl, nxr
	2047	DO j = nys, nyn
	2048	q(nzt+1,j,i) = interpolate_in_time( nest_offl%q_top(0,j,i), &
	2049	nest_offl%q_top(1,j,i), &
	2050	fac_dt )
	2051	q_ref_l(nzt+1) = q_ref_l(nzt+1) + q(nzt+1,j,i)
	2052	ENDDO
	2053	ENDDO
	2054	ELSE
	2055	var_1d(nzt+1) = interpolate_in_time( nest_offl%q_top(0,1,1), &
	2056	nest_offl%q_top(1,1,1), &
	2057	fac_dt )
	2058	q(nzt+1,nys:nyn,nxl:nxr) = var_1d(nzt+1)
	2059	q_ref_l(nzt+1) = q_ref_l(nzt+1) + SUM( q(nzt+1,nys:nyn,nxl:nxr) )
	2060	ENDIF
	2061	ENDIF
	2062
	2063	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	2064	DO n = 1, UBOUND( chem_species, 1 )
	2065	IF ( nest_offl%chem_from_file_t(n) ) THEN
	2066	IF ( lod == 2 ) THEN
	2067	DO i = nxl, nxr
	2068	DO j = nys, nyn
	2069	chem_species(n)%conc(nzt+1,j,i) = interpolate_in_time( &
	2070	nest_offl%chem_top(0,j,i,n), &
	2071	nest_offl%chem_top(1,j,i,n), &
	2072	fac_dt )
	2073	ref_chem_l(nzt+1,n) = ref_chem_l(nzt+1,n) + chem_species(n)%conc(nzt+1,j,i)
	2074	ENDDO
	2075	ENDDO
	2076	ELSE
	2077	var_1d(nzt+1) = interpolate_in_time( nest_offl%chem_top(0,1,1,n), &
	2078	nest_offl%chem_top(1,1,1,n), &
	2079	fac_dt )
	2080	chem_species(n)%conc(nzt+1,nys:nyn,nxl:nxr) = var_1d(nzt+1)
	2081	ref_chem_l(nzt+1,n) = ref_chem_l(nzt+1,n) + &
	2082	SUM( chem_species(n)%conc(nzt+1,nys:nyn,nxl:nxr) )
	2083	ENDIF
	2084	ENDIF
	2085	ENDDO
	2086	ENDIF
	2087	!
	2088	!-- Moreover, set Neumann boundary condition for subgrid-scale TKE, passive scalar, dissipation, and
	2089	!-- chemical species if required.
	2090	IF ( rans_mode .AND. rans_tke_e ) THEN
	2091	IF ( bc_dirichlet_l ) diss(:,:,nxl-1) = diss(:,:,nxl)
	2092	IF ( bc_dirichlet_r ) diss(:,:,nxr+1) = diss(:,:,nxr)
	2093	IF ( bc_dirichlet_s ) diss(:,nys-1,:) = diss(:,nys,:)
	2094	IF ( bc_dirichlet_n ) diss(:,nyn+1,:) = diss(:,nyn,:)
	2095	ENDIF
[4079]	2096	! IF ( .NOT. constant_diffusion ) THEN
	2097	! IF ( bc_dirichlet_l ) e(:,:,nxl-1) = e(:,:,nxl)
	2098	! IF ( bc_dirichlet_r ) e(:,:,nxr+1) = e(:,:,nxr)
	2099	! IF ( bc_dirichlet_s ) e(:,nys-1,:) = e(:,nys,:)
	2100	! IF ( bc_dirichlet_n ) e(:,nyn+1,:) = e(:,nyn,:)
	2101	! e(nzt+1,:,:) = e(nzt,:,:)
	2102	! ENDIF
	2103	! IF ( passive_scalar ) THEN
	2104	! IF ( bc_dirichlet_l ) s(:,:,nxl-1) = s(:,:,nxl)
	2105	! IF ( bc_dirichlet_r ) s(:,:,nxr+1) = s(:,:,nxr)
	2106	! IF ( bc_dirichlet_s ) s(:,nys-1,:) = s(:,nys,:)
	2107	! IF ( bc_dirichlet_n ) s(:,nyn+1,:) = s(:,nyn,:)
	2108	! ENDIF
[3347]	2109
[4834]	2110	CALL exchange_horiz( u, nbgp )
	2111	CALL exchange_horiz( v, nbgp )
	2112	CALL exchange_horiz( w, nbgp )
	2113	IF ( .NOT. neutral ) CALL exchange_horiz( pt, nbgp )
	2114	IF ( humidity ) CALL exchange_horiz( q, nbgp )
	2115	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	2116	DO n = 1, UBOUND( chem_species, 1 )
[3858]	2117	!
[4834]	2118	!-- Do local exchange only when necessary, i.e. when data is coming from dynamic file.
	2119	IF ( nest_offl%chem_from_file_t(n) ) CALL exchange_horiz( chem_species(n)%conc, nbgp )
	2120	ENDDO
	2121	ENDIF
[3347]	2122	!
[4834]	2123	!-- Set top boundary condition at all horizontal grid points, also at the lateral boundary grid
	2124	!-- points.
	2125	w(nzt+1,:,:) = w(nzt,:,:)
[4079]	2126	!
[4834]	2127	!-- Offline nesting for salsa
	2128	IF ( salsa ) CALL salsa_nesting_offl_bc
[4270]	2129	!
[4834]	2130	!-- Calculate the mean profiles. These are later stored on u_init, v_init, etc., in order to adjust
	2131	!-- the Rayleigh damping under time-evolving atmospheric conditions accordingly - damping against
	2132	!-- the representative mean profiles, not against the initial profiles. Note, in LOD = 1 case no
	2133	!-- averaging is required.
[3347]	2134	#if defined( __parallel )
[4834]	2135	CALL MPI_ALLREDUCE( u_ref_l, u_ref, nzt+1-nzb+1, MPI_REAL, MPI_SUM, comm2d, ierr )
	2136	CALL MPI_ALLREDUCE( v_ref_l, v_ref, nzt+1-nzb+1, MPI_REAL, MPI_SUM, comm2d, ierr )
	2137	IF ( humidity ) THEN
	2138	CALL MPI_ALLREDUCE( q_ref_l, q_ref, nzt+1-nzb+1, MPI_REAL, MPI_SUM, comm2d, ierr )
	2139	ENDIF
	2140	IF ( .NOT. neutral ) THEN
	2141	CALL MPI_ALLREDUCE( pt_ref_l, pt_ref, nzt+1-nzb+1, MPI_REAL, MPI_SUM, comm2d, ierr )
	2142	ENDIF
	2143	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	2144	CALL MPI_ALLREDUCE( ref_chem_l, ref_chem, ( nzt+1-nzb+1 ) * SIZE( ref_chem(nzb,:) ), &
	2145	MPI_REAL, MPI_SUM, comm2d, ierr )
	2146	ENDIF
[3347]	2147	#else
[4834]	2148	u_ref = u_ref_l
	2149	v_ref = v_ref_l
	2150	IF ( humidity ) q_ref = q_ref_l
	2151	IF ( .NOT. neutral ) pt_ref = pt_ref_l
	2152	IF ( air_chemistry .AND. nesting_offline_chem ) ref_chem = ref_chem_l
[3347]	2153	#endif
	2154	!
[4834]	2155	!-- Average data. Note, reference profiles up to nzt are derived from lateral boundaries, at the
	2156	!-- model top it is derived from the top boundary. Thus, number of input data is different from
	2157	!-- nzb:nzt compared to nzt+1.
	2158	!-- Derived from lateral boundaries.
	2159	u_ref(nzb:nzt) = u_ref(nzb:nzt) / REAL( 2.0_wp * ( ny + 1 + nx ), KIND = wp )
	2160	v_ref(nzb:nzt) = v_ref(nzb:nzt) / REAL( 2.0_wp * ( ny + nx + 1 ), KIND = wp )
	2161	IF ( humidity ) &
	2162	q_ref(nzb:nzt) = q_ref(nzb:nzt) / REAL( 2.0_wp * ( ny + 1 + nx + 1 ), KIND = wp )
	2163	IF ( .NOT. neutral ) &
	2164	pt_ref(nzb:nzt) = pt_ref(nzb:nzt) / REAL( 2.0_wp * ( ny + 1 + nx + 1 ), KIND = wp )
	2165	IF ( air_chemistry .AND. nesting_offline_chem ) &
	2166	ref_chem(nzb:nzt,:) = ref_chem(nzb:nzt,:) / REAL( 2.0_wp * ( ny + 1 + nx + 1 ), KIND = wp )
[3347]	2167	!
[4834]	2168	!-- Derived from top boundary.
	2169	u_ref(nzt+1) = u_ref(nzt+1) / REAL( ( ny + 1 ) * ( nx ), KIND = wp )
	2170	v_ref(nzt+1) = v_ref(nzt+1) / REAL( ( ny ) * ( nx + 1 ), KIND = wp )
	2171	IF ( humidity ) &
	2172	q_ref(nzt+1) = q_ref(nzt+1) / REAL( ( ny + 1 ) * ( nx + 1 ), KIND = wp )
	2173	IF ( .NOT. neutral ) &
	2174	pt_ref(nzt+1) = pt_ref(nzt+1) / REAL( ( ny + 1 ) * ( nx + 1 ), KIND = wp )
	2175	IF ( air_chemistry .AND. nesting_offline_chem ) &
	2176	ref_chem(nzt+1,:) = ref_chem(nzt+1,:) / REAL( ( ny + 1 ) * ( nx + 1 ),KIND = wp )
[3347]	2177	!
[4834]	2178	!-- Write onto init profiles, which are used for damping. Also set lower boundary condition for
	2179	!-- scalars (not required for u and v as these are zero at k=nzb).
	2180	u_init = u_ref
	2181	v_init = v_ref
	2182	IF ( humidity ) THEN
	2183	q_init = q_ref
	2184	q_init(nzb) = q_init(nzb+1)
	2185	ENDIF
	2186	IF ( .NOT. neutral ) THEN
	2187	pt_init = pt_ref
	2188	pt_init(nzb) = pt_init(nzb+1)
	2189	ENDIF
[4230]	2190
[4834]	2191	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	2192	DO n = 1, UBOUND( chem_species, 1 )
	2193	IF ( nest_offl%chem_from_file_t(n) ) THEN
	2194	chem_species(n)%conc_pr_init(:) = ref_chem(:,n)
	2195	chem_species(n)%conc_pr_init(nzb) = chem_species(n)%conc_pr_init(nzb+1)
	2196	ENDIF
	2197	ENDDO
	2198	ENDIF
	2199	IF ( ALLOCATED( ref_chem ) ) DEALLOCATE( ref_chem )
	2200	IF ( ALLOCATED( ref_chem_l ) ) DEALLOCATE( ref_chem_l )
[3347]	2201	!
[4834]	2202	!-- Further, adjust Rayleigh damping height in case of time-changing conditions.
	2203	!-- Therefore, calculate boundary-layer depth first.
	2204	CALL nesting_offl_calc_zi
	2205	CALL adjust_sponge_layer
[4273]	2206
[4834]	2207	CALL cpu_log( log_point(58), 'offline nesting', 'stop' )
[4226]	2208
[4834]	2209	IF ( debug_output_timestep ) CALL debug_message( 'nesting_offl_bc', 'end' )
[4226]	2210
	2211
[4834]	2212	END SUBROUTINE nesting_offl_bc
[4226]	2213
[4834]	2214	!--------------------------------------------------------------------------------------------------!
[4226]	2215	! Description:
[4834]	2216	!--------------------------------------------------------------------------------------------------!
	2217	!> Update of the geostrophic wind components. Note, currently this routine is not invoked.
	2218	!--------------------------------------------------------------------------------------------------!
	2219	SUBROUTINE nesting_offl_geostrophic_wind
[4226]	2220
[4834]	2221	INTEGER(iwp) :: k
[3347]	2222	!
[4834]	2223	!-- Update geostrophic wind components from dynamic input file.
	2224	DO k = nzb+1, nzt
	2225	ug(k) = interpolate_in_time( nest_offl%ug(0,k), nest_offl%ug(1,k), fac_dt )
	2226	vg(k) = interpolate_in_time( nest_offl%vg(0,k), nest_offl%vg(1,k), fac_dt )
	2227	ENDDO
	2228	ug(nzt+1) = ug(nzt)
	2229	vg(nzt+1) = vg(nzt)
[3347]	2230
[4834]	2231	END SUBROUTINE nesting_offl_geostrophic_wind
[3987]	2232
[4834]	2233	!--------------------------------------------------------------------------------------------------!
[4226]	2234	! Description:
[4834]	2235	!--------------------------------------------------------------------------------------------------!
	2236	!> Determine the interpolation constant for time interpolation. The calculation is separated from
	2237	!> the nesting_offl_bc and nesting_offl_geostrophic_wind in order to be independent on the order of
	2238	!> calls.
	2239	!--------------------------------------------------------------------------------------------------!
	2240	SUBROUTINE nesting_offl_interpolation_factor
[4226]	2241	!
[4834]	2242	!-- Determine interpolation factor and limit it to 1. This is because t+dt can slightly exceed
	2243	!-- time(tind_p) before boundary data is updated again.
	2244	fac_dt = ( time_since_reference_point - nest_offl%time(nest_offl%tind) + dt_3d ) / &
	2245	( nest_offl%time(nest_offl%tind_p) - nest_offl%time(nest_offl%tind) )
[3987]	2246
[4834]	2247	fac_dt = MIN( 1.0_wp, fac_dt )
[3347]	2248
[4834]	2249	END SUBROUTINE nesting_offl_interpolation_factor
[4226]	2250
[4834]	2251	!--------------------------------------------------------------------------------------------------!
[3347]	2252	! Description:
[4834]	2253	!--------------------------------------------------------------------------------------------------!
	2254	!> Calculates the boundary-layer depth from the boundary data, according to bulk-Richardson
	2255	!> criterion.
	2256	!--------------------------------------------------------------------------------------------------!
	2257	SUBROUTINE nesting_offl_calc_zi
[3347]	2258
[4834]	2259	INTEGER(iwp) :: i !< loop index in x-direction
	2260	INTEGER(iwp) :: j !< loop index in y-direction
	2261	INTEGER(iwp) :: k !< loop index in z-direction
	2262	INTEGER(iwp) :: k_max_loc !< index of maximum wind speed along z-direction
	2263	INTEGER(iwp) :: k_surface !< topography top index in z-direction
	2264	INTEGER(iwp) :: num_boundary_gp_non_cyclic !< number of non-cyclic boundaries, used for averaging ABL depth
	2265	INTEGER(iwp) :: num_boundary_gp_non_cyclic_l !< number of non-cyclic boundaries, used for averaging ABL depth
[3347]	2266
[4834]	2267	REAL(wp) :: ri_bulk !< bulk Richardson number
	2268	REAL(wp) :: ri_bulk_crit = 0.25_wp !< critical bulk Richardson number
	2269	REAL(wp) :: topo_max !< maximum topography level in model domain
	2270	REAL(wp) :: topo_max_l !< maximum topography level in subdomain
	2271	REAL(wp) :: vpt_surface !< near-surface virtual potential temperature
	2272	REAL(wp) :: zi_l !< mean boundary-layer depth on subdomain
	2273	REAL(wp) :: zi_local !< local boundary-layer depth
	2274
	2275	REAL(wp), DIMENSION(nzb:nzt+1) :: vpt_col !< vertical profile of virtual potential temperature at (j,i)-grid point
	2276	REAL(wp), DIMENSION(nzb:nzt+1) :: uv_abs !< vertical profile of horizontal wind speed at (j,i)-grid point
	2277
	2278
[3347]	2279	!
[4834]	2280	!-- Calculate mean boundary-layer height from boundary data.
	2281	!-- Start with the left and right boundaries.
	2282	zi_l = 0.0_wp
	2283	num_boundary_gp_non_cyclic_l = 0
	2284	IF ( bc_dirichlet_l .OR. bc_dirichlet_r ) THEN
[3347]	2285	!
[4834]	2286	!-- Sum-up and store number of boundary grid points used for averaging ABL depth
	2287	num_boundary_gp_non_cyclic_l = num_boundary_gp_non_cyclic_l + nxr - nxl + 1
[4022]	2288	!
[4834]	2289	!-- Determine index along x. Please note, index indicates boundary grid point for scalars.
	2290	i = MERGE( -1, nxr + 1, bc_dirichlet_l )
	2291
	2292	DO j = nys, nyn
[3347]	2293	!
[4834]	2294	!-- Determine topography top index at current (j,i) index
	2295	k_surface = topo_top_ind(j,i,0)
[3347]	2296	!
[4834]	2297	!-- Pre-compute surface virtual temperature. Therefore, use 2nd prognostic level according to
	2298	!-- Heinze et al. (2017).
	2299	IF ( humidity ) THEN
	2300	vpt_surface = pt(k_surface+2,j,i) * ( 1.0_wp + 0.61_wp * q(k_surface+2,j,i) )
	2301	vpt_col = pt(:,j,i) * ( 1.0_wp + 0.61_wp * q(:,j,i) )
	2302	ELSE
	2303	vpt_surface = pt(k_surface+2,j,i)
	2304	vpt_col = pt(:,j,i)
	2305	ENDIF
[3347]	2306	!
[4834]	2307	!-- Calculate local boundary layer height from bulk Richardson number, i.e. the height where
	2308	!-- the bulk Richardson number exceeds its critical value of 0.25
	2309	!-- (according to Heinze et al., 2017).
	2310	!-- Note, no interpolation of u- and v-component is made, as both are mainly mean inflow
	2311	!-- profiles with very small spatial variation.
	2312	!-- Add a safety factor in case the velocity term becomes zero. This may happen if overhanging
	2313	!-- 3D structures are directly located at the boundary, where velocity inside the building is
	2314	!-- zero (k_surface is the index of the lowest upward-facing surface).
	2315	uv_abs(:) = SQRT( MERGE( u(:,j,i+1), u(:,j,i), bc_dirichlet_l )2 + v(:,j,i)2 )
[4022]	2316	!
[4834]	2317	!-- Determine index of the maximum wind speed
	2318	k_max_loc = MAXLOC( uv_abs(:), DIM = 1 ) - 1
[4022]	2319
[4834]	2320	zi_local = 0.0_wp
	2321	DO k = k_surface+1, nzt
	2322	ri_bulk = zu(k) * g / vpt_surface * &
	2323	( vpt_col(k) - vpt_surface ) / ( uv_abs(k) + 1E-5_wp )
[4022]	2324	!
[4834]	2325	!-- Check if critical Richardson number is exceeded. Further, check if there is a maxium in
	2326	!-- the wind profile in order to detect also ABL heights in the stable boundary layer.
	2327	IF ( zi_local == 0.0_wp .AND. ( ri_bulk > ri_bulk_crit .OR. k == k_max_loc ) ) &
	2328	zi_local = zu(k)
[3347]	2329	ENDDO
	2330	!
[4834]	2331	!-- Assure that the minimum local boundary-layer depth is at least at the second vertical grid
	2332	!-- level.
	2333	zi_l = zi_l + MAX( zi_local, zu(k_surface+2) )
	2334
	2335	ENDDO
	2336
	2337	ENDIF
[4022]	2338	!
[4834]	2339	!-- Do the same at the north and south boundaries.
	2340	IF ( bc_dirichlet_s .OR. bc_dirichlet_n ) THEN
	2341
	2342	num_boundary_gp_non_cyclic_l = num_boundary_gp_non_cyclic_l + nxr - nxl + 1
	2343
	2344	j = MERGE( -1, nyn + 1, bc_dirichlet_s )
	2345
	2346	DO i = nxl, nxr
	2347	k_surface = topo_top_ind(j,i,0)
	2348
	2349	IF ( humidity ) THEN
	2350	vpt_surface = pt(k_surface+2,j,i) * ( 1.0_wp + 0.61_wp * q(k_surface+2,j,i) )
	2351	vpt_col = pt(:,j,i) * ( 1.0_wp + 0.61_wp * q(:,j,i) )
	2352	ELSE
	2353	vpt_surface = pt(k_surface+2,j,i)
	2354	vpt_col = pt(:,j,i)
	2355	ENDIF
	2356
	2357	uv_abs(:) = SQRT( u(:,j,i)2 + MERGE( v(:,j+1,i), v(:,j,i), bc_dirichlet_s )2 )
[4022]	2358	!
[4834]	2359	!-- Determine index of the maximum wind speed
	2360	k_max_loc = MAXLOC( uv_abs(:), DIM = 1 ) - 1
	2361
	2362	zi_local = 0.0_wp
	2363	DO k = k_surface+1, nzt
	2364	ri_bulk = zu(k) * g / vpt_surface * &
	2365	( vpt_col(k) - vpt_surface ) / ( uv_abs(k) + 1E-5_wp )
	2366	!
	2367	!-- Check if critical Richardson number is exceeded. Further, check if there is a maxium in
	2368	!-- the wind profile in order to detect also ABL heights in the stable boundary layer.
	2369	IF ( zi_local == 0.0_wp .AND. ( ri_bulk > ri_bulk_crit .OR. k == k_max_loc ) ) &
	2370	zi_local = zu(k)
[3347]	2371	ENDDO
[4834]	2372	zi_l = zi_l + MAX( zi_local, zu(k_surface+2) )
	2373
	2374	ENDDO
	2375
	2376	ENDIF
	2377
[3347]	2378	#if defined( __parallel )
[4834]	2379	CALL MPI_ALLREDUCE( zi_l, zi_ribulk, 1, MPI_REAL, MPI_SUM, comm2d, ierr )
	2380	CALL MPI_ALLREDUCE( num_boundary_gp_non_cyclic_l, num_boundary_gp_non_cyclic, &
	2381	1, MPI_INTEGER, MPI_SUM, comm2d, ierr )
[3347]	2382	#else
[4834]	2383	zi_ribulk = zi_l
	2384	num_boundary_gp_non_cyclic = num_boundary_gp_non_cyclic_l
[3347]	2385	#endif
[4834]	2386	zi_ribulk = zi_ribulk / REAL( num_boundary_gp_non_cyclic, KIND = wp )
[3347]	2387	!
[4834]	2388	!-- Finally, check if boundary layer depth is not below the any topography.
	2389	!-- zi_ribulk will be used to adjust rayleigh damping height, i.e. the lower level of the sponge
	2390	!-- layer, as well as to adjust the synthetic turbulence generator accordingly. If Rayleigh damping
	2391	!-- would be applied near buildings, etc., this would spoil the simulation results.
	2392	topo_max_l = zw(MAXVAL( topo_top_ind(nys:nyn,nxl:nxr,0) ))
	2393
[3347]	2394	#if defined( __parallel )
[4834]	2395	CALL MPI_ALLREDUCE( topo_max_l, topo_max, 1, MPI_REAL, MPI_MAX, comm2d, ierr )
[3347]	2396	#else
[4834]	2397	topo_max = topo_max_l
[3347]	2398	#endif
[4022]	2399	! zi_ribulk = MAX( zi_ribulk, topo_max )
[4834]	2400
	2401	END SUBROUTINE nesting_offl_calc_zi
	2402
	2403
	2404	!--------------------------------------------------------------------------------------------------!
[3347]	2405	! Description:
[4834]	2406	!--------------------------------------------------------------------------------------------------!
	2407	!> Adjust the height where the rayleigh damping starts, i.e. the lower level of the sponge layer.
	2408	!--------------------------------------------------------------------------------------------------!
	2409	SUBROUTINE adjust_sponge_layer
[3347]	2410
[4834]	2411	INTEGER(iwp) :: k !< loop index in z-direction
	2412
	2413	REAL(wp) :: rdh !< updated Rayleigh damping height
	2414
	2415
	2416	IF ( rayleigh_damping_height > 0.0_wp .AND. rayleigh_damping_factor > 0.0_wp ) THEN
[3347]	2417	!
[4834]	2418	!-- Update Rayleigh-damping height and re-calculate height-depending damping coefficients.
	2419	!-- Assure that rayleigh damping starts well above the boundary layer.
	2420	rdh = MIN( MAX( zi_ribulk * 1.3_wp, 10.0_wp * dz(1) ), &
	2421	0.8_wp * zu(nzt), rayleigh_damping_height )
[3347]	2422	!
	2423	!-- Update Rayleigh damping factor
[4834]	2424	DO k = nzb+1, nzt
	2425	IF ( zu(k) >= rdh ) THEN
	2426	rdf(k) = rayleigh_damping_factor * &
	2427	( SIN( pi * 0.5_wp * ( zu(k) - rdh ) / ( zu(nzt) - rdh ) ) )**2
	2428	ENDIF
	2429	ENDDO
	2430	rdf_sc = rdf
[3347]	2431
[4834]	2432	ENDIF
	2433
	2434	END SUBROUTINE adjust_sponge_layer
	2435
	2436	!--------------------------------------------------------------------------------------------------!
[3347]	2437	! Description:
	2438	! ------------
	2439	!> Performs consistency checks
[4834]	2440	!--------------------------------------------------------------------------------------------------!
	2441	SUBROUTINE nesting_offl_check_parameters
[3347]	2442	!
[4834]	2443	!-- Check if offline nesting is applied in nested child domain.
	2444	IF ( nesting_offline .AND. child_domain ) THEN
	2445	message_string = 'Offline nesting is only applicable in root model.'
	2446	CALL message( 'offline_nesting_check_parameters', 'PA0622', 1, 2, 0, 6, 0 )
	2447	ENDIF
[3347]	2448
[4834]	2449	END SUBROUTINE nesting_offl_check_parameters
	2450
	2451	!--------------------------------------------------------------------------------------------------!
[3347]	2452	! Description:
	2453	! ------------
	2454	!> Reads the parameter list nesting_offl_parameters
[4834]	2455	!--------------------------------------------------------------------------------------------------!
	2456	SUBROUTINE nesting_offl_parin
[3347]	2457
[4842]	2458	CHARACTER(LEN=100) :: line !< dummy string that contains the current line of the parameter file
[4843]	2459
[4842]	2460	INTEGER(iwp) :: io_status !< status after reading the namelist file
[3347]	2461
[4843]	2462	LOGICAL :: switch_off_module = .FALSE. !< local namelist parameter to switch off the module
	2463	!< although the respective module namelist appears in
	2464	!< the namelist file
[3347]	2465
[4843]	2466	NAMELIST /nesting_offl_parameters/ switch_off_module
[3347]	2467
[4843]	2468
[3347]	2469	!
[4842]	2470	!-- Move to the beginning of the namelist file and try to find and read the namelist.
	2471	REWIND( 11 )
	2472	READ( 11, nesting_offl_parameters, IOSTAT=io_status )
[3347]	2473
	2474	!
[4842]	2475	!-- Action depending on the READ status
[4843]	2476	IF ( io_status == 0 ) THEN
[4842]	2477	!
[4843]	2478	!-- nesting_offl_parameters namelist was found and read correctly. Enable the
	2479	!-- offline nesting.
	2480	IF ( .NOT. switch_off_module ) nesting_offline = .TRUE.
	2481
	2482	ELSEIF ( io_status > 0 ) THEN
	2483	!
[4842]	2484	!-- nesting_offl_parameters namelist was found but contained errors. Print an error message
	2485	!-- including the line that caused the problem.
	2486	BACKSPACE( 11 )
	2487	READ( 11 , '(A)' ) line
	2488	CALL parin_fail_message( 'nesting_offl_parameters', line )
[3347]	2489
[4842]	2490	ENDIF
[3347]	2491
[4842]	2492	END SUBROUTINE nesting_offl_parin
[3347]	2493
	2494
[4834]	2495	!--------------------------------------------------------------------------------------------------!
[3347]	2496	! Description:
	2497	! ------------
	2498	!> Writes information about offline nesting into HEADER file
[4834]	2499	!--------------------------------------------------------------------------------------------------!
	2500	SUBROUTINE nesting_offl_header ( io )
[3347]	2501
[4834]	2502	INTEGER(iwp), INTENT(IN) :: io !< Unit of the output file
[3347]	2503
[4834]	2504	WRITE ( io, 1 )
	2505	IF ( nesting_offline ) THEN
	2506	WRITE ( io, 3 )
	2507	ELSE
	2508	WRITE ( io, 2 )
	2509	ENDIF
[3347]	2510
[4834]	2511	1 FORMAT (//' Offline nesting in COSMO model:'/' -------------------------------'/)
[3347]	2512	2 FORMAT (' --> No offlince nesting is used (default) ')
	2513	3 FORMAT (' --> Offlince nesting is used. Boundary data is read from dynamic input file ')
	2514
[4834]	2515	END SUBROUTINE nesting_offl_header
[3347]	2516
[4842]	2517
[4834]	2518	!--------------------------------------------------------------------------------------------------!
[3347]	2519	! Description:
	2520	! ------------
[4834]	2521	!> Allocate arrays used to read boundary data from NetCDF file and initialize boundary data.
	2522	!--------------------------------------------------------------------------------------------------!
	2523	SUBROUTINE nesting_offl_init
[3347]	2524
[4834]	2525	INTEGER(iwp) :: i !< loop index for x-direction
	2526	INTEGER(iwp) :: j !< loop index for y-direction
	2527	INTEGER(iwp) :: n !< running index for chemical species
[4724]	2528
	2529	!
[4834]	2530	!-- Before arrays for the boundary data are allocated, the LOD of the dynamic input data at the
	2531	!-- boundaries is read.
[4724]	2532	#if defined ( __netcdf )
	2533	!
[4834]	2534	!-- Open file in read-only mode
	2535	CALL open_read_file( TRIM( input_file_dynamic ) // TRIM( coupling_char ), pids_id )
[4724]	2536	!
[4834]	2537	!-- Read attributes for LOD. In order to gurantee that also older drivers, where attribute is not
	2538	!-- given, are working, do not abort the run but assume LOD2 forcing.
	2539	CALL get_attribute( pids_id, char_lod, nest_offl%lod_east_pt, .FALSE., 'ls_forcing_left_pt', .FALSE. )
	2540	CALL get_attribute( pids_id, char_lod, nest_offl%lod_east_qv, .FALSE., 'ls_forcing_left_qv', .FALSE. )
	2541	CALL get_attribute( pids_id, char_lod, nest_offl%lod_east_u, .FALSE., 'ls_forcing_left_u', .FALSE. )
	2542	CALL get_attribute( pids_id, char_lod, nest_offl%lod_east_v, .FALSE., 'ls_forcing_left_v', .FALSE. )
	2543	CALL get_attribute( pids_id, char_lod, nest_offl%lod_east_w, .FALSE., 'ls_forcing_left_w', .FALSE. )
[4724]	2544
[4834]	2545	CALL get_attribute( pids_id, char_lod, nest_offl%lod_north_pt, .FALSE., 'ls_forcing_north_pt', .FALSE. )
	2546	CALL get_attribute( pids_id, char_lod, nest_offl%lod_north_qv, .FALSE., 'ls_forcing_north_qv', .FALSE. )
	2547	CALL get_attribute( pids_id, char_lod, nest_offl%lod_north_u, .FALSE., 'ls_forcing_north_u', .FALSE. )
	2548	CALL get_attribute( pids_id, char_lod, nest_offl%lod_north_v, .FALSE., 'ls_forcing_north_v', .FALSE. )
	2549	CALL get_attribute( pids_id, char_lod, nest_offl%lod_north_w, .FALSE., 'ls_forcing_north_w', .FALSE. )
[4724]	2550
[4834]	2551	CALL get_attribute( pids_id, char_lod, nest_offl%lod_south_pt, .FALSE., 'ls_forcing_south_pt', .FALSE. )
	2552	CALL get_attribute( pids_id, char_lod, nest_offl%lod_south_qv, .FALSE., 'ls_forcing_south_qv', .FALSE. )
	2553	CALL get_attribute( pids_id, char_lod, nest_offl%lod_south_u, .FALSE., 'ls_forcing_south_u', .FALSE. )
	2554	CALL get_attribute( pids_id, char_lod, nest_offl%lod_south_v, .FALSE., 'ls_forcing_south_v', .FALSE. )
	2555	CALL get_attribute( pids_id, char_lod, nest_offl%lod_south_w, .FALSE., 'ls_forcing_south_w', .FALSE. )
[4724]	2556
[4834]	2557	CALL get_attribute( pids_id, char_lod, nest_offl%lod_west_pt, .FALSE., 'ls_forcing_right_pt', .FALSE. )
	2558	CALL get_attribute( pids_id, char_lod, nest_offl%lod_west_qv, .FALSE., 'ls_forcing_right_qv', .FALSE. )
	2559	CALL get_attribute( pids_id, char_lod, nest_offl%lod_west_u, .FALSE., 'ls_forcing_right_u', .FALSE. )
	2560	CALL get_attribute( pids_id, char_lod, nest_offl%lod_west_v, .FALSE., 'ls_forcing_right_v', .FALSE. )
	2561	CALL get_attribute( pids_id, char_lod, nest_offl%lod_west_w, .FALSE., 'ls_forcing_right_w', .FALSE. )
[4724]	2562
[4834]	2563	CALL get_attribute( pids_id, char_lod, nest_offl%lod_top_pt, .FALSE., 'ls_forcing_top_pt', .FALSE. )
	2564	CALL get_attribute( pids_id, char_lod, nest_offl%lod_top_qv, .FALSE., 'ls_forcing_top_qv', .FALSE. )
	2565	CALL get_attribute( pids_id, char_lod, nest_offl%lod_top_u, .FALSE., 'ls_forcing_top_u', .FALSE. )
	2566	CALL get_attribute( pids_id, char_lod, nest_offl%lod_top_v, .FALSE., 'ls_forcing_top_v', .FALSE. )
	2567	CALL get_attribute( pids_id, char_lod, nest_offl%lod_top_w, .FALSE., 'ls_forcing_top_w', .FALSE. )
[4724]	2568
[4834]	2569	CALL close_input_file( pids_id )
[4724]	2570	#endif
	2571	!
[4834]	2572	!-- Temporary workaround until most of the dynamic drivers contain a LOD attribute. So far INIFOR
	2573	!-- did not provide the LOD attribute. In order to still use these older dynamic drivers, provide a
	2574	!-- temporary workaround. If the LOD is not given, a NetCDF interal error will occur but the
	2575	!-- simulation will not be aborted since the no_abort flag is passed. However, the respective
	2576	!-- attribute value might be given an arbitrary number. Hence, check for valid LOD's and manually
	2577	!-- set them to LOD 2 (as assumed so far). Note, this workaround should be removed later (date of
	2578	!-- reference: 6. Oct. 2020).
	2579	IF ( nest_offl%lod_east_pt /= 1 .AND. nest_offl%lod_east_pt /= 2 ) nest_offl%lod_east_pt = 2
	2580	IF ( nest_offl%lod_east_qv /= 1 .AND. nest_offl%lod_east_qv /= 2 ) nest_offl%lod_east_qv = 2
	2581	IF ( nest_offl%lod_east_u /= 1 .AND. nest_offl%lod_east_u /= 2 ) nest_offl%lod_east_u = 2
	2582	IF ( nest_offl%lod_east_v /= 1 .AND. nest_offl%lod_east_v /= 2 ) nest_offl%lod_east_v = 2
	2583	IF ( nest_offl%lod_east_w /= 1 .AND. nest_offl%lod_east_w /= 2 ) nest_offl%lod_east_w = 2
[4724]	2584
[4834]	2585	IF ( nest_offl%lod_north_pt /= 1 .AND. nest_offl%lod_north_pt /= 2 ) nest_offl%lod_north_pt = 2
	2586	IF ( nest_offl%lod_north_qv /= 1 .AND. nest_offl%lod_north_qv /= 2 ) nest_offl%lod_north_qv = 2
	2587	IF ( nest_offl%lod_north_u /= 1 .AND. nest_offl%lod_north_u /= 2 ) nest_offl%lod_north_u = 2
	2588	IF ( nest_offl%lod_north_v /= 1 .AND. nest_offl%lod_north_v /= 2 ) nest_offl%lod_north_v = 2
	2589	IF ( nest_offl%lod_north_w /= 1 .AND. nest_offl%lod_north_w /= 2 ) nest_offl%lod_north_w = 2
[4724]	2590
[4834]	2591	IF ( nest_offl%lod_south_pt /= 1 .AND. nest_offl%lod_south_pt /= 2 ) nest_offl%lod_south_pt = 2
	2592	IF ( nest_offl%lod_south_qv /= 1 .AND. nest_offl%lod_south_qv /= 2 ) nest_offl%lod_south_qv = 2
	2593	IF ( nest_offl%lod_south_u /= 1 .AND. nest_offl%lod_south_u /= 2 ) nest_offl%lod_south_u = 2
	2594	IF ( nest_offl%lod_south_v /= 1 .AND. nest_offl%lod_south_v /= 2 ) nest_offl%lod_south_v = 2
	2595	IF ( nest_offl%lod_south_w /= 1 .AND. nest_offl%lod_south_w /= 2 ) nest_offl%lod_south_w = 2
[4724]	2596
[4834]	2597	IF ( nest_offl%lod_west_pt /= 1 .AND. nest_offl%lod_west_pt /= 2 ) nest_offl%lod_west_pt = 2
	2598	IF ( nest_offl%lod_west_qv /= 1 .AND. nest_offl%lod_west_qv /= 2 ) nest_offl%lod_west_qv = 2
	2599	IF ( nest_offl%lod_west_u /= 1 .AND. nest_offl%lod_west_u /= 2 ) nest_offl%lod_west_u = 2
	2600	IF ( nest_offl%lod_west_v /= 1 .AND. nest_offl%lod_west_v /= 2 ) nest_offl%lod_west_v = 2
	2601	IF ( nest_offl%lod_west_w /= 1 .AND. nest_offl%lod_west_w /= 2 ) nest_offl%lod_west_w = 2
[4724]	2602
[4834]	2603	IF ( nest_offl%lod_top_pt /= 1 .AND. nest_offl%lod_top_pt /= 2 ) nest_offl%lod_top_pt = 2
	2604	IF ( nest_offl%lod_top_qv /= 1 .AND. nest_offl%lod_top_qv /= 2 ) nest_offl%lod_top_qv = 2
	2605	IF ( nest_offl%lod_top_u /= 1 .AND. nest_offl%lod_top_u /= 2 ) nest_offl%lod_top_u = 2
	2606	IF ( nest_offl%lod_top_v /= 1 .AND. nest_offl%lod_top_v /= 2 ) nest_offl%lod_top_v = 2
	2607	IF ( nest_offl%lod_top_w /= 1 .AND. nest_offl%lod_top_w /= 2 ) nest_offl%lod_top_w = 2
[4724]	2608	!
[4834]	2609	!-- For consistency, check if all boundary input variables have the same LOD.
	2610	IF ( MAX( nest_offl%lod_east_pt, nest_offl%lod_east_qv, nest_offl%lod_east_u, &
	2611	nest_offl%lod_east_v, nest_offl%lod_east_w, &
	2612	nest_offl%lod_north_pt, nest_offl%lod_north_qv, nest_offl%lod_north_u, &
	2613	nest_offl%lod_north_v, nest_offl%lod_north_w, &
	2614	nest_offl%lod_south_pt, nest_offl%lod_south_qv, nest_offl%lod_south_u, &
	2615	nest_offl%lod_south_v, nest_offl%lod_south_w, &
	2616	nest_offl%lod_north_pt, nest_offl%lod_north_qv, nest_offl%lod_north_u, &
	2617	nest_offl%lod_north_v, nest_offl%lod_north_w, &
	2618	nest_offl%lod_top_pt, nest_offl%lod_top_qv, nest_offl%lod_top_u, &
	2619	nest_offl%lod_top_v, nest_offl%lod_top_w ) &
	2620	/= &
	2621	MIN( nest_offl%lod_east_pt, nest_offl%lod_east_qv, nest_offl%lod_east_u, &
	2622	nest_offl%lod_east_v, nest_offl%lod_east_w, &
	2623	nest_offl%lod_north_pt, nest_offl%lod_north_qv, nest_offl%lod_north_u, &
	2624	nest_offl%lod_north_v, nest_offl%lod_north_w, &
	2625	nest_offl%lod_south_pt, nest_offl%lod_south_qv, nest_offl%lod_south_u, &
	2626	nest_offl%lod_south_v, nest_offl%lod_south_w, &
	2627	nest_offl%lod_north_pt, nest_offl%lod_north_qv, nest_offl%lod_north_u, &
	2628	nest_offl%lod_north_v, nest_offl%lod_north_w, &
	2629	nest_offl%lod_top_pt, nest_offl%lod_top_qv, nest_offl%lod_top_u, &
	2630	nest_offl%lod_top_v, nest_offl%lod_top_w ) ) THEN
	2631	message_string = 'A mixture of different LOD for the provided boundary data is not ' // &
	2632	'possible.'
	2633	CALL message( 'nesting_offl_init', 'PA0504', 1, 2, 0, 6, 0 )
	2634	ENDIF
[4724]	2635	!
[4834]	2636	!-- As all LODs are the same, store it.
	2637	lod = nest_offl%lod_east_u
[4724]	2638	!
[4834]	2639	!-- Allocate arrays for geostrophic wind components. Arrays will incorporate 2 time levels in order
	2640	!-- to interpolate in between.
	2641	ALLOCATE( nest_offl%ug(0:1,1:nzt) )
	2642	ALLOCATE( nest_offl%vg(0:1,1:nzt) )
[3347]	2643	!
[4834]	2644	!-- Set index range according to the given LOD in order to allocate the input arrays.
	2645	IF ( bc_dirichlet_l .OR. bc_dirichlet_r ) THEN
	2646	IF ( lod == 2 ) THEN
	2647	j_start = nys
	2648	j_start_v = nysv
	2649	j_end = nyn
	2650	ELSE
	2651	j_start = 1
	2652	j_start_v = 1
	2653	j_end = 1
[4724]	2654	ENDIF
[4834]	2655	ENDIF
[4724]	2656
[4834]	2657	IF ( bc_dirichlet_n .OR. bc_dirichlet_s ) THEN
	2658	IF( lod == 2 ) THEN
	2659	i_start = nxl
	2660	i_start_u = nxlu
	2661	i_end = nxr
[4125]	2662	ELSE
[4834]	2663	i_start = 1
	2664	i_start_u = 1
	2665	i_end = 1
[3347]	2666	ENDIF
[4834]	2667	ENDIF
[4125]	2668	!
[4834]	2669	!-- Allocate arrays for reading left/right boundary values. Arrays will incorporate 2 time levels in
	2670	!-- order to interpolate in between. Depending on the given LOD, the x-, or y-dimension will be
	2671	!-- either nxl:nxr, or nys:nyn (for LOD=2), or it reduces to one element for LOD=1. If the core has
	2672	!-- no lateral boundary, allocate a dummy array as well, in order to enable netcdf parallel access.
	2673	!-- Dummy arrays will be allocated with dimension length zero.
	2674	IF ( bc_dirichlet_l ) THEN
	2675	ALLOCATE( nest_offl%u_l(0:1,nzb+1:nzt,j_start:j_end) )
	2676	ALLOCATE( nest_offl%v_l(0:1,nzb+1:nzt,j_start_v:j_end) )
	2677	ALLOCATE( nest_offl%w_l(0:1,nzb+1:nzt-1,j_start:j_end) )
	2678	IF ( humidity ) ALLOCATE( nest_offl%q_l(0:1,nzb+1:nzt,j_start:j_end) )
	2679	IF ( .NOT. neutral ) ALLOCATE( nest_offl%pt_l(0:1,nzb+1:nzt,j_start:j_end) )
	2680	IF ( air_chemistry .AND. nesting_offline_chem ) &
	2681	ALLOCATE( nest_offl%chem_l(0:1,nzb+1:nzt,j_start:j_end,1:UBOUND( chem_species, 1 )) )
	2682	ELSE
	2683	ALLOCATE( nest_offl%u_l(1:1,1:1,1:1) )
	2684	ALLOCATE( nest_offl%v_l(1:1,1:1,1:1) )
	2685	ALLOCATE( nest_offl%w_l(1:1,1:1,1:1) )
	2686	IF ( humidity ) ALLOCATE( nest_offl%q_l(1:1,1:1,1:1) )
	2687	IF ( .NOT. neutral ) ALLOCATE( nest_offl%pt_l(1:1,1:1,1:1) )
	2688	IF ( air_chemistry .AND. nesting_offline_chem ) &
	2689	ALLOCATE( nest_offl%chem_l(1:1,1:1,1:1,1:UBOUND( chem_species, 1 )) )
	2690	ENDIF
	2691	IF ( bc_dirichlet_r ) THEN
	2692	ALLOCATE( nest_offl%u_r(0:1,nzb+1:nzt,j_start:j_end) )
	2693	ALLOCATE( nest_offl%v_r(0:1,nzb+1:nzt,j_start_v:j_end) )
	2694	ALLOCATE( nest_offl%w_r(0:1,nzb+1:nzt-1,j_start:j_end) )
	2695	IF ( humidity ) ALLOCATE( nest_offl%q_r(0:1,nzb+1:nzt,j_start:j_end) )
	2696	IF ( .NOT. neutral ) ALLOCATE( nest_offl%pt_r(0:1,nzb+1:nzt,j_start:j_end) )
	2697	IF ( air_chemistry .AND. nesting_offline_chem ) &
	2698	ALLOCATE( nest_offl%chem_r(0:1,nzb+1:nzt,j_start:j_end,1:UBOUND( chem_species, 1 )) )
	2699	ELSE
	2700	ALLOCATE( nest_offl%u_r(1:1,1:1,1:1) )
	2701	ALLOCATE( nest_offl%v_r(1:1,1:1,1:1) )
	2702	ALLOCATE( nest_offl%w_r(1:1,1:1,1:1) )
	2703	IF ( humidity ) ALLOCATE( nest_offl%q_r(1:1,1:1,1:1) )
	2704	IF ( .NOT. neutral ) ALLOCATE( nest_offl%pt_r(1:1,1:1,1:1) )
	2705	IF ( air_chemistry .AND. nesting_offline_chem ) &
	2706	ALLOCATE( nest_offl%chem_r(1:1,1:1,1:1,1:UBOUND( chem_species, 1 )) )
	2707	ENDIF
[4125]	2708	!
[4834]	2709	!-- Allocate arrays for reading north/south boundary values. Arrays will incorporate 2 time levels
	2710	!-- in order to interpolate in between. If the core has no boundary, allocate a dummy array, in
	2711	!-- order to enable netcdf parallel access. Dummy arrays will be allocated with dimension length
	2712	!-- zero.
	2713	IF ( bc_dirichlet_n ) THEN
	2714	ALLOCATE( nest_offl%u_n(0:1,nzb+1:nzt,i_start_u:i_end) )
	2715	ALLOCATE( nest_offl%v_n(0:1,nzb+1:nzt,i_start:i_end) )
	2716	ALLOCATE( nest_offl%w_n(0:1,nzb+1:nzt-1,i_start:i_end) )
	2717	IF ( humidity ) ALLOCATE( nest_offl%q_n(0:1,nzb+1:nzt,i_start:i_end) )
	2718	IF ( .NOT. neutral ) ALLOCATE( nest_offl%pt_n(0:1,nzb+1:nzt,i_start:i_end) )
	2719	IF ( air_chemistry .AND. nesting_offline_chem ) &
	2720	ALLOCATE( nest_offl%chem_n(0:1,nzb+1:nzt,i_start:i_end,1:UBOUND( chem_species, 1 )) )
	2721	ELSE
	2722	ALLOCATE( nest_offl%u_n(1:1,1:1,1:1) )
	2723	ALLOCATE( nest_offl%v_n(1:1,1:1,1:1) )
	2724	ALLOCATE( nest_offl%w_n(1:1,1:1,1:1) )
	2725	IF ( humidity ) ALLOCATE( nest_offl%q_n(1:1,1:1,1:1) )
	2726	IF ( .NOT. neutral ) ALLOCATE( nest_offl%pt_n(1:1,1:1,1:1) )
	2727	IF ( air_chemistry .AND. nesting_offline_chem ) &
	2728	ALLOCATE( nest_offl%chem_n(1:1,1:1,1:1,1:UBOUND( chem_species, 1 )) )
	2729	ENDIF
	2730	IF ( bc_dirichlet_s ) THEN
	2731	ALLOCATE( nest_offl%u_s(0:1,nzb+1:nzt,i_start_u:i_end) )
	2732	ALLOCATE( nest_offl%v_s(0:1,nzb+1:nzt,i_start:i_end) )
	2733	ALLOCATE( nest_offl%w_s(0:1,nzb+1:nzt-1,i_start:i_end) )
	2734	IF ( humidity ) ALLOCATE( nest_offl%q_s(0:1,nzb+1:nzt,i_start:i_end) )
	2735	IF ( .NOT. neutral ) ALLOCATE( nest_offl%pt_s(0:1,nzb+1:nzt,i_start:i_end) )
	2736	IF ( air_chemistry .AND. nesting_offline_chem ) &
	2737	ALLOCATE( nest_offl%chem_s(0:1,nzb+1:nzt,i_start:i_end,1:UBOUND( chem_species, 1 )) )
	2738	ELSE
	2739	ALLOCATE( nest_offl%u_s(1:1,1:1,1:1) )
	2740	ALLOCATE( nest_offl%v_s(1:1,1:1,1:1) )
	2741	ALLOCATE( nest_offl%w_s(1:1,1:1,1:1) )
	2742	IF ( humidity ) ALLOCATE( nest_offl%q_s(1:1,1:1,1:1) )
	2743	IF ( .NOT. neutral ) ALLOCATE( nest_offl%pt_s(1:1,1:1,1:1) )
	2744	IF ( air_chemistry .AND. nesting_offline_chem ) &
	2745	ALLOCATE( nest_offl%chem_s(1:1,1:1,1:1,1:UBOUND( chem_species, 1 )) )
	2746	ENDIF
[3347]	2747	!
[4834]	2748	!-- Allocate arrays for reading data at the top boundary. In contrast to the lateral boundaries,
	2749	!-- each core reads these data so that no dummy arrays need to be allocated.
	2750	IF ( lod == 2 ) THEN
	2751	ALLOCATE( nest_offl%u_top(0:1,nys:nyn,nxlu:nxr) )
	2752	ALLOCATE( nest_offl%v_top(0:1,nysv:nyn,nxl:nxr) )
	2753	ALLOCATE( nest_offl%w_top(0:1,nys:nyn,nxl:nxr) )
	2754	IF ( humidity ) ALLOCATE( nest_offl%q_top(0:1,nys:nyn,nxl:nxr) )
	2755	IF ( .NOT. neutral ) ALLOCATE( nest_offl%pt_top(0:1,nys:nyn,nxl:nxr) )
	2756	IF ( air_chemistry .AND. nesting_offline_chem ) &
	2757	ALLOCATE( nest_offl%chem_top(0:1,nys:nyn,nxl:nxr,1:UBOUND( chem_species, 1 )) )
	2758	ELSE
	2759	ALLOCATE( nest_offl%u_top(0:1,1:1,1:1) )
	2760	ALLOCATE( nest_offl%v_top(0:1,1:1,1:1) )
	2761	ALLOCATE( nest_offl%w_top(0:1,1:1,1:1) )
	2762	IF ( humidity ) ALLOCATE( nest_offl%q_top(0:1,1:1,1:1) )
	2763	IF ( .NOT. neutral ) ALLOCATE( nest_offl%pt_top(0:1,1:1,1:1) )
	2764	IF ( air_chemistry .AND. nesting_offline_chem ) &
	2765	ALLOCATE( nest_offl%chem_top(0:1,1:1,1:1,1:UBOUND( chem_species, 1 )) )
	2766	ENDIF
	2767	!
	2768	!-- For chemical species, create the names of the variables. This is necessary to identify the
	2769	!-- respective variable and write it onto the correct array in the chem_species datatype.
	2770	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	2771	ALLOCATE( nest_offl%chem_from_file_l(1:UBOUND( chem_species, 1 )) )
	2772	ALLOCATE( nest_offl%chem_from_file_n(1:UBOUND( chem_species, 1 )) )
	2773	ALLOCATE( nest_offl%chem_from_file_r(1:UBOUND( chem_species, 1 )) )
	2774	ALLOCATE( nest_offl%chem_from_file_s(1:UBOUND( chem_species, 1 )) )
	2775	ALLOCATE( nest_offl%chem_from_file_t(1:UBOUND( chem_species, 1 )) )
[4273]	2776
[4834]	2777	ALLOCATE( nest_offl%var_names_chem_l(1:UBOUND( chem_species, 1 )) )
	2778	ALLOCATE( nest_offl%var_names_chem_n(1:UBOUND( chem_species, 1 )) )
	2779	ALLOCATE( nest_offl%var_names_chem_r(1:UBOUND( chem_species, 1 )) )
	2780	ALLOCATE( nest_offl%var_names_chem_s(1:UBOUND( chem_species, 1 )) )
	2781	ALLOCATE( nest_offl%var_names_chem_t(1:UBOUND( chem_species, 1 )) )
[3737]	2782	!
[4834]	2783	!-- Initialize flags that indicate whether the variable is on file or not. Please note, this is
	2784	!-- only necessary for chemistry variables.
	2785	nest_offl%chem_from_file_l(:) = .FALSE.
	2786	nest_offl%chem_from_file_n(:) = .FALSE.
	2787	nest_offl%chem_from_file_r(:) = .FALSE.
	2788	nest_offl%chem_from_file_s(:) = .FALSE.
	2789	nest_offl%chem_from_file_t(:) = .FALSE.
[4273]	2790
[4834]	2791	DO n = 1, UBOUND( chem_species, 1 )
	2792	nest_offl%var_names_chem_l(n) = nest_offl%char_l // TRIM(chem_species(n)%name)
	2793	nest_offl%var_names_chem_n(n) = nest_offl%char_n // TRIM(chem_species(n)%name)
	2794	nest_offl%var_names_chem_r(n) = nest_offl%char_r // TRIM(chem_species(n)%name)
	2795	nest_offl%var_names_chem_s(n) = nest_offl%char_s // TRIM(chem_species(n)%name)
	2796	nest_offl%var_names_chem_t(n) = nest_offl%char_t // TRIM(chem_species(n)%name)
	2797	ENDDO
	2798	ENDIF
[3737]	2799	!
[4834]	2800	!-- Offline nesting for salsa
	2801	IF ( salsa ) CALL salsa_nesting_offl_init
[4270]	2802	!
[4834]	2803	!-- Before initial data input is initiated, check if dynamic input file is present.
	2804	IF ( .NOT. input_pids_dynamic ) THEN
	2805	message_string = 'nesting_offline = .TRUE. requires dynamic ' // &
	2806	'input file ' // TRIM( input_file_dynamic ) // TRIM( coupling_char )
	2807	CALL message( 'nesting_offl_init', 'PA0546', 1, 2, 0, 6, 0 )
	2808	ENDIF
[4226]	2809	!
[4834]	2810	!-- Read COSMO data at lateral and top boundaries
	2811	CALL nesting_offl_input
[3347]	2812	!
[4834]	2813	!-- Check if sufficient time steps are provided to cover the entire simulation. Note, dynamic input
	2814	!-- is only required for the 3D simulation, not for the soil/wall spinup. However, as the spinup
	2815	!-- time is added to the end_time, this must be considered here.
	2816	IF ( end_time - spinup_time > nest_offl%time(nest_offl%nt-1) ) THEN
	2817	message_string = 'end_time of the simulation exceeds the ' // &
	2818	'time dimension in the dynamic input file.'
	2819	CALL message( 'nesting_offl_init', 'PA0183', 1, 2, 0, 6, 0 )
	2820	ENDIF
[4169]	2821	!
[4834]	2822	!-- Set indicies for boundary grid points
	2823	IF ( bc_dirichlet_l .OR. bc_dirichlet_r ) THEN
	2824	i_bound = MERGE( nxl - 1, nxr + 1, bc_dirichlet_l )
	2825	i_bound_u = MERGE( nxlu - 1, nxr + 1, bc_dirichlet_l )
	2826	ENDIF
	2827	IF ( bc_dirichlet_n .OR. bc_dirichlet_s ) THEN
	2828	j_bound = MERGE( nys - 1, nyn + 1, bc_dirichlet_s )
	2829	j_bound_v = MERGE( nysv - 1, nyn + 1, bc_dirichlet_s )
	2830	ENDIF
[4724]	2831	!
[4834]	2832	!-- Initialize boundary data. Please note, do not initialize boundaries in case of restart runs.
	2833	!-- This case the boundaries are already initialized and the boundary data from file would be on the
	2834	!-- wrong time level.
	2835	IF ( TRIM( initializing_actions ) /= 'read_restart_data' ) THEN
[4724]	2836	!
[4834]	2837	!-- Distinguish between LOD = 1 and LOD = 2 inititialization
	2838	IF ( lod == 2 ) THEN
	2839	IF ( bc_dirichlet_l ) THEN
	2840	u(nzb+1:nzt,nys:nyn,i_bound_u) = nest_offl%u_l(0,nzb+1:nzt,nys:nyn)
	2841	v(nzb+1:nzt,nysv:nyn,i_bound) = nest_offl%v_l(0,nzb+1:nzt,nysv:nyn)
	2842	w(nzb+1:nzt-1,nys:nyn,i_bound) = nest_offl%w_l(0,nzb+1:nzt-1,nys:nyn)
	2843	IF ( .NOT. neutral ) pt(nzb+1:nzt,nys:nyn,i_bound) = nest_offl%pt_l(0,nzb+1:nzt,nys:nyn)
	2844	IF ( humidity ) q(nzb+1:nzt,nys:nyn,i_bound) = nest_offl%q_l(0,nzb+1:nzt,nys:nyn)
	2845	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	2846	DO n = 1, UBOUND( chem_species, 1 )
	2847	IF( nest_offl%chem_from_file_l(n) ) THEN
	2848	chem_species(n)%conc(nzb+1:nzt,nys:nyn,i_bound) = &
	2849	nest_offl%chem_l(0,nzb+1:nzt,nys:nyn,n)
	2850	ENDIF
	2851	ENDDO
[4724]	2852	ENDIF
[4834]	2853	ENDIF
	2854	IF ( bc_dirichlet_r ) THEN
	2855	u(nzb+1:nzt,nys:nyn,i_bound_u) = nest_offl%u_r(0,nzb+1:nzt,nys:nyn)
	2856	v(nzb+1:nzt,nysv:nyn,i_bound) = nest_offl%v_r(0,nzb+1:nzt,nysv:nyn)
	2857	w(nzb+1:nzt-1,nys:nyn,i_bound) = nest_offl%w_r(0,nzb+1:nzt-1,nys:nyn)
	2858	IF ( .NOT. neutral ) pt(nzb+1:nzt,nys:nyn,i_bound) = nest_offl%pt_r(0,nzb+1:nzt,nys:nyn)
	2859	IF ( humidity ) q(nzb+1:nzt,nys:nyn,i_bound) = nest_offl%q_r(0,nzb+1:nzt,nys:nyn)
	2860	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	2861	DO n = 1, UBOUND( chem_species, 1 )
	2862	IF( nest_offl%chem_from_file_r(n) ) THEN
	2863	chem_species(n)%conc(nzb+1:nzt,nys:nyn,i_bound) = &
	2864	nest_offl%chem_r(0,nzb+1:nzt,nys:nyn,n)
	2865	ENDIF
	2866	ENDDO
[4724]	2867	ENDIF
[4834]	2868	ENDIF
[4724]	2869
[4834]	2870	IF ( bc_dirichlet_n) THEN
	2871	u(nzb+1:nzt,j_bound,nxlu:nxr) = nest_offl%u_n(0,nzb+1:nzt,nxlu:nxr)
	2872	v(nzb+1:nzt,j_bound_v,nxl:nxr) = nest_offl%v_n(0,nzb+1:nzt,nxl:nxr)
	2873	w(nzb+1:nzt-1,j_bound,nxl:nxr) = nest_offl%w_n(0,nzb+1:nzt-1,nxl:nxr)
	2874	IF ( .NOT. neutral ) pt(nzb+1:nzt,j_bound,nxl:nxr) = nest_offl%pt_n(0,nzb+1:nzt,nxl:nxr)
	2875	IF ( humidity ) q(nzb+1:nzt,j_bound,nxl:nxr) = nest_offl%q_n(0,nzb+1:nzt,nxl:nxr)
	2876	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	2877	DO n = 1, UBOUND( chem_species, 1 )
	2878	IF( nest_offl%chem_from_file_n(n) ) THEN
	2879	chem_species(n)%conc(nzb+1:nzt,j_bound,nxl:nxr) = &
	2880	nest_offl%chem_n(0,nzb+1:nzt,nxl:nxr,n)
	2881	ENDIF
	2882	ENDDO
[4724]	2883	ENDIF
[4834]	2884	ENDIF
	2885	IF ( bc_dirichlet_s) THEN
	2886	u(nzb+1:nzt,j_bound,nxlu:nxr) = nest_offl%u_s(0,nzb+1:nzt,nxlu:nxr)
	2887	v(nzb+1:nzt,j_bound_v,nxl:nxr) = nest_offl%v_s(0,nzb+1:nzt,nxl:nxr)
	2888	w(nzb+1:nzt-1,j_bound,nxl:nxr) = nest_offl%w_s(0,nzb+1:nzt-1,nxl:nxr)
	2889	IF ( .NOT. neutral ) pt(nzb+1:nzt,j_bound,nxl:nxr) = nest_offl%pt_s(0,nzb+1:nzt,nxl:nxr)
	2890	IF ( humidity ) q(nzb+1:nzt,j_bound,nxl:nxr) = nest_offl%q_s(0,nzb+1:nzt,nxl:nxr)
[4273]	2891	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
[3891]	2892	DO n = 1, UBOUND( chem_species, 1 )
[4834]	2893	IF( nest_offl%chem_from_file_s(n) ) THEN
	2894	chem_species(n)%conc(nzb+1:nzt,j_bound,nxl:nxr) = &
	2895	nest_offl%chem_s(0,nzb+1:nzt,nxl:nxr,n)
[3891]	2896	ENDIF
	2897	ENDDO
	2898	ENDIF
[4834]	2899	ENDIF
	2900
	2901	u(nzt+1,nys:nyn,nxlu:nxr) = nest_offl%u_top(0,nys:nyn,nxlu:nxr)
	2902	v(nzt+1,nysv:nyn,nxl:nxr) = nest_offl%v_top(0,nysv:nyn,nxl:nxr)
	2903	w(nzt,nys:nyn,nxl:nxr) = nest_offl%w_top(0,nys:nyn,nxl:nxr)
	2904	w(nzt+1,nys:nyn,nxl:nxr) = nest_offl%w_top(0,nys:nyn,nxl:nxr)
	2905	IF ( .NOT. neutral ) pt(nzt+1,nys:nyn,nxl:nxr) = nest_offl%pt_top(0,nys:nyn,nxl:nxr)
	2906	IF ( humidity ) q(nzt+1,nys:nyn,nxl:nxr) = nest_offl%q_top(0,nys:nyn,nxl:nxr)
	2907	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	2908	DO n = 1, UBOUND( chem_species, 1 )
	2909	IF( nest_offl%chem_from_file_t(n) ) THEN
	2910	chem_species(n)%conc(nzt+1,nys:nyn,nxl:nxr) = &
	2911	nest_offl%chem_top(0,nys:nyn,nxl:nxr,n)
	2912	ENDIF
	2913	ENDDO
	2914	ENDIF
[4724]	2915	!
[4834]	2916	!-- LOD 1
	2917	ELSE
	2918	IF ( bc_dirichlet_l ) THEN
	2919	DO j = nys, nyn
	2920	u(nzb+1:nzt,j,i_bound_u) = nest_offl%u_l(0,nzb+1:nzt,1)
	2921	w(nzb+1:nzt-1,j,i_bound) = nest_offl%w_l(0,nzb+1:nzt-1,1)
	2922	ENDDO
	2923	DO j = nysv, nyn
	2924	v(nzb+1:nzt,j,i_bound) = nest_offl%v_l(0,nzb+1:nzt,1)
	2925	ENDDO
	2926	IF ( .NOT. neutral ) THEN
[4724]	2927	DO j = nys, nyn
[4834]	2928	pt(nzb+1:nzt,j,i_bound) = nest_offl%pt_l(0,nzb+1:nzt,1)
[3891]	2929	ENDDO
[4834]	2930	ENDIF
	2931	IF ( humidity ) THEN
	2932	DO j = nys, nyn
	2933	q(nzb+1:nzt,j,i_bound) = nest_offl%q_l(0,nzb+1:nzt,1)
[4724]	2934	ENDDO
[3891]	2935	ENDIF
[4834]	2936	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	2937	DO n = 1, UBOUND( chem_species, 1 )
	2938	IF( nest_offl%chem_from_file_l(n) ) THEN
	2939	DO j = nys, nyn
	2940	chem_species(n)%conc(nzb+1:nzt,j,i_bound) = &
	2941	nest_offl%chem_l(0,nzb+1:nzt,1,n)
	2942	ENDDO
	2943	ENDIF
	2944	ENDDO
	2945	ENDIF
	2946	ENDIF
	2947	IF ( bc_dirichlet_r ) THEN
	2948	DO j = nys, nyn
	2949	u(nzb+1:nzt,j,i_bound_u) = nest_offl%u_r(0,nzb+1:nzt,1)
	2950	w(nzb+1:nzt-1,j,i_bound) = nest_offl%w_r(0,nzb+1:nzt-1,1)
	2951	ENDDO
	2952	DO j = nysv, nyn
	2953	v(nzb+1:nzt,j,i_bound) = nest_offl%v_r(0,nzb+1:nzt,1)
	2954	ENDDO
	2955	IF ( .NOT. neutral ) THEN
[4724]	2956	DO j = nys, nyn
[4834]	2957	pt(nzb+1:nzt,j,i_bound) = nest_offl%pt_r(0,nzb+1:nzt,1)
[3891]	2958	ENDDO
[4834]	2959	ENDIF
	2960	IF ( humidity ) THEN
	2961	DO j = nys, nyn
	2962	q(nzb+1:nzt,j,i_bound) = nest_offl%q_r(0,nzb+1:nzt,1)
[4724]	2963	ENDDO
[3891]	2964	ENDIF
[4834]	2965	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	2966	DO n = 1, UBOUND( chem_species, 1 )
	2967	IF( nest_offl%chem_from_file_r(n) ) THEN
	2968	DO j = nys, nyn
	2969	chem_species(n)%conc(nzb+1:nzt,j,i_bound) = &
	2970	nest_offl%chem_r(0,nzb+1:nzt,1,n)
	2971	ENDDO
	2972	ENDIF
[4724]	2973	ENDDO
[4834]	2974	ENDIF
	2975	ENDIF
	2976	IF ( bc_dirichlet_n ) THEN
	2977	DO i = nxlu, nxr
	2978	u(nzb+1:nzt,j_bound,i) = nest_offl%u_n(0,nzb+1:nzt,1)
	2979	ENDDO
	2980	DO i = nxl, nxr
	2981	v(nzb+1:nzt,j_bound_v,i) = nest_offl%v_n(0,nzb+1:nzt,1)
	2982	w(nzb+1:nzt-1,j_bound,i) = nest_offl%w_n(0,nzb+1:nzt-1,1)
	2983	ENDDO
	2984	IF ( .NOT. neutral ) THEN
[4724]	2985	DO i = nxl, nxr
[4834]	2986	pt(nzb+1:nzt,j_bound,i) = nest_offl%pt_n(0,nzb+1:nzt,1)
[4724]	2987	ENDDO
	2988	ENDIF
[4834]	2989	IF ( humidity ) THEN
	2990	DO i = nxl, nxr
	2991	q(nzb+1:nzt,j_bound,i) = nest_offl%q_n(0,nzb+1:nzt,1)
[4724]	2992	ENDDO
[4834]	2993	ENDIF
	2994	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	2995	DO n = 1, UBOUND( chem_species, 1 )
	2996	IF( nest_offl%chem_from_file_n(n) ) THEN
	2997	DO i = nxl, nxr
	2998	chem_species(n)%conc(nzb+1:nzt,j_bound,i) = &
	2999	nest_offl%chem_n(0,nzb+1:nzt,1,n)
	3000	ENDDO
	3001	ENDIF
	3002	ENDDO
	3003	ENDIF
	3004	ENDIF
	3005	IF ( bc_dirichlet_s ) THEN
	3006	DO i = nxlu, nxr
	3007	u(nzb+1:nzt,j_bound,i) = nest_offl%u_s(0,nzb+1:nzt,1)
	3008	ENDDO
	3009	DO i = nxl, nxr
	3010	v(nzb+1:nzt,j_bound_v,i) = nest_offl%v_s(0,nzb+1:nzt,1)
	3011	w(nzb+1:nzt-1,j_bound,i) = nest_offl%w_s(0,nzb+1:nzt-1,1)
	3012	ENDDO
	3013	IF ( .NOT. neutral ) THEN
[4724]	3014	DO i = nxl, nxr
[4834]	3015	pt(nzb+1:nzt,j_bound,i) = nest_offl%pt_s(0,nzb+1:nzt,1)
[4724]	3016	ENDDO
	3017	ENDIF
[4834]	3018	IF ( humidity ) THEN
	3019	DO i = nxl, nxr
	3020	q(nzb+1:nzt,j_bound,i) = nest_offl%q_s(0,nzb+1:nzt,1)
	3021	ENDDO
	3022	ENDIF
[4273]	3023	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
[3891]	3024	DO n = 1, UBOUND( chem_species, 1 )
[4834]	3025	IF( nest_offl%chem_from_file_s(n) ) THEN
	3026	DO i = nxl, nxr
	3027	chem_species(n)%conc(nzb+1:nzt,j_bound,i) = &
	3028	nest_offl%chem_s(0,nzb+1:nzt,1,n)
	3029	ENDDO
[3891]	3030	ENDIF
	3031	ENDDO
	3032	ENDIF
[3737]	3033	ENDIF
[4581]	3034
[4834]	3035	u(nzt+1,nys:nyn,nxlu:nxr) = nest_offl%u_top(0,1,1)
	3036	v(nzt+1,nysv:nyn,nxl:nxr) = nest_offl%v_top(0,1,1)
	3037	w(nzt,nys:nyn,nxl:nxr) = nest_offl%w_top(0,1,1)
	3038	w(nzt+1,nys:nyn,nxl:nxr) = nest_offl%w_top(0,1,1)
	3039	IF ( .NOT. neutral ) pt(nzt+1,nys:nyn,nxl:nxr) = nest_offl%pt_top(0,1,1)
	3040	IF ( humidity ) q(nzt+1,nys:nyn,nxl:nxr) = nest_offl%q_top(0,1,1)
	3041	IF ( air_chemistry .AND. nesting_offline_chem ) THEN
	3042	DO n = 1, UBOUND( chem_species, 1 )
	3043	IF( nest_offl%chem_from_file_t(n) ) THEN
	3044	chem_species(n)%conc(nzt+1,nys:nyn,nxl:nxr) = nest_offl%chem_top(0,1,1,n)
	3045	ENDIF
	3046	ENDDO
	3047	ENDIF
[4273]	3048	ENDIF
[3347]	3049	!
[4834]	3050	!-- In case of offline nesting the pressure forms itself based on the prescribed lateral
	3051	!-- boundary conditions. Hence, explicit forcing by pressure gradients via geostrophic wind
	3052	!-- components is not necessary and would be canceled out by the perturbation pressure otherwise.
	3053	!-- For this reason, set geostrophic wind components to zero.
	3054	ug(nzb+1:nzt) = 0.0_wp
	3055	vg(nzb+1:nzt) = 0.0_wp
	3056
	3057	ENDIF
[3891]	3058	!
[4834]	3059	!-- After boundary data is initialized, mask topography at the boundaries for the velocity
	3060	!-- components.
	3061	u = MERGE( u, 0.0_wp, BTEST( wall_flags_total_0, 1 ) )
	3062	v = MERGE( v, 0.0_wp, BTEST( wall_flags_total_0, 2 ) )
	3063	w = MERGE( w, 0.0_wp, BTEST( wall_flags_total_0, 3 ) )
[3347]	3064	!
[4834]	3065	!-- Initial calculation of the boundary layer depth from the prescribed boundary data. This is
	3066	!-- required for initialize the synthetic turbulence generator correctly.
	3067	CALL nesting_offl_calc_zi
	3068	!
	3069	!-- After boundary data is initialized, ensure mass conservation. Not necessary in restart runs.
	3070	IF ( TRIM( initializing_actions ) /= 'read_restart_data' ) THEN
	3071	CALL nesting_offl_mass_conservation
	3072	ENDIF
[3347]	3073
[4834]	3074	END SUBROUTINE nesting_offl_init
	3075
	3076	!--------------------------------------------------------------------------------------------------!
[3347]	3077	! Description:
[4834]	3078	!--------------------------------------------------------------------------------------------------!
	3079	!> Interpolation function, used to interpolate boundary data in time.
	3080	!--------------------------------------------------------------------------------------------------!
	3081	FUNCTION interpolate_in_time( var_t1, var_t2, fac )
[3347]	3082
[4834]	3083	REAL(wp) :: fac !< interpolation factor
	3084	REAL(wp) :: interpolate_in_time !< time-interpolated boundary value
	3085	REAL(wp) :: var_t1 !< boundary value at t1
	3086	REAL(wp) :: var_t2 !< boundary value at t2
[3347]	3087
[4834]	3088	interpolate_in_time = ( 1.0_wp - fac ) * var_t1 + fac * var_t2
[3347]	3089
[4834]	3090	END FUNCTION interpolate_in_time
[3347]	3091
	3092
	3093
	3094	END MODULE nesting_offl_mod

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

Context Navigation

source: palm/trunk/SOURCE/nesting_offl_mod.f90 @ 4875

Download in other formats: