Update SinglePhaseBase.cpp

GEOS-DEV · Jul 5, 2024 · 3839b4a · 3839b4a
1 parent fa97b6a
commit 3839b4a
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/coreComponents/physicsSolvers/fluidFlow/SinglePhaseBase.cpp b/src/coreComponents/physicsSolvers/fluidFlow/SinglePhaseBase.cpp
@@ -462,12 +462,12 @@ void SinglePhaseBase::initializePostInitialConditionsPreSubGroups()
     {
       region.forElementSubRegions< FaceElementSubRegion >( [&]( FaceElementSubRegion & subRegion )
       {
+        ConstitutiveBase & fluid = getConstitutiveModel( subRegion, subRegion.getReference< string >( viewKeyStruct::fluidNamesString() )  );
+        real64 const defaultDensity = getFluidProperties( fluid ).defaultDensity;
+
         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( defaultDensity * region.getDefaultAperture() );
       } );