24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2870  |  回复: 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 ]
回复此楼
已阅   回复此楼   关注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 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程302分求调剂 +18 zyx上岸! 2026-04-04 18/900 2026-04-10 10:07 by 314126402
[考研] 283求调剂 +20 那个噜子 2026-04-09 20/1000 2026-04-10 09:50 by liuhuiying09
[考研] 求调剂 +13 MAX怅惘 2026-04-09 14/700 2026-04-10 09:14 by Delta2012
[考研] 求调剂希望还是希望在山河四省附近 +3 快乐的小白鸽 2026-04-05 3/150 2026-04-09 17:36 by wp06
[硕博家园] 新一代电子信息294求调剂 不挑学校 +5 Ytyt11 2026-04-09 6/300 2026-04-09 14:40 by Ytyt11
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[考研] 材料考研求调剂总分280 +30 mkjlz1 2026-04-06 35/1750 2026-04-08 21:25 by cyh—315
[考研] 机械工程313分找工科调剂 +3 双一流本科机械 2026-04-08 3/150 2026-04-08 20:41 by 土木硕士招生
[考研] 324求调剂 +17 想上学求调 2026-04-03 17/850 2026-04-08 20:04 by 我减肥1
[考研] 298求调剂 +4 manman511 2026-04-05 4/200 2026-04-08 16:50 by tjzhao
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 22:00 by hemengdong
[考博] 博士申请 +3 IQwQl 2026-04-05 3/150 2026-04-07 20:31 by greychen00
[考研] 生物学363调剂求助 +7 fanzhang6666 2026-04-06 9/450 2026-04-07 17:37 by lijunpoly
[考研] 求调剂 +10 chenxrlkx 2026-04-05 10/500 2026-04-06 11:31 by 猪会飞
[考研] 求调剂到0856材料工程 +3 程9915 2026-04-05 3/150 2026-04-05 18:15 by 蓝云思雨
[考研] 一志愿西北农林畜牧专硕336分求调剂 +3 5ourr 2026-04-03 3/150 2026-04-05 10:40 by JOKER0401
[考研] 求调剂 +4 晟功? 2026-04-03 4/200 2026-04-04 21:58 by hemengdong
[考研] 085400电子信息319求调剂(接受跨专业调剂) +5 星星不眨眼喽 2026-04-03 6/300 2026-04-04 21:50 by hemengdong
[考研] 085701求调剂 +7 龚禹铭 2026-04-04 8/400 2026-04-04 13:49 by 小小树2024
[考研] 268求调剂 +8 你好tg 2026-04-03 9/450 2026-04-04 05:08 by gswylq
信息提示
请填处理意见