diff --git a/convert.py b/convert.py
old mode 100644
new mode 100755
diff --git a/convert_gui.py b/convert_gui.py
old mode 100644
new mode 100755
diff --git a/spc/spc.py b/spc/spc.py
index 0baadca..ad185d3 100644
--- a/spc/spc.py
+++ b/spc/spc.py
@@ -4,396 +4,44 @@
 
 author: Rohan Isaac
 """
+# pylint: disable=invalid-name
 
 from __future__ import division, absolute_import, unicode_literals, print_function
 import struct
 import numpy as np
 
-from .sub import subFile, subFileOld
-from .global_fun import read_subheader, flag_bits
+from .sub import subFile, subFileOld, subFileShimadzu
+from .util import read_subheader, flag_bits
 
-
-class File:
-    """
-    Starts loading the data from a .SPC spectral file using data from the
-    header. Stores all the attributes of a spectral file:
-
-    Data
-    ----
-    content: Full raw data
-    sub[i]: sub file object for each subfileFor each subfile
-        sub[i].y: y data for each subfile
-    x: x-data, global, or for the first subheader
-
-    Examples
-    --------
-    >>> import spc
-    >>> ftir_1 = spc.File('/path/to/ftir.spc')
-    """
-
-    # Format strings for various parts of the file
-    # calculate size of strings using `struct.calcsize(string)`
-    head_str = "<cccciddicccci9s9sh32s130s30siicchf48sfifc187s"
-    old_head_str = "<cchfffcchcccc8shh28s130s30s32s"
+class FileFormat:
+    # Format string for various parts of the file
+    # Calculate size of strings using `struct.calcsize(string)`
     logstc_str = "<iiiii44s"
 
-    # byte positon of various parts of the file
-    head_siz = 512
-    old_head_siz = 256
+    # byte position of various parts of the file
     subhead_siz = 32
     log_siz = 64
 
-    subhead1_pos = head_siz + subhead_siz
-
-    # ------------------------------------------------------------------------
-    # CONSTRUCTOR
-    # ------------------------------------------------------------------------
-
-    def __init__(self, filename):
-        # load entire into memory temporarly
-        with open(filename, "rb") as fin:
-            content = fin.read()
-            # print "Read raw data"
-
-        self.length = len(content)
-        # extract first two bytes to determine file type version
-        self.ftflg, self.fversn = struct.unpack('<cc'.encode('utf8'), content[:2])
-        # --------------------------------------------
-        # NEW FORMAT (LSB)
-        # --------------------------------------------
-        if self.fversn == b'\x4b':
-            # format: new LSB 1st
-            # -------------
-            # unpack header
-            # -------------
-            # use little-endian format with standard sizes
-            # use naming scheme in SPC.H header file
-            self.ftflg, \
-                self.fversn, \
-                self.fexper, \
-                self.fexp, \
-                self.fnpts, \
-                self.ffirst, \
-                self.flast, \
-                self.fnsub, \
-                self.fxtype, \
-                self.fytype, \
-                self.fztype, \
-                self.fpost, \
-                self.fdate, \
-                self.fres, \
-                self.fsource, \
-                self.fpeakpt, \
-                self.fspare, \
-                self.fcmnt, \
-                self.fcatxt, \
-                self.flogoff, \
-                self.fmods, \
-                self.fprocs, \
-                self.flevel, \
-                self.fsampin, \
-                self.ffactor, \
-                self.fmethod, \
-                self.fzinc, \
-                self.fwplanes, \
-                self.fwinc, \
-                self.fwtype, \
-                self.freserv \
-                = struct.unpack(self.head_str.encode('utf8'), content[:self.head_siz])
-
-            # Flag bits
-            self.tsprec, \
-                self.tcgram, \
-                self.tmulti, \
-                self.trandm, \
-                self.tordrd, \
-                self.talabs, \
-                self.txyxys, \
-                self.txvals = flag_bits(self.ftflg)[::-1]
-
-            # fix data types if necessary
-            self.fnpts = int(self.fnpts)  # of points should be int
-            self.fexp = ord(self.fexp)
-
-            self.ffirst = float(self.ffirst)
-            self.flast = float(self.flast)
-
-            self.flogoff = int(self.flogoff)  # byte; should be int
-
-            self.fxtype = ord(self.fxtype)
-            self.fytype = ord(self.fytype)
-            self.fztype = ord(self.fztype)
-
-            self.fexper = ord(self.fexper)
-            self.fcmnt = str(self.fcmnt)
-
-            # Convert date time to appropriate format
-            d = self.fdate
-            self.year = d >> 20
-            self.month = (d >> 16) % (2**4)
-            self.day = (d >> 11) % (2**5)
-            self.hour = (d >> 6) % (2**5)
-            self.minute = d % (2**6)
-
-            # null terminated string, replace null characters with spaces
-            # split and join to remove multiple spaces
-            try:
-                self.cmnt = ' '.join((self.fcmnt.replace('\x00', ' ')).split())
-            except:
-                self.cmnt = self.fcmnt
-
-            # figure out type of file
-            if self.fnsub > 1:
-                self.dat_multi = True
-
-            if self.txyxys:
-                # x values are given
-                self.dat_fmt = '-xy'
-            elif self.txvals:
-                # only one subfile, which contains the x data
-                self.dat_fmt = 'x-y'
-            else:
-                # no x values are given, but they can be generated
-                self.dat_fmt = 'gx-y'
-
-            print('{}({})'.format(self.dat_fmt, self.fnsub))
-
-            sub_pos = self.head_siz
-
-            if not self.txyxys:
-                # txyxys don't have global x data
-                if self.txvals:
-                    # if global x data is given
-                    x_dat_pos = self.head_siz
-                    x_dat_end = self.head_siz + (4 * self.fnpts)
-                    self.x = np.array(
-                        [struct.unpack_from(
-                            'f', content[x_dat_pos:x_dat_end], 4 * i)[0]
-                            for i in range(0, self.fnpts)])
-                    sub_pos = x_dat_end
-                else:
-                    # otherwise generate them
-                    self.x = np.linspace(self.ffirst, self.flast, num=self.fnpts)
-
-            # make a list of subfiles
-            self.sub = []
-
-            # if subfile directory is given
-            if self.dat_fmt == '-xy' and self.fnpts > 0:
-                self.directory = True
-                # loop over entries in directory
-                for i in range(0, self.fnsub):
-                    ssfposn, ssfsize, ssftime = struct.unpack(
-                        '<iif'.encode('utf8'), content[self.fnpts + (i * 12):self.fnpts + ((i + 1) * 12)])
-                    # add sufile, load defaults for npts and exp
-                    self.sub.append(subFile(content[ssfposn:ssfposn + ssfsize], 0, 0, True, self.tsprec, self.tmulti))
-
-            else:
-                # don't have directory, for each subfile
-                for i in range(self.fnsub):
-                    # figure out its size
-                    if self.txyxys:
-                        # use points in subfile
-                        subhead_lst = read_subheader(content[sub_pos:(sub_pos + 32)])
-                        pts = subhead_lst[6]
-                        # 4 bytes each for x and y, and 32 for subheader
-                        dat_siz = (8 * pts) + 32
-                    else:
-                        # use global points
-                        pts = self.fnpts
-                        dat_siz = (4 * pts) + 32
-
-                    sub_end = sub_pos + dat_siz
-                    # read into object, add to list
-                    self.sub.append(subFile(content[sub_pos:sub_end],
-                                            self.fnpts, self.fexp, self.txyxys, self.tsprec, self.tmulti))
-                    # update positions
-                    sub_pos = sub_end
-
-            # if log data exists
-            # flog offset to log data offset not zero (bytes)
-            if self.flogoff:
-                log_head_end = self.flogoff + self.log_siz
-                self.logsizd, \
-                    self.logsizm, \
-                    self.logtxto, \
-                    self.logbins, \
-                    self.logdsks, \
-                    self.logspar \
-                    = struct.unpack(self.logstc_str.encode('utf8'),
-                                    content[self.flogoff:log_head_end])
-                log_pos = self.flogoff + self.logtxto
-
-                log_end_pos = log_pos + self.logsizd
-
-                # line endings: get rid of any '\r' and then split on '\n'
-                self.log_content = content[log_pos:log_end_pos].replace(b'\r', b'').split(b'\n')
-
-                # split log data into dictionary based on =
-                self.log_dict = dict()
-                self.log_other = []  # put the rest into a list
-                for x in self.log_content:
-                    if x.find(b'=') >= 0:
-                        # stop it from breaking if there is more than 1 =
-                        key, value = x.split(b'=')[:2]
-                        self.log_dict[key] = value
-                    else:
-                        self.log_other.append(x)
-
-            # spacing between data
-            self.spacing = (self.flast - self.ffirst) / (self.fnpts - 1)
-
-            # call functions
-            self.set_labels()
-            self.set_exp_type()
-
-        # --------------------------------------------
-        # NEW FORMAT (MSB)
-        # --------------------------------------------
-        elif self.fversn == b'\x4c':
-            # new MSB 1st
-            print("New MSB 1st, yet to be implemented")
-            pass  # To be implemented
-
-        # --------------------------------------------
-        # OLD FORMAT
-        # --------------------------------------------
-        elif self.fversn == b'\x4d':
-            # old format
-            # oxtype -> fxtype
-            # oytype -> fytype
-            self.oftflgs, \
-                self.oversn, \
-                self.oexp, \
-                self.onpts, \
-                self.ofirst, \
-                self.olast, \
-                self.fxtype, \
-                self.fytype, \
-                self.oyear, \
-                self.omonth, \
-                self.oday, \
-                self.ohour, \
-                self.ominute, \
-                self.ores, \
-                self.opeakpt, \
-                self.onscans, \
-                self.ospare, \
-                self.ocmnt, \
-                self.ocatxt, \
-                self.osubh1 = struct.unpack(self.old_head_str.encode('utf8'),
-                                            content[:self.old_head_siz])
-
-            # Flag bits (assuming same)
-            self.tsprec, \
-                self.tcgram, \
-                self.tmulti, \
-                self.trandm, \
-                self.tordrd, \
-                self.talabs, \
-                self.txyxys, \
-                self.txvals = flag_bits(self.oftflgs)[::-1]
-
-            # fix data types
-            self.oexp = int(self.oexp)
-            self.onpts = int(self.onpts)  # can't have floating num of pts
-            self.ofirst = float(self.ofirst)
-            self.olast = float(self.olast)
-
-            # Date information
-            # !! to fix !!
-            # Year collected (0=no date/time) - MSB 4 bits are Z type
-
-            # extracted as characters, using ord
-            self.omonth = ord(self.omonth)
-            self.oday = ord(self.oday)
-            self.ohour = ord(self.ohour)
-            self.ominute = ord(self.ominute)
-
-            # number of scans (? subfiles sometimes ?)
-            self.onscans = int(self.onscans)
-
-            # null terminated strings
-            self.ores = self.ores.split(b'\x00')[0]
-            self.ocmnt = self.ocmnt.split(b'\x00')[0]
-
-            # can it have separate x values ?
-            self.x = np.linspace(self.ofirst, self.olast, num=self.onpts)
-
-            # make a list of subfiles
-            self.sub = []
-
-            # already have subheader from main header, retrace steps
-            sub_pos = self.old_head_siz - self.subhead_siz
-
-            # for each subfile
-            # in the old format we don't know how many subfiles to expect,
-            # just looping till we run out
-            i = 0
-            while True:
-                try:
-                    # read in subheader
-                    subhead_lst = read_subheader(content[sub_pos:sub_pos + self.subhead_siz])
-
-                    if subhead_lst[6] > 0:
-                        # default to subfile points, unless it is zero
-                        pts = subhead_lst[6]
-                    else:
-                        pts = self.onpts
-
-                    # figure out size of subheader
-                    dat_siz = (4 * pts)
-                    sub_end = sub_pos + self.subhead_siz + dat_siz
-
-                    # read into object, add to list
-                    # send it pts since we have already figured that out
-                    self.sub.append(subFileOld(
-                        content[sub_pos:sub_end], pts, self.oexp, self.txyxys))
-                    # update next subfile postion, and index
-                    sub_pos = sub_end
-
-                    i += 1
-                except:
-                    # zero indexed, set the total number of subfile
-                    self.fnsub = i + 1
-                    break
-
-            # assuming it can't have separate x values
-            self.dat_fmt = 'gx-y'
-            print('{}({})'.format(self.dat_fmt, self.fnsub))
-
-            self.fxtype = ord(self.fxtype)
-            self.fytype = ord(self.fytype)
-            # need to find from year apparently
-            self.fztype = 0
-            self.set_labels()
-
-        # --------------------------------------------
-        # SHIMADZU
-        # --------------------------------------------
-        elif self.fversn == b'\xcf':
-            print("Highly experimental format, may not work ")
-            raw_data = content[10240:]  # data starts here (maybe every time)
-            # spacing between y and x data is atleast 0 bytes
-            s_32 = chr(int('0', 2)) * 32
-            s_8 = chr(int('0', 2)) * 8  # zero double
-            dat_len = raw_data.find(s_32)
-            for i in range(dat_len, len(raw_data), 8):
-                # find first non zero double
-                if raw_data[i:i + 8] != s_8:
-                    break
-            dat_siz = int(dat_len / 8)
-            self.y = struct.unpack(('<' + dat_siz * 'd').encode('utf8'), raw_data[:dat_len])
-            self.x = struct.unpack(('<' + dat_siz * 'd').encode('utf8'), raw_data[i:i + dat_len])
+    def unpack_flag_bits(self):
+        # Flag bits (assuming same)
+        self.tsprec, \
+            self.tcgram, \
+            self.tmulti, \
+            self.trandm, \
+            self.tordrd, \
+            self.talabs, \
+            self.txyxys, \
+            self.txvals = flag_bits(self.tflgs)[::-1]
 
+        if self.txyxys:
+            # x values are given
+            self.dat_fmt = '-xy'
+        elif self.txvals:
+            # only one subfile, which contains the x data
+            self.dat_fmt = 'x-y'
         else:
-            print("File type %s not supported yet. Please add issue. "
-                  % hex(ord(self.fversn)))
-            self.content = content
-
-    # ------------------------------------------------------------------------
-    # Process other data
-    # ------------------------------------------------------------------------
+            # no x values are given, but they can be generated
+            self.dat_fmt = 'gx-y'
 
     def set_labels(self):
         """
