OPM
Restart Error
Input.zip Hello Happy New Year.
Please why am I getting this restart error
Total number of active cells: 3600 / total pore volume: 512599201 RB
===============Saturation Functions Diagnostics===============
System: Oil-Gas system.
Relative permeability input format: Saturation Family I.
Number of saturation regions: 1
Warning: PVT Table evaluation:Attempt to get tabulated value for (0.000000, 24209168.403533) on a table of extent 280.000000 to 400.000000 times 100000.000000 to 100000000.000000. Will use extrapolation
Warning: PVT Table evaluation:Attempt to get tabulated value for (0.000000, 24209168.403533) on a table of extent 280.000000 to 400.000000 times 100000.000000 to 100000000.000000. Will use extrapolation
Warning: PVT Table evaluation:Attempt to get tabulated value for (275.000000, 24209168.403533) on a table of extent 280.000000 to 400.000000 times 100000.000000 to 100000000.000000. Will use extrapolation
flow: /usr/include/opm/models/blackoil/blackoilintensivequantities.hh:482: void Opm::BlackOilIntensiveQuantities<TypeTag>::update(const ElementContext&, unsigned int, unsigned int) [with TypeTag = Opm::Properties::TTag::EclFlowGasOilProblem; Opm::BlackOilIntensiveQuantities<TypeTag>::ElementContext = Opm::FvBaseElementContext<Opm::Properties::TTag::EclFlowGasOilProblem>]: Assertion `isfinite(fluidState_.density(phaseIdx))' failed.
[LAPTOP-H5DAMARN:12564] *** Process received signal ***
[LAPTOP-H5DAMARN:12564] Signal: Aborted (6)
[LAPTOP-H5DAMARN:12564] Signal code: (-6)
Warning: PVT Table evaluation:Attempt to get tabulated value for (275.000000, 24209168.403533) on a table of extent 280.000000 to 400.000000 times 100000.000000 to 100000000.000000. Will use extrapolation
flow: /usr/include/opm/models/blackoil/blackoilintensivequantities.hh:482: void Opm::BlackOilIntensiveQuantities<TypeTag>::update(const ElementContext&, unsigned int, unsigned int) [with TypeTag = Opm::Properties::TTag::EclFlowGasOilProblem; Opm::BlackOilIntensiveQuantities<TypeTag>::ElementContext = Opm::FvBaseElementContext<Opm::Properties::TTag::EclFlowGasOilProblem>]: Assertion `isfinite(fluidState_.density(phaseIdx))' failed.
[LAPTOP-H5DAMARN:12564] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x42520)[0x7fdd9b40b520]
[LAPTOP-H5DAMARN:12564] [ 1] /lib/x86_64-linux-gnu/libc.so.6(pthread_kill+0x12c)[0x7fdd9b45f9fc]
[LAPTOP-H5DAMARN:12564] [ 2] /lib/x86_64-linux-gnu/libc.so.6(raise+0x16)[0x7fdd9b40b476]
[LAPTOP-H5DAMARN:12564] Aborted
I guess there has been a division by zero somewhere. This might be hard to find. If you only experience this when you restart and the model runs otherwise this might point to a bug.
I am not sure if or when somebody finds the time, but to help we might need the version of OPM you are using and the models to replicate as a starting point.
Thank you for bringing the issue of OPM version to our attention. It appears that in the 2023.10 release notes of opm (https://opm-project.org/?p=1621), the 'RPTRST' keyword now requires an additional parameter "TEMP" to function properly. This means that to ensure compatibility with newer versions, the correct syntax for the input file should include the 'TEMP' parameter, as shown below:
RPTRST 'BASIC=1' 'TEMP'/
We don't understand why the input file works without the "TEMP" keyword on older versions, but we are happy!!.