求助：非平衡冲击模拟冲击波速如何计算？

ddddgc

对长程无序，短程有序的二维体系，以一定的粒子速度，在NVE系综下，采用动量镜法，做了非平衡冲击模拟。

请教下，如何计算冲击波速？

尝试过，输出分块温度，以温度差异判断冲击波面位置。但是因为冲击波后原子温度有波动，所以判断出的波面位置不是很准确，计算出的冲击波速有些不稳定。想请问一下有没有更好的办法？