Lammps Shock Package

Hi Lammps users

I have simulated shock wave in Lammps with fix wall/atoms and fix append/atoms.
I took the length of thermostated region of fix append/atoms equal to the length of appended region and my fix append/atoms, fix wall/piston and fix nve commands were:
fix walls all wall/piston zlo vel 0.0 0.0 ${VLammps} temp 1.0 10 589999 $l units box

fix 1 all nve

fix 2 all append/atoms zhi size y freq {FREQ} temp $T 100 579946 $y units lattice

But when I saw the temperature profile it seems that the shock on the border of appended region was splited into two waves: one in direction of original wave and the other on the opposite direction. And when I visualized my results in AtomEye when one region is appended, the previously appended region cross section shrinks and it’s seems that the newly appended region have zero temperature.

Can you help me find out what the problem is?

Any help is appreciated in advance

Mahdi Tavakol

Aidan may be able to suggest something, but
this is the kind of thing you have to debug. I am
guessing the commands are doing exactly
what you are asking them to do.