def get_vl(data, vdir):
t = data[:,0]
# lazy for now
vpl = (data[:,3] - (n_h2o * v_w[wm][str(temp)]))/n_lip
plt_save(t, vpl, vdir + 'vl')
vl_deets = mean_stderr(vpl)
return vpl, vl_deets
def kap_ts(time, alts, vdir):
# correct units
al_ts = alts * 1e-18
kap_ts = [0]
for t in range(len(al_ts)):
if t == 0:
continue
else:
kap_ts.append(get_kap(al_ts[:t]))
kap = np.array(kap_ts)
# flucs are too crazy to plot beginning
s0 = len(kap)*.25
plt_save(time[s0:], kap[s0:], vdir + 'kap')
kap_deets = mean_stderr(kap)
return kap[-1], kap_deets
def get_al(data, vdir):
t = data[:,0]
apl = data[:,1]*data[:,2]/(n_lip/2.0)
plt_save(t, apl, vdir + 'al')
al_deets = mean_stderr(apl)
return apl, al_deets
def main():
data = np.genfromtxt(args.ts)
print mean_stderr(data[:,1][::100])
