Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |