diff --git a/Source/TimeIntegration/ERF_SlowRhsPost.cpp b/Source/TimeIntegration/ERF_SlowRhsPost.cpp index 7f0ddc0e7..269740ade 100644 --- a/Source/TimeIntegration/ERF_SlowRhsPost.cpp +++ b/Source/TimeIntegration/ERF_SlowRhsPost.cpp @@ -106,7 +106,7 @@ void erf_slow_rhs_post (int level, int finest_level, TurbChoice tc = solverChoice.turbChoice[level]; const MultiFab* t_mean_mf = nullptr; - if (most) t_mean_mf = most->get_mac_avg(0,2); + if (most) t_mean_mf = most->get_mac_avg(level,2); const bool l_use_terrain = (solverChoice.mesh_type != MeshType::ConstantDz); const bool l_moving_terrain = (solverChoice.terrain_type == TerrainType::MovingFittedMesh); diff --git a/Source/TimeIntegration/ERF_SlowRhsPre.cpp b/Source/TimeIntegration/ERF_SlowRhsPre.cpp index aa05a41e3..c28cf84d9 100644 --- a/Source/TimeIntegration/ERF_SlowRhsPre.cpp +++ b/Source/TimeIntegration/ERF_SlowRhsPre.cpp @@ -122,7 +122,7 @@ void erf_slow_rhs_pre (int level, int finest_level, TurbChoice tc = solverChoice.turbChoice[level]; const MultiFab* t_mean_mf = nullptr; - if (most) t_mean_mf = most->get_mac_avg(0,2); + if (most) t_mean_mf = most->get_mac_avg(level,2); int start_comp = 0; int num_comp = 2;