finite_abundances.py

gistfile1.txt

In [30]: ok = image[1].data["physical"]

In [31]: elements = [_ for _ in image[1].data.dtype.names if _.endswith("_H") and _ not in ("SRC_H", "PARAM_M_H")]

In [32]: for element in elements:
    print(element, np.sum(np.isfinite(image[1].data[element][ok])))
   ....:     
('AL_H', 99127)
('CA_H', 99127)
('C_H', 99127)
('FE_H', 99127)
('K_H', 99127)
('MG_H', 99127)
('MN_H', 99127)
('NA_H', 99127)
('NI_H', 99127)
('N_H', 99127)
('O_H', 99127)
('SI_H', 99127)
('S_H', 99127)
('TI_H', 99127)
('V_H', 99127)