thelrgbird

金蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 1250.1
散金: 25
帖子: 150
在線: 89.8小時
蟲號: 597343
注冊: 2008-09-09
性別: GG
專業(yè): 理論和計算化學(xué)

[求助] 關(guān)于腳本書寫求助

在提交作業(yè)的時候?qū)懥诉@樣一個腳本：
#!/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這個方式簡單的表示算數(shù)關(guān)系，第二個是什么原因我也看不明白了…… 好像說我這樣寫的POSCAR有問題，麻煩各位幫忙看看，應(yīng)該怎么修改一下可以呢？謝謝啦！
PS：我想建的模型就是簡單的立方面心的Pt

回復(fù)此樓

» 猜你喜歡

309求調(diào)劑已經(jīng)有8人回復(fù)
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷已經(jīng)有8人回復(fù)
求調(diào)劑已經(jīng)有5人回復(fù)
求化學(xué)調(diào)劑已經(jīng)有5人回復(fù)
085600，材料與化工321分求調(diào)劑已經(jīng)有10人回復(fù)
0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分已經(jīng)有5人回復(fù)
一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑已經(jīng)有3人回復(fù)
085600，專業(yè)課化工原理，320分求調(diào)劑已經(jīng)有4人回復(fù)
材料與化工272求調(diào)劑已經(jīng)有16人回復(fù)
290求調(diào)劑已經(jīng)有3人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

求助一個簡單的腳本已經(jīng)有4人回復(fù)
求助高斯腳本同時提交兩個任務(wù)，一個結(jié)束，另一個也結(jié)束了已經(jīng)有8人回復(fù)
shell腳本語言代碼求助已經(jīng)有10人回復(fù)
【求助】求助差分電荷密度的腳本！已經(jīng)有45人回復(fù)
【求助成功】vasp測試K點收斂性的腳本問題已經(jīng)有20人回復(fù)
【求助】求將vasp的輸出文件轉(zhuǎn)化成xyz坐標(biāo)的腳本已經(jīng)有11人回復(fù)
【求助】comsol 3.5a 如何寫腳本文件？已經(jīng)有6人回復(fù)

1樓 2012-02-16 13:50:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

賀儀

鐵桿木蟲 (著名寫手)

1ST強(qiáng)帖: 5
應(yīng)助: 93 (初中生)
貴賓: 0.636
金幣: 5937
散金: 1063
紅花: 38
沙發(fā): 2
帖子: 2499
在線: 588.7小時
蟲號: 995125
注冊: 2010-04-12
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +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

二，把你提示的錯誤貼上來

贊一下

回復(fù)此樓

嘟啦啦嘟啦啦嘟

2樓2012-02-17 09:31:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

souledge

專家顧問 (著名寫手)

