read_restart confined.restart 

variable    nb_points equal 1000 
variable    eps_max equal 0.1 
variable    shearrate equal 1e-1
variable    dt equal 0.01

# Don't touch
variable    eps equal xy/lx
variable    rate equal v_shearrate*lx*v_dt
variable    nb_step equal $(round(v_eps_max/(v_rate*v_dt)))
variable    print_freq equal $(round(v_nb_step/v_nb_points))

# ===== GEOMETRY ====  

change_box all triclinic

# ==== DUMP ====

thermo_style       	custom step atoms lx xy pxy v_eps v_nb_step
thermo_modify      	norm no flush yes
thermo             	5000

# ==== FIXES ====

timestep ${dt}

fix 2 all deform 1 xy erate ${rate}
fix FIX5 all enforce2d

fix output_file all print ${print_freq} "${eps} $(pxy)" file StrainStress.txt screen no

run ${nb_step}

write_data box_sheared.data