binflow

Hello every body.
I would be thankful if you could help me.
I am trying to run Binflow ,However there are some errors :(LIGGGHTS (Version LIGGGHTS-PUBLIC 3.3.1, compiled 2016-01-14-14:02:34 by kloss, git commit 9b1fc57b77c483deaf2c0a6c69241c5d4f91e3c2 based on LAMMPS 23 Nov 2013)
Created orthogonal box = (-0.138 -0.138 -0.0045) to (0.138 0.138 0.43)
1 by 1 by 1 MPI processor grid
ERROR: Fix particletemplate/sphere (id pts): unrecognized keyword (/build/liggghts-QCz4in/liggghts-3.3.1+repack1/src/fix_template_sphere.cpp:277)

***********************************
#binflow

## Initialization
# Preliminaries

units si
atom_style sphere
atom_modify map array
boundary f f f
newton off
communicate single vel yes

#processors 2 2 3
# Declare domain

region domain block -0.138 0.138 -0.138 0.138 -0.0045 0.43 units box
create_box 2 domain

# Neighbor listing

neighbor 0.003 bin
neigh_modify every 1 check no

### Setup

#Material and interaction properties required

fix m1 all property/global youngsModulus peratomtype 2.5e7 2.5e7
fix m2 all property/global poissonsRatio peratomtype 0.25 0.25
fix m3 all property/global coefficientRestitution peratomtypepair 2 0.5 0.5 0.5 0.5
fix m4 all property/global coefficientFriction peratomtypepair 2 0.2 0.175 0.175 0.5

#Particle insertion

fix pts all particletemplate/sphere 1 atom_type 1 density constant 1000 & radius constant 0.0015

fix pdd all particledistribution/discrete 63243 1 pts 1.0
region factory cylinder z 0 0 0.10 0.4 0.41 units box

fix ins all insert/rate/region seed 123481 distributiontemplate pdd &
nparticles 150000 particlerate 150000 insert_every 1000 &
overlapcheck yes vel constant 0. 0. -1.0 regions factory ntry_mc 10000

# Import mesh from cad:

fix cad1 all mesh/surface file hopper.stl type 2 scale 0.001
# Use the imported mesh as granular wall

fix geometry all wall/gran/hertz/history mesh n_meshes 1 meshes cad1

# Create stopper for funnel
fix stopper all wall/gran/hertz/history primitive type 2
zplane 0.0

#Define the physics
pair_style gran/hertz/history
pair_coeff * *

### Detailed settings

# Integrator

fix integrate all nve/sphere
# Gravity

fix grav all gravity 9.81 vector 0.0 0.0 -1.0
# Time step

timestep 0.00001
#Thermo dynamic output settings

thermo style custom step atoms ke cpu
thermo 10000
thermo_modify norm no lost ignore
# Check time step and initialize dump file

fix ctg all check/timestep/gran 1 0.01 0.01
run 1
unfix ctg
# Create imaging information

dump dumpstl all stl 10000 dump*.stl
dump dmp all custom 10000 dump.1 id type type x y z ix iy iz vx vy vz
& fx fy fz omegax omegay omegaz radius

### Execution and further settings

#Fill the hopper
run 150000 upto
# Pull the stopper
unfix stopper

run 300000