Merge pull request #374 from thetisproject/fix-nonhydrostatic-adaptivity

Fix NonHydrostaticTimeIntegrator2D for projection of solutions with adaptivity

thetis/coupled_timeintegrator_2d.py

@@ -199,11 +199,12 @@ def initialize(self, solution2d):
         self.timesteppers.swe2d.initialize(self.fields.solution_2d)
         if self.nh_options.update_free_surface:
             self.timesteppers.fs2d.initialize(self.fields.elev_2d)
-            self.elev_old.assign(self.fields.elev_2d)
 
     @PETSc.Log.EventDecorator("thetis.NonHydrostaticTimeIntegrator2D.advance")
     def advance(self, t, update_forcings=None):
         """Advances equations for one time step."""
+        if self.nh_options.update_free_surface:
+            self.elev_old.assign(self.fields.elev_2d)
         if self.serial_advancing:
             # --- advance in serial ---
             self.timesteppers.swe2d.advance(t, update_forcings=update_forcings)
@@ -213,7 +214,6 @@ def advance(self, t, update_forcings=None):
             if self.nh_options.update_free_surface:
                 self.fields.elev_2d.assign(self.elev_old)
                 self.timesteppers.fs2d.advance(t, update_forcings=update_forcings)
-                self.elev_old.assign(self.fields.elev_2d)
             # update old solution
             if self.options.swe_timestepper_type == 'SSPIMEX':
                 self.timesteppers.swe2d.erk.solution_old.assign(self.fields.solution_2d)
@@ -234,4 +234,3 @@ def advance(self, t, update_forcings=None):
                     elif last_stage:
                         self.fields.elev_2d.assign(self.elev_old)
                         self.timesteppers.fs2d.advance(t, update_forcings=update_forcings)
-                        self.elev_old.assign(self.fields.elev_2d)