24小时热门版块排行榜    

查看: 624  |  回复: 0

wqg1806

木虫 (著名写手)

科研农民工

[求助] 求助弹性常数并行计算脚本

本人属于菜鸟,用这样的脚本,能算,总感觉不对,没有看到输出SUMMARY.求高人指点!
#PBS -o script.out
#PBS -j oe
#PBS -N akak
#PBS -q batch
#PBS -l nodes=1:ppn=6
#PBS -l walltime=2400:00:00
ulimit -s unlimited
TMPFILE=`whoami`_mpich_gm.tmp
GM_NODEFILE=/tmp/$TMPFILE
cat $PBS_NODEFILE > /tmp/$TMPFILE
echo Working directory is $PBS_O_WORKDIR
cd $PBS_O_WORKDIR
echo Runing on host `hostname`
echo Starting Time is `date`
echo Directory is `pwd`
echo This jobs runs on the following processors:
echo `cat $GM_NODEFILE`
NPROCS=`wc -l < $GM_NODEFILE`
echo This job has allocated $NPROCS nodes
## User Parallel Program ###########

mpirun -v -machinefile $GM_NODEFILE -np $NPROCS  /share/apps/vasp/vasp &> output


##User Parallel Program ###########
rm -f /tmp/$TMPFILE
echo Ending Time
echo Ending Time is `date`
#!/bin/sh
for i in -0.018 -0.015 -0.012 -0.09 -0.06 -0.03 0.00 \
0.03 0.06 0.09 0.012 0.015 0.018
do
echo $i | defvector.x
cp fort.3 POSCAR
####
cat > INCAR < SYSTEM = Nb
ENCUT = 400
ISTART = 0
ICHARG = 2
ISMEAR = 1; SIGMA = 0.2
NSW = 160; IBRION = 2
EDIFF = 1E-5
EDIFFG = -1E-2
ISIF = 2
POTIM = 0.2
PREC = Accurate
LWAVE = .FALSE.
LCHARG =.FALSE.
!
echo "delta = $i "; vasp
cp CONTCAR pos.$i
cp CONTCAR POSCAR
cat > INCAR < SYSTEM = Nb
ENCUT = 400
ISTART = 0
ICHARG = 2
ISMEAR = -5
EDIFF = 1E-5
PREC = Accurate
LWAVE = .FALSE.
LCHARG =.FALSE.
!
echo "delta = $i "; vasp
E=‘grep "TOTEN" out.$i | tail -1 | awk ’{printf "%12.6f \n", $5 }’‘
echo $i $E >>SUMMARY
done
回复此楼
一路走到黑
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wqg1806 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +53 医学老男孩 2026-08-13 118/5900 2026-08-19 18:14 by 医学老男孩
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 PnssoZglMj3v 2026-08-19 3/150 2026-08-19 18:04 by oCMhZAsgNFoO
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 PnssoZglMj3v 2026-08-19 4/200 2026-08-19 17:54 by oCMhZAsgNFoO
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 PnssoZglMj3v 2026-08-19 5/250 2026-08-19 17:42 by oCMhZAsgNFoO
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +4 瞬息宇宙 2026-08-15 4/200 2026-08-19 16:23 by allmydream
[基金申请] 今天基金会出结果吗?20260819 +8 kkkl_v 2026-08-19 8/400 2026-08-19 12:34 by Vivilian
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 时间戳变了,能看出什么问题? +10 基诺咪客 2026-08-17 14/700 2026-08-19 10:34 by 基诺咪客
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 重要消息,中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
信息提示
请填处理意见