In [93]:

    

import numpy as np
from pydelft.grd import grd
from pydelft.bnd import bnd
from pydelft.bct import bct
from PyQt4 import QtGui
import sys
import os
import pandas as pd
import linecache


    

In [2]:

    

b = bct(fname=r'i:\Japan\Data\Sendai\T3\create_riemann\testing\run7\natori_10_v2.bct')
b.filename


    

    
Out[2]:





'i:\\Japan\\Data\\Sendai\\T3\\create_riemann\\testing\\run7\\natori_10_v2.bct'

In [110]:

    

f = open(b.filename,'r')
names = []
contents = []
locations = []
timefunctions = []
referencetimes = []
timeunits = []
interpolations = []
parameters = []
records = []
data = []

i = 0
while
for i, line in enumerate(f):
    if 'table-name' in line:
        names.append(line.split("'")[1])
    elif 'contents' in line:
        contents.append(line.split("'")[1])
    elif 'location' in line:
        locations.append(line.split("'")[1])
    elif 'time-function' in line:
        timefunctions.append(line.split("'")[1])
    elif 'reference-time' in line:
        referencetimes.append(int(line.split()[1]))
    elif 'time-unit' in line:
        timeunits.append(line.split("'")[1])
    elif 'interpolation' in line:
        interpolations.append(line.split("'")[1])
    elif 'parameter' in line:
        parameters.append([line.split("'")[1], line.split("'")[3]])
    elif 'records-in-table' in line:
        records.append(int(line.split()[1]))
    else:
        break   
        
bcts = []    
for line in range(i+1,i+1+records[0]):
        time_1_2 = linecache.getline(b.filename, line).split()
        bcts.append((float(time_1_2[0]),float(time_1_2[1]),
                     float(time_1_2[2])))
data.append(np.array(bcts).T)

i = line+1

f.close()


    

In [112]:

    

i


    

    
Out[112]:





14

In [92]:

    

linecache.getline(b.filename, i+1).split()


    

    
Out[92]:





['0.00', '0.00000e+00', '0.00000e+00']

In [ ]: