Generate LineEstimate Spectrum for a given radial range
It can be useful to just plot a line estimate spectrum for a given radial bin ... e.g. from 0 to 1 au.
One could use the flux column from the rInfo for that. This is for one line
linest=m1.getLineEstimate("CO",1300)
# this is the flux for each x column
fcolumn = [x.Fcolumn for x in linest.rInfo]
# simply summing them up gives the total flux
print(linest.flux,np.sum(fcolumn))
# this should then be the flax between 1 and 10 au, taking the closest x values to 1 and 10 au
print(np.sum(fcolumn[np.argmin(np.abs(m1.x[:,0]-1)):np.argmin(np.abs(m1.x[:,0]-10))]))It could look something like this (from Jelke, he made a kind of hack for that).
