版塊導航: 正在加載中...

調(diào)劑小程序

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

>論壇更新日志 (2294)
>考研 (134)
>導師招生 (133)
>蟲友互識 (83)
>基金申請 (50)
>休閑灌水 (32)
>碩博家園 (18)
>博后之家 (17)
>文獻求助 (14)
>微米和納米 (13)
>教師之家 (13)
>論文投稿 (12)
>考博 (10)
>找工作 (10)
>公派出國 (8)
>論文道賀祈福 (7)

返回列表

thelrgbird

金蟲 (小有名氣)

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

[求助] 關于腳本書寫求助

在提交作業(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ù)關系，第二個是什么原因我也看不明白了…… 好像說我這樣寫的POSCAR有問題，麻煩各位幫忙看看，應該怎么修改一下可以呢？謝謝啦！
PS：我想建的模型就是簡單的立方面心的Pt

回復此樓

» 猜你喜歡

一志愿華中科技大學071000，求調(diào)劑已經(jīng)有4人回復
生物學一志愿985，分數(shù)349求調(diào)劑已經(jīng)有4人回復
生物學調(diào)劑已經(jīng)有4人回復
求調(diào)劑院校信息已經(jīng)有4人回復
085600材料與化工306 已經(jīng)有4人回復
286求調(diào)劑已經(jīng)有10人回復
328求調(diào)劑，英語六級551，有科研經(jīng)歷已經(jīng)有9人回復
一志愿北京化工大學070300 學碩336求調(diào)劑已經(jīng)有4人回復
286分人工智能專業(yè)請求調(diào)劑愿意跨考！已經(jīng)有8人回復
資源與環(huán)境調(diào)劑申請(333分) 已經(jīng)有5人回復

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

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

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

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

賀儀

鐵桿木蟲 (著名寫手)

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

【答案】應助回帖

感謝參與，應助指數(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

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

贊一下

回復此樓

嘟啦啦嘟啦啦嘟

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

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

souledge

專家顧問 (著名寫手)

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

【答案】應助回帖

感謝參與，應助指數(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就可以了～

贊一下

回復此樓

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

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

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

dxcharlary

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

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

【答案】應助回帖

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

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

贊一下

回復此樓

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

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

dxcharlary

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

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

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

贊一下

回復此樓

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

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

thelrgbird

金蟲 (小有名氣)

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

引用回帖:

樓: 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不認……不能繼續(xù)計算……

贊一下

回復此樓

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

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

thelrgbird

金蟲 (小有名氣)

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

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  {一個運算公式} {一個運算公式}

贊一下

回復此樓

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

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

thelrgbird

金蟲 (小有名氣)

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

引用回帖:

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

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

贊一下

回復此樓

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

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

thelrgbird

金蟲 (小有名氣)

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

引用回帖:

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

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

贊一下

回復此樓

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

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

dxcharlary

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

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

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吧。

贊一下

回復此樓

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

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

相關版塊跳轉(zhuǎn) 我要訂閱樓主 thelrgbird 的主題更新

返回列表

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

[考研] 生物學調(diào)劑 +4	Surekei 2026-03-21	4/200	2026-03-22 07:18 by ilovexiaobin
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3	邱gl 2026-03-20	7/350	2026-03-21 19:05 by 15709483992
[考研] 297求調(diào)劑 +11	戲精丹丹丹 2026-03-17	12/600	2026-03-21 17:47 by ColorlessPI
[考研] 299求調(diào)劑 +5	shxchem 2026-03-20	7/350	2026-03-21 17:09 by ColorlessPI
[考研] 0856材料專碩353求調(diào)劑 +3	NIFFFfff 2026-03-20	3/150	2026-03-21 10:23 by luoyongfeng
[考研] 南昌大學材料專碩311分求調(diào)劑 +6	77chaselx 2026-03-20	6/300	2026-03-21 07:24 by JourneyLucky
[考研] 265求調(diào)劑 +3	Jack?k?y 2026-03-17	3/150	2026-03-21 03:17 by JourneyLucky
[考研] 303求調(diào)劑 +5	睿08 2026-03-17	7/350	2026-03-21 03:11 by JourneyLucky
[考研] 299求調(diào)劑 +6	△小透明* 2026-03-17	6/300	2026-03-21 02:42 by JourneyLucky
[考研] 265求調(diào)劑 +9	梁梁校校 2026-03-17	9/450	2026-03-21 02:17 by JourneyLucky
[考研] 南京大學化學376求調(diào)劑 +3	hisfailed 2026-03-19	6/300	2026-03-20 23:43 by hisfailed
[考研] 294求調(diào)劑材料與化工專碩 +15	陌の森林 2026-03-18	15/750	2026-03-20 23:28 by JourneyLucky
[考研] 0817 化學工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +22	rare12345 2026-03-18	22/1100	2026-03-20 20:39 by zhukairuo
[考研] 材料與化工專碩調(diào)劑 +7	heming3743 2026-03-16	7/350	2026-03-20 19:31 by zhukairuo
[考研] 0703化學336分求調(diào)劑 +6	zbzihdhd 2026-03-15	7/350	2026-03-18 09:53 by zhukairuo
[考研] 334求調(diào)劑 +3	志存高遠意在機?/a> 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 283求調(diào)劑 +3	聽風就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 0856求調(diào)劑 +3	劉夢微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu