gemini3d / pygemini

Python interface for Gemini3D
Apache License 2.0
8 stars 4 forks source link

Eq with gridflag=2 (inverted grid) causes all simulations using it to not work due to negative density #22

Open jdiazpena opened 1 year ago

jdiazpena commented 1 year ago

What happened?

For the SAID simulations, I have the following Eq config.nml

If I set gridflag=2 for an inverted grid instead of 0 (close dipole) running a full simulation does not work with pygemini

Using the following config.nml to do the run, and selecting the same gridflag as the EQ, it errors out with negative density check.

If the EQ is run with a gridflag=0 (close dipole), and then the full simulation is run with gridflag=2, it pygemini does not start the simulation since it detects negative density.

When using 0 grid flag on both of them it works just fine.