You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
**************************************************************************************************************
process.plasma_profiles - WARNING - Parabolic plasma profiles is used for an L-Mode plasma, but the physics variables do not describe an L-Mode plasma. 'rhopedt', 'rhopedn', 'teped', 'tesep', 'neped', 'nesep', and 'tbeta' have all been reset to L-Mode appropriate values
/Users/timothynunn/process/process/impurity_radiation.py:398: RuntimeWarning: divide by zero encountered in log
zav_of_te = yi + c * (numpy.log(tprofile) - xi)
/Users/timothynunn/process/process/impurity_radiation.py:456: RuntimeWarning: divide by zero encountered in log
pimpden = numpy.exp(yi + c * (numpy.log(tprofile) - xi))
6 | Convergence Parameter: 1.037E-07
************************************* PROCESS found a feasible solution **************************************
/Users/timothynunn/process/process/physics.py:2258: RuntimeWarning: divide by zero encountered in scalar divide
1
/Users/timothynunn/process/process/physics.py:1743: RuntimeWarning: invalid value encountered in divide
] = (
/Users/timothynunn/process/process/physics.py:1753: RuntimeWarning: invalid value encountered in divide
physics_variables.fhe3 * physics_variables.deni / physics_variables.dene
/Users/timothynunn/process/process/physics.py:1807: RuntimeWarning: divide by zero encountered in log
- (np.log(physics_variables.dene) / 2.0)
/Users/timothynunn/process/process/physics.py:1808: RuntimeWarning: divide by zero encountered in log
+ np.log(physics_variables.te * 1000.0)
/Users/timothynunn/process/process/physics.py:1806: RuntimeWarning: invalid value encountered in scalar add
31.0
/Users/timothynunn/process/process/physics.py:1812: RuntimeWarning: divide by zero encountered in log
- (np.log(physics_variables.dene) / 2.0)
/Users/timothynunn/process/process/physics.py:1813: RuntimeWarning: divide by zero encountered in log
+ np.log(physics_variables.te * 1000.0)
/Users/timothynunn/process/process/physics.py:1811: RuntimeWarning: invalid value encountered in scalar add
31.3
/Users/timothynunn/process/process/physics.py:1863: RuntimeWarning: invalid value encountered in divide
physics_variables.aion = physics_variables.aion / physics_variables.dnitot
/Users/timothynunn/process/process/physics.py:1866: RuntimeWarning: invalid value encountered in divide
physics_variables.zeffai = (
/Users/timothynunn/process/process/profiles.py:310: RuntimeWarning: invalid value encountered in divide
physics_variables.ti / physics_variables.te * physics_variables.te0
/Users/timothynunn/process/process/profiles.py:178: RuntimeWarning: invalid value encountered in divide
physics_variables.dnitot / physics_variables.dene * physics_variables.ne0
/Users/timothynunn/process/process/plasma_profiles.py:281: RuntimeWarning: invalid value encountered in scalar divide
-dtdrho_max * physics_variables.rminor * rho_te_max / te_max
/Users/timothynunn/process/process/plasma_profiles.py:285: RuntimeWarning: invalid value encountered in scalar divide
-dndrho_max * physics_variables.rminor * rho_ne_max / ne_max
Traceback (most recent call last):
File "/Users/timothynunn/process/.venv/bin/process", line 33, in <module>
sys.exit(load_entry_point('process', 'console_scripts', 'process')())
File "/Users/timothynunn/process/process/main.py", line 657, in main
Process(args)
File "/Users/timothynunn/process/process/main.py", line 127, in __init__
self.run_mode()
File "/Users/timothynunn/process/process/main.py", line 204, in run_mode
self.run.run()
File "/Users/timothynunn/process/process/main.py", line 362, in run
self.run_scan(self.solver)
File "/Users/timothynunn/process/process/main.py", line 465, in run_scan
self.scan = Scan(self.models, solver)
File "/Users/timothynunn/process/process/scan.py", line 22, in __init__
self.run_scan()
File "/Users/timothynunn/process/process/scan.py", line 37, in run_scan
final.finalise(self.models, ifail)
File "/Users/timothynunn/process/process/final.py", line 23, in finalise
op.write(models, ft.constants.nout)
File "/Users/timothynunn/process/process/output.py", line 22, in write
models.stellarator.run(output=True)
File "/Users/timothynunn/process/process/stellarator.py", line 121, in run
self.stopt(True)
File "/Users/timothynunn/process/process/stellarator.py", line 354, in stopt
self.stopt_output(
File "/Users/timothynunn/process/process/stellarator.py", line 3235, in stopt_output
powerht_local, pscalingmw_local = self.power_at_ignition_point(
File "/Users/timothynunn/process/process/stellarator.py", line 3285, in power_at_ignition_point
self.stphys(False)
File "/Users/timothynunn/process/process/stellarator.py", line 4065, in stphys
physics_variables.piepv = rether(
ZeroDivisionError: division by zero
Environment
OS: MacOS Ventura 13.0.1
Version: 3.1.0
The text was updated successfully, but these errors were encountered:
Describe the bug
Running the HELIAS-5B input file results in a divide by 0 error after a feasible solution is found.
Steps to reproduce
FELIX.IN.DAT.txt
Expected behaviour
MFile is written once feasible solution is found.
Evidence
Environment
The text was updated successfully, but these errors were encountered: