Missed some Get() fixes. Was not a bug, but nice if you want timestep…

…s on the fly

kharma/grmhd/grmhd.cpp

@@ -297,7 +297,7 @@ Real EstimateTimestep(MeshData<Real> *md)
     // TODO maybe split normal, ISMR timesteps? Excised pole/recalculated ctop too?
     double min_ndt = std::numeric_limits<double>::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<Real> *md)
     auto pmesh = md->GetMeshPointer();
     auto& params = pmesh->packages.Get<KHARMAPackage>("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);