24小时热门版块排行榜    

查看: 442  |  回复: 0

vallen

银虫 (正式写手)

[求助] Discover中怎样用BTCL编程计算两组分子之间的相互作用能

我现在有一段轨迹,想分析里面的两组分子之间作用能随时间的变化

网上有个例子是这样的
readFile returns the total energy of the ith step
# This loop will read all available MD records in the his file
for {$i = 1} {[set energy [readFile history frame = $i]] != ""} {incr i} {
# We can control which records will be calculated
if { $i < 3000 } {
# Output the step number to stdout
echo $i
# Write the calculation results to the output text file
puts $file "Step = $i $energy"
# Perform the nonbond energy calculation here,
# the subsets' names should be changed for a different system
EnergyAnalysis pv19,ureum -fid $file
# Output a newline
puts $file ""
} else {
break
}
}

这样的写法不知道对不对?
我不知道怎么从我的体系中选择两组分子  即  subset,用BTCL语句是怎么写的?
是不是要用subset define  命令啊?具体怎么写呢?
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 vallen 的主题更新
信息提示
请填处理意见