initialize water saturation for gas dissolved in water

This commit is contained in:
jakobtorben 2024-07-05 10:23:13 +02:00
parent db0b4e1e6d
commit 5eac0880a9

View File

@ -584,6 +584,8 @@ public:
sw = (*this)[waterSwitchIdx];
if (primaryVarsMeaningGas() == GasMeaning::Sg)
sg = (*this)[compositionSwitchIdx];
if (primaryVarsMeaningWater() == WaterMeaning::Rsw)
sw = 1.0;
if (primaryVarsMeaningGas() == GasMeaning::Disabled && gasEnabled)
sg = 1.0 - sw; // water + gas case