Addressing stability issues with introducing refinement levels on a r…

…estart (Exawind#942)

putting postsolve into PDE initialize and regrid
-- will remove redundancy of fillpatch operations in future PR

amr-wind/equation_systems/PDE.H

@@ -71,6 +71,9 @@ public:
             m_fields, m_sim.has_overset(), variable_density,
             m_sim.has_mesh_mapping()));
         m_src_op.init_source_terms(m_sim);
+
+        // Post-solve operations should also apply after initialization
+        m_post_solve_op(m_time.current_time());
     }
 
     //! Perform update actions after a regrid is performed
@@ -90,6 +93,9 @@ public:
         m_adv_op.reset(new AdvectionOp<PDE, Scheme>(
             m_fields, m_sim.has_overset(), variable_density,
             m_sim.has_mesh_mapping()));
+
+        // Post-solve operations should also apply after a regrid
+        m_post_solve_op(m_time.current_time());
     }
 
     //! Return the object holding the fields necessary for solving this PDE