@@ -434,13 +82,13 @@ def set_labels(self):
                      "Millimeters (mm)",
                      "Hours"]
 
-        if self.fxtype < 30:
-            self.xlabel = fxtype_op[self.fxtype]
+        if self.xtype < 30:
+            self.xlabel = fxtype_op[self.xtype]
         else:
             self.xlabel = "Unknown"
 
-        if self.fztype < 30:
-            self.zlabel = fxtype_op[self.fztype]
+        if self.ztype < 30:
+            self.zlabel = fxtype_op[self.ztype]
         else:
             self.zlabel = "Unknown"
 
@@ -481,10 +129,10 @@ def set_labels(self):
                       "Arbitrary or Single Beam with Valley Peaks",
                       "Emission"]
 
-        if self.fytype < 27:
-            self.ylabel = fytype_op[self.fytype]
-        elif self.fytype > 127 and self.fytype < 132:
-            self.ylabel = fytype_op2[self.fytype - 128]
+        if self.ytype < 27:
+            self.ylabel = fytype_op[self.ytype]
+        elif 127 < self.ytype < 132:
+            self.ylabel = fytype_op2[self.ytype - 128]
         else:
             self.ylabel = "Unknown"
 
@@ -495,7 +143,7 @@ def set_labels(self):
         # split it based on 00 string
         # format x, y, z
         if self.talabs:
-            ll = self.fcatxt.split(b'\x00')
+            ll = self.catxt.split(b'\x00')
             if len(ll) > 2:
                 # make sure there are enough items to extract from
                 xl, yl, zl = ll[:3]
@@ -526,7 +174,393 @@ def set_exp_type(self):
                      "Atomic Spectrum",
                      "Chromatography Diode Array Spectra"]
 
