sdf

icartt/icartt.py

@@ -174,6 +174,26 @@ def read(ict_file):
 
     return(out)
 
+def read_first_and_last(ict_file):
+
+    out = read_header(ict_file)
+
+    f = open( ict_file, 'r' )
+
+    nul = [f.readline() for i in xrange(out['nheader'])]
+
+    def nan_miss_float(raw, vmiss):
+        return [ float(x.replace(vmiss[i], 'NaN')) for i, x in enumerate(raw) ]
+
+    out['data']  = [ nan_miss_float(f.readline().split(','), out['vmiss']) ]
+    for line in f:
+        pass
+    last = line
+    out['data'] += [ nan_miss_float(last.split(','), out['vmiss']) ]
+
+    f.close()
+
+    return(out)