Wave_100 = linspace(Wave[ind3], Wave[ind4], 100, endpoint=True)
continuum_g = Dm.Continuum_Gradient * Wave_g + Dm.Continuum_n
continuum_100 = Dm.Continuum_Gradient * Wave_100 + Dm.Continuum_n
#------------------------------------Calculations------------------------------------------------------
Dm.Num = 1
Dm.Deblend_Check = False
Dm.ind1, Dm.ind2, Dm.ind3, Dm.ind4, Dm.ind5, Dm.ind6 = ind1, ind2, ind3, ind4, ind5, ind6
Dm.Selections = [Wave_1, Wave_2, Wave_3, Wave_4, Wave_5, Wave_6]
#Calculation and plotting of continuum properties
Dm.LocalMedian, Dm.ContinuumFlux, Dm.SigmaContinuum, Dm.Continuum_Gradient, Dm.Continuum_n, WPoint1, Wpoint2, FPoint1, FPoint2 = Dm.ContinuumRegions(Int, Wave)
print 'continuum', Dm.ContinuumFlux
#Check if emission or absorption line as well as line mixture
Dm.Check_GaussianMixture(Wave)
#This comes from the Measure_LineIntensity method
if Dm.Deblend_Check == True:
Wave_Gaussian, Int_Gaussian = Dm.Measure_LineIntensity(Wave, Int, SaveData=False)
