Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |