fix restart

erf-model · Oct 25, 2024 · 5bf6ed5 · 5bf6ed5
1 parent dce97f7
commit 5bf6ed5
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 1 deletion.
diff --git 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
@@ -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)  {
-Original file line number
+Diff line change
@@ Expand Up / @@ -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)  {
@@ Expand Down @@