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)