From 3fa93e6c05b52ca1979ac62c55d7014873115d9c Mon Sep 17 00:00:00 2001
From: Liwei Ji <jiliwei.phys@gmail.com>
Date: Tue, 10 Mar 2026 21:15:02 -0400
Subject: [PATCH] TestNuPcsArdBH: add par of containing_bx

---
 TestNuPcsArdBH/param.ccl              | 36 +++++++++++++++++++++++----
 TestNuPcsArdBH/src/testnupcsardbh.cxx |  7 +++---
 TestNuPcsArdBH/test/nupcsardbh.par    |  6 +++++
 3 files changed, 41 insertions(+), 8 deletions(-)

diff --git a/TestNuPcsArdBH/param.ccl b/TestNuPcsArdBH/param.ccl
index c5b1992..d347dc0 100644
--- a/TestNuPcsArdBH/param.ccl
+++ b/TestNuPcsArdBH/param.ccl
@@ -1,10 +1,5 @@
 # Parameter definitions for thorn TestNuPcsArdBH
 
-CCTK_REAL BHmass "BH mass" STEERABLE=always
-{
-  0:* :: ""
-} 1.0
-
 CCTK_INT num_particles "Number of particles to create" STEERABLE=never
 {
   1:* :: "positive"
@@ -14,3 +9,34 @@ CCTK_INT random_seed "Random seed for particle initialization" STEERABLE=never
 {
   1:* :: "positive"
 } 42
+
+CCTK_REAL containing_bx_xmin "Spatial bounds for particle generation"
+{
+  *:* :: ""
+} -1.0
+
+CCTK_REAL containing_bx_xmax "Spatial bounds for particle generation"
+{
+  *:* :: ""
+} 1.0
+
+CCTK_REAL containing_bx_ymin "Spatial bounds for particle generation"
+{
+  *:* :: ""
+} -1.0
+
+CCTK_REAL containing_bx_ymax "Spatial bounds for particle generation"
+{
+  *:* :: ""
+} 1.0
+
+CCTK_REAL containing_bx_zmin "Spatial bounds for particle generation"
+{
+  *:* :: ""
+} -1.0
+
+CCTK_REAL containing_bx_zmax "Spatial bounds for particle generation"
+{
+  *:* :: ""
+} 1.0
+
diff --git a/TestNuPcsArdBH/src/testnupcsardbh.cxx b/TestNuPcsArdBH/src/testnupcsardbh.cxx
index 3813657..2ecc811 100644
--- a/TestNuPcsArdBH/src/testnupcsardbh.cxx
+++ b/TestNuPcsArdBH/src/testnupcsardbh.cxx
@@ -33,10 +33,11 @@ extern "C" void TestNuPcsArdBH_InitParticles(CCTK_ARGUMENTS) {
     NuParticleContainer::ParticleInitData pdata = {
         {}, {}, {0.0, 0.0, 0.0, 0.0, 0.0, 0.0}, {}};
 
-    // Constrain positions to the target cell near (5.05, 5.05, 0)
+    // Constrain positions to the target cell
     const auto dx = pc->Geom(lev).CellSizeArray();
-    RealBox containing_bx({5.0, 5.0, -0.5 * dx[2]},
-                          {5.0 + dx[0], 5.0 + dx[1], 0.5 * dx[2]});
+    RealBox containing_bx(
+        {containing_bx_xmin, containing_bx_ymin, containing_bx_zmin},
+        {containing_bx_xmax, containing_bx_ymax, containing_bx_zmax});
 
     pc->InitRandom(num_particles, random_seed, pdata, true, containing_bx);
 
diff --git a/TestNuPcsArdBH/test/nupcsardbh.par b/TestNuPcsArdBH/test/nupcsardbh.par
index aadec7e..7fd7b6a 100644
--- a/TestNuPcsArdBH/test/nupcsardbh.par
+++ b/TestNuPcsArdBH/test/nupcsardbh.par
@@ -53,6 +53,12 @@ SingleBH::mass = 0.0
 
 TestNuPcsArdBH::num_particles = 512
 TestNuPcsArdBH::random_seed = 42
+TestNuPcsArdBH::containing_bx_xmin = +5.0
+TestNuPcsArdBH::containing_bx_xmax = +5.0 + 0.1
+TestNuPcsArdBH::containing_bx_ymin = +5.0
+TestNuPcsArdBH::containing_bx_ymax = +5.0 + 0.1
+TestNuPcsArdBH::containing_bx_zmin = -0.05
+TestNuPcsArdBH::containing_bx_zmax = +0.05
 
 IO::out_dir = $parfile
 IO::out_every = 32