move override sync for velocities to before ubar/vbar are calculated

hklion · Aug 19, 2024 · 4c189a6 · 4c189a6
1 parent b7f6019
commit 4c189a6
Showing 1 changed file with 6 additions and 5 deletions.
diff --git a/Source/REMORA.cpp b/Source/REMORA.cpp
@@ -668,18 +668,19 @@ REMORA::init_only (int lev, Real time)
         amrex::Abort("Need to specify T init procedure");
     }
 
+    // Ensure that the face-based data are the same on both sides of a periodic domain.
+    // The data associated with the lower grid ID is considered the correct value.
+    xvel_new[lev]->OverrideSync(geom[lev].periodicity());
+    yvel_new[lev]->OverrideSync(geom[lev].periodicity());
+    zvel_new[lev]->OverrideSync(geom[lev].periodicity());
+
     set_2darrays(lev);
 
     init_set_vmix(lev);
     set_hmixcoef(lev);
     set_coriolis(lev);
     init_custom_smflux(geom[lev], time, *vec_sustr[lev], *vec_svstr[lev], solverChoice);
 
-    // Ensure that the face-based data are the same on both sides of a periodic domain.
-    // The data associated with the lower grid ID is considered the correct value.
-    xvel_new[lev]->OverrideSync(geom[lev].periodicity());
-    yvel_new[lev]->OverrideSync(geom[lev].periodicity());
-    zvel_new[lev]->OverrideSync(geom[lev].periodicity());
 }
 
 // read in some parameters from inputs file