Finding spectrum line flux via scipy trapz numerical integration returns wrong flux

Question

I am working on extracting spectral activity indices from stars. To do this, I look at the spectra of the star and calculate necessary line fluxes. An example of a spectrum