| 6 | 1/1 | 返回列表 |
| 查看: 1081 | 回復(fù): 5 | ||
linda19918銀蟲 (小有名氣)
|
[求助]
第一次用VASP,學(xué)了幾天了,怎么用腳本來找合適的ENCUT,按照說明說總是出錯 已有2人參與
|
|
#!/bin/sh rm WAVECAR for i in 150 200 250 300 350 400 do cat > INCAR <<! SYSTEM = Si-Diamond ENCUT = $i ISTART = 0 ; ICHARG = 2 ISMEAR = -5 PREC = Accurate ! echo "ENCUT = $i eV" ; time vasp E='grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'' echo $i $E >>comment done 然后運行該腳本,該腳本的名字為new,輸入的sh new過后出現(xiàn)一下提示: : No such file or directory rm: cannot remove `WAVECAR\r': No such file or directory 'ew: line 4: syntax error near unexpected token `do 'ew: line 4: `do 我該怎么去修改或者我的操作是錯誤的, PS:INCAR,POSCAR,POTCAR,KPONTS都有 |
|
我的腳本,希望對你有用。 #!/bin/sh APP_NAME=intelg_mid NP=12 NP_PER_NODE=12 RUN="RAW" CURDIR=$PWD VASP=~/bin/vasp.5.2.11-mp source /home-gg/compiler/intel/composer_xe_2011_sp1/bin/compilervars.sh intel64 source /home-gg/compiler/mpi/openmpi-1.4.4-intel.sh cd $CURDIR #start creating .nodelist rm -rf $CURDIR/nodelist >& /dev/null for i in `echo $LSB_HOSTS` do echo "$i" >> $CURDIR/nodelist done nodelist=$(cat $CURDIR/nodelist | uniq | awk '{print $1}' | tr '\n' ',') rm WAVECAR for i in 400 450 500 550 600 do cat > INCAR <<! SYSTEM=SiC ENCUT=$i ISTART=0;ICHARG=2 ISMEAR=-5 PREC=high ! echo "ENCUT= $i eV" mpirun -np $NP -machinefile $PBS_NODEFILE vasp E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5}'` KP=`grep "irreducible" OUTCAR | tail -1 | awk '{printf"%12.4f \n", $5}'` echo $i $E >> comment # # echo -n "start time " > time date >> time mpirun -np $NP -machinefile $CURDIR/nodelist $VASP &> vasp.out echo -n "end time " >> time ; date >> time cat time >> vasp.out done # # exit |

銀蟲 (小有名氣)
銀蟲 (小有名氣)
至尊木蟲 (著名寫手)

| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分,一志愿四川大學(xué),誠心求調(diào)劑 +9 | 吃吃吃才有意義 2026-03-19 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考! +3 | lemonzzn 2026-03-17 | 5/250 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +6 | 墨墨漠 2026-03-20 | 6/300 |
|
|
[考研] 278求調(diào)劑 +6 | 煙火先于春 2026-03-17 | 6/300 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
[考研] 354求調(diào)劑 +5 | Tyoumou 2026-03-18 | 8/400 |
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研] 一志愿中南化學(xué)(0703)總分337求調(diào)劑 +8 | niko- 2026-03-19 | 9/450 |
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 ,六級已過,有科研經(jīng)歷 +13 | 曦熙兮 2026-03-15 | 13/650 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|