run_casa.sh: #!/bin/sh #Don't put any commands before the #SBATCH options or they will not work
#SBATCH --export ALL # Export all environment variables to the job.
#SBATCH --mem=64G # Amount of memory needed by the whole job.
#SBATCH -D /lustre/aoc/users/mlacy/srdp/selfcal/req1236/working # Working directory set to your Lustre area
#SBATCH --mail-type=END,FAIL # Send email when Jobs end or fail # casa's python requires a DISPLAY for matplot, so create a virtual X server
xvfb-run /home/casa/packages/pipeline/casa-6.5.34-289-pipeline-2023.01.0.36124/bin/casa --pipeline --nogui -c selfcal_script.py >& casa.out
selfcal_script.py: # This CASA pipescript is meant for use with CASA 6.5.3 and pipeline 2023.1.0.42
context = h_init()
context.set_state('ProjectSummary', 'observatory', 'Karl G. Jansky Very Large Array')
context.set_state('ProjectSummary', 'telescope', 'EVLA')
try:
hifv_importdata(vis=['23A-157.sb44369379.eb44370721.60147.8327970949.ms'], datacolumns={'data': 'raw','corrected': 'regcal_contline_all'})
hifv_flagtargetsdata(pipelinemode="automatic")
hif_mstransform(pipelinemode="automatic")
hif_checkproductsize(maximsize=16384)
hif_makeimlist(specmode='cont',datatype='regcal')
hif_makeimages(hm_cyclefactor=3.0)
hif_selfcal(pipelinemode="automatic")
hif_makeimlist(specmode='cont',datatype='selfcal')
hif_makeimages(hm_cyclefactor=3.0)
hifv_pbcor(pipelinemode="automatic")
hifv_exportdata(imaging_products_only=True)
finally:
h_save()
| 
