Hi,
I am trying to use the `-conserveLin` in krome. However, when I run the code, sometimes the `krome_conserveLin_x` routine returned a singular matrix and crashed because the A(idx_atom_E, idx_atom_E) = 0. I am wondering whether this term should be removed from the matrix because you use another way to fix the electron abundance. Or if the term must be retained, should the value be set to 1 to work around the singularity?
The option also generated some variables like idx_atom__dust or idx_atom__ortho, I think this is a bug.
I will appreciate any suggestions.
Best regards,
Chia-Jung