|
gromacs获取每帧的原子受力与速度时间：2025-04-28 #!/bin/bash #本次运算dt单位为mdp的fs，frame为nstfout或nstvout的值，n为nsteps/frame+1，最大为200，所以小于201 dt=2 frame=500 for ((n=0; n<201; n++)) do # 计算当前时间点（整数运算），因为win下的linux无法使用bc current_time=$(($n*$dt*$frame/1000))
# 执行轨迹分析命令，echo后面不需要加""，否则只是打印这个字符串，并不是进行调用 echo -e system | gmx traj -f nvt1.trr -s nvt1.tpr -cv veloc$current_time.pdb -cf force$current_time.pdb -b $current_time -e $current_time #必须删除all_veloc.xvg all_force.xvg，因为超过99个重复的###会报错，并且不需要这个如果需要在上面的命令后面加上去就好 rm -r all_veloc.xvg all_force.xvg done |