dump forcechain

Hello, everyone,
I have been bothering by a small problem for a long time. That is the local/gran/vtk. I have run the case by CFDEM successfully. The command "dump dmp2 all local/gran/vtk 100000 post/forcechain/forcechain*.vtk cpg1 " can output what I need. However, when run the case in LIGGGHTS, it doesn't work any more. The output is ERROR: Invalid dump style (/build/liggghts-QCz4in/liggghts-3.3.1+repack1/src/output.cpp:591). I don‘t know why this happens. I found that the screen output that LIGGGHTS version is 3.3.1. What I installed is CFDEM 3.8, that is, the LIGGGHTS is 3.8. Why does it say the LIGGGHTS version is 3.3.1? Is this the problem? Could anyone help me? Thank you very much.
I have upload the in.liggghts_run file. The following is what I get after run the case in LIGGGHTS.

lcraul@ubuntu:~/Mywork/testcase$ mpirun -np 4 liggghts < in.liggghts_run
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)
log ../log.liggghts
thermo_log post/thermo.txt

### Initialization

# Preliminaries
units si
atom_style granular
atom_modify map array
boundary f f f
newton off
communicate single vel yes
processors 4 1 1

# Declare domain
region reg block -0.004 0.107 -0.001 0.026 -0.001 0.002 units box
create_box 2 reg
Created orthogonal box = (-0.004 -0.001 -0.001) to (0.107 0.026 0.002)
4 by 1 by 1 MPI processor grid

# Neighbor listing
neighbor 0.00033 bin
neigh_modify delay 0

### Setup

# Material properties required for new pair styles
fix m1 all property/global youngsModulus peratomtype 1e9 1e9
fix m2 all property/global poissonsRatio peratomtype 0.3 0.3
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.4 0.6 0.6 0.8

# partciles insertion
region bc1 block -0.0033 0 0.0 0.025 0.0 0.001 units box
fix pts1 all particletemplate/sphere 15485863 atom_type 1 density constant 1360 radius constant 0.000165
fix pdd1 all particledistribution/discrete/numberbased 49979687 1 pts1 1
Fix particledistribution/discrete (id pdd1): distribution based on number%:
pts1: d=3.300000e-04 (max. bounding sphere) number%=100.000000%
fix ins2 all insert/rate/region seed 67867967 distributiontemplate pdd1 mass 0.00225 massrate 0.00675 insert_every 33000 overlapcheck yes all_in yes vel constant 1 0 0 region bc1 ntry_mc 10000

# Import mesh from cad:
fix cad1 all mesh/surface file Wedgefracture.stl type 2

Reading STL file 'Wedgefracture.stl'

# Use the imported mesh as granular wall
fix geometry all wall/gran model hertz tangential history mesh n_meshes 1 meshes cad1

# Define the physics
pair_style gran model hertz tangential history #Hertzian without cohesion
pair_coeff * *

# cfd coupling
# fix cfd all couple/cfd couple_every 100 mpi
# fix cfd2 all couple/cfd/force/implicit

### Detailed settings

# Integrator
fix integrate all nve/sphere

# Gravity
fix gravi all gravity 9.81 vector 0.0 -1.0 0.0

# Time step (keep < 20% T_Rayleigh)
timestep 0.0000001

# Thermodynamic output settings
compute 1 all erotate/sphere
compute cpg1 all pair/gran/local pos vel force force_normal force_tangential delta
thermo_style custom step atoms ke c_1 vol
thermo 1000
thermo_modify lost ignore norm no
compute_modify thermo_temp dynamic yes

# Check time step and initialize dump file
fix ctg all check/timestep/gran 10000 0.2 0.2
run 1
Setting up run ...
Import and parallelization of mesh cad1 containing 6720 triangle(s) successful
INFO: Resetting random generator for region bc1
INFO: Particle insertion ins2: 870.439042 particles every 33000 steps - particle rate 263769.406724 (mass rate 0.006750)
87923 particles (mass 0.002250) within 3333000 steps
Memory usage per processor = 7.05982 Mbytes
Step Atoms KinEng 1 Volume
0 0 -0 0 8.991e-06
insertion: proc 0 at 0 %
insertion: proc 0 at 10 %
insertion: proc 0 at 20 %
insertion: proc 0 at 30 %
insertion: proc 0 at 40 %
insertion: proc 0 at 50 %
insertion: proc 0 at 60 %
insertion: proc 0 at 70 %
insertion: proc 0 at 80 %
insertion: proc 0 at 90 %
INFO: Particle insertion ins2: inserted 870 particle templates (mass 0.000022) at step 1
- a total of 870 particle templates (mass 0.000022) inserted so far.
1 870 1.1131882e-05 0 8.991e-06
Loop time of 0.0981231 on 4 procs for 1 steps with 870 atoms

Pair time (%) = 4.37498e-05 (0.0445867)
Neigh time (%) = 0.00629884 (6.41933)
Comm time (%) = 0.00139028 (1.41687)
Outpt time (%) = 0.0332977 (33.9346)
Other time (%) = 0.0570925 (58.1846)

Nlocal: 217.5 ave 870 max 0 min
Histogram: 3 0 0 0 0 0 0 0 0 1
Nghost: 0 ave 0 max 0 min
Histogram: 4 0 0 0 0 0 0 0 0 0
Neighs: 1111.25 ave 4445 max 0 min
Histogram: 3 0 0 0 0 0 0 0 0 1

Total # of neighbors = 4445
Ave neighs/atom = 5.1092
Neighbor list builds = 1
Dangerous builds = 0
unfix ctg

# Create imaging information
dump stl1 all mesh/stl 1000000 post/dump*.stl
dump dmp all custom 10000 post/dump.liggghts_run id type x y z vx vy vz fx fy fz omegax omegay omegaz radius
dump dmp2 all local/gran/vtk 100000 post/forcechain/forcechain*.vtk cpg1
ERROR: Invalid dump style (/build/liggghts-QCz4in/liggghts-3.3.1+repack1/src/output.cpp:591)
-------------------------------------------------------
Primary job terminated normally, but 1 process returned
a non-zero exit code.. Per user-direction, the job has been aborted.
-------------------------------------------------------
--------------------------------------------------------------------------
mpirun detected that one or more processes exited with non-zero status, thus causing
the job to be terminated. The first process to do so was:

Process name: [[4500,1],2]
Exit code: 1