import sys, os
from astrodata import AstroData
from astrodata import DescriptorUnits as Units
from astrodata.DescriptorUnits import iunits
import GemCalcUtil
## Intergrating the Units module test
ad = AstroData("../../../test_data/gmosspect/N20020810S0117.fits")
print "\n============= CURRENT setup using gemCalcUtil.py in DV =============\n"
centwave = ad.central_wavelength()
print("CASE 1: Get central_wavelength as Meter (DEFAULT):")
print(" centwave = ad.central_wavelength()")
print(" centwave = %s" % centwave.as_str())
print("\nCASE 2: Get central_wavelength as Micrometer:")
print(" um_centwave = ad.central_wavelength(asMicrometers=True)")
umcentwave = ad.central_wavelength(asMicrometers=True)
print(" um_centwave = %s" % umcentwave.as_str())
print("\nCASE 3: Get central_wavelength as Nanometer:")
print(" nm_centwave = ad.central_wavelength(asNanometers=True)")
nmcentwave = ad.central_wavelength(asNanometers=True)
print(" nm_centwave = %s" % nmcentwave.as_str())
print("\nCASE 4: Get central_wavelength as Angstroms:")
print(" ang_centwave = ad.central_wavelength(asAngstroms=True)")
from astrodata import AstroData
gmos = AstroData("../../../test_data/gmosspect/N20020810S0117.fits")
niri = AstroData("../../../test_data/niri/spectra/N20100104S0385.fits")
nifs = AstroData("../../../test_data/recipedata/stackable/set001/N20090902S0079.fits")
print gmos.central_wavelength()
print nifs.central_wavelength()
print gmos.central_wavelength()
v = 1.0/xbin
y = 1<<xbin # all operators are overloaded
#print v,y
data = ad[0].data
#print data
d2 = data/xbin.for_numpy() # see craig for what numpy does if you leave of .for_numpy())
#print d2
print "\n==Old setup results for units==\n"
print "centwave m (default) = 7.5e-07 <type 'float'>"
print "centwave um = 0.75"
print "centwave nm = 750.0:"
print "centwave Angstroms = 7500.0"
ad = AstroData("../../../test_data/gmosspect/N20020810S0117.fits")
centwave = ad.central_wavelength()
print "centwave m (default) = ", centwave, str(type(centwave))
centwave = ad.central_wavelength(asMicrometers=True)
print "centwave um = ", centwave
centwave = ad.central_wavelength(asNanometers=True)
print "centwave nm = ", centwave
centwave = ad.central_wavelength(asAngstroms=True)
print "centwave Angstroms = ", centwave
