Non-numeric pressure - simulation unstable

Dear All,

I run a simulation and encounter following error: Non-numeric pressure -
simulation unstable. I use FeNiCr.eam.alloy potential and read the model
from a input file.

The detail information as following:

# ------------------------ INITIALIZATION ----------------------------
units metal
dimension 3
boundary p p p
atom_style atomic

# ----------------------- ATOM DEFINITION ----------------------------

read_data "Three_elements.data" group all

# ------------------------ FORCE FIELDS ------------------------------

pair_style eam/alloy
pair_coeff * * FeNiCr_2013.eam.alloy Fe Ni Cr

# ------------------------- SETTINGS ---------------------------------

compute peratom all pe/atom
compute virial all stress/atom NULL

#---------------------- EQUILIBRATION ----------------------------

# ---------------------set timestep
reset_timestep 0
timestep 0.001

# ---------------set temperature

velocity all create 300 12345

# Assign fixes

fix 1 all npt temp 300 300 0.5 iso 0 0 5 drag 1

# Set thermo output

thermo 400
thermo_style custom step lx ly lz press pxx pyy pzz pe temp

# ------------Use custom dump for Ovito output

# Run for at least 20 picosecond (assuming 1 fs timestep)

run 20000

Dear All,

I run a simulation and encounter following error: Non-numeric pressure -
simulation unstable. I use FeNiCr.eam.alloy potential and read the model
from a input file.

‚Äčthis is a frequently discussed error message. please consult the LAMMPS
manual and the (many) previous discussions in the LAMMPS mailing‚Äč archives
for suggestions of how to address this and report back only in case none of
the suggested remedies apply.

axel.