freegs-plasma / freegs

Free boundary Grad-Shafranov solver

http://freegs.readthedocs.io/en/latest/

`core_mask` fails for `psi_bndr = 0.0`

kripnerl opened this issue 4 years ago · comments

Lukáš Kripner commented 4 years ago

This condition is satisfied for psi_bndr = 0, or 0.0 ...

https://github.com/bendudson/freegs/blob/87bc845a8c1f4f48249c6fd930bc97b436468dd2/freegs/critical.py#L268

I would replace it by
if psi_bndry is None: