Skip to content

[Bug]: Tokhoku[2,3]d: cannot generate python inputs #11

@scivision

Description

@scivision

What happened?

equilibrium simulations tohoku[2,3]d fail to setup input with PyGemini.
by copy pasting the command, this is also seen on a laptop.

python3  -m gemini3d.model cfg/equilibrium/tohoku2d_eq /tmp/tohoku2d_eq

Relevant log output

23/154 Test  #73: setup:python:tohoku3d_eq ............................***Failed   19.62 sec
Traceback (most recent call last):
  File "/miniconda3/lib/python3.10/runpy.py", line 196, in _run_module_as_main
    return _run_code(code, main_globals, None,
  File "/miniconda3/lib/python3.10/runpy.py", line 86, in _run_code
    exec(code, run_globals)
  File "/code/pygemini/src/gemini3d/model.py", line 220, in <module>
    cli()
  File "/code/pygemini/src/gemini3d/model.py", line 216, in cli
    setup(P.config_file, P.out_dir, P.root)
  File "/code/pygemini/src/gemini3d/model.py", line 143, in setup
    equilibrium(cfg)
  File "/code/pygemini/src/gemini3d/model.py", line 159, in equilibrium
    dat = equilibrium_state(cfg, xg)
  File "/code/pygemini/src/gemini3d/plasma.py", line 423, in equilibrium_state
    ns = np.concatenate((ns, ns[:, lx1, :, :], ns[:, ::-1, :, :]), 1)
IndexError: index 95 is out of bounds for axis 1 with size 95
 11/154 Test  #64: setup:python:tohoku2d_eq ............................***Failed    2.60 sec
Traceback (most recent call last):
  File "/miniconda3/lib/python3.10/runpy.py", line 196, in _run_module_as_main
    return _run_code(code, main_globals, None,
  File "/miniconda3/lib/python3.10/runpy.py", line 86, in _run_code
    exec(code, run_globals)
  File "/code/pygemini/src/gemini3d/model.py", line 220, in <module>
    cli()
  File "/code/pygemini/src/gemini3d/model.py", line 216, in cli
    setup(P.config_file, P.out_dir, P.root)
  File "/code/pygemini/src/gemini3d/model.py", line 143, in setup
    equilibrium(cfg)
  File "/code/pygemini/src/gemini3d/model.py", line 159, in equilibrium
    dat = equilibrium_state(cfg, xg)
  File "/code/pygemini/src/gemini3d/plasma.py", line 423, in equilibrium_state
    ns = np.concatenate((ns, ns[:, lx1, :, :], ns[:, ::-1, :, :]), 1)
IndexError: index 127 is out of bounds for axis 1 with size 127

Metadata

Metadata

Assignees

No one assigned

    Labels

    bugSomething isn't working

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions