24小时热门版块排行榜    

查看: 1559  |  回复: 3

little_lion

木虫 (小有名气)

[求助] trjconv命令等间隔输出pdb

求教如何用trjconv命令等间隔地输出pdb,比如输出,0ps,2ps,4ps,6ps的构型呢?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
zh1987hs: 金币+3, 谢谢 2013-03-03 09:30:50
仔细读手册:
仔细看参数:
体会:
-b
-e
-dt
2楼2013-02-27 16:57:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

469187892

银虫 (小有名气)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
zh1987hs: 金币+3, 谢谢 2013-03-03 09:30:56
假设想输出2ps的构型:trjconv -f  .trr/xtc -s .tpr -o -dump 2 同理可输出其他的结构,要是一些命令不懂的话,就可以输入trjconv -h
加油!
3楼2013-03-02 15:06:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

旺仔小馒头01

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
little_lion: 金币+15, ★★★★★最佳答案, 非常感谢你的帮助,祝你科研顺利 2013-03-09 14:16:16
使用-sep选项就可以
可以根据你的.mdp文件中设置的输出步长,例如:
nstxout                = 1000                ; save coordinates every 2 ps
那你就可以直接使用
trjconv -s topol.tpr -f traj.xtc -o conf.pdb -sep
这样就可以直接每隔2ps输出一个pdb,即输出,0ps, 2ps, 4ps, 6ps, ..., 构型,分别自动命名为conf0.pdb,conf1.pdb,conf2.pdb,conf3.pdb,...。
如果.mdp文件中设置的输出步长小,如为100
就可以再加上-skip 9选项,即:
trjconv -s topol.tpr -f traj.xtc -skip 9 -o conf.pdb -sep
就可以达到相同的效果。

希望能帮到你~~~
我也来这学习啦,希望大家多多指教。
4楼2013-03-09 11:22:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 little_lion 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 【长期招聘】复旦大学 AI 制药实验室 博士后招聘 +3 第三种态度 2026-10-28 9/450 2026-11-02 13:06 by Grace_PepLib
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 YXM1FDCSE8Bn 2026-10-31 6/300 2026-11-02 12:51 by n0P4fnsXKbBK
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 YXM1FDCSE8Bn 2026-10-31 5/250 2026-11-02 09:42 by dcc2yOIRGGTT
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-02 09:29 by dcc2yOIRGGTT
[基金申请] 各位基础老师们,看完中央《从严治党》的决定后,对学术界更加有信心了吧 +3 iaeyuan 2026-11-01 3/150 2026-11-02 07:49 by kudofaye
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-02 06:13 by xzcv0Tehd8lQ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-02 06:07 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-02 06:07 by xzcv0Tehd8lQ
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 xvJS1BYDxvBS 2026-10-31 10/500 2026-11-02 06:00 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 xvJS1BYDxvBS 2026-10-31 11/550 2026-11-02 06:00 by xzcv0Tehd8lQ
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 05:54 by xzcv0Tehd8lQ
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-02 05:41 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-02 02:34 by xzcv0Tehd8lQ
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 22:24 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 21:35 by xzcv0Tehd8lQ
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 19:50 by xzcv0Tehd8lQ
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 X5bHxOULxFqb 2026-10-30 22/1100 2026-11-01 16:35 by BCdvBAEpt3zm
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 12:40 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
信息提示
请填处理意见