In [1]:

    

from ZarcfitCalculations import ZarcfitCalculations, SetDefaultParameters
import numpy as np
import pandas as pd
import matplotlib

%pylab inline


    

    




Populating the interactive namespace from numpy and matplotlib

In [5]:

    

path = "../data/HVC2014_10Grenon/"
fnameobs = "BC13867-A 2014-10-23.z"
fnamepred = "HVC 2014-10-27simplified.csv"
pathobs = path+fnameobs
pathpred = path+fnamepred
temp = np.loadtxt(pathobs, skiprows=11, delimiter=",")
data = pd.read_csv(pathpred)
data_active = data.loc[np.logical_and(data["Fl"]<20., data["Fm"]>1.)]


    

In [17]:

    

obs = temp[:,4]+1j*temp[:,5]
frequency = temp[:,0].copy()
nfreq = len(frequency)
zarc = ZarcfitCalculations(obs, frequency)
Linf, Rinf, Rh, Fh, Ph, Rl, Fl, Pl, Rm, Fm, Pm, Re, Qe, Pef, Pei  = SetDefaultParameters()
zarc.SetParametersSeries(Linf, Rinf, Rh, Fh, Ph, Rl, Fl, Pl, Rm, Fm, Pm, Re, Qe, Pef, Pei)


    

In [ ]:

    




    

In [ ]:

    

zarc = ZarcfitCalculations()