Home

Context Navigation

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

Last change on this file since 4805 was 4795, checked in by suehring, 4 years ago
mesoscale nesting: bugfix in obtaining the correct timestamp in case of restart runs; virtual measurements: add missing control flags
Property svn:keywords set to `Id`
File size: 177.9 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |