sobereva 发表于 2023-6-10 07:14 尝试按照老师说的原则修改了一下,如下,一直没有成功,想请教老师具体还需要怎么修改,谢谢! proc showtrj {fps1 fps2 space selection output} { set selnow [atomselect top $selection frame $fps1] set selnext [atomselect top $selection frame $fps1] set num [$selnow num] set outfile [open $output w] for {set fps $fps1} {$fps < $fps2} {incr fps $space} { $selnow frame $fps $selnext frame [expr $fps + $space] $selnow update $selnext update puts $outfile "MODEL $fps" for {set i 0} {$i < $num} {incr i 1} { set x [lindex [$selnext get x] $i] set y [lindex [$selnext get y] $i] set z [lindex [$selnext get z] $i] set atomname [lindex [$selnext get name] $i] set resid [lindex [$selnext get resid] $i] set segname [lindex [$selnext get segname] $i] draw line [lindex [$selnow get {x y z}] $i] [lindex [$selnext get {x y z}] $i] draw sphere [lindex [$selnext get {x y z}] $i] radius 0.12 puts $outfile "ATOM $i $atomname $segname $resid $x $y $z 1.00 0.00" } puts $outfile "ENDMDL" puts "Frame $fps done" } close $outfile } |
landian666 发表于 2023-6-9 09:16 修改脚本,把各帧里某个原子的坐标作为导出的xyz文件里的不同原子坐标,则载入xyz文件后直接能看到所有红色原子 |
sobereva 发表于 2023-6-9 06:22 不好意思,没有表达清楚,我的意思是,在利用这个方法绘制出一个原子的轨迹后,如下图,有办法通过修改脚本或者别的方法,将图中红色轨迹以及原始的分子导出一个pdb文件吗?以便之后可以直接导入这个pdb观看整个三维状态。 |
20150609214054_20186.gif (11.91 KB, 下载次数 Times of downloads: 33)
landian666 发表于 2023-6-8 16:44 看不懂你的意思 如果你想用pdb格式记录轨迹,save coordinate的时候选上所有帧,格式选pdb就行了。但pdb格式记录轨迹尺寸很大,很不推荐 |
sobereva 发表于 2022-4-12 12:54 想请教sob老师,我根据这个方法做出了轨迹图,然后需要怎么修改脚本能够将生成的轨迹和原始分子一起导出一个pdb的构型文件那?以便用来估计粒子轨迹的可及运动半径,也方便后期作图。 |
sobereva 发表于 2022-8-30 18:09 非常感谢sob老师的回复! |
tangzi_via 发表于 2022-8-30 16:37 多帧叠加显示,并且coloring method用timestep |
sobereva 发表于 2022-4-12 12:54 sob 老师您好!请问一下显示粒子运动轨迹的这个VMD脚本,能否把时间元素加上,使柱子是一个渐变的颜色。才开始学,希望得到您的帮助 |
显示粒子运动轨迹的VMD脚本 http://sobereva.com/65 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-26 14:52 , Processed in 1.637407 second(s), 25 queries , Gzip On.