Skip to content

Commit

Permalink
move override sync for velocities to before ubar/vbar are calculated
Browse files Browse the repository at this point in the history
  • Loading branch information
hklion committed Aug 19, 2024
1 parent b7f6019 commit 4c189a6
Showing 1 changed file with 6 additions and 5 deletions.
11 changes: 6 additions & 5 deletions Source/REMORA.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down

0 comments on commit 4c189a6

Please sign in to comment.