예제 #1
0
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
예제 #2
0
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
예제 #3
0
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
예제 #4
0
def main():
    data = np.genfromtxt(args.ts)
    print mean_stderr(data[:,1][::100])