dimension 3 units lj boundary p p p pair_style morse 4.0 variable a0 equal 4.0 variable T equal 0.25 variable T2 equal 0.3 variable R equal 20 variable NE equal 100 variable NS equal 500 read_data resume.data pair_coeff * * 1.0 ${a0} 1.0 pair_modify shift yes mass 1 1.0 # Dynamics: fix step all nve fix temp all langevin ${T2} ${T2} 1.0 1238914 gjf yes fix noro all momentum 1 angular rescale # Commenting this line removes issue... velocity all create $T 1249 thermo 10 run ${NE} print "DONE!"