def main_func():
conn = general.open_conn()
cmps = conn.execute("select comp_key from comps where function = 'phi6'").fetchall();
for c in cmps:
(fname,) = conn.execute("select fout from comps where comp_key = ? ;",c).fetchone()
print fname
print c
f = h5py.File(fname,'r+')
nf = f.attrs["number-of-planes"]
del f["/parameters/scaler_order_parameter_%(#)07d"%{"#":c[0]}]
del f["/parameters/neighborhood_size_%(#)07d"%{"#":c[0]}]
for fr in range(0,nf):
del f["/frame%(#)06d"%{"#":fr}+"/neighborhood_size_%(#)07d"%{"#":c[0]}]
del f["/frame%(#)06d"%{"#":fr}+"/scaler_order_parameter_%(#)07d"%{"#":c[0]}]
pass
import lib.general as gen
import lib.msd as lm
import lib.plots as plts
conn = gen.open_conn("/home/tcaswell/colloids/proc_db.db")
ce = conn.execute
steps = [5,10,15,20,25]
styles = ['x','^','o','<','>']
msds = [ce("select comp_key from msd where date = '2010-12-08' and msd_steps = ?",(c,)).fetchall()
for c in steps]
fits = [ [lm.fit_msd(c[0],conn) for c in m]
for m in msds]
ax = set_up_axis('T [C]',r'r [$\mu$m]','NIPAM radius vs Temperature')
for fts,sty,lab in zip(fits,styles,steps):
ax.plot([f[1] for f in fts],[f[2] for f in fts],sty,label=str(lab) + ' steps')
fig.save('test')
def main_loop():
conn = gen.open_conn();
gofr_date(conn,'2010-05-24')
conn.close()
def main_loop():
conn = gen.open_conn();
iden_date(conn,'2010-05-24')
def main_loop():
conn = gen.open_conn();
gofr_group(conn,'2010-04-26-2')
gofr_group(conn,'2010-04-26-6-thin')
gofr_group(conn,'2010-04-26-6-q')
