ABACUS/scripts/build_LiebLin_c_scan_k_fixe...

#! /bin/zsh

# Runs `build_LiebLin_c_scan_k_fixed_dsfs` and `_plots` for nk = 1 up to nkmax.

if [[ $# -ne 6 ]]; then
    echo "Arguments needed: whichDSF, kBT, target_sumrule, N, nkmax (max momentum in units of kF/4), width."
    exit 1
fi

if [[ $1 != 'd' && $1 != 'g' && $1 != 'o' ]]; then
    echo "Only the d, g and o scanning options are implemented."
    exit 1
fi

if [[ $2 -lt 0 ]]; then
    echo "Temperature kBT must be > 0."
    exit 1
fi

if [[ $3 -lt 0 || $3 -gt 1 ]]; then
    echo "Requirement: 0 < target_sumrule < 1."
    exit 1
fi

if [[ $4 -lt 0 ]]; then
    echo "N must be > 0."
    exit 1
fi

if [[ $5 -lt 0 ]]; then
    echo "nkmax must be > 0."
    exit 1
fi

if [[ $6 -lt 0 ]]; then
    echo "width must be > 0."
    exit 1
fi


whichDSF=$1
kBT=$2
target_sumrule=$3
N=$4
nkmax=$5
width=$6

for nk in {1..$nkmax}
do
    # produce dsfs files
    build_LiebLin_c_scan_k_fixed_dsfs.sh $whichDSF $kBT $target_sumrule $N $nk $width
    # produce plots
    build_LiebLin_c_scan_k_fixed_dsfs_plot.sh $whichDSF $kBT $target_sumrule $N $width
done