计算化学公社

标题: lammps的密度分布如何像rdf和msd一样输出时间的平均值？ [打印本页]

作者
Author: sigusigu 时间: 2023-3-2 12:20
标题: lammps的密度分布如何像rdf和msd一样输出时间的平均值？
我在处理lammps数据时发现密度分布数据是根据规定的时间间隔输出的，大概有几千组，需要取平均数，感觉处理起来很麻烦。请问大家可以像rdf和msd计算时一样把fix ave/time命令融合进来嘛？这样也许就能输出平均值了。compute H2Ozchunk OH2O chunk/atom bin/1d z lower 0.003 nchunk once discard yes units reduced
fix H2Ozaxisd OH2O ave/chunk ${Nevery} ${Nrepeat} ${Nfreq} H2Ozchunk density/mass ave running file H2Oz_density.log

或者大家有直接用代码处理这种平均值的吗？还请不吝赐教

作者
Author: sigusigu 时间: 2023-3-2 14:16
解决了，最后一组就是总的平均值。事实上用了ave running从第二组开始每组都是之前所有组的平均数

欢迎光临计算化学公社 (http://bbs.keinsci.com/)