revert and give up

GEOS-DEV · Jul 5, 2024 · fa97b6a · fa97b6a
1 parent 90923b9
commit fa97b6a
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 20 deletions.
diff --git a/src/coreComponents/physicsSolvers/fluidFlow/SinglePhaseBase.cpp b/src/coreComponents/physicsSolvers/fluidFlow/SinglePhaseBase.cpp
@@ -465,8 +465,11 @@ void SinglePhaseBase::initializePostInitialConditionsPreSubGroups()
         subRegion.getWrapper< real64_array >( fields::flow::hydraulicAperture::key() ).
           setApplyDefaultValue( region.getDefaultAperture() );
 
+        SingleFluidBase & fluid =
+          getConstitutiveModel< SingleFluidBase >( subRegion, subRegion.getReference< string >( FlowSolverBase::viewKeyStruct::fluidNamesString() ) );
+        real64 const defaultDensity = fluid.defaultDensity();
         subRegion.getWrapper< real64_array >( FaceElementSubRegion::viewKeyStruct::creationMassString() ).
-          setApplyDefaultValue( 0.0 );
+          setApplyDefaultValue( defaultDensity * region.getDefaultAperture() );
       } );
     } );
 

diff --git a/src/coreComponents/physicsSolvers/multiphysics/HydrofractureSolver.cpp b/src/coreComponents/physicsSolvers/multiphysics/HydrofractureSolver.cpp
@@ -129,25 +129,6 @@ void HydrofractureSolver< POROMECHANICS_SOLVER >::initializePostInitialCondition
 {
   Base::initializePostInitialConditionsPreSubGroups();
 
-  DomainPartition & domain = this->template getGroupByPath< DomainPartition >( "/Problem/domain" );
-  forDiscretizationOnMeshTargets( domain.getMeshBodies(), [&]( string const &,
-                                                               MeshLevel & mesh,
-                                                               arrayView1d< string const > const & regionNames )
-  {
-    mesh.getElemManager().forElementRegions< SurfaceElementRegion >( regionNames,
-                                                                     [&]( localIndex const,
-                                                                          SurfaceElementRegion & region )
-    {
-      region.forElementSubRegions< FaceElementSubRegion >( [&]( FaceElementSubRegion & subRegion )
-      {
-        SingleFluidBase & fluid =
-          this->template getConstitutiveModel< SingleFluidBase >( subRegion, subRegion.getReference< string >( FlowSolverBase::viewKeyStruct::fluidNamesString() ) );
-        real64 const defaultDensity = fluid.defaultDensity();
-        subRegion.getWrapper< real64_array >( FaceElementSubRegion::viewKeyStruct::creationMassString() ).
-          setApplyDefaultValue( defaultDensity * region.getDefaultAperture() );
-      } );
-    } );
-  } );
 }
 
 template< typename POROMECHANICS_SOLVER >