-        self.exp_type = fexper_op[self.fexper]
+        self.exp_type = fexper_op[self.exper]
+
+
+class OldFormat(FileFormat):
+    # Format string for the header
+    # Calculate size of strings using `struct.calcsize(string)`
+    head_str = "<cchfffcchcccc8shh28s130s30s32s"
+
+    # byte position
+    head_siz = 256
+
+    def __init__(self, content):
+        self.unpack_header(content)
+        self.unpack_flag_bits()
+
+        # can it have separate x values ?
+        self.x = np.linspace(self.first, self.last, num=self.npts)
+
+        # make a list of subfiles
+        self.sub = []
+
+        # already have subheader from main header, retrace steps
+        sub_pos = self.head_siz - self.subhead_siz
+
+        # for each subfile
+        # in the old format we don't know how many subfiles to expect,
+        # just looping till we run out
+        i = 0
+        while True:
+            try:
+                # read in subheader
+                subhead_lst = read_subheader(content[sub_pos:sub_pos + self.subhead_siz])
+
+                if subhead_lst[6] > 0:
+                    # default to subfile points, unless it is zero
+                    pts = subhead_lst[6]
+                else:
+                    pts = self.npts
+
+                # figure out size of subheader
+                dat_siz = (4 * pts)
+                sub_end = sub_pos + self.subhead_siz + dat_siz
+
+                # read into object, add to list
+                # send it pts since we have already figured that out
+                self.sub.append(subFileOld(
+                    content[sub_pos:sub_end], pts, self.exp, self.txyxys))
+                # update next subfile postion, and index
+                sub_pos = sub_end
+
+                i += 1
+            except:
+                # zero indexed, set the total number of subfile
+                self.nsub = i + 1
+                break
+
+        print('{}({})'.format(self.dat_fmt, self.nsub))
+
+        self.set_labels()
+
+    def unpack_header(self, content):
+        self.tflgs, \
+            self.versn, \
+            self.exp, \
+            self.npts, \
+            self.first, \
+            self.last, \
+            self.xtype, \
+            self.ytype, \
+            self.year, \
+            self.month, \
+            self.day, \
+            self.hour, \
+            self.minute, \
+            self.res, \
+            self.peakpt, \
+            self.nscans, \
+            self.spare, \
+            self.cmnt, \
+            self.catxt, \
+            self.subh1 = struct.unpack(self.head_str.encode('utf8'),
+                                        content[:self.head_siz])
+
+        # fix data types
+        self.exp = int(self.exp)
+        self.npts = int(self.npts)  # can't have floating num of pts
+        self.first = float(self.first)
+        self.last = float(self.last)
+
+        # Date information
+        # !! to fix !!
+        # Year collected (0=no date/time) - MSB 4 bits are Z type
+
+        # extracted as characters, using ord
+        self.month = ord(self.month)
+        self.day = ord(self.day)
+        self.hour = ord(self.hour)
+        self.minute = ord(self.minute)
+
+        # number of scans (? subfiles sometimes ?)
+        self.nscans = int(self.nscans)
+
+        # null terminated strings
+        self.res = self.res.split(b'\x00')[0]
+        self.cmnt = self.cmnt.split(b'\x00')[0]
+
+        self.xtype = ord(self.xtype)
+        self.ytype = ord(self.ytype)
+        # need to find from year apparently
+        self.ztype = 0
+
+class NewFormat(FileFormat):
+    # Format string for the header
+    # Calculate size of strings using `struct.calcsize(string)`
+    head_str = "<cccciddicccci9s9sh32s130s30siicchf48sfifc187s"
+
+    # byte position
+    head_siz = 512
+
+    subhead1_pos = head_siz + FileFormat.subhead_siz
+
+
+class NewFormatLSB(NewFormat):
+
+    def __init__(self, content):
+        # use naming scheme in SPC.H header file
+        self.unpack_header(content)
+        self.unpack_flag_bits()
+
+        # TODO use __str__ instead?
+        print('{}({})'.format(self.dat_fmt, self.nsub))
+
+        sub_pos = self.head_siz
+
+        if not self.txyxys:
+            # txyxys don't have global x data
+            if self.txvals:
+                # if global x data is given
+                x_dat_pos = self.head_siz
+                x_dat_end = self.head_siz + (4 * self.npts)
+                self.x = np.array(
+                    [struct.unpack_from(
+                        'f', content[x_dat_pos:x_dat_end], 4 * i)[0]
+                        for i in range(0, self.npts)])
+                sub_pos = x_dat_end
+            else:
+                # otherwise generate them
+                self.x = np.linspace(self.first, self.last, num=self.npts)
+
+        # make a list of subfiles
+        self.sub = []
+
+        # if subfile directory is given
+        if self.dat_fmt == '-xy' and self.npts > 0:
+            self.directory = True
+            # loop over entries in directory
+            for i in range(0, self.nsub):
+                ssfposn, ssfsize, ssftime = struct.unpack(
+                    '<iif'.encode('utf8'), content[self.npts + (i * 12):self.npts + ((i + 1) * 12)])
+                # add sufile, load defaults for npts and exp
+                self.sub.append(subFile(content[ssfposn:ssfposn + ssfsize], 0, 0, True, self.tsprec, self.tmulti))
+
+        else:
+            # don't have directory, for each subfile
+            for i in range(self.nsub):
+                # figure out its size
+                if self.txyxys:
+                    # use points in subfile
+                    subhead_lst = read_subheader(content[sub_pos:(sub_pos + 32)])
+                    pts = subhead_lst[6]
+                    # 4 bytes each for x and y, and 32 for subheader
+                    dat_siz = (8 * pts) + 32
+                else:
+                    # use global points
+                    pts = self.npts
+                    dat_siz = (4 * pts) + 32
+
+                sub_end = sub_pos + dat_siz
+                # read into object, add to list
+                self.sub.append(subFile(content[sub_pos:sub_end],
+                                        self.npts, self.exp, self.txyxys, self.tsprec, self.tmulti))
+                # update positions
+                sub_pos = sub_end
+
+        # if log data exists
+        # flog offset to log data offset not zero (bytes)
+        if self.logoff:
+            log_head_end = self.logoff + self.log_siz
+            self.logsizd, \
+                self.logsizm, \
+                self.logtxto, \
+                self.logbins, \
+                self.logdsks, \
+                self.logspar \
+                = struct.unpack(FileFormat.logstc_str.encode('utf8'),
+                                content[self.logoff:log_head_end])
+            log_pos = self.logoff + self.logtxto
+
+            log_end_pos = log_pos + self.logsizd
+
+            # line endings: get rid of any '\r' and then split on '\n'
+            self.log_content = content[log_pos:log_end_pos].replace(b'\r', b'').split(b'\n')
+
+            # split log data into dictionary based on =
+            self.log_dict = dict()
+            self.log_other = []  # put the rest into a list
+            for x in self.log_content:
+                if x.find(b'=') >= 0:
+                    # stop it from breaking if there is more than 1 =
+                    key, value = x.split(b'=')[:2]
+                    self.log_dict[key] = value
+                else:
+                    self.log_other.append(x)
+
+        # call functions
+        self.set_labels()
+        self.set_exp_type()
+
+    def unpack_header(self, content):
+        # use little-endian format with standard sizes
+        # use naming scheme in SPC.H header file
+        self.tflgs, \
+            self.versn, \
+            self.exper, \
+            self.exp, \
+            self.npts, \
+            self.first, \
+            self.last, \
+            self.nsub, \
+            self.xtype, \
+            self.ytype, \
+            self.ztype, \
+            self.post, \
+            self.date, \
+            self.res, \
+            self.source, \
+            self.peakpt, \
+            self.spare, \
+            self.cmnt, \
+            self.catxt, \
+            self.logoff, \
+            self.mods, \
+            self.procs, \
+            self.level, \
+            self.sampin, \
+            self.factor, \
+            self.method, \
+            self.zinc, \
+            self.wplanes, \
+            self.winc, \
+            self.wtype, \
+            self.reserv \
+            = struct.unpack(self.head_str.encode('utf8'), content[:self.head_siz])
+
+        # fix data types if necessary
+        self.npts = int(self.npts)
+        self.exp = ord(self.exp)
+
+        self.first = float(self.first)
+        self.last = float(self.last)
+
+        # spacing between data
+        self.spacing = (self.last - self.first) / (self.npts - 1)
+
+        self.logoff = int(self.logoff)  # byte; should be int
+
+        self.xtype = ord(self.xtype)
+        self.ytype = ord(self.ytype)
+        self.ztype = ord(self.ztype)
+
+        self.exper = ord(self.exper)
+        self.cmnt = str(self.cmnt)
+
+        # Convert date time to appropriate format
+        d = self.date
+        self.year = d >> 20
+        self.month = (d >> 16) % (2**4)
+        self.day = (d >> 11) % (2**5)
+        self.hour = (d >> 6) % (2**5)
+        self.minute = d % (2**6)
+
+        # null terminated string, replace null characters with spaces
+        # split and join to remove multiple spaces
+        try:
+            self.cmnt = ' '.join((self.cmnt.replace('\x00', ' ')).split())
+        except:
+            pass
+
+        # figure out type of file
+        self.dat_multi = self.nsub > 1
+
+class NewFormatMSB(NewFormat):
+    def __init__(self, content):
+        print("New MSB 1st, yet to be implemented")
+        pass  # To be implemented
+
+class ShimadzuFormat(FileFormat):
+    def __init__(self, content):
+        print("Highly experimental format, may not work ")
+        raw_data = content[10240:]  # data starts here (maybe every time)
+        # spacing between y and x data is atleast 0 bytes
+        s_32 = chr(int('0', 2)) * 32
+        s_8 = chr(int('0', 2)) * 8  # zero double
+        dat_len = raw_data.find(s_32)
+        for i in range(dat_len, len(raw_data), 8):
+            # find first non zero double
+            if raw_data[i:i + 8] != s_8:
+                break
+        dat_siz = int(dat_len / 8)
+
+        self.dat_fmt = 'x-y'
+        self.nsub = 1
+
+        self.y = struct.unpack(('<' + dat_siz * 'd').encode('utf8'), raw_data[:dat_len])
+        self.x = struct.unpack(('<' + dat_siz * 'd').encode('utf8'), raw_data[i:i + dat_len])
+
+        self.ylabel = ''
+        self.xlabel = ''
+
+        # creating a `sub` member to maintain consistency
+        self.sub = []
+        self.sub.append(subFileShimadzu(self.y))
+
+class File:
+    """
+    Starts loading the data from a .SPC spectral file using data from the
+    header. Stores all the attributes of a spectral file:
+
+    Data
+    ----
+    content: Full raw data
+    sub[i]: sub file object for each subfileFor each subfile
+        sub[i].y: y data for each subfile
+    x: x-data, global, or for the first subheader
+
+    Examples
+    --------
+    >>> import spc
+    >>> ftir_1 = spc.File('/path/to/ftir.spc')
+    """
+
+    # ------------------------------------------------------------------------
+    # CONSTRUCTOR
+    # ------------------------------------------------------------------------
+
+    def __init__(self, filename):
+        # load entire into memory temporarly
+        with open(filename, "rb") as fin:
+            content = fin.read()
+            # print "Read raw data"
+
+        self.length = len(content)
+        # extract first two bytes to determine file type version
+        self.tflg, self.versn = struct.unpack('<cc'.encode('utf8'), content[:2])
+
+        # --------------------------------------------
+        # NEW FORMAT (LSB)
+        # --------------------------------------------
+        if self.versn == b'\x4b':
+            # format: new LSB 1st
+            self.format = NewFormatLSB(content)
+
+        # --------------------------------------------
+        # NEW FORMAT (MSB)
+        # --------------------------------------------
+        elif self.versn == b'\x4c':
+            # new MSB 1st
+            self.format = NewFormatMSB(content)
+
+        # --------------------------------------------
+        # OLD FORMAT
+        # --------------------------------------------
+        elif self.versn == b'\x4d':
+            # old format
+            self.format = OldFormat(content)
+
+        # --------------------------------------------
+        # SHIMADZU
+        # --------------------------------------------
+        elif self.versn == b'\xcf':
+            self.format = ShimadzuFormat(content)
+            self.content = content
+
+        else:
+            print("File type %s not supported yet. Please add issue. "
+                  % hex(ord(self.versn)))
+            self.content = content
 
     # ------------------------------------------------------------------------
     # output
