Skip to content

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

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
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions