| 24小時(shí)熱門(mén)版塊排行榜 |
| 10 | 1/1 | 返回列表 |
| 查看: 3394 | 回復(fù): 9 | |||
[交流]
求教vasp中的PBS腳本如何實(shí)現(xiàn)循環(huán)? 已有4人參與
|
|
比如我想實(shí)現(xiàn)以下的循環(huán),用PBS腳本提交任務(wù)的話,應(yīng)該怎么編寫(xiě)? 新人,沒(méi)金幣,不好意思了 ![]() #!/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 |
|
本帖內(nèi)容被屏蔽 |
鐵蟲(chóng) (正式寫(xiě)手)
|
大致就是你寫(xiě)的這樣,但是要加vasp運(yùn)行的命令 手機(jī)回復(fù),不方便添加我自己的腳本,如果有需要的話回頭給你貼一個(gè) 發(fā)自小木蟲(chóng)IOS客戶端 |
鐵蟲(chóng) (正式寫(xiě)手)
|
PBS腳本不是shell腳本么? 這是我們用來(lái)測(cè)試K點(diǎn)的腳本,題主可以參考一下 #!/bin/bash #PBS -N kpoints #PBS -j oe #PBS -l nodes=2:ppn=20 #PBS -l walltime=800:00:00 #PBS -q lhj nprocs=`wc -l < $PBS_NODEFILE` cd $PBS_O_WORKDIR for i in 3 5 7 9 11 13 15 do cat >KPOINTS <<! auto generate 0 M $i $i 1 0 0 0 ! /public/intel/impi/4.1.0.024/intel64/bin/mpirun -np $nprocs ~/bin/vasp.5.3 E=`grep "entropy =" OUTCAR |tail -1| awk '{printf "%12.6f \n",$5}'` echo $i $E >>newkpoints rm WAVECAR CHG* done |
鐵蟲(chóng) (正式寫(xiě)手)
鐵蟲(chóng) (正式寫(xiě)手)
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|