#!/bin/bash

#
PARATEC="$HOME/paratecSGL/bin/paratec.mpi"
VOLUME="$HOME/BerkeleyGW/Visual/volume.py"
SURFACE="$HOME/BerkeleyGW/bin/surface.x"
AVERAGE="$HOME/BerkeleyGW/Visual/average.py"
MPIRUN=mpirun
#
cd ./1-scf
mkdir ./vis
$MPIRUN $PARATEC
$VOLUME ./input paratec ./vis/CHARGE_.a3Dr a3dr ./rho.cube cube true re false
$VOLUME ./input paratec ./vis/VEFFAVG.a3Dr a3dr ./vtot.cube cube true re false
$VOLUME ./input paratec ./vis/VELECAVG.a3Dr a3dr ./velec.cube cube true re false
$SURFACE ./surface.par &> ./surface.out
$AVERAGE ./velec.cube cube 0.5 0.5 0.5 &> ./velec_avg_face.out
$AVERAGE ./vtot.cube cube &> ./vtot_avg_vol.out
cd ..
#
cd ./2-nscf
$MPIRUN $PARATEC
cd ..
#
