Back to home page

sPhenix code displayed by LXR

 
 

    


File indexing completed on 2026-07-16 08:16:52

0001 #! /bin/bash
0002 
0003 set -euo pipefail
0004 
0005 # Electric field grid
0006 emin=360
0007 emax=440
0008 ne=21
0009 
0010 # Magnetic field grid
0011 bmin=1.15
0012 bmax=1.45
0013 nb=16
0014 
0015 # E-B angle grid
0016 amin=0.0
0017 amax=0.2
0018 na=20
0019 
0020 jobid=$1
0021 
0022 ie=$(( jobid / nb ))
0023 ib=$(( jobid % nb ))
0024 
0025 if (( jobid < 0 || jobid >= ne * nb )); then
0026   echo "ERROR: jobid=$jobid outside allowed range 0 to $((ne * nb - 1))"
0027   exit 1
0028 fi
0029 
0030 enow=$(( emin + 4 * ie ))
0031 #enow=$(awk -v emin="$emin" \
0032 #             -v emax="$emax" \
0033 #             -v ne="$ne" \
0034 #             -v i="$ie" \
0035 #'BEGIN {
0036 #   if (ne == 1) print emin;
0037 #   else print emin + (emax - emin) * i / (ne - 1);
0038 #}')
0039 
0040 bnow=$(awk -v ib="$ib" 'BEGIN { printf "%.2f", 1.15 + 0.02 * ib }')
0041 #bnow=$(awk -v bmin="$bmin" \
0042 #             -v bmax="$bmax" \
0043 #             -v nb="$nb" \
0044 #             -v i="$ib" \
0045 #'BEGIN {
0046 #   if (nb == 1) print bmin;
0047 #   else print bmin + (bmax - bmin) * i / (nb - 1);
0048 #}')
0049 
0050 echo "jobid=$jobid ie=$ie ib=$ib enow=$enow bnow=$bnow"
0051 
0052 script_dir="$(cd "$(dirname "${BASH_SOURCE[0]}")" && pwd)"
0053 
0054 if [[ ! -d "${script_dir}/gasfiles" ]]; then
0055   echo "ERROR: ${script_dir}/gasfiles does not exist."
0056   echo "Create it or make it a symbolic link before launching jobs."
0057   exit 1
0058 fi
0059 
0060 efile=$(printf "E%03d_B%03d.gas" "$ie" "$ib")
0061 output="${script_dir}/gasfiles/${efile}"
0062 
0063 echo GasModel "$enow" "$enow" 1 "$bnow" "$bnow" 1 "$amin" "$amax" "$na" "$output"
0064 GasModel "$enow" "$enow" 1 "$bnow" "$bnow" 1 "$amin" "$amax" "$na" "$output"
0065 
0066 echo "all done"