@@ -549,27 +583,27 @@ def data_txt(self, delimiter='\t', newline='\n'):
         """
 
         dat = ''
-        if self.fnsub == 1:
-            if self.dat_fmt.endswith('-xy'):
-                x = self.sub[0].x
+        if self.format.nsub == 1:
+            if self.format.dat_fmt.endswith('-xy'):
+                x = self.format.sub[0].x
             else:
-                x = self.x
-            y = self.sub[0].y
+                x = self.format.x
+            y = self.format.sub[0].y
 
             for x1, y1 in zip(x, y):
                 dat += '{}{}{}{}'.format(x1, delimiter, y1, newline)
         else:
-            if not self.dat_fmt.endswith('-xy'):
+            if not self.format.dat_fmt.endswith('-xy'):
                 # does not have separate x data
-                for i in range(len(self.x)):
-                    dat += '{}'.format(self.x[i])
-                    for s in self.sub:
+                for i in range(len(self.format.x)):
+                    dat += '{}'.format(self.format.x[i])
+                    for s in self.format.sub:
                         dat += '{}{}'.format(delimiter, s.y[i])
                     dat += newline
             else:
                 # txyxy format, return one long xy file with subfiles
                 # separated by blank lines
-                for i in self.sub:
+                for i in self.format.sub:
                     for x1, y1 in zip(i.x, i.y):
                         dat += '{}{}{}{}'.format(x1, delimiter, y1, newline)
                     dat += newline
@@ -577,27 +611,27 @@ def data_txt(self, delimiter='\t', newline='\n'):
 
     # Writes out data to a stream (significantly faster than appending to a string)
     def stream_data_txt(self, stream, delimiter='\t', newline='\n'):
-        if self.fnsub == 1:
-            if self.dat_fmt.endswith('-xy'):
-                x = self.sub[0].x
+        if self.format.nsub == 1:
+            if self.format.dat_fmt.endswith('-xy'):
+                x = self.format.sub[0].x
             else:
-                x = self.x
-            y = self.sub[0].y
+                x = self.format.x
+            y = self.format.sub[0].y
 
             for x1, y1 in zip(x, y):
                 stream.write('{}{}{}{}'.format(x1, delimiter, y1, newline))
         else:
-            if not self.dat_fmt.endswith('-xy'):
+            if not self.format.dat_fmt.endswith('-xy'):
                 # does not have separate x data
-                for i in range(len(self.x)):
-                    stream.write('{}'.format(self.x[i]))
-                    for s in self.sub:
+                for i in range(len(self.format.x)):
+                    stream.write('{}'.format(self.format.x[i]))
+                    for s in self.format.sub:
                         stream.write('{}{}'.format(delimiter, s.y[i]))
                     stream.write(newline)
             else:
                 # txyxy format, return one long xy file with subfiles
                 # separated by blank lines
-                for i in self.sub:
+                for i in self.format.sub:
                     for x1, y1 in zip(i.x, i.y):
                         stream.write('{}{}{}{}'.format(x1, delimiter, y1, newline))
                     stream.write(newline)
@@ -624,11 +658,11 @@ def write_file(self, path, delimiter='\t', newline='\n'):
 
     def print_metadata(self):
         """ Print out select metadata"""
-        print("Scan: ", self.log_dict['Comment'], "\n",
-              float(self.log_dict['Start']), "to ",
-              float(self.log_dict['End']), "; ",
-              float(self.log_dict['Increment']), "cm-1;",
-              float(self.log_dict['Integration Time']), "s integration time")
+        print("Scan: ", self.format.log_dict['Comment'], "\n",
+              float(self.format.log_dict['Start']), "to ",
+              float(self.format.log_dict['End']), "; ",
+              float(self.format.log_dict['Increment']), "cm-1;",
+              float(self.format.log_dict['Integration Time']), "s integration time")
 
     def plot(self):
         """ Plots data, and use column headers, returns figure object plotted
@@ -641,15 +675,15 @@ def plot(self):
 
         """
         import matplotlib.pyplot as plt
-        if self.dat_fmt.endswith('-xy'):
-            for s in self.sub:
+        if self.format.dat_fmt.endswith('-xy'):
+            for s in self.format.sub:
                 plt.plot(s.x, s.y)
         else:
-            x = self.x
-            for s in self.sub:
+            x = self.format.x
+            for s in self.format.sub:
                 plt.plot(x, s.y)
-        plt.xlabel(self.xlabel)
-        plt.ylabel(self.ylabel)
+        plt.xlabel(self.format.xlabel)
+        plt.ylabel(self.format.ylabel)
         return plt.gcf()
 
     def debug_info(self):
@@ -664,49 +698,49 @@ def debug_info(self):
         """
         print("\nDEBUG INFO\nFlags:\n")
         # Flag bits
-        if self.tsprec:
+        if self.format.tsprec:
             print("16-bit y data")
-        if self.tcgram:
+        if self.format.tcgram:
             print("enable fexper")
-        if self.tmulti:
+        if self.format.tmulti:
             print("multiple traces")
-        if self.trandm:
+        if self.format.trandm:
             print("arb time (z) values")
-        if self.tordrd:
+        if self.format.tordrd:
             print("ordered but uneven subtimes")
-        if self.talabs:
+        if self.format.talabs:
             print("use fcatxt axis not fxtype")
-        if self.txyxys:
+        if self.format.txyxys:
             print("each subfile has own x's")
-        if self.txvals:
+        if self.format.txvals:
             print("floating x-value array preceeds y's")
 
         print('----\n')
         # spc format version
-        if self.fversn == chr(0x4b):
-            self.pr_versn = "new LSB 1st"
-        elif self.fversn == chr(0x4c):
-            self.pr_versn = "new MSB 1st"
-        elif self.fversn == chr(0x4d):
-            self.pr_versn = "old format"
+        if self.format.versn == chr(0x4b):
+            self.format.pr_versn = "new LSB 1st"
+        elif self.format.versn == chr(0x4c):
+            self.format.pr_versn = "new MSB 1st"
+        elif self.format.versn == chr(0x4d):
+            self.format.pr_versn = "old format"
         else:
-            self.pr_versn = "unknown version"
+            self.format.pr_versn = "unknown version"
 
-        print("Version:", self.pr_versn)
+        print("Version:", self.format.pr_versn)
 
         # subfiles
-        if self.fnsub == 1:
+        if self.format.nsub == 1:
             print("Single file only")
         else:
-            print("Multiple subfiles:", self.fnsub)
+            print("Multiple subfiles:", self.format.nsub)
 
         # multiple y values
-        if self.tmulti:
+        if self.format.tmulti:
             print("Multiple y-values")
         else:
             print("Single set of y-values")
 
-        # print "There are ", self.fnpts, \
-        #    " points between ", self.ffirst, \
-        #    " and ", self.flast, \
+        # print "There are ", self.npts, \
+        #    " points between ", self.first, \
+        #    " and ", self.last, \
         #    " in steps of ", self.pr_spacing
diff --git a/spc/sub.py b/spc/sub.py
index a20b38a..888b69e 100644
--- a/spc/sub.py
+++ b/spc/sub.py
@@ -9,7 +9,7 @@
 import struct
 import numpy as np
 
-from .global_fun import read_subheader
+from .util import read_subheader
 
 
 class subFile:
@@ -26,6 +26,8 @@ class subFile:
     """
 
     def __init__(self, data, fnpts, fexp, txyxy, tsprec, tmulti):
+        # fixed header size, 32 bytes
+        y_dat_pos = 32
 
         # extract subheader info
         self.subflgs, \
@@ -38,16 +40,10 @@ def __init__(self, data, fnpts, fexp, txyxy, tsprec, tmulti):
             self.subscan, \
             self.subwlevel, \
             self.subresv \
-            = read_subheader(data[:32])
-
-        # header is 32 bytes
-        y_dat_pos = 32
+            = read_subheader(data[:y_dat_pos])
 
-        if txyxy:
-            # only reason to use subnpts if x data is here
-            pts = self.subnpts
-        else:
-            pts = fnpts
+        # only reason to use subnpts if x data is here
+        pts = self.subnpts if txyxy else fnpts
 
         # Choosing exponent
         # -----------------
@@ -118,7 +114,7 @@ class subFileOld:
     """
 
     def __init__(self, data, pts, fexp, txyxy):
-        # fixed header size
+        # fixed header size, 32 bytes
         y_dat_pos = 32
 
         # extract subheader info
@@ -192,3 +188,10 @@ def __init__(self, data, pts, fexp, txyxy):
         # if 1 subfile changed
         # if 8 if peak table should not be used
         # if 128 if subfile modified by arithmetic
+
+class subFileShimadzu:
+    """A "subFile" class for Shimadzu files
+    This is a dummy class, currently used only to maintain consistency
+    """
+    def __init__(self, y):
+        self.y = y
diff --git a/spc/global_fun.py b/spc/util.py
similarity index 100%
rename from spc/global_fun.py
rename to spc/util.py
index dbd21da..f3f72d9 100644
--- a/spc/global_fun.py
+++ b/spc/util.py
@@ -3,14 +3,15 @@
 
 author: Rohan Isaac
 """
+
 from __future__ import division, absolute_import, unicode_literals, print_function
 import struct
 
+
 # ------------------------------------------------------------------------
 # Process subfile data
 # ------------------------------------------------------------------------
 
-
 def read_subheader(subheader):
     """
     Return the subheader as a list
@@ -44,14 +45,13 @@ def read_subheader(subheader):
 
     return item_cpy
 
+
 # ------------------------------------------------------------------------
 # Decode a character to boolean array
 # ------------------------------------------------------------------------
 
-
 def flag_bits(n):
     """Return the bits of a byte as a boolean array:
-
     n (charater):
         8-bit character passed
 
diff --git a/test_data/Renishaw Ramascope Raman.spc b/test_data/Renishaw Ramascope Raman.spc
new file mode 100644
index 0000000..b5db07c
Binary files /dev/null and b/test_data/Renishaw Ramascope Raman.spc differ
diff --git a/test_data/Shimadzu UV-Vis.spc b/test_data/Shimadzu UV-Vis.spc
new file mode 100644
index 0000000..facc7cd
Binary files /dev/null and b/test_data/Shimadzu UV-Vis.spc differ
diff --git a/test_data/txt2/Renishaw Ramascope Raman.spc.txt b/test_data/txt2/Renishaw Ramascope Raman.spc.txt
new file mode 100644
index 0000000..a50ba86
--- /dev/null
+++ b/test_data/txt2/Renishaw Ramascope Raman.spc.txt	
@@ -0,0 +1,43 @@
+368.1315625	1739.97631836
+369.204176897	1670.51464844
+370.276791295	1772.1072998
+371.349405692	1729.35314941
+372.422020089	1637.3717041
+373.494634487	1650.23095703
+374.567248884	1707.82946777
+375.639863281	1652.10180664
+376.712477679	1703.54992676
+377.785092076	1768.54199219
+378.857706473	1679.75842285
+379.930320871	1747.78295898
+381.002935268	1659.63220215
+382.075549665	1673.01416016
+383.148164062	1700.8236084
+384.22077846	1699.73352051
+385.293392857	1661.68554688
+386.366007254	1687.59521484
+387.438621652	1603.41943359
+388.511236049	1688.19128418
+389.583850446	1628.76904297
+390.656464844	1626.72216797
+391.729079241	1617.99267578
+392.801693638	1552.7010498
+393.874308036	1615.81506348
+394.946922433	1611.68823242
+396.01953683	1625.23168945
+397.092151228	1644.00622559
+398.164765625	1622.30163574
+399.237380022	1646.70056152
+400.30999442	1661.34460449
+401.382608817	1725.60058594
+402.455223214	1588.95654297
+403.527837612	1690.34729004
+404.600452009	1583.52783203
+405.673066406	1707.29284668
+406.745680804	1627.45410156
+407.818295201	1573.32958984
+408.890909598	1660.47509766
+409.963523996	1703.47851562
+411.036138393	1625.2565918
+412.10875279	1587.6986084
+413.181367187	1587.55078125
diff --git a/test_data/txt2/Shimadzu UV-Vis.spc.txt b/test_data/txt2/Shimadzu UV-Vis.spc.txt
new file mode 100644
index 0000000..13a8ab9
--- /dev/null
+++ b/test_data/txt2/Shimadzu UV-Vis.spc.txt	
@@ -0,0 +1,1001 @@
+800.0	-0.0126953125
+799.5	-0.0126342773438
+799.0	-0.0126037597656
+798.5	-0.0124969482422
+798.0	-0.0123596191406
+797.5	-0.0124359130859
+797.0	-0.0125885009766
+796.5	-0.0127563476562
+796.0	-0.0127410888672
+795.5	-0.0126953125
+795.0	-0.0128021240234
+794.5	-0.0129089355469
+794.0	-0.0129547119141
+793.5	-0.0130310058594
+793.0	-0.0132141113281
+792.5	-0.0131378173828
+792.0	-0.0132751464844
+791.5	-0.0135803222656
+791.0	-0.0135955810547
+790.5	-0.0134887695312
+790.0	-0.0133209228516
+789.5	-0.0132751464844
+789.0	-0.0132751464844
+788.5	-0.0133209228516
+788.0	-0.0135192871094
+787.5	-0.0136260986328
+787.0	-0.0136108398438
+786.5	-0.0134887695312
+786.0	-0.013427734375
+785.5	-0.0133972167969
+785.0	-0.0132751464844
+784.5	-0.0134887695312
+784.0	-0.0135192871094
+783.5	-0.0133056640625
+783.0	-0.0133666992188
+782.5	-0.0135040283203
+782.0	-0.0135955810547
+781.5	-0.0136260986328
+781.0	-0.0134735107422
+780.5	-0.0135650634766
+780.0	-0.0136260986328
+779.5	-0.0132293701172
+779.0	-0.01318359375
+778.5	-0.0135955810547
+778.0	-0.0137023925781
+777.5	-0.0136108398438
+777.0	-0.0138244628906
+776.5	-0.0140838623047
+776.0	-0.0142211914062
+775.5	-0.0143585205078
+775.0	-0.014404296875
+774.5	-0.0142822265625
+774.0	-0.0142364501953
+773.5	-0.0141906738281
+773.0	-0.0142517089844
+772.5	-0.0142517089844
+772.0	-0.0141754150391
+771.5	-0.0144805908203
+771.0	-0.0146484375
+770.5	-0.0146026611328
+770.0	-0.0147094726562
+769.5	-0.0147399902344
+769.0	-0.0146942138672
+768.5	-0.0145721435547
+768.0	-0.0145263671875
+767.5	-0.0146789550781
+767.0	-0.0146942138672
+766.5	-0.0146789550781
+766.0	-0.0146331787109
+765.5	-0.0145568847656
+765.0	-0.0147399902344
+764.5	-0.014892578125
+764.0	-0.0150146484375
+763.5	-0.0152893066406
+763.0	-0.0153045654297
+762.5	-0.0151214599609
+762.0	-0.0149993896484
+761.5	-0.0150756835938
+761.0	-0.0152435302734
+760.5	-0.0152130126953
+760.0	-0.0152130126953
+759.5	-0.0152587890625
+759.0	-0.0153350830078
+758.5	-0.015380859375
+758.0	-0.0153350830078
+757.5	-0.015380859375
+757.0	-0.0153045654297
+756.5	-0.0153350830078
+756.0	-0.0154113769531
+755.5	-0.0154876708984
+755.0	-0.0155792236328
+754.5	-0.0154724121094
+754.0	-0.0155944824219
+753.5	-0.0156707763672
+753.0	-0.0155639648438
+752.5	-0.0156860351562
+752.0	-0.0158538818359
+751.5	-0.0157470703125
+751.0	-0.015625
+750.5	-0.0157012939453
+750.0	-0.0158996582031
+749.5	-0.0158996582031
+749.0	-0.0157775878906
+748.5	-0.0157012939453
+748.0	-0.0157318115234
+747.5	-0.0157165527344
+747.0	-0.0156402587891
+746.5	-0.0159301757812
+746.0	-0.0160980224609
+745.5	-0.0160064697266
+745.0	-0.0160827636719
+744.5	-0.0162200927734
+744.0	-0.0162811279297
+743.5	-0.0161590576172
+743.0	-0.0162200927734
+742.5	-0.0165100097656
+742.0	-0.0168304443359
+741.5	-0.0167083740234
+741.0	-0.0164642333984
+740.5	-0.0165405273438
+740.0	-0.0167388916016
+739.5	-0.0167999267578
+739.0	-0.0168609619141
+738.5	-0.01708984375
+738.0	-0.0170593261719
+737.5	-0.0169525146484
+737.0	-0.0169219970703
+736.5	-0.0168609619141
+736.0	-0.0168151855469
+735.5	-0.0169372558594
+735.0	-0.0169525146484
+734.5	-0.0170593261719
+734.0	-0.0172271728516
+733.5	-0.0173034667969
+733.0	-0.0173950195312
+732.5	-0.0175476074219
+732.0	-0.0177764892578
+731.5	-0.0179443359375
+731.0	-0.0180816650391
+730.5	-0.0181427001953
+730.0	-0.0180206298828
+729.5	-0.0179443359375
+729.0	-0.0179443359375
+728.5	-0.0180053710938
+728.0	-0.0181121826172
+727.5	-0.0182495117188
+727.0	-0.0185852050781
+726.5	-0.0186309814453
+726.0	-0.0182800292969
+725.5	-0.0182800292969
+725.0	-0.0185241699219
+724.5	-0.0184478759766
+724.0	-0.0184326171875
+723.5	-0.0186462402344
+723.0	-0.0188446044922
+722.5	-0.0188446044922
+722.0	-0.0188903808594
+721.5	-0.0190124511719
+721.0	-0.0188446044922
+720.5	-0.0189361572266
+720.0	-0.0190582275391
+719.5	-0.0191345214844
+719.0	-0.0193023681641
+718.5	-0.0192718505859
+718.0	-0.0194854736328
+717.5	-0.0198059082031
+717.0	-0.019775390625
+716.5	-0.0198211669922
+716.0	-0.02001953125
+715.5	-0.0201416015625
+715.0	-0.0202789306641
+714.5	-0.0203552246094
+714.0	-0.0203857421875
+713.5	-0.0204772949219
+713.0	-0.0205078125
+712.5	-0.0202941894531
+712.0	-0.020263671875
+711.5	-0.0204620361328
+711.0	-0.0206451416016
+710.5	-0.0208435058594
+710.0	-0.0209655761719
+709.5	-0.0210723876953
+709.0	-0.0212249755859
+708.5	-0.021240234375
+708.0	-0.0212554931641
+707.5	-0.0215301513672
+707.0	-0.0217895507812
+706.5	-0.0216979980469
+706.0	-0.0217132568359
+705.5	-0.0220184326172
+705.0	-0.0222015380859
+704.5	-0.0223388671875
+704.0	-0.0224304199219
+703.5	-0.0225677490234
+703.0	-0.0226593017578
+702.5	-0.0226593017578
+702.0	-0.0227203369141
+701.5	-0.0228576660156
+701.0	-0.0231323242188
+700.5	-0.0233764648438
+700.0	-0.0235900878906
+699.5	-0.0238952636719
+699.0	-0.0241394042969
+698.5	-0.0242462158203
+698.0	-0.0243530273438
+697.5	-0.0244445800781
+697.0	-0.0245971679688
+696.5	-0.0248413085938
+696.0	-0.0249938964844
+695.5	-0.0250396728516
+695.0	-0.0252838134766
+694.5	-0.0256958007812
+694.0	-0.0260772705078
+693.5	-0.0263824462891
+693.0	-0.0265045166016
+692.5	-0.0264129638672
+692.0	-0.0265197753906
+691.5	-0.0267181396484
+691.0	-0.0268402099609
+690.5	-0.0271759033203
+690.0	-0.0274200439453
+689.5	-0.0276184082031
+689.0	-0.0278015136719
+688.5	-0.0281066894531
+688.0	-0.0285186767578
+687.5	-0.0289916992188
+687.0	-0.0293884277344
+686.5	-0.0294342041016
+686.0	-0.0294189453125
+685.5	-0.0294952392578
+685.0	-0.0297393798828
+684.5	-0.0300903320312
+684.0	-0.0304565429688
+683.5	-0.0308227539062
+683.0	-0.0314483642578
+682.5	-0.0321197509766
+682.0	-0.0326080322266
+681.5	-0.0327606201172
+681.0	-0.0327453613281
+680.5	-0.0330810546875
+680.0	-0.0337371826172
+679.5	-0.0344390869141
+679.0	-0.0349426269531
+678.5	-0.0350494384766
+678.0	-0.0349884033203
+677.5	-0.0353393554688
+677.0	-0.0360260009766
+676.5	-0.0368347167969
+676.0	-0.0373229980469
+675.5	-0.0375061035156
+675.0	-0.0375518798828
+674.5	-0.0378265380859
+674.0	-0.0382537841797
+673.5	-0.0385437011719
+673.0	-0.0388336181641
+672.5	-0.0390472412109
+672.0	-0.0391693115234
+671.5	-0.0392913818359
+671.0	-0.0394439697266
+670.5	-0.0395965576172
+670.0	-0.0398254394531
+669.5	-0.0401000976562
+669.0	-0.0403442382812
+668.5	-0.0406341552734
+668.0	-0.040771484375
+667.5	-0.0404205322266
+667.0	-0.0402069091797
+666.5	-0.0404205322266
+666.0	-0.0404510498047
+665.5	-0.0404205322266
+665.0	-0.0406799316406
+664.5	-0.0408782958984
+664.0	-0.040771484375
+663.5	-0.0405426025391
+663.0	-0.0403900146484
+662.5	-0.0403442382812
+662.0	-0.0406036376953
+661.5	-0.0407257080078
+661.0	-0.040771484375
+660.5	-0.0408935546875
+660.0	-0.0409240722656
+659.5	-0.0409545898438
+659.0	-0.0409240722656
+658.5	-0.0407409667969
+658.0	-0.0406341552734
+657.5	-0.0405883789062
+657.0	-0.040283203125
+656.5	-0.0399780273438
+656.0	-0.0399932861328
+655.5	-0.0400848388672
+655.0	-0.0401153564453
+654.5	-0.0403442382812
+654.0	-0.0406036376953
+653.5	-0.0408935546875
+653.0	-0.0411071777344
+652.5	-0.0412292480469
+652.0	-0.0412292480469
+651.5	-0.0409851074219
+651.0	-0.0409240722656
+650.5	-0.040771484375
+650.0	-0.0407104492188
+649.5	-0.0410003662109
+649.0	-0.0412750244141
+648.5	-0.0413513183594
+648.0	-0.0413665771484
+647.5	-0.0415954589844
+647.0	-0.041748046875
+646.5	-0.0416870117188
+646.0	-0.0417022705078
+645.5	-0.0418701171875
+645.0	-0.0419921875
+644.5	-0.0419464111328
+644.0	-0.0419616699219
+643.5	-0.0421600341797
+643.0	-0.0423278808594
+642.5	-0.0424194335938
+642.0	-0.0425262451172
+641.5	-0.0426483154297
+641.0	-0.0427703857422
+640.5	-0.0428619384766
+640.0	-0.04296875
+639.5	-0.0430450439453
+639.0	-0.0430755615234
+638.5	-0.0431976318359
+638.0	-0.0434417724609
+637.5	-0.0435485839844
+637.0	-0.0434112548828
+636.5	-0.0435180664062
+636.0	-0.0437469482422
+635.5	-0.0439300537109
+635.0	-0.0440673828125
+634.5	-0.0442199707031
+634.0	-0.0443725585938
+633.5	-0.0446014404297
+633.0	-0.0447845458984
+632.5	-0.0447845458984
+632.0	-0.0449676513672
+631.5	-0.0453033447266
+631.0	-0.0455169677734
+630.5	-0.0456390380859
+630.0	-0.0459289550781
+629.5	-0.0463409423828
+629.0	-0.0467224121094
+628.5	-0.0469207763672
+628.0	-0.047119140625
+627.5	-0.0473175048828
+627.0	-0.0475616455078
+626.5	-0.0479736328125
+626.0	-0.0483245849609
+625.5	-0.0486297607422
+625.0	-0.0487976074219
+624.5	-0.0489044189453
+624.0	-0.0491180419922
+623.5	-0.0493621826172
+623.0	-0.0496520996094
+622.5	-0.0497894287109
+622.0	-0.0499267578125
+621.5	-0.0502166748047
+621.0	-0.0505218505859
+620.5	-0.0509185791016
+620.0	-0.0512084960938
+619.5	-0.0514068603516
+619.0	-0.0516662597656
+618.5	-0.0517883300781
+618.0	-0.0519409179688
+617.5	-0.0520477294922
+617.0	-0.0523376464844
+616.5	-0.052734375
+616.0	-0.0528564453125
+615.5	-0.052978515625
+615.0	-0.0531921386719
+614.5	-0.0533447265625
+614.0	-0.0535125732422
+613.5	-0.0539093017578
+613.0	-0.0541839599609
+612.5	-0.0542907714844
+612.0	-0.0542755126953
+611.5	-0.0542297363281
+611.0	-0.0544586181641
+610.5	-0.0547332763672
+610.0	-0.0550079345703
+609.5	-0.05517578125
+609.0	-0.0550384521484
+608.5	-0.0550689697266
+608.0	-0.0551452636719
+607.5	-0.05517578125
+607.0	-0.0553283691406
+606.5	-0.0553894042969
+606.0	-0.0554809570312
+605.5	-0.0555267333984
+605.0	-0.0554809570312
+604.5	-0.0554656982422
+604.0	-0.0554809570312
+603.5	-0.0554504394531
+603.0	-0.0554809570312
+602.5	-0.0554809570312
+602.0	-0.0554504394531
+601.5	-0.0554504394531
+601.0	-0.055419921875
+600.5	-0.0553894042969
+600.0	-0.0553131103516
+599.5	-0.0550994873047
+599.0	-0.0549926757812
+598.5	-0.0550842285156
+598.0	-0.0549468994141
+597.5	-0.0548095703125
+597.0	-0.0549163818359
+596.5	-0.0549774169922
+596.0	-0.0549621582031
+595.5	-0.0548858642578
+595.0	-0.0548400878906
+594.5	-0.0549011230469
+594.0	-0.0547637939453
+593.5	-0.0545654296875
+593.0	-0.0544891357422
+592.5	-0.0546112060547
+592.0	-0.0548095703125
+591.5	-0.0548400878906
+591.0	-0.0547637939453
+590.5	-0.0547332763672
+590.0	-0.0547180175781
+589.5	-0.0545043945312
+589.0	-0.0545043945312
+588.5	-0.0546112060547
+588.0	-0.0546722412109
+587.5	-0.0546112060547
+587.0	-0.0544891357422
+586.5	-0.0545959472656
+586.0	-0.0546264648438
+585.5	-0.0545349121094
+585.0	-0.0544281005859
+584.5	-0.0543670654297
+584.0	-0.0544128417969
+583.5	-0.0544738769531
+583.0	-0.0543975830078
+582.5	-0.0544738769531
+582.0	-0.0547027587891
+581.5	-0.0548248291016
+581.0	-0.0548858642578
+580.5	-0.0550231933594
+580.0	-0.0549926757812
+579.5	-0.054931640625
+579.0	-0.0549774169922
+578.5	-0.0548858642578
+578.0	-0.0549621582031
+577.5	-0.0551147460938
+577.0	-0.0550689697266
+576.5	-0.0550689697266
+576.0	-0.0552215576172
+575.5	-0.0552368164062
+575.0	-0.0553894042969
+574.5	-0.055419921875
+574.0	-0.0554656982422
+573.5	-0.0555572509766
+573.0	-0.0556030273438
+572.5	-0.0557708740234
+572.0	-0.0557861328125
+571.5	-0.0559387207031
+571.0	-0.0561981201172
+570.5	-0.0562591552734
+570.0	-0.0562286376953
+569.5	-0.0562438964844
+569.0	-0.0563049316406
+568.5	-0.0565185546875
+568.0	-0.0567626953125
+567.5	-0.0567932128906
+567.0	-0.0569152832031
+566.5	-0.0571136474609
+566.0	-0.0572509765625
+565.5	-0.0573120117188
+565.0	-0.057373046875
+564.5	-0.0575256347656
+564.0	-0.0576019287109
+563.5	-0.0576629638672
+563.0	-0.0577545166016
+562.5	-0.0579223632812
+562.0	-0.05810546875
+561.5	-0.0581970214844
+561.0	-0.0583343505859
+560.5	-0.05859375
+560.0	-0.0587768554688
+559.5	-0.0587768554688
+559.0	-0.0588836669922
+558.5	-0.0590515136719
+558.0	-0.0592041015625
+557.5	-0.0593566894531
+557.0	-0.0593872070312
+556.5	-0.0593872070312
+556.0	-0.0596313476562
+555.5	-0.0598907470703
+555.0	-0.0599670410156
+554.5	-0.0599975585938
+554.0	-0.0602416992188
+553.5	-0.0605621337891
+553.0	-0.0606384277344
+552.5	-0.0606384277344
+552.0	-0.0608978271484
+551.5	-0.06103515625
+551.0	-0.0611724853516
+550.5	-0.0613098144531
+550.0	-0.0613555908203
+549.5	-0.0615692138672
+549.0	-0.0617370605469
+548.5	-0.0620880126953
+548.0	-0.0623168945312
+547.5	-0.0623779296875
+547.0	-0.0625610351562
+546.5	-0.0626831054688
+546.0	-0.0627899169922
+545.5	-0.0629577636719
+545.0	-0.0629577636719
+544.5	-0.0628356933594
+544.0	-0.0630035400391
+543.5	-0.0634002685547
+543.0	-0.0637664794922
+542.5	-0.0638580322266
+542.0	-0.0639190673828
+541.5	-0.0641937255859
+541.0	-0.064453125
+540.5	-0.0646209716797
+540.0	-0.064697265625
+539.5	-0.0648345947266
+539.0	-0.0650939941406
+538.5	-0.0653533935547
+538.0	-0.0657196044922
+537.5	-0.0659637451172
+537.0	-0.0661773681641
+536.5	-0.0663909912109
+536.0	-0.0666351318359
+535.5	-0.0668182373047
+535.0	-0.0668182373047
+534.5	-0.0670013427734
+534.0	-0.0672912597656
+533.5	-0.0674438476562
+533.0	-0.0676422119141
+532.5	-0.0679168701172
+532.0	-0.0682678222656
+531.5	-0.068603515625
+531.0	-0.0687866210938
+530.5	-0.0689392089844
+530.0	-0.0691375732422
+529.5	-0.0694274902344
+529.0	-0.0698394775391
+528.5	-0.0701599121094
+528.0	-0.0704803466797
+527.5	-0.0706634521484
+527.0	-0.0707244873047
+526.5	-0.0707855224609
+526.0	-0.0707855224609
+525.5	-0.0710754394531
+525.0	-0.0714263916016
+524.5	-0.0717010498047
+524.0	-0.0720672607422
+523.5	-0.0724182128906
+523.0	-0.0727081298828
+522.5	-0.0728454589844
+522.0	-0.0729370117188
+521.5	-0.0732116699219
+521.0	-0.0736083984375
+520.5	-0.0739898681641
+520.0	-0.07421875
+519.5	-0.0745239257812
+519.0	-0.0749206542969
+518.5	-0.0751647949219
+518.0	-0.0755462646484
+517.5	-0.0760345458984
+517.0	-0.0762939453125
+516.5	-0.0765991210938
+516.0	-0.0769500732422
+515.5	-0.0771179199219
+515.0	-0.0774230957031
+514.5	-0.0778503417969
+514.0	-0.0781097412109
+513.5	-0.0782775878906
+513.0	-0.0787353515625
+512.5	-0.0791473388672
+512.0	-0.0795593261719
+511.5	-0.0801086425781
+511.0	-0.0804595947266
+510.5	-0.0807037353516
+510.0	-0.0811004638672
+509.5	-0.0814666748047
+509.0	-0.0818023681641
+508.5	-0.0823669433594
+508.0	-0.0829467773438
+507.5	-0.0834655761719
+507.0	-0.0838623046875
+506.5	-0.0840606689453
+506.0	-0.0844268798828
+505.5	-0.0850219726562
+505.0	-0.0853118896484
+504.5	-0.0855255126953
+504.0	-0.0860137939453
+503.5	-0.0865783691406
+503.0	-0.0870971679688
+502.5	-0.0875854492188
+502.0	-0.0880889892578
+501.5	-0.0884552001953
+501.0	-0.0889892578125
+500.5	-0.0895385742188
+500.0	-0.0898895263672
+499.5	-0.0904235839844
+499.0	-0.0909576416016
+498.5	-0.0914916992188
+498.0	-0.0919952392578
+497.5	-0.0923156738281
+497.0	-0.0927581787109
+496.5	-0.0934448242188
+496.0	-0.0940704345703
+495.5	-0.0945281982422
+495.0	-0.0949249267578
+494.5	-0.0954132080078
+494.0	-0.0959930419922
+493.5	-0.0966033935547
+493.0	-0.0971832275391
+492.5	-0.0976867675781
+492.0	-0.0980987548828
+491.5	-0.0986480712891
+491.0	-0.0993194580078
+490.5	-0.0999298095703
+490.0	-0.100341796875
+489.5	-0.100692749023
+489.0	-0.101272583008
+488.5	-0.101989746094
+488.0	-0.10270690918
+487.5	-0.103302001953
+487.0	-0.103652954102
+486.5	-0.103851318359
+486.0	-0.104293823242
+485.5	-0.10498046875
+485.0	-0.105590820312
+484.5	-0.10612487793
+484.0	-0.106811523438
+483.5	-0.107482910156
+483.0	-0.108047485352
+482.5	-0.108627319336
+482.0	-0.109115600586
+481.5	-0.109512329102
+481.0	-0.110214233398
+480.5	-0.110946655273
+480.0	-0.111618041992
+479.5	-0.112213134766
+479.0	-0.112686157227
+478.5	-0.113342285156
+478.0	-0.113754272461
+477.5	-0.114242553711
+477.0	-0.114898681641
+476.5	-0.115386962891
+476.0	-0.115966796875
+475.5	-0.116577148438
+475.0	-0.117141723633
+474.5	-0.117889404297
+474.0	-0.118774414062
+473.5	-0.119400024414
+473.0	-0.120071411133
+472.5	-0.120819091797
+472.0	-0.121215820312
+471.5	-0.121612548828
+471.0	-0.122283935547
+470.5	-0.122787475586
+470.0	-0.123580932617
+469.5	-0.124481201172
+469.0	-0.125076293945
+468.5	-0.125732421875
+468.0	-0.126358032227
+467.5	-0.127075195312
+467.0	-0.127731323242
+466.5	-0.128448486328
+466.0	-0.129043579102
+465.5	-0.129669189453
+465.0	-0.130386352539
+464.5	-0.131088256836
+464.0	-0.131942749023
+463.5	-0.13249206543
+463.0	-0.133026123047
+462.5	-0.133743286133
+462.0	-0.134613037109
+461.5	-0.135406494141
+461.0	-0.135757446289
+460.5	-0.136077880859
+460.0	-0.136749267578
+459.5	-0.137390136719
+459.0	-0.137954711914
+458.5	-0.138717651367
+458.0	-0.139450073242
+457.5	-0.140075683594
+457.0	-0.140533447266
+456.5	-0.14111328125
+456.0	-0.141677856445
+455.5	-0.142105102539
+455.0	-0.142639160156
+454.5	-0.143127441406
+454.0	-0.143630981445
+453.5	-0.143981933594
+453.0	-0.14453125
+452.5	-0.145278930664
+452.0	-0.145660400391
+451.5	-0.145874023438
+451.0	-0.146133422852
+450.5	-0.14665222168
+450.0	-0.147125244141
+449.5	-0.147766113281
+449.0	-0.148391723633
+448.5	-0.148849487305
+448.0	-0.14909362793
+447.5	-0.149276733398
+447.0	-0.149612426758
+446.5	-0.149688720703
+446.0	-0.150100708008
+445.5	-0.150772094727
+445.0	-0.1513671875
+444.5	-0.151641845703
+444.0	-0.151840209961
+443.5	-0.152008056641
+443.0	-0.152221679688
+442.5	-0.152725219727
+442.0	-0.153091430664
+441.5	-0.153274536133
+441.0	-0.153625488281
+440.5	-0.154159545898
+440.0	-0.154571533203
+439.5	-0.154846191406
+439.0	-0.155059814453
+438.5	-0.155334472656
+438.0	-0.155456542969
+437.5	-0.155578613281
+437.0	-0.155838012695
+436.5	-0.155990600586
+436.0	-0.156372070312
+435.5	-0.156677246094
+435.0	-0.15657043457
+434.5	-0.156814575195
+434.0	-0.157196044922
+433.5	-0.157592773438
+433.0	-0.158081054688
+432.5	-0.158493041992
+432.0	-0.158874511719
+431.5	-0.158981323242
+431.0	-0.159210205078
+430.5	-0.159515380859
+430.0	-0.159774780273
+429.5	-0.160049438477
+429.0	-0.160491943359
+428.5	-0.160720825195
+428.0	-0.160842895508
+427.5	-0.161346435547
+427.0	-0.161544799805
+426.5	-0.161727905273
+426.0	-0.162002563477
+425.5	-0.162094116211
+425.0	-0.162338256836
+424.5	-0.162780761719
+424.0	-0.163131713867
+423.5	-0.163162231445
+423.0	-0.163360595703
+422.5	-0.163543701172
+422.0	-0.163528442383
+421.5	-0.163955688477
+421.0	-0.164520263672
+420.5	-0.164794921875
+420.0	-0.164978027344
+419.5	-0.165176391602
+419.0	-0.165161132812
+418.5	-0.165496826172
+418.0	-0.166229248047
+417.5	-0.16650390625
+417.0	-0.16650390625
+416.5	-0.166702270508
+416.0	-0.167053222656
+415.5	-0.167221069336
+415.0	-0.167419433594
+414.5	-0.167907714844
+414.0	-0.168426513672
+413.5	-0.168518066406
+413.0	-0.168701171875
+412.5	-0.168899536133
+412.0	-0.168930053711
+411.5	-0.169418334961
+411.0	-0.169723510742
+410.5	-0.169677734375
+410.0	-0.170059204102
+409.5	-0.170501708984
+409.0	-0.1708984375
+408.5	-0.171157836914
+408.0	-0.171035766602
+407.5	-0.171188354492
+407.0	-0.171401977539
+406.5	-0.172073364258
+406.0	-0.172592163086
+405.5	-0.172302246094
+405.0	-0.17253112793
+404.5	-0.173065185547
+404.0	-0.173446655273
+403.5	-0.173568725586
+403.0	-0.173843383789
+402.5	-0.174255371094
+402.0	-0.174118041992
+401.5	-0.174118041992
+401.0	-0.174423217773
+400.5	-0.174514770508
+400.0	-0.174560546875
+399.5	-0.175003051758
+399.0	-0.175216674805
+398.5	-0.175018310547
+398.0	-0.175094604492
+397.5	-0.17529296875
+397.0	-0.175354003906
+396.5	-0.175598144531
+396.0	-0.176025390625
+395.5	-0.176147460938
+395.0	-0.17610168457
+394.5	-0.175888061523
+394.0	-0.175598144531
+393.5	-0.175750732422
+393.0	-0.175689697266
+392.5	-0.17578125
+392.0	-0.176055908203
+391.5	-0.17610168457
+391.0	-0.175888061523
+390.5	-0.175598144531
+390.0	-0.175445556641
+389.5	-0.17512512207
+389.0	-0.174865722656
+388.5	-0.174591064453
+388.0	-0.174728393555
+387.5	-0.174713134766
+387.0	-0.174575805664
+386.5	-0.174758911133
+386.0	-0.174362182617
+385.5	-0.173950195312
+385.0	-0.173934936523
+384.5	-0.173614501953
+384.0	-0.173446655273
+383.5	-0.173309326172
+383.0	-0.172821044922
+382.5	-0.172698974609
+382.0	-0.172302246094
+381.5	-0.171524047852
+381.0	-0.171249389648
+380.5	-0.170944213867
+380.0	-0.17057800293
+379.5	-0.170608520508
+379.0	-0.170394897461
+378.5	-0.169998168945
+378.0	-0.169448852539
+377.5	-0.169067382812
+377.0	-0.168701171875
+376.5	-0.167953491211
+376.0	-0.167663574219
+375.5	-0.167053222656
+375.0	-0.166381835938
+374.5	-0.166320800781
+374.0	-0.166152954102
+373.5	-0.165313720703
+373.0	-0.164703369141
+372.5	-0.164443969727
+372.0	-0.163848876953
+371.5	-0.163345336914
+371.0	-0.163314819336
+370.5	-0.162902832031
+370.0	-0.161773681641
+369.5	-0.161178588867
+369.0	-0.160690307617
+368.5	-0.15998840332
+368.0	-0.159683227539
+367.5	-0.159378051758
+367.0	-0.158752441406
+366.5	-0.158599853516
+366.0	-0.157974243164
+365.5	-0.157135009766
+365.0	-0.157073974609
+364.5	-0.156356811523
+364.0	-0.156295776367
+363.5	-0.155807495117
+363.0	-0.15478515625
+362.5	-0.154479980469
+362.0	-0.153701782227
+361.5	-0.153381347656
+361.0	-0.152908325195
+360.5	-0.152008056641
+360.0	-0.151153564453
+359.5	-0.150650024414
+359.0	-0.150665283203
+358.5	-0.150283813477
+358.0	-0.149703979492
+357.5	-0.149353027344
+357.0	-0.149215698242
+356.5	-0.148742675781
+356.0	-0.148086547852
+355.5	-0.147888183594
+355.0	-0.147842407227
+354.5	-0.147186279297
+354.0	-0.146728515625
+353.5	-0.146865844727
+353.0	-0.146514892578
+352.5	-0.145935058594
+352.0	-0.145614624023
+351.5	-0.145431518555
+351.0	-0.145126342773
+350.5	-0.145172119141
+350.0	-0.144927978516
+349.5	-0.144561767578
+349.0	-0.145294189453
+348.5	-0.144836425781
+348.0	-0.14421081543
+347.5	-0.144195556641
+347.0	-0.144134521484
+346.5	-0.144302368164
+346.0	-0.143905639648
+345.5	-0.144073486328
+345.0	-0.144485473633
+344.5	-0.144927978516
+344.0	-0.144989013672
+343.5	-0.144760131836
+343.0	-0.14518737793
+342.5	-0.144180297852
+342.0	-0.143493652344
+341.5	-0.144622802734
+341.0	-0.144500732422
+340.5	-0.144424438477
+340.0	-0.145004272461
+339.5	-0.145431518555
+339.0	-0.144989013672
+338.5	-0.144927978516
+338.0	-0.145538330078
+337.5	-0.144897460938
+337.0	-0.145202636719
+336.5	-0.146041870117
+336.0	-0.145812988281
+335.5	-0.14582824707
+335.0	-0.146545410156
+334.5	-0.146347045898
+334.0	-0.145935058594
+333.5	-0.146713256836
+333.0	-0.146759033203
+332.5	-0.146926879883
+332.0	-0.147964477539
+331.5	-0.148300170898
+331.0	-0.149078369141
+330.5	-0.149490356445
+330.0	-0.149627685547
+329.5	-0.149703979492
+329.0	-0.149398803711
+328.5	-0.149475097656
+328.0	-0.149688720703
+327.5	-0.149826049805
+327.0	-0.150329589844
+326.5	-0.151077270508
+326.0	-0.150863647461
+325.5	-0.150405883789
+325.0	-0.151123046875
+324.5	-0.152542114258
+324.0	-0.153549194336
+323.5	-0.152984619141
+323.0	-0.152328491211
+322.5	-0.152984619141
+322.0	-0.1533203125
+321.5	-0.153076171875
+321.0	-0.153396606445
+320.5	-0.154220581055
+320.0	-0.153854370117
+319.5	-0.15412902832
+319.0	-0.154754638672
+318.5	-0.154876708984
+318.0	-0.155548095703
+317.5	-0.156066894531
+317.0	-0.156341552734
+316.5	-0.156860351562
+316.0	-0.15705871582
+315.5	-0.157608032227
+315.0	-0.15803527832
+314.5	-0.158050537109
+314.0	-0.159103393555
+313.5	-0.159805297852
+313.0	-0.160079956055
+312.5	-0.159576416016
+312.0	-0.159225463867
+311.5	-0.160278320312
+311.0	-0.16162109375
+310.5	-0.162170410156
+310.0	-0.16194152832
+309.5	-0.162094116211
+309.0	-0.162704467773
+308.5	-0.162963867188
+308.0	-0.162612915039
+307.5	-0.163131713867
+307.0	-0.164779663086
+306.5	-0.165588378906
+306.0	-0.165725708008
+305.5	-0.166290283203
+305.0	-0.166442871094
+304.5	-0.16667175293
+304.0	-0.167343139648
+303.5	-0.167465209961
+303.0	-0.168090820312
+302.5	-0.168838500977
+302.0	-0.168807983398
+301.5	-0.16926574707
+301.0	-0.169967651367
+300.5	-0.170684814453
+300.0	-0.171188354492