It looks like in the `psc_flatfoil_yz_small` case, both `rho_e` and `rho_i` are negative -- obviously, `rho_i` should not.