| 5 | 1/1 | 返回列表 |
| 查看: 1344 | 回復(fù): 4 | |||
| 本帖產(chǎn)生 2 個(gè) 1ST強(qiáng)帖 ,點(diǎn)擊這里進(jìn)行查看 | |||
dingfengbo新蟲 (小有名氣)
|
[求助]
encut測試是否需要修改要提交的作業(yè)腳本
|
||
|
我想做encut測試這是我的encut腳本: #!/bin/sh rm WAVECAR for a in 400 450 500 550 600 650 700 750 do cat >INCAR < SYSTEM=CGa LPLANE=.TRUE. NPAR=8 LREAL = .TRUE. PREC = Accurate EDIFF =1E-6 EDIFFG = -0.03 IALGO = 48 NELMIN = 4 ISYM = 0 ISPIN = 1 NBANDS = 240 LCHARG = .TRUE. LWAVE = .TRUE. LVTOT = .FALSE. NBLOCK = 1 NSW = 10 IBRION = 2 POTIM=0.3 ISMEAR = 0 SIGMA = 0.05 ENCUT=$a ISTART=0;ICHARG=2 ! echo "ENCUT=$a eV" vasp E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'` echo $a $E >>comment done 測試的時(shí)候我想將之提交,下面是我的作業(yè)腳本: #!/bin/sh #PBS -l nodes=2:ppn=8 #PBS -u XX #PBS -N YY JOBINFOR=$PBS_JOBID MASTERNODE=`hostname` SCRATCHDIR=$PBS_JOBID NCPU=`wc -l < $PBS_NODEFILE` SERVER=$PBS_O_HOST WORKDIR=$PBS_O_WORKDIR MKDIR=/bin/mkdir RSH=/usr/bin/rsh CP=/bin/cp cat $PBS_NODEFILE > $WORKDIR/mpihosts LAUNCH="mpirun --hostfile $WORKDIR/mpihosts -np $NCPU" PROGRAMEXEC="/opt/vasp/vasp" ${RSH} ${MASTERNODE} cd ${WORKDIR} ${MKDIR} ${SCRATCHDIR} cd ${SCRATCHDIR} ${CP} ${WORKDIR}/* . ${LAUNCH} ${PROGRAMEXEC} > out 2>&1 請問我怎樣修改作業(yè)腳本才能進(jìn)行encut測試的計(jì)算呢?謝謝 |
Photochemistry |
金蟲 (著名寫手)
木蟲 (著名寫手)

新蟲 (小有名氣)
木蟲 (正式寫手)
宅男甲
|
個(gè)人意見:PROGRAMEXEC="/opt/vasp/vasp" 改為 PROGRAMEXEC='腳本位置”是不行的,最后一行其實(shí)就是一個(gè)運(yùn)行并行的vasp程序的指令。如果改成腳本位置,那運(yùn)行的程序也就是vasp都不能調(diào)用。只是樓主注意下變量PROGRAMEXEC="/opt/vasp/vasp"設(shè)置在你的機(jī)器上是正確的就行。還有就是我覺得這個(gè)腳本寫得羅里羅嗦, 我整合個(gè)簡單點(diǎn)的 看著清楚的 供樓主參考下: #!/bin/bash #PBS -N 5-zmos2nrbig #PBS -j oe #PBS -q vasp #PBS -l walltime=9999:00:00 #PBS -l nodes=2:ppn=8 rm WAVECAR for a in 400 450 500 550 600 650 700 750 do cat >INCAR < SYSTEM=CGa LPLANE=.TRUE. NPAR=8 LREAL = .TRUE. PREC = Accurate EDIFF =1E-6 EDIFFG = -0.03 IALGO = 48 NELMIN = 4 ISYM = 0 ISPIN = 1 NBANDS = 240 LCHARG = .TRUE. LWAVE = .TRUE. LVTOT = .FALSE. NBLOCK = 1 NSW = 10 IBRION = 2 POTIM=0.3 ISMEAR = 0 SIGMA = 0.05 ENCUT=$a ISTART=0;ICHARG=2 ! echo "ENCUT=$a eV" mpirun XXX(此處是mpi的設(shè)置參數(shù),根據(jù)自己的機(jī)器設(shè)置) ./vasp > ./out(此處將vasp程序 vasp必須的文件和腳本放同一目錄,要嫌麻煩就在腳本里加幾個(gè)變量把路徑寫進(jìn)去) E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'` echo $a $E >>comment done |

| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +4 | vv迷 2026-03-22 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|
[考研] 263求調(diào)劑 +4 | yqdszhdap- 2026-03-22 | 4/200 |
|
|
[考研] 310求調(diào)劑 +4 | baibai1314 2026-03-16 | 4/200 |
|
|
[考研] 材料與化工085600,總分304,本科有兩篇sci參與,求調(diào)劑 +4 | 幸運(yùn)的醬醬 2026-03-22 | 5/250 |
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +5 | 我愛生物生物愛?/a> 2026-03-17 | 5/250 |
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 4/200 |
|
|
[考研] 一志愿南大,0703化學(xué),分?jǐn)?shù)336,求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +5 | 紙魚ly 2026-03-21 | 5/250 |
|
|
[考研] 279分求調(diào)劑 一志愿211 +14 | chaojifeixia 2026-03-19 | 15/750 |
|
|
[考研] 化學(xué)求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
[考研] 085600材料與化工 +8 | 安全上岸! 2026-03-16 | 8/400 |
|
|
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
[考研] 環(huán)境工程調(diào)劑 +9 | 大可digkids 2026-03-16 | 9/450 |
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +9 | 手機(jī)用戶 2026-03-17 | 9/450 |
|
|
[考研] 【同濟(jì)軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|