Skip to content

fix: single phase statistics CSV generation #3928

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
arng40 wants to merge 9 commits into
base: develop
Choose a base branch
from
+16 −13
Open

fix: single phase statistics CSV generation #3928

Show file tree
Hide file tree
Changes from 5 commits
Commits
Show all changes
9 commits
Select commit Hold shift + click to select a range
2e8e411
fix singlePhaseStats
arng40 Dec 22, 2025
850946c
Merge branch 'develop' into fix/dudes/csv-generation
arng40 Dec 22, 2025
c73a68e
xsd
arng40 Dec 22, 2025
3ce7eaa
Merge branch 'fix/dudes/csv-generation' of https://github.com/GEOS-DE…
arng40 Dec 22, 2025
8014180
use table fonctions
arng40 Dec 22, 2025
d36de97
Merge branch 'develop' into fix/dudes/csv-generation
MelReyCG Dec 30, 2025
73fefb5
Merge remote-tracking branch 'origin/develop' into fix/dudes/csv-gene…
arng40 Jan 7, 2026
c09be58
Merge branch 'fix/dudes/csv-generation' of https://github.com/GEOS-DE…
arng40 Jan 7, 2026
e51d89a
Merge branch 'develop' into fix/dudes/csv-generation
MelReyCG Jan 26, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
29 changes: 16 additions & 13 deletions src/coreComponents/physicsSolvers/fluidFlow/SinglePhaseStatistics.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -70,20 +70,23 @@ void SinglePhaseStatistics::registerDataOnMesh( Group & meshBodies )
// write output header
if( m_writeCSV > 0 && MpiWrapper::commRank() == 0 )
{
TableLayout tableLayout( {
TableLayout::Column().setName( GEOS_FMT( "Time [{}]", units::getSymbol( units::Unit::Time ))),
TableLayout::Column().setName( GEOS_FMT( "Min pressure [{}]", units::getSymbol( units::Unit::Pressure ))),
TableLayout::Column().setName( GEOS_FMT( "Average pressure [{}]", units::getSymbol( units::Unit::Pressure ))),
TableLayout::Column().setName( GEOS_FMT( "Max pressure [{}]", units::getSymbol( units::Unit::Pressure ))),
TableLayout::Column().setName( GEOS_FMT( "Min delta pressure [{}]", units::getSymbol( units::Unit::Pressure ))),
TableLayout::Column().setName( GEOS_FMT( "Max delta pressure [{}]", units::getSymbol( units::Unit::Pressure )) ),
TableLayout::Column().setName( GEOS_FMT( "Min temperature [{}]", units::getSymbol( units::Unit::Temperature ))),
TableLayout::Column().setName( GEOS_FMT( "Average temperature [{}]", units::getSymbol( units::Unit::Temperature ))),
TableLayout::Column().setName( GEOS_FMT( "Max temperature [{}]", units::getSymbol( units::Unit::Temperature ))),
TableLayout::Column().setName( GEOS_FMT( "Total dynamic pore volume [{}]", units::getSymbol( units::Unit::ReservoirVolume ) )),
TableLayout::Column().setName( GEOS_FMT( "Total fluid mass [{}]", units::getSymbol( units::Unit::Mass )))
} );

TableCSVFormatter csvFormatter( tableLayout );
std::ofstream outputFile( m_outputDir + "/" + regionNames[i] + ".csv" );
outputFile <<
GEOS_FMT( "Time [{}]", units::getSymbol( units::Unit::Time ))<<
GEOS_FMT( "Average pressure [{}]", units::getSymbol( units::Unit::Pressure ))<<
GEOS_FMT( "Min pressure [{}]", units::getSymbol( units::Unit::Pressure ))<<
GEOS_FMT( "Max pressure [{}]", units::getSymbol( units::Unit::Pressure ))<<
GEOS_FMT( "Min delta pressure [{}]", units::getSymbol( units::Unit::Pressure ))<<
GEOS_FMT( "Max delta pressure [{}]", units::getSymbol( units::Unit::Pressure )) <<
GEOS_FMT( "Min temperature [{}]", units::getSymbol( units::Unit::Temperature ))<<
GEOS_FMT( "Average temperature [{}]", units::getSymbol( units::Unit::Temperature ))<<
GEOS_FMT( "Max temperature [{}]", units::getSymbol( units::Unit::Temperature ))<<
GEOS_FMT( "Total dynamic pore volume [{}]", units::getSymbol( units::Unit::ReservoirVolume ) ) <<
GEOS_FMT( "Total fluid mass [{}]", units::getSymbol( units::Unit::Mass ));
outputFile << std::endl;
outputFile << csvFormatter.headerToString();
outputFile.close();
}
}
Expand Down
7 changes: 4 additions & 3 deletions src/coreComponents/schema/schema.xsd
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3655,8 +3655,9 @@ Information output from lower logLevels is added with the desired log level
<!--referenceElevation => Reference elevation where BHP control is enforced [m]-->
<xsd:attribute name="referenceElevation" type="real64" use="required" />
<!--referenceReservoirRegion => Name of reservoir region used for obtaining average region pressure used in volume rate constraint calculations.
Frequency of pressure update is set in Single/CompositionalMultiPhaseStatistics definition.
Setting cycleFrequency='1' will update the pressure every timestep, note that is a lagged property in constraint properties-->
Frequency of pressure update is set in SinglePhase/CompositionalMultiphaseStatistics definition.
Setting cycleFrequency='1' will update the pressure every timestep, note that is a lagged property in constraint propertiesNote the event associated with the statists task must be entered before the solver event.
-->
<xsd:attribute name="referenceReservoirRegion" type="groupNameRef" default="" />
<!--statusTableName => Name of the well status table when the status of the well is a time dependent function.
If the status function evaluates to a positive value at the current time, the well will be open otherwise the well will be shut.-->
Expand Down Expand Up @@ -8403,7 +8404,7 @@ The expected format is "{ waterMax, oilMax }", in that order-->
<xsd:attribute name="phaseCapPressureExponentInv" type="real64_array" default="{0.5}" />
<!--phaseCapPressureMultiplier => Entry pressure value for each phase-->
<xsd:attribute name="phaseCapPressureMultiplier" type="real64_array" default="{1}" />
<!--phaseMinVolumeFraction => Minimum volume fraction value for each phase (between 0 and 1) used to compute the capillary pressure-->
<!--phaseMinVolumeFraction => Minimum volume fraction value for each phase-->
<xsd:attribute name="phaseMinVolumeFraction" type="real64_array" default="{0}" />
<!--phaseNames => List of fluid phases-->
<xsd:attribute name="phaseNames" type="groupNameRef_array" use="required" />
Expand Down
Loading