24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2872  |  回复: 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 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 药学专硕调剂 +5 ? 一路生?花? 2026-04-10 6/300 2026-04-10 11:54 by pengliang8036
[考研] 296求调剂 +6 汪!?! 2026-04-08 6/300 2026-04-10 11:02 by mattzhming
[考研] 332求调剂 +8 蕉蕉123 2026-04-10 8/400 2026-04-10 10:49 by 紫翼精灵
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +20 慕绝cc 2026-04-09 24/1200 2026-04-10 10:22 by xujun0624
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +16 努力奋斗112 2026-04-06 18/900 2026-04-10 10:06 by may_新宇
[考研] 材料专硕(0856) 339分求调剂 +11 哈哈哈鹅哈哈哈 2026-04-04 11/550 2026-04-10 09:37 by 690616278
[考研] 087100初试311求调剂 +3 任雅琴 2026-04-09 3/150 2026-04-09 22:42 by lbsjt
[考研] 调剂 +19 2261744733 2026-04-08 19/950 2026-04-09 19:11 by vgtyfty
[考研] 材料299专硕求调剂 +10 +21 2026-04-09 10/500 2026-04-09 17:34 by 1753564080
[考研] 322求调剂,08工科 +3 今天是个小号 2026-04-08 3/150 2026-04-09 15:53 by wp06
[考研] 085501机械专硕 302分 不挑专业求调剂 +5 汪某. 2026-04-09 5/250 2026-04-09 15:38 by 蒋皓禹
[考研] 070300化学学硕311分求调剂 +18 梁富贵险中求 2026-04-04 20/1000 2026-04-09 11:18 by 哒哒哒呱呱呱
[考研] 本科南方医科大学 一志愿985 药学学硕284分 求调剂 +3 弱水听文 2026-04-09 3/150 2026-04-09 09:06 by susuqq
[考研] 材料调剂 +14 一样YWY 2026-04-06 14/700 2026-04-08 23:00 by 猪会飞
[考研] 一志愿南京航空航天大学 材料与化工329分求调剂 +11 Mr. Z 2026-04-05 12/600 2026-04-08 16:15 by luoyongfeng
[考研] 生物医药调剂|SCI中科院三区一作+多项科研成果 +8 likangxing 2026-04-07 11/550 2026-04-08 00:02 by lys0704
[考研] 求调剂到材料 +5 程9915 2026-04-06 5/250 2026-04-06 15:21 by yulian1987
[考研] 考研调剂 +5 美丽的youth_ 2026-04-04 6/300 2026-04-06 06:57 by houyaoxu
[考研] 085600调剂 +9 东照照照 2026-04-04 9/450 2026-04-05 13:44 by ujn_zhuj
[考研] 283分求调剂 +7 小聂爱学习 2026-04-03 7/350 2026-04-04 21:51 by hemengdong
信息提示
请填处理意见