| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 1050 | 回復(fù): 2 | |||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話(huà)題的所有回帖 | |||
mlbiao新蟲(chóng) (初入文壇)
|
[交流]
請(qǐng)教一個(gè)VASP腳本問(wèn)題,謝謝~ 已有2人參與
|
||
|
剛接觸VASP,對(duì)linux SHELL也不是很熟,故在此請(qǐng)教: 如下有兩文件,第一個(gè)是常見(jiàn)的vasp運(yùn)行PBS ; 第二個(gè)是好像是用來(lái)進(jìn)行ENCUT 收斂情況的測(cè)試(或類(lèi)似其他測(cè)試)。 請(qǐng)問(wèn)怎樣這兩個(gè)合成一個(gè)PBS呢?也就是要仍以PBS運(yùn)行VASP。謝謝~ PS: 沒(méi)金幣只好發(fā)成資源帖 #PBS -N vasp #PBS -j oe #PBS -l nodes=5:ppn=8 source /public/software/profile.d/intel-env.sh source /public/software/profile.d/openmpi-intel-env.sh NP='cat $PBS_NODEFILE | wc -l' cd $PBS_O_WORKDIR mpirun -np -$NP -machinefile $PBS_NODEFILE /public/software/vasp5.2/vasp ---------------------------------------------------------- #!/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 |
鐵桿木蟲(chóng) (正式寫(xiě)手)
|
正規(guī)的做法是用-W選項(xiàng),見(jiàn)下: It is possible to start a job on the condition that another one completes beforehand; this may be necessary for instance if the input to one job is generated by another job. Job dependency is defined in PBS using the -W flag. To illustrate with an example, suppose you need to start a job using the script second_job.sh after another job finished successfully. Assume the first job is started using script first_job.sh and the command to start the first job $ qsub first_job.sh returns the job ID 7777. Then, the command to start the second job is $ qsub -W depend=afterok:7777 second_job.sh This job dependency can be further automated (possibly to be included in a bash script) using environment variables: JOB_ID_1=`qsub first_job.sh` JOB_ID-2=`qsub -W depend=afterok:$JOB_ID_1 second_job.sh` |
鐵桿木蟲(chóng) (正式寫(xiě)手)
|
合并會(huì)造成一個(gè)問(wèn)題就是提交作業(yè)之后,還沒(méi)有算,就開(kāi)始執(zhí)行讀取total energy 命令。除非你加入反饋代碼,VASP每計(jì)算完畢一次反饋給shell,之后再執(zhí)行讀取操作。建議你提交和讀取分兩次進(jìn)行,代碼還需修改一下。 我寫(xiě)一個(gè),供參考,些許有小錯(cuò),根據(jù)情況修改。 (1) 提交腳本 #!/bin/sh cd=`pwd` for i in 150 200 250 300 350 400 do mkdir $i cd $i cp ../POSCAR . cp ../KPOINTS . cp ../POTCAR . cat > INCAR <<! SYSTEM = Si-Diamond ENCUT = $i ISTART = 0 ; ICHARG = 2 ISMEAR = -5 PREC = Accurate ! bash ../vasp_pbs.sh cd .. done (2)讀取腳本 #!/bin/sh cd=`pwd` for i in 150 200 250 300 350 400 do cd $i echo "ENCUT = $i eV" ; time vasp E='grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'' echo $i $E >>../comment cd .. done (3)提交作業(yè)腳本vasp_pbs.sh #PBS -N vasp #PBS -j oe #PBS -l nodes=5:ppn=8 source /public/software/profile.d/intel-env.sh source /public/software/profile.d/openmpi-intel-env.sh NP='cat $PBS_NODEFILE | wc -l' cd $PBS_O_WORKDIR mpirun -np -$NP -machinefile $PBS_NODEFILE /public/software/vasp5.2/vasp |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工272求調(diào)劑 +16 | 阿斯蒂芬2004 2026-03-28 | 16/800 |
|
|---|---|---|---|---|
|
[考研] 300求調(diào)劑,材料科學(xué)英一數(shù)二 +10 | leaflight 2026-03-24 | 10/500 |
|
|
[考研] 考研調(diào)劑 +7 | 小蠟新筆 2026-03-29 | 7/350 |
|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|
[考研] 0856求調(diào)劑 +7 | 楒桉 2026-03-28 | 7/350 |
|
|
[考研] 321求調(diào)劑 +7 | 璞玉~~ 2026-03-25 | 8/400 |
|
|
[考研] 一志愿211 初試270分 求調(diào)劑 +8 | 谷雨上岸 2026-03-23 | 9/450 |
|
|
[考研] 0703化學(xué)調(diào)劑,求導(dǎo)師收 +9 | 天天好運(yùn)來(lái)上岸?/a> 2026-03-24 | 10/500 |
|
|
[考研] 壓國(guó)家一區(qū)線(xiàn),求導(dǎo)師收留,有恩必謝! +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
[考研] 求調(diào)劑 +4 | 零八# 2026-03-27 | 4/200 |
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 12/600 |
|
|
[考研] 一志愿華東理工大學(xué)081700,初試分?jǐn)?shù)271 +6 | kotoko_ik 2026-03-23 | 7/350 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
[考研] 0703化學(xué)338求調(diào)劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|
|
[考研] 342求調(diào)劑 +3 | 加油a李zs 2026-03-26 | 3/150 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
[考研] 打過(guò)很多競(jìng)賽,085406控制工程300分,求調(diào)劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
[考研] 296求調(diào)劑 +4 | 汪?! 2026-03-25 | 7/350 |
|