#################################
variable up equal 0.5e-2
units           metal
boundary        p p p
variable        L equal 8
#lattice         diamond 5.431 orient x 1 -1 0 orient y 1 1 -2 orient z 1 1 1
lattice         diamond 5.431
region          regbox prism 0 ${L} 0 ${L} 0 ${L} 0.0 0.0 0.0
create_box      1 regbox
create_atoms    1 box
thermo   2000
thermo_modify flush yes
################## Potential Definition ######################################################
mass 1 1
pair_style meam
pair_coeff * * library.meam Si NULL Si
neighbor        2 bin
neigh_modify once no every 1 delay 0 check yes
variable T2 equal 2*$T
################### Relax the box #############################################################
fix r1 all box/relax  aniso 0.0
minimize 0.0 1.0e-11 100000 10000000
################### Temperature Control #######################################################
compute           mobile_temp all temp
velocity          all create ${T2} 5812775 temp mobile_temp
fix               1 all npt temp ${T} ${T} 0.1 iso 0 0 3 drag 1
fix               2 all temp/rescale 10 $T $T .01 1.0 
fix_modify        2 temp mobile_temp
################### Display and Output control ############################################### 
thermo_style custom step pxx pyy pzz pxy pxz pyz lx ly lz xy xz yz
run 400000
unfix 1


#################### Strain Loop Test #########################################################
variable        loop1  equal 5
variable        n loop ${loop1}
label strain
variable        lambda equal ${up}*$L
displace_box    all yz delta ${lambda}
fix 1 all nvt temp $T $T 0.1 drag 1 
run 400000
next n
variable l delete 
jump in.elastic_NVT_Shear strain
###############################################################################################