adiabatic InitialProfile will only run once now

geodynamics · Jun 1, 2014 · d6e745c · d6e745c
1 parent b872a6c
commit d6e745c
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/source/adiabatic_conditions/initial_profile.cc b/source/adiabatic_conditions/initial_profile.cc
@@ -44,6 +44,8 @@ namespace aspect
     void
     InitialProfile<dim>::initialize()
     {
+      if (initialized)
+        return;
       delta_z = this->get_geometry_model().maximal_depth() / (n_points-1);
 
       const unsigned int n_compositional_fields = this->n_compositional_fields();