
echo | awk 'END{n=1000;L=20; for(cop=1; cop<5;cop++){print n; print ""; for(i=1;i<=n;i++){print "C",L*rand(),L*rand(),L*rand()}}}' > random_4.xyz
#less NaCl_basis.xyz
xyz_gofr.py random_4.xyz --dr 0.4 -o -v 8000 > uniform.gofr &

xyz_lattice_translate -a1 1.5 1.5 0 -a2 1.5 0 1.5 -a3 0 1.5 1.5 NaCl_basis.xyz -n 40 > NaCl_chunk.xyz
cut_cube_in_xyz NaCl_chunk.xyz -s 24 -x 0.3 -y 0.3 -z 0.3 > NaCl.xyz
xyz_add_random.py -n 5 -s 0.02 NaCl.xyz > NaCl_5.xyz
xyz_gofr.py NaCl_5.xyz --dr 0.02 -o -v 13824 > obc.gofr &
#xyz_gofr.py NaCl_5.xyz --dr 0.02 --a1 24 0 0 --a2 0 24 0 --a3 0 0 24 > pbc.gofr &
