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.3-28-pipeline-2023.0.0.36/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_contl inecontline_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()
|