Ejemplo n.º 1
0
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)
Ejemplo n.º 2
0
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)