fix 'pure virtual method called' error (#774)

### Description
Move `DeviceVector` inside `setInitialConditionsOnGrid` to avoid the
`pure virtual method called` error, first seen in #761 .

This branch is based on dependabot/submodules/extern/amrex-6d9c25b . 

### Related issues
Fix issues in PR #761 

### Checklist
_Before this pull request can be reviewed, all of these tasks should be
completed. Denote completed tasks with an `x` inside the square brackets
`[ ]` in the Markdown source below:_
- [ ] I have added a description (see above).
- [ ] I have added a link to any related issues see (see above).
- [ ] I have read the [Contributing
Guide](https://github.com/quokka-astro/quokka/blob/development/CONTRIBUTING.md).
- [ ] I have added tests for any new physics that this PR adds to the
code.
- [ ] I have tested this PR on my local computer and all tests pass.
- [ ] I have manually triggered the GPU tests with the magic comment
`/azp run`.
- [ ] I have requested a reviewer for this PR.

---------

Signed-off-by: dependabot[bot] <support@github.com>
Co-authored-by: dependabot[bot] <49699333+dependabot[bot]@users.noreply.github.com>

src/problems/RadTube/test_radiation_tube.cpp

@@ -85,17 +85,12 @@ RadSystem<TubeProblem>::DefineOpacityExponentsAndLowerValues(amrex::GpuArray<dou
 	return exponents_and_values;
 }
 
-// declare global variables
-// initial conditions read from file
-amrex::Gpu::HostVector<double> x_arr;
-amrex::Gpu::HostVector<double> rho_arr;
-amrex::Gpu::HostVector<double> Pgas_arr;
-amrex::Gpu::HostVector<double> Erad_arr;
-
-amrex::Gpu::DeviceVector<double> x_arr_g;
-amrex::Gpu::DeviceVector<double> rho_arr_g;
-amrex::Gpu::DeviceVector<double> Pgas_arr_g;
-amrex::Gpu::DeviceVector<double> Erad_arr_g;
+template <> struct SimulationData<TubeProblem> {
+	amrex::Gpu::DeviceVector<double> x_arr_g;
+	amrex::Gpu::DeviceVector<double> rho_arr_g;
+	amrex::Gpu::DeviceVector<double> Pgas_arr_g;
+	amrex::Gpu::DeviceVector<double> Erad_arr_g;
+};
 
 template <> void QuokkaSimulation<TubeProblem>::preCalculateInitialConditions()
 {
@@ -106,6 +101,11 @@ template <> void QuokkaSimulation<TubeProblem>::preCalculateInitialConditions()
 	std::string header;
 	std::getline(fstream, header);
 
+	std::vector<double> x_arr;
+	std::vector<double> rho_arr;
+	std::vector<double> Pgas_arr;
+	std::vector<double> Erad_arr;
+
 	for (std::string line; std::getline(fstream, line);) {
 		std::istringstream iss(line);
 		std::vector<double> values;
@@ -123,16 +123,16 @@ template <> void QuokkaSimulation<TubeProblem>::preCalculateInitialConditions()
 		Erad_arr.push_back(Erad);
 	}
 
-	x_arr_g.resize(x_arr.size());
-	rho_arr_g.resize(rho_arr.size());
-	Pgas_arr_g.resize(Pgas_arr.size());
-	Erad_arr_g.resize(Erad_arr.size());
+	userData_.x_arr_g.resize(x_arr.size());
+	userData_.rho_arr_g.resize(rho_arr.size());
+	userData_.Pgas_arr_g.resize(Pgas_arr.size());
+	userData_.Erad_arr_g.resize(Erad_arr.size());
 
 	// copy to device
-	amrex::Gpu::copyAsync(amrex::Gpu::hostToDevice, x_arr.begin(), x_arr.end(), x_arr_g.begin());
-	amrex::Gpu::copyAsync(amrex::Gpu::hostToDevice, rho_arr.begin(), rho_arr.end(), rho_arr_g.begin());
-	amrex::Gpu::copyAsync(amrex::Gpu::hostToDevice, Pgas_arr.begin(), Pgas_arr.end(), Pgas_arr_g.begin());
-	amrex::Gpu::copyAsync(amrex::Gpu::hostToDevice, Erad_arr.begin(), Erad_arr.end(), Erad_arr_g.begin());
+	amrex::Gpu::copyAsync(amrex::Gpu::hostToDevice, x_arr.begin(), x_arr.end(), userData_.x_arr_g.begin());
+	amrex::Gpu::copyAsync(amrex::Gpu::hostToDevice, rho_arr.begin(), rho_arr.end(), userData_.rho_arr_g.begin());
+	amrex::Gpu::copyAsync(amrex::Gpu::hostToDevice, Pgas_arr.begin(), Pgas_arr.end(), userData_.Pgas_arr_g.begin());
+	amrex::Gpu::copyAsync(amrex::Gpu::hostToDevice, Erad_arr.begin(), Erad_arr.end(), userData_.Erad_arr_g.begin());
 	amrex::Gpu::streamSynchronizeAll();
 }
 
@@ -144,11 +144,11 @@ template <> void QuokkaSimulation<TubeProblem>::setInitialConditionsOnGrid(quokk
 	const amrex::Box &indexRange = grid_elem.indexRange_;
 	const amrex::Array4<double> &state_cc = grid_elem.array_;
 
-	auto const &x_ptr = x_arr_g.dataPtr();
-	auto const &rho_ptr = rho_arr_g.dataPtr();
-	auto const &Pgas_ptr = Pgas_arr_g.dataPtr();
-	auto const &Erad_ptr = Erad_arr_g.dataPtr();
-	int x_size = static_cast<int>(x_arr_g.size());
+	auto const &x_ptr = userData_.x_arr_g.dataPtr();
+	auto const &rho_ptr = userData_.rho_arr_g.dataPtr();
+	auto const &Pgas_ptr = userData_.Pgas_arr_g.dataPtr();
+	auto const &Erad_ptr = userData_.Erad_arr_g.dataPtr();
+	int x_size = static_cast<int>(userData_.x_arr_g.size());
 
 	const auto radBoundaries_g = RadSystem_Traits<TubeProblem>::radBoundaries;