variable t index 300 350 400 450 500 550 600 label loop read_restart restart.minimized.lammps log log.equillib.lammps fix 1 all box/relax aniso 1.0 include mod.neigh-output-minsets.lammps #dump mov all movie 1000 "equillib-$t.mpg" type type bond type 0.1 adiam 0.5 view 45 45 velocity all create $t 123 fix 2 all nvt temp $t $t 100.0 #timestep 0.1 run 50000 unfix 2 fix 2 all npt temp $t $t 100.0 aniso 1.0 1.0 350.0 run 400000 minimize ${etol} ${ftol} ${maxiter} ${maxeval} write_restart restart.equilibd-minimizd-$t.lammps #write_data data.equilibd-minimizd-$t.lammps clear next t jump in.equillib.lammps loop