24小时热门版块排行榜    

查看: 2857  |  回复: 12
【奖励】 本帖被评价10次,作者小虫009增加金币 9.2

小虫009

金虫 (初入文坛)


[资源] 【分享】一个将CASTEP的.cell文件转换为POSCAR的简单脚本

在CASTEP下建立好结构模型,保存文件,找到*.cell文件(存有结构信息)
拷贝到LINUX下,输入perl cell2poscar.pl .cell, 生成POSCAR文件;
或者将 cell2poscar.pl放在~/bin下,输入chmod +x cell2poscar.pl转为可执行形式,直接输入cell2poscar.cell .cell即可。

目前只是半成品,生成的POSCAR没有告诉原子数,可以自己手动加上去,不支持Selective模式
不过基本达到转换目的。
哎偶的perl水平太太菜了,希望各位高手多提宝贵意见!



cell2poscar.pl脚本内容:

eval '(exit $?0)' && eval 'exec perl -S $0 ${1+"$@"}' && eval 'exec perl -S $0 $argv:q' if 0;
#;-*- Perl -*-

#  used to convert .cell(Castep) to POSCAR(vasp)
#
#  usage: cell2poscar.pl
#
@args=@ARGV;
@args==1 || die "usage: cell2poscar.pl \n";

open (IN,@args[0]) || die ("Can't open the file $!";
open (OUT,">POSCAR";
print " Please input the first line of your POSCAR\n";

$filename=;

print " Plesse tell me what you want for (Direct) or (Selective dynamic)   \n";

$Direct_or_selectivedynamic="Direct";
chomp($yourchose=;
if ($yourchose=~ m/^(d|D)/)
   {$Direct_or_selectivedynamic="Direct";}
      elsif($yourchose=~m/^(s|S)/)
           {$Direct_or_selectivedynamic="Selective dynamic";}
            else{$Direct_or_selectivedynamic="Direct";}
$i=1;
while(
{
      @header[0]=$filename;
      @header[1]="1\n";
      if (/\%BLOCK LATTICE_CART/)
      { $_="\n";}
        elsif (/\%ENDBLOCK POSITIONS_FRAC/)
           { $_="\n"; last;}
           elsif (/\%ENDBLOCK LATTICE_CART/)
           {$_=" ";}
           elsif (/\%BLOCK POSITIONS_FRAC/)
           {@header[6]="$Direct_or_selectivedynamic\n";  $i++;}
           else
              {  if(/[a-z]/i){ s/[a-z][a-z]//i;@header[$i+1]=$_; $i++;}
                 else{@header[$i+1]=$_; $i++;}
              }
}

print OUT @header;
close OUT;
close IN;

[ Last edited by zxzj05 on 2009-12-4 at 21:42 ]
回复此楼

» 收录本帖的淘帖专辑推荐

VASP资料

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

★★★★★ 五星级,优秀推荐

呵呵,不错!!!
2楼2009-12-04 21:37:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanchen021

铁杆木虫 (正式写手)


★★★★★ 五星级,优秀推荐

呵呵,我发过一个类似的帖子,LZ可以去看一下
http://muchong.com/bbs/viewthread.php?tid=1663483&fpage=1
4楼2009-12-05 00:08:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wjhwoods

铜虫 (正式写手)


t

对不起,本贴要求 金币高于 500 个 才可浏览,请返回
9楼2010-02-05 04:52:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cjinf0614

金虫 (小有名气)


★★★★★ 五星级,优秀推荐

顶一个
12楼2012-05-03 14:20:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

★★★★★ 五星级,优秀推荐

★★★★★ 五星级
13楼2014-10-20 09:58:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
zxzj053楼
2009-12-04 21:42   回复  
 
2009-12-05 11:19   回复  
 谢谢分享
akakcolin6楼
2009-12-05 13:12   回复  
2009-12-06 09:34   回复  
 
0303deer8楼
2009-12-07 10:16   回复  
 谢谢分享!
2011-04-01 11:09   回复  
五星好评  
huyuxiav11楼
2011-04-15 11:28   回复  
五星好评  谢谢分享!
相关版块跳转 我要订阅楼主 小虫009 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿武汉理工材料工程专硕调剂 +7 Doleres 2026-03-19 7/350 2026-03-20 10:39 by EBSD
[考研] 329求调剂 +6 想上学吖吖 2026-03-19 6/300 2026-03-20 09:00 by 每天只摆一小会
[考研] 材料专硕英一数二306 +6 z1z2z3879 2026-03-18 6/300 2026-03-20 08:49 by xingguangj
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 生物学调剂招人!!! +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考博] 东华理工大学化材专业26届硕士博士申请 +8 zlingli 2026-03-13 8/400 2026-03-19 16:32 by 轻松不少随
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +5 枫桥ZL 2026-03-18 7/350 2026-03-19 14:52 by 功夫疯狂
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-19 3/150 2026-03-19 13:36 by houyaoxu
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[考研] 280求调剂 +6 咕噜晓晓 2026-03-18 7/350 2026-03-18 11:25 by 无际的草原
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[基金申请] 被我言中:新模板不强调格式了,假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 278求调剂 +3 Yy7400 2026-03-13 3/150 2026-03-17 08:24 by laoshidan
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
信息提示
请填处理意见