Processing stops when usine macro 'co2pres'
Nobody has claimed this yet.
- Dominant language
- Python
- Stars
- 8
- Forks
- 10
- PR merge metrics
- No merged PRs in 30d
Description
Hello,
I have an issue using the carbon module when I want to have the co2pres macro initialized : the FEHM.exe stops after the first time step (without any error message, the python code then keeps running).
This is in a code trying to make a small co2 injection model inspired from the tutorial.
If I did adapt the library to my newest 3.11 Python version (Windows) by solving the vocabulary issues, the carbon module not working, this problem doesn't have an explicit way to debug.
Because of the early stop in the .exe, the error reported by Python is about the self.times variable that is empty in this case when trying to do a cont.slice_plot(........).
Here you can see a few parts of my code after the initializations that I skip here.
dat.add(fcarb(iprtype = 4))
dat.carb.on(iprtype = 4)
dat.files.co2in='C:\\Users\\xxxxx\\Documents\\xxxxxxxxxxxxxxxx\\PyFEHM\\co2_interp_table.txt
dat.add(fmacro('pres',zone=0,param=(('pressure',19),('temperature',(50)),('saturation',0.))))
dat.add(fmacro('co2pres',zone=0,param=(('pressure',19),('temperature',(50)),('phase',4))))
dat.add (fmacro('co2frac', zone='reservoir', param=(('water_rich_sat',0.0), ('co2_rich_sat',1.0 ), ('co2_mass_frac',1.0), ('init_salt_conc',0.), ('override_flag',1))))
dat.add (fmacro('co2flow', zone='wellpoint', param=(('rate',1), ('energy',-40), ('impedance',1.e-2), ('bc_flag',1))))
In general, wether I comment or uncoment the other macros doesn't change anything, it really is happening when using co2pres macro, which is necessary to have the P & T initialized (I think)
I don't know if this issue is well known, but thank you in advance for your help.
Contributor guide
No contributing guide indexed for this repository
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
Research direction
The report names the co2pres macro, FEHM.exe, cont.slice_plot, and co2_interp_table.txt, but omits repository files or a complete reproducer. Start by reducing the setup to the carbon initialization and co2pres call, then inspect the FEHM.exe run and Python state after the first time step. Done means the failure is reproducible with a clear diagnostic or the co2pres case completes without leaving self.times empty.
Written by the indexing model from the issue text.
Assessment
- Tech stack
- python
- Domain
- backend
- Issue type
- Bug
- Difficulty
- 4/5
- Estimated time
- 3-5 days
- Activity status
- Stale
- Clarity
- Needs clarification
- Newbie friendliness
- 30/100