# -*- coding: utf-8 -*-
"""
Created on Sat Dec 10 16:42:34 2016

@author: grav
"""

try:
   import pyfits
except:
   from astropy.io import fits as pyfits
import numpy as np
import matplotlib.pyplot as plt
import os


#file= "/data/1/2016-12-10/GRAVITY.2016-12-10T03-43-40.fits"
file= "/data/1/2016-12-10/GRAVITY.2016-12-10T03-32-11.fits" # SC MED COMB
#file= "/data/1/2016-12-10/GRAVITY.2016-12-10T05-59-26.fits"
hdulist = pyfits.open(file)
data_sc  = hdulist['IMAGING_DATA_SC'].data

x=95

plt.figure(23)
plt.clf()
plt.plot(data_sc[:,138,x])

#file= "/data/1/2016-12-10/GRAVITY.2016-12-10T03-37-54.fits"

#hdulist = pyfits.open(file)
#data_sc  = hdulist['IMAGING_DATA_SC'].data

plt.figure(13)
plt.clf()
plt.plot(data_sc[:,283,x])

plt.figure(12)
plt.clf()
plt.plot(data_sc[:,177,x])


#file= "/data/1/2016-12-10/reduced_reference/GRAVITY.2016-12-10T03-32-11_debug.fits"
##file= "/data/1/2016-12-09/reduced/GRAVITY.2016-12-09T15-49-03_debug.fits"
#
#hdulist = pyfits.open(file)
#spectrum_SC  = hdulist['SPECTRUM_DATA_SC']
#
#
#plt.figure(12)
#plt.clf()
#plt.plot(spectrum_SC.data["DATA12"][:,x])
#
#plt.figure(13)
#plt.clf()
#plt.plot(spectrum_SC.data["DATA24"][:,x])
#
#plt.figure(14)
#plt.clf()
#plt.plot(spectrum_SC.data["DATA8"][:,x])
#
#plt.figure(23)
#plt.clf()
#plt.plot(spectrum_SC.data["DATA20"][:,x])
#
#plt.figure(24)
#plt.clf()
#plt.plot(spectrum_SC.data["DATA4"][:,x])
#
#plt.figure(34)
#plt.clf()
#plt.plot(spectrum_SC.data["DATA16"][:,x])
#
##file= "/data/1/2016-12-10/reduced/GRAVITY.2016-12-10T03-18-47_preproc.fits"
##hdulist = pyfits.open(file)
#spectrum_FT  = hdulist['SPECTRUM_DATA_FT']
#plt.figure("FT 13")
#plt.clf()
#ft=spectrum_FT.data["DATA16"][:,1]
#plt.plot(ft[::100])
#
##
#plt.figure("FT 13")
#plt.clf()
#fft = np.fft.fft(OPD_FT-OPD_FT.mean())
##
##
##
##OPD_FT  = (hdulist['OPD_FT'].data["OPD"])
##plt.plot(OPD_FT[:])
##
##plt.figure("OPD_FT")
##plt.clf()
##fft = np.fft.fft(OPD_FT-OPD_FT.mean())
##freq = np.fft.fftfreq(OPD_FT.shape[-1], d=1/300.)
##plt.plot(freq, fft)