这段命令是用来获取分子中与溶剂中的锂离子(resname LI)相邻且距离小于3.5的氢原子(noh same resid as (resname SOL and within 3.5 of resname LI)),并将每个氢原子的数量写入一个名为"ry.txt"的文件中。该命令会遍历分子中的每一帧(numFrame),在每一帧中更新氢原子的选择($hnumber update),并将每一帧中符合条件的氢原子数量写入文件中。最后,关闭文件。

set numFrame molinfo top get numframesset out open rytxt w for set i 0 $i $numFrame incr i set hnumber atomselect top noh same resid as resname SOL and within 35 of resname LI$hnumber frame $i$hnumbe

原文地址: https://www.cveoy.top/t/topic/ibQ4 著作权归作者所有。请勿转载和采集!

免费AI点我,无需注册和登录