def plot_crd_z_hexd(data1,dtc,outdir,zpbc): from functionsdistance import plot_vs_time print data1.shape coor = data1[:,::2] toplot = coor plot_vs_time(toplot,dtc,outdir+"/fig_crd-z.png") coor -= zpbc*np.floor(data1[:,::2]/zpbc+0.5) toplot = add_constant_toarray( coor[:,::2], [zpbc/2.,-zpbc/2.]) kwargs = {"linewidth":0,"marker":"o","markersize":1,"markeredgewidth":0} plot_vs_time(toplot,dtc,outdir+"/fig_pbccrd-z.png",**kwargs)
def plot_crd_z(data1,data2,dtc,outdir,zpbc): print "data", data1.shape, data2.shape from functionsdistance import plot_vs_time toplot = add_vec_toarray(data1[:,::2],[data2[:,0],data2[:,0]+zpbc/2.,data2[:,0]-zpbc/2.]) plot_vs_time(toplot,dtc,outdir+"/fig_crd-z.png") coor = shift_wrt_layer(data1,data2) toplot = add_constant_toarray(coor[:,::2], [zpbc/2.,-zpbc/2.]) plot_vs_time(toplot,dtc,outdir+"/fig_crd-z.shift.png") coor -= zpbc*np.floor(coor/zpbc+0.5) toplot = add_constant_toarray( coor[:,::2], [zpbc/2.,-zpbc/2.]) kwargs = {"linewidth":0,"marker":"o","markersize":1,"markeredgewidth":0} plot_vs_time(toplot,dtc,outdir+"/fig_pbccrd-z.shift.png",**kwargs)