24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 1468  |  回复: 12
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

thelrgbird

金虫 (小有名气)

[求助] 关于脚本书写求助

在提交作业的时候写了这样一个脚本:
#!/bin/sh
rm WAVECAR
rm fcc_mesh
for i in 2.6 2.65 2.7 2.75 2.8 2.85 2.9 2.95  3.0 3.05
do
cat > POSCAR << !
Pt
1.000
$i/2    $i/2     0.000
$i/2    0.000  $i/2
0.000 $i/2     $i/2
1
d
0.000 0.000 0.000
!
~/vasp
grep F= OS* >> fcc_mesh
done

运行出了不少错误…… 第一个就是不能用$i/2这个方式简单的表示算数关系,第二个是什么原因我也看不明白了…… 好像说我这样写的POSCAR有问题,麻烦各位帮忙看看,应该怎么修改一下可以呢?谢谢啦!
PS:我想建的模型就是简单的立方面心的Pt
回复此楼

» 猜你喜欢

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

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

thelrgbird

金虫 (小有名气)

引用回帖:
: Originally posted by dxcharlary at 2012-02-18 21:24:30:
另外你的grep F= OS* >> fcc_mesh是想往fcc_mesh里面写什么?

这个就是想手动检测晶格参数的意思,OS*其实就是OSZICAR,懒得继续输入了
9楼2012-02-20 09:36:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

贺仪

铁杆木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
uuv2010(金币+1): 多谢提示 2012-02-17 09:43:37
thelrgbird(金币+2): 有帮助 谢谢啦~ 2012-02-20 16:03:32
thelrgbird(金币+2): 有帮助 谢谢帮助 2012-02-20 16:04:58
一,既然不可以运算,那么你可以先算一下呀
for i in 2.6 2.65 2.7 2.75 2.8 2.85 2.9 2.95  3.0 3.05(在这里算一下不就行了?)
do
cat > POSCAR << !
Pt
1.000
$i/2    $i/2     0.000
$i/2    0.000  $i/2
0.000 $i/2     $i/2

二,把你提示的错误贴上来
嘟啦啦嘟啦啦嘟
2楼2012-02-17 09:31:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
WDD880227(金币+1): 感谢交流 2012-02-18 19:45:57
thelrgbird(金币+2): 有帮助 谢谢~ 2012-02-20 16:03:46
thelrgbird(金币+2): 有帮助 谢谢帮助 2012-02-20 16:05:19
另外,可以用这样的方法计算:
HalfI=`echo "scale=5; $i / 2" | bc -l`
然后在晶格常数位置用$HalfI就可以了~
思想重于技巧,内涵重于表象
3楼2012-02-18 18:04:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
liliangfang(金币+2): 多谢指教 2012-02-19 08:45:34
thelrgbird(金币+2): 有帮助 谢谢 2012-02-20 16:03:57
应该是你的第一个错误导致$i/2没被执行,所以可能晶格参数没写POSCAR里面去。所以vasp根本就不认识你的POSCAR。
4楼2012-02-18 21:22:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 335求调剂 +19 想上岸呀!! 2026-04-12 21/1050 2026-04-14 16:23 by Art1977
[考研] 求调剂学校 +9 不会吃肉 2026-04-13 9/450 2026-04-14 16:12 by bljnqdcc
[考研] 085400电子信息类(川大控制工程)求调剂可跨专业 求老师联系 +4 626776879 2026-04-08 4/200 2026-04-14 14:49 by htcscq
[考研] 284求调剂 +17 让我上岸吧阿西 2026-04-09 17/850 2026-04-14 14:44 by 不我拉绿卡
[考研] 085600材料与化工349分求调剂 +16 李木子啊哈哈 2026-04-12 17/850 2026-04-14 09:11 by fenglj492
[考研] 调剂求收留 +32 果然有我 2026-04-10 33/1650 2026-04-14 08:49 by 木木mumu~
[考研] 考研英一数一338分 +9 长江大学东校区 2026-04-13 10/500 2026-04-14 00:41 by 王珺璞
[考研] 307求调剂 +10 tzq94092 2026-04-10 10/500 2026-04-12 08:18 by wise999
[考研] 一志愿厦大0856,306求调剂 +15 Bblinging 2026-04-11 15/750 2026-04-11 22:53 by 314126402
[考研] 材料工程日语考生求调剂 +7 0856?调剂 2026-04-10 7/350 2026-04-11 21:33 by 蓝云思雨
[考研] 调剂 +4 电气300求调剂不 2026-04-08 7/350 2026-04-11 10:44 by 紫曦紫棋
[考研] 农业管理302分求调剂 +3 xuening1 2026-04-10 3/150 2026-04-11 10:18 by zhq0425
[考研] 337求调剂 +4 研s. 2026-04-10 4/200 2026-04-11 08:57 by zhq0425
[考研] 本科211 工科085400 280分求调剂 可跨专业 +11 LZH(等待调剂中 2026-04-10 11/550 2026-04-11 08:39 by zhq0425
[考研] 263能源动力专硕求调剂 +3 加大号饭盒袋 2026-04-10 3/150 2026-04-10 22:23 by 286640313
[考研] 初试261 +3 Asht少 2026-04-10 6/300 2026-04-10 16:38 by Asht少
[考研] 266求调剂 +29 阳阳哇塞 2026-04-07 29/1450 2026-04-10 16:20 by 高维春
[考博] 博士自荐 +7 可可小胖 2026-04-08 7/350 2026-04-10 08:28 by kimhero
[考研] 材料专硕调剂 +16 哈哈哈吼吼吼哈 2026-04-07 17/850 2026-04-09 21:16 by wutongshun
[考研] 求调剂 +13 柒luck 2026-04-07 13/650 2026-04-08 22:46 by 猪会飞
信息提示
请填处理意见