下图是我模拟的问题：两侧是固壁，中间是气体原子。我想把体系划分成5个区域（下图中6条白色的线条之间的部分），然后在模拟的过程中实时输出每一部分内所有原子在每个时刻的温度和压力。
但由于模拟的物质是气体，所以每时每刻各个region里面的原子都不同，这导致我需要每算一步就重新定义一遍group，再重新定义 compute temp和compute stress/atom，然后通过控制thermo_style将计算结果输出到log文件中。这样非常麻烦，而且输出的数据还混杂在log文件里，也不方便提取。不知道大家有没有什么简便的办法呢？
我现在的方法是：
label ......定义循环
# --
group ......重新定义group
compute ... temp ...对每个group重新定义计算温度
compute ... stress/atom ...对每个group重新定义计算压力张量
therom 1
therom_style custom ...将计算的温度和压力输出到log文件中
run 1计算1步就输出
# --
jump ......循环结束

我记得以前有人讨论过，有个命令可以实现，固定空间并更新空间里的原子，你可以搜搜看

有啥线索没？。。。这大海捞针的。。。

compute reduce/region

compute temp
compute stress/atom
fix ave/spatial

这个改代码，自己加一个命令，判断原子属于region直接用domain->regions[iregion]->match(x[i][0],x[i][1],x[i][2])。