from pylab import *
import lamost
files=lamost.get_files("gband.list")
indexfilename="gband.dat"
result="gband.result"
points=lamost.get_points(indexfilename)
lamost.calculate(files, points,area=True,subregion=[4260.,4310.],result=result)
lamost.plotscatter(result=result1,legendy=False)
from pylab import *
import lamost
indexfilename="lithium_B5606304_sp02.dat"
#indexfilename="lithium_set1.dat"
#files=lamost.get_files("lithium.list")
files=lamost.get_files("test_1.list")
result1="lithium1.result"
result2="lithium2.result"
points=lamost.get_points(indexfilename)
#region=[[6680.000, 6700.000 ],[ 6720.000, 6735.000]]
lamost.calculate(files, points, line_strength=3,subregion=[6695.,6720.],result=result1,makemovie=False)
lamost.plotscatter(result=result1,legendy=False)
#lamost.calculate(files, points, line_strength=2,subregion=[6703.,6710.],result=result1,plotting=True)
#lamost.calculate(files, points, line_strength=2,subregion=[6715.,6720.],result=result2,plotting=True)
#lamost.calculate(files, points, wavelength, flux,line_strength=3,result=result)
from pylab import *
import lamost
files=lamost.get_files("calcium.list")
indexfilename="calcium.dat"
result="calcium.result"
points=lamost.get_points(indexfilename)
lamost.calculate(files, points, core_depth=3,result=result)
lamost.plotscatter(result=result,legendy=False)
