lammps的密度分布如何像rdf和msd一样输出时间的平均值？

sigusigu · 发表于 Post on 2023-3-2 12:20:25

我在处理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

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

sigusigu · 发表于 Post on 2023-3-2 14:16:50

解决了，最后一组就是总的平均值。事实上用了ave running从第二组开始每组都是之前所有组的平均数

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Lammps] lammps的密度分布如何像rdf和msd一样输出时间的平均值？

评分 Rate

浏览过的版块