專家經(jīng)驗: +174
1ST強(qiáng)帖: 19
應(yīng)助: 303 (大學(xué)生)
金幣: 11306.4
散金: 1123
紅花: 108
帖子: 2814
在線: 916.1小時
蟲號: 941486
注冊: 2010-01-12
專業(yè): 結(jié)構(gòu)陶瓷
管轄: 第一性原理

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +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`
然后在晶格常數(shù)位置用$HalfI就可以了～

贊一下

回復(fù)此樓

思想重于技巧，內(nèi)涵重于表象

3樓2012-02-18 18:04:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dxcharlary

專家顧問 (職業(yè)作家)

專家經(jīng)驗: +83
1ST強(qiáng)帖: 2
應(yīng)助: 251 (大學(xué)生)
金幣: 10421.6
散金: 155
紅花: 65
帖子: 3137
在線: 757.6小時
蟲號: 936718
注冊: 2010-01-04
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
管轄: 第一性原理

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1
liliangfang(金幣+2): 多謝指教 2012-02-19 08:45:34
thelrgbird(金幣+2): ★有幫助謝謝 2012-02-20 16:03:57

應(yīng)該是你的第一個錯誤導(dǎo)致$i/2沒被執(zhí)行，所以可能晶格參數(shù)沒寫POSCAR里面去。所以vasp根本就不認(rèn)識你的POSCAR。

贊一下

回復(fù)此樓

4樓2012-02-18 21:22:48

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dxcharlary

專家顧問 (職業(yè)作家)

專家經(jīng)驗: +83
1ST強(qiáng)帖: 2
應(yīng)助: 251 (大學(xué)生)
金幣: 10421.6
散金: 155
紅花: 65
帖子: 3137
在線: 757.6小時
蟲號: 936718
注冊: 2010-01-04
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
管轄: 第一性原理

另外你的grep F= OS* >> fcc_mesh是想往fcc_mesh里面寫什么？

贊一下

回復(fù)此樓

5樓2012-02-18 21:24:30

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

thelrgbird

金蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 1250.1
散金: 25
帖子: 150
在線: 89.8小時
蟲號: 597343
注冊: 2008-09-09
性別: GG
專業(yè): 理論和計算化學(xué)

引用回帖:

樓: Originally posted by 賀儀 at 2012-02-17 09:31:08:
一，既然不可以運算，那么你可以先算一下呀
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 ...

這樣還是不行的，$i/2在這里出現(xiàn)的話，最后反映在POSCAR上就是2.6/2這樣的東西，貌似VASP不認(rèn)……不能繼續(xù)計算……

贊一下

回復(fù)此樓

6樓2012-02-20 09:29:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

thelrgbird

金蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 1250.1
散金: 25
帖子: 150
在線: 89.8小時
蟲號: 597343
注冊: 2008-09-09
性別: GG
專業(yè): 理論和計算化學(xué)

WDD880227(金幣+1): 感謝交流 2012-02-20 09:39:12

引用回帖:

樓: Originally posted by souledge at 2012-02-18 18:04:45:
另外，可以用這樣的方法計算：
HalfI=`echo "scale=5; $i / 2" | bc -l`
然后在晶格常數(shù)位置用$HalfI就可以了～

我最后也是在木蟲上查的，用的是a= python -c "print $i/2" 但是這樣的話所有的簡單計算都會要重新命名一個新的變量，覺得可麻煩了…… 不知道能不能直接用某個命令讓其直接算出呢？例如：
cat > POSCAR << !
Pt
1.000
{一個運算公式} {一個運算公式}    0.000
{一個運算公式} 0.000  {一個運算公式}
0.000  {一個運算公式} {一個運算公式}

贊一下

回復(fù)此樓

7樓2012-02-20 09:33:18

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

thelrgbird

金蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 1250.1
散金: 25
帖子: 150
在線: 89.8小時
蟲號: 597343
注冊: 2008-09-09
性別: GG
專業(yè): 理論和計算化學(xué)

引用回帖:

樓: Originally posted by dxcharlary at 2012-02-18 21:22:48:
應(yīng)該是你的第一個錯誤導(dǎo)致$i/2沒被執(zhí)行，所以可能晶格參數(shù)沒寫POSCAR里面去。所以vasp根本就不認(rèn)識你的POSCAR。

您說的第一個錯誤是哪里？

贊一下

回復(fù)此樓

8樓2012-02-20 09:34:52

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

thelrgbird

金蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 1250.1
散金: 25
帖子: 150
在線: 89.8小時
蟲號: 597343
注冊: 2008-09-09
性別: GG
專業(yè): 理論和計算化學(xué)

引用回帖:

樓: Originally posted by dxcharlary at 2012-02-18 21:24:30:
另外你的grep F= OS* >> fcc_mesh是想往fcc_mesh里面寫什么？

這個就是想手動檢測晶格參數(shù)的意思，OS*其實就是OSZICAR，懶得繼續(xù)輸入了

贊一下

回復(fù)此樓

9樓2012-02-20 09:36:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dxcharlary

專家顧問 (職業(yè)作家)

專家經(jīng)驗: +83
1ST強(qiáng)帖: 2
應(yīng)助: 251 (大學(xué)生)
金幣: 10421.6
散金: 155
紅花: 65
帖子: 3137
在線: 757.6小時
蟲號: 936718
注冊: 2010-01-04
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
管轄: 第一性原理

WDD880227(金幣+1): 感謝提示 2012-02-20 15:09:02

引用回帖:

7樓: Originally posted by thelrgbird at 2012-02-20 09:33:18:
我最后也是在木蟲上查的，用的是a= python -c "print $i/2" 但是這樣的話所有的簡單計算都會要重新命名一個新的變量，覺得可麻煩了…… 不知道能不能直接用某個命令讓其直接算出呢？例如：
cat > ...

我記得shell根本就不支持浮點運算。用python或者是matlab吧。

贊一下

回復(fù)此樓

10樓2012-02-20 09:51:02

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 085600，材料與化工321分求調(diào)劑 +10	大饞小子 2026-03-28	10/500	2026-03-29 23:35 by 飛行日記西
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +8	柳溪邊 2026-03-26	8/400	2026-03-29 20:42 by 唐沐兒
[考研] 求調(diào)劑 +4	青春裁為三截 2026-03-29	4/200	2026-03-29 15:01 by 唐沐兒
[考研] 一志愿雙一流機(jī)械285分求調(diào)劑 +4	幸運的三木 2026-03-29	5/250	2026-03-29 14:49 by Miko19
[考研] 調(diào)劑考研 +3	王杰一 2026-03-29	3/150	2026-03-29 08:09 by fmesaito
[考研] 329求調(diào)劑 +7	星野? 2026-03-26	7/350	2026-03-29 06:43 by 544594351
[考研] 一志愿211 初試270分求調(diào)劑 +8	谷雨上岸 2026-03-23	9/450	2026-03-29 01:00 by 我是小康
[考研] 數(shù)一英一271專碩（085401）求調(diào)劑，可跨 +7	前行必有光 2026-03-28	8/400	2026-03-28 23:22 by 小木蟲tim
[考研] 295材料工程專碩求調(diào)劑 +7	1428151015 2026-03-27	7/350	2026-03-28 19:58 by S240
[考研] 一志愿北化085600材料專碩275\|有文章專利｜求調(diào)劑 +7	Micky11223 2026-03-25	7/350	2026-03-28 18:34 by 無際的草原
[考研] 一志愿廈門大學(xué)化學(xué)學(xué)碩307求調(diào)劑 +10	y7czhao 2026-03-26	10/500	2026-03-28 14:23 by 唐沐兒
[考研] 347求調(diào)劑 +3	山頂見α 2026-03-25	3/150	2026-03-28 14:13 by 唐沐兒
[考研] 0703化學(xué)求調(diào)劑，各位老師看看我�。。� +5	祁祺祺 2026-03-25	5/250	2026-03-27 21:44 by 東方豬豬
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4	哇呼哼呼哼 2026-03-26	4/200	2026-03-27 17:02 by zhyzzh
[考研] 0856調(diào)劑 +5	求求讓我有書讀?/a> 2026-03-26	6/300	2026-03-27 15:12 by caszguilin
[考研] 一志愿河工大 081700 276求調(diào)劑 +4	地球繞著太陽轉(zhuǎn) 2026-03-23	4/200	2026-03-26 14:27 by zzll406
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +4	gdlf9999 2026-03-24	4/200	2026-03-25 23:01 by boxking200
[考研] 302求調(diào)劑 +4	錦衣衛(wèi)藤椒 2026-03-25	4/200	2026-03-25 16:29 by 功夫瘋狂
[考研] 085404電子信息284分求調(diào)劑 +4	13659058978 2026-03-24	4/200	2026-03-24 12:15 by syl20081243
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4	allen-yin 2026-03-23	6/300	2026-03-23 15:04 by 汪�。�！