From 4850465f8602f133d820cca2ad799dc5decce785 Mon Sep 17 00:00:00 2001 From: Ben Prather Date: Tue, 7 Jan 2025 17:32:43 -0500 Subject: [PATCH] Missed some Get() fixes. Was not a bug, but nice if you want timesteps on the fly --- kharma/grmhd/grmhd.cpp | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/kharma/grmhd/grmhd.cpp b/kharma/grmhd/grmhd.cpp index 831070ec..81e5060d 100644 --- a/kharma/grmhd/grmhd.cpp +++ b/kharma/grmhd/grmhd.cpp @@ -297,7 +297,7 @@ Real EstimateTimestep(MeshData *md) // TODO maybe split normal, ISMR timesteps? Excised pole/recalculated ctop too? double min_ndt = std::numeric_limits::max(); for (auto &pmb : pmesh->block_list) { - auto rc = pmb->meshblock_data.Get().get(); + auto rc = pmb->meshblock_data.Get(md->StageName()).get(); // We only need this block-wise to check boundary flags for ISMR, could special-case that const bool polar_inner_x2 = pmb->boundary_flag[BoundaryFace::inner_x2] == BoundaryFlag::user; const bool polar_outer_x2 = pmb->boundary_flag[BoundaryFace::outer_x2] == BoundaryFlag::user; @@ -629,7 +629,7 @@ void UpdateAveragedCtop(MeshData *md) auto pmesh = md->GetMeshPointer(); auto& params = pmesh->packages.Get("Boundaries")->AllParams(); for (auto &pmb : pmesh->block_list) { - auto &rc = pmb->meshblock_data.Get(); + auto &rc = pmb->meshblock_data.Get(md->StageName()); for (int i = 0; i < BOUNDARY_NFACES; i++) { BoundaryFace bface = (BoundaryFace)i; auto bname = KBoundaries::BoundaryName(bface);