Wave_Center, Wave_Center, ind3, ind4 = ExtractSubRegion(Wave, Int, Wave_3, Wave_4)
Wave_Right, Flux_Right, ind5, ind6 = ExtractSubRegion(Wave, Int, Wave_5, Wave_6)
Cont_Error = Determine_Sides_Continuum(Int, Wave, ind1, ind2, ind3, ind4, ind5, ind6)
Continuum_Error = [Cont_Error] * len(Wave)
print
Wave_g = Wave[ind3:ind4]
flux_g = Int[ind3:ind4]
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)
