请选择 进入手机版 | 继续访问电脑版

分子模拟论坛 Molecular Simulation Forums

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 460|回复: 4

轨迹档转成PDB檔的SCRIPT

[复制链接]

30

主题

146

帖子

328

积分

中级会员

Rank: 3Rank: 3

积分
328
发表于 2008-7-16 15:25:00 | 显示全部楼层 |阅读模式
[ol]
  • set FILE 400
  • set N 1
  • set START 100 !starting frame no.
  • set END 200000 !ending frame no.
  • set INCR 100 !change in frame no. from oneframe to the next
  • calc LOOP= @END-@START
  • read rtf card name top_all27_prot_na.rtf
  • read para card namepar_all27_prot_na_gbsw.prm
  • ! Get starting psf and coordinates
  • open read formatted unit 3 namestep1_pdbreader.psf
  • read psf card unit 3
  • !open unit 4 card read name n1.crd
  • !read coor card unit 4
  • open read unit 4 card namestep1_pdbreader.pdb
  • read coor unit 4 pdb
  • resid
  • open unit 67 write form name@FILE-@N-cor.pdb !output pdb trajectory
  • open unit 31 read unform name run.dcd!input cor trajectory
  • traj iread 31 nrea 1 begin @START stop @ENDskip @INCR
  • set tr @INCR
  • LABEL ILP
  • traj read
  • write coor pdb unit 67
  • incr tr by @INCR
  • if tr .le. @LOOP goto ILP
  • stop[/ol]复制代码
  • 回复

    使用道具 举报

    24

    主题

    93

    帖子

    216

    积分

    中级会员

    Rank: 3Rank: 3

    积分
    216
    发表于 2008-9-4 14:20:00 | 显示全部楼层
    thank you !
    回复

    使用道具 举报

    8

    主题

    50

    帖子

    114

    积分

    注册会员

    Rank: 2

    积分
    114
    发表于 2009-3-11 10:48:00 | 显示全部楼层
    多谢共享!

    回复

    使用道具 举报

    18

    主题

    216

    帖子

    456

    积分

    中级会员

    Rank: 3Rank: 3

    积分
    456
    发表于 2009-3-19 22:43:00 | 显示全部楼层
    thanks.....
    回复

    使用道具 举报

    4

    主题

    21

    帖子

    48

    积分

    新手上路

    Rank: 1

    积分
    48
    发表于 2009-4-14 07:29:00 | 显示全部楼层
    很有用的script,thanks
    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 立即注册

    本版积分规则

    Archiver|手机版|小黑屋|分子模拟论坛

    GMT+8, 2020-8-4 22:28 , Processed in 0.068149 second(s), 26 queries .

    Powered by Discuz! X3.4

    © 2001-2017 Comsenz Inc.

    快速回复 返回顶部 返回列表