Setonix MWA Stack 2024-11-05 EoR.sh

# request 1/8 of a GPU node for an hour
salloc --nodes=1 --partition=gpu --account=mwaeor-gpu -t 01:00:00 --gres=gpu:1

# load all the modules
module load mwalib/1.5.0 birli/0.15.1 hyperdrive-amd-gfx90a/0.4.1 wsclean/3.4-idg-everybeam

export metafits="1090701368.metafits"
wget 'https://projects.pawsey.org.au/birli-test/'$metafits
export raw="1090701368_20140729203555_gpubox01_00.fits"
wget 'https://projects.pawsey.org.au/birli-test/'$raw
export srclist="srclist_pumav3_EoR0LoBES_EoR1pietro_CenA-GP_2023-11-07.fits"
wget 'https://github.com/JLBLine/srclists/raw/master/'$srclist

# mwalib
wget https://raw.githubusercontent.com/MWATelescope/mwalib/main/examples/mwalib-print-context.py
python mwalib-print-context.py -m $metafits $raw

# birli
export prep_uvfits=birli_${metafits%.*}.uvfits
birli -m $metafits --uvfits-out=$prep_uvfits $raw

# hyperdrive di-cal
export MWA_BEAM_FILE="/scratch/references/mwa/beam-models/mwa_full_embedded_element_pattern.h5"
export solns=hyp_${metafits%.*}.fits
hyperdrive di-calibrate --data $metafits $prep_uvfits --source-list $srclist --outputs $solns

# hyperdrive apply
export cal_ms=hyp_${metafits%.*}.ms
hyperdrive solutions-apply --data $metafits $prep_uvfits --solutions $solns --outputs $cal_ms

# wsclean
wsclean \
    -name "wsclean_${metafits%.*}" \
    -size 4096 4096 \
    -scale 20asec \
    -weight briggs 0 \
    -nmiter 2 \
    -niter 10000 \
    -mgain 0.85 \
    -gain 0.1 \
    -gridder idg -grid-with-beam \
    -mwa-path /scratch/references/mwa/beam-models/ \
    -temp-dir /tmp \
    $cal_ms