diff --git a/Source/ERF.cpp b/Source/ERF.cpp index b61ef8cba..10bbaf566 100644 --- a/Source/ERF.cpp +++ b/Source/ERF.cpp @@ -722,8 +722,10 @@ ERF::InitData_post () restart(); // Create the physbc objects for {cons, u, v, w, base state} + // We fill the additional base state ghost cells just in case we have read the old format for (int lev(0); lev <= max_level; ++lev) { make_physbcs(lev); + (*physbcs_base[lev])(base_state[lev],0,base_state[lev].nComp(),base_state[lev].nGrowVect()); } } diff --git a/Source/IO/ERF_Checkpoint.cpp b/Source/IO/ERF_Checkpoint.cpp index 55b1ab710..92b8b6438 100644 --- a/Source/IO/ERF_Checkpoint.cpp +++ b/Source/IO/ERF_Checkpoint.cpp @@ -415,7 +415,6 @@ ERF::ReadCheckpointFile () }); } } - (*physbcs_base[lev])(base_state[lev],0,base_state[lev].nComp(),base_state[lev].nGrowVect()); base_state[lev].FillBoundary(geom[lev].periodicity()); if (solverChoice.use_terrain) {