24小时热门版块排行榜    

查看: 300  |  回复: 0

490105179

木虫 (正式写手)

[求助] 测试晶格参数出现问题

我根据侯博士手册上的教程测试晶格参数,因为我这提交任务需要脚本,所以将run_a0添加在我的脚本里,但是测试的结果总是少一个,请大家帮忙分析一下。先贴上脚本文件。#!/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 j in 3.5  3.9  4.3

do
cat > POSCAR <<!
Pt
$j
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
1   
Direct
0.0000000000000000  0.0000000000000000  0.0000000000000000   
!
echo "a = $j";
mpirun -np $NP -machinefile $PBS_NODEFILE vasp
E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'`
V=`grep "volume" OUTCAR | tail -1 | awk '{printf "%12.4f \n" , $5}'`
echo $V $E >> EtVo.dat
#
#
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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 490105179 的主题更新
信息提示
请填处理意见