check should not be called for hybrid

src/coreComponents/physicsSolvers/fluidFlow/CompositionalMultiphaseFVM.cpp

@@ -162,6 +162,8 @@ void CompositionalMultiphaseFVM::initializePreSubGroups()
   m_linearSolverParameters.get().mgr.strategy = m_isThermal
                                                 ? LinearSolverParameters::MGR::StrategyType::thermalCompositionalMultiphaseFVM
                                                 : LinearSolverParameters::MGR::StrategyType::compositionalMultiphaseFVM;
+
+  checkDiscretizationName();
 }
 
 void CompositionalMultiphaseFVM::setupDofs( DomainPartition const & domain,

src/coreComponents/physicsSolvers/fluidFlow/FlowSolverBase.cpp

@@ -365,10 +365,17 @@ void FlowSolverBase::initializePreSubGroups()
       }
     } );
   }
-  else
+}
+
+void FlowSolverBase::checkDiscretizationName()
+{
+  DomainPartition & domain = this->getGroupByPath< DomainPartition >( "/Problem/domain" );
+  NumericalMethodsManager const & numericalMethodManager = domain.getNumericalMethodManager();
+  FiniteVolumeManager const & fvManager = numericalMethodManager.getFiniteVolumeManager();
+  if( !fvManager.hasGroup< FluxApproximationBase >( m_discretizationName ) )
   {
     GEOS_ERROR( GEOS_FMT( "{}: can not find discretization named '{}' (a discretization deriving from FluxApproximationBase must be selected for {} solver '{}' )",
-                          getDataContext(), m_discretizationName, getCatalogName(), getName() ) );
+                          getDataContext(), m_discretizationName, getCatalogName(), getName()));
   }
 }
 

src/coreComponents/physicsSolvers/fluidFlow/FlowSolverBase.hpp

@@ -208,6 +208,8 @@ class FlowSolverBase : public PhysicsSolverBase
 
   virtual void initializePreSubGroups() override;
 
+  void checkDiscretizationName();
+
   virtual void initializePostInitialConditionsPreSubGroups() override;
 
   void initializeState( DomainPartition & domain );

src/coreComponents/physicsSolvers/fluidFlow/ReactiveCompositionalMultiphaseOBL.cpp

@@ -121,6 +121,11 @@ ReactiveCompositionalMultiphaseOBL::ReactiveCompositionalMultiphaseOBL( const st
   m_linearSolverParameters.get().mgr.strategy = LinearSolverParameters::MGR::StrategyType::reactiveCompositionalMultiphaseOBL;
 }
 
+void ReactiveCompositionalMultiphaseOBL::initializePreSubGroups()
+{
+  checkDiscretizationName();
+}
+
 void ReactiveCompositionalMultiphaseOBL::setupDofs( DomainPartition const & domain,
                                                     DofManager & dofManager ) const
 {

src/coreComponents/physicsSolvers/fluidFlow/ReactiveCompositionalMultiphaseOBL.hpp

@@ -306,6 +306,8 @@ class ReactiveCompositionalMultiphaseOBL : public FlowSolverBase
 
   virtual void initializePostInitialConditionsPreSubGroups() override;
 
+  virtual void initializePreSubGroups() override;
+
 
 private:
 

src/coreComponents/physicsSolvers/fluidFlow/SinglePhaseFVM.cpp

@@ -68,6 +68,8 @@ void SinglePhaseFVM< BASE >::initializePreSubGroups()
 {
   BASE::initializePreSubGroups();
 
+  checkDiscretizationName();
+
   if( m_isThermal )
   {
     // For thermal simulations 2 pdes are considered so we let AMG know.