FillPatch before tagging cells for regridding (#2086)

erf-model · Jan 24, 2025 · ce96f80 · ce96f80
1 parent 202be36
commit ce96f80
Show file tree

Hide file tree

Showing 2 changed files with 20 additions and 1 deletion.
diff --git a/Source/ERF_Tagging.cpp b/Source/ERF_Tagging.cpp
@@ -17,6 +17,25 @@ ERF::ErrorEst (int levc, TagBoxArray& tags, Real time, int /*ngrow*/)
     const int clearval = TagBox::CLEAR;
     const int   tagval = TagBox::SET;
 
+    //
+    // Make sure the ghost cells of the level we are tagging at are filled
+    //    in case we take differences that require them
+    // NOTE: We are Fillpatching only the cell-centered variables here
+    //
+    MultiFab& S_new = vars_new[levc][Vars::cons];
+    MultiFab& U_new = vars_new[levc][Vars::xvel];
+    MultiFab& V_new = vars_new[levc][Vars::yvel];
+    MultiFab& W_new = vars_new[levc][Vars::zvel];
+    //
+    if (levc == 0) {
+        FillPatch(levc, time, {&S_new, &U_new, &V_new, &W_new});
+    } else {
+        FillPatch(levc, time, {&S_new, &U_new, &V_new, &W_new},
+                              {&S_new, &rU_new[levc], &rV_new[levc], &rW_new[levc]},
+                              base_state[levc], base_state[levc],
+                              false, true);
+    }
+
     for (int j=0; j < ref_tags.size(); ++j)
     {
         //

diff --git a/Source/SourceTerms/ERF_MakeSources.cpp b/Source/SourceTerms/ERF_MakeSources.cpp
@@ -424,7 +424,7 @@ void make_sources (int level,
             {
                 const Real t_blank = t_blank_arr(i, j, k);
                 cell_src(i, j, k, RhoTheta_comp) -= t_blank * CdT * U_s * (cell_data(i,j,k,RhoTheta_comp) - dptr_t_plane(k));
-                cell_src(i, j, k, Rho_comp)      -= t_blank * CdT * U_s * (cell_data(i,j,k,Rho_comp) - dptr_r_plane(k));                 
+                cell_src(i, j, k, Rho_comp)      -= t_blank * CdT * U_s * (cell_data(i,j,k,Rho_comp) - dptr_r_plane(k));
             });
         }