...
| Code Block |
|---|
| language | bash |
|---|
| theme | Midnight |
|---|
|
# On Garrawarla, test vcsbeam version
$ module use /pawsey/mwa/software/python3/modulefiles
$ module load vcsbeam
$ make_mwa_tied_array_beam -V
MWA Beamformer v2.18.8_99df2a4 |
...
| Code Block |
|---|
| language | bash |
|---|
| theme | Midnight |
|---|
|
$ cd /astro/mwavcs/asvo/251861
$ wget -O 1318345216.fits http://ws.mwatelescope.org/metadata/fits?obs_id=1318345216
$ cd /astro/mwavcs/vcs/1318345216/combined
$ for freq in {109..120}
do
sed "s/FREQ/${freq}/" < make_legacy.sbatch > make_legacy_ch${freq}.sbatch
sbatch make_legacy_ch${freq}.sbatch
done |
...
| Code Block |
|---|
| language | bash |
|---|
| theme | Midnight |
|---|
|
$ cd /astro/mwavcs/vcs/1318345216/vcsbeam/incoh
$ sbatch vcsbeam_incoh.sbatch
# ... After beamforming finishes...
$ module load singularity
$ pawsey/mwa/singularity/presto/presto.sif prepfold -psr J0034-0721 -nosearch -nodmsearch incoh*.fits |
| Code Block |
|---|
| language | bash |
|---|
| theme | Midnight |
|---|
| title | /astro/mwavcs/vcs/1318345216/vcsbeam/incoh/vcsbeam_incoh.sbatch |
|---|
| collapse | true |
|---|
|
#!/bin/bash -l
#SBATCH --nodes=12
#SBATCH --ntasks-per-node=1
#SBATCH --cpus-per-task=1
#SBATCH --mem=370gb
#SBATCH --partition=gpuq
#SBATCH --gres=gpu:1
#SBATCH --time=00:30:00
#SBATCH --account=mwavcs
#SBATCH --nice=0
#SBATCH --export=NONE
module use /pawsey/mwa/software/python3/modulefiles
module load vcsbeam
srun -N 12 -n 12 make_mwa_incoh_beam \
-m /astro/mwavcs/vcs/1318345216/combined/1318345216_legacy.metafits \
-b 1318345216 \
-T 600 \
-f 109 \
-d /astro/mwavcs/vcs/1318345216/combined \
-o incoh |
...
| Code Block |
|---|
| language | bash |
|---|
| theme | Midnight |
|---|
|
$ module load singularity
$ /pawsey/mwa/singularity/presto/presto.sif prepfold -psr J0034-0721 -nosearch -nodmsearch incoh*.fits |
...
| Code Block |
|---|
| language | bash |
|---|
| theme | Midnight |
|---|
| title | /astro/mwavcs/vcs/1318345216/cal/1318345216/rts/flagged_channels.txt |
|---|
| 0
1
16
30
31 |
0
1
16
30
31 |
Selecting only the last 38% of the observation for calibration
One of the parameters passed to the RTS is the "basename" of the files containing the visibilities. In order to distinguish between the RFI-ridden visibilities and the clean ones, I will create soft links to the clean gpubox files (created in an earlier step) and change the case of the filenames (gpubox → GPUBOX). The first and last filenames are 1318345216_20111011150016_gpubox01_00.fits and 1318345216_20111011151015_gpubox12_00.fits, with the timestamp in the format YYYYMMDDHHMMSS. To get the last 38%, I need to start 6:12 (mm:ss) into the observation, which means that my first clean file will be 1318345216_20111011150628_gpubox01_00.fits.
| Code Block |
|---|
| language | bash |
|---|
| theme | Midnight |
|---|
| title | Create source list |
|---|
|
$ cd /astro/mwavcs/vcs/1318345216/cal/1318345216/vis
$ ls 1318345216_2011101115062[89]_gpubox*.fits 1318345216_201110111506[3-5]?_gpubox*.fits 1318345216_20111011150[7-9]??_gpubox*.fits 1318345216_20111011151???_gpubox*.fits | while read f; do ln -s ${f} $(echo ${f} | tr 'gpubox' 'GPUBOX'); done |
After doing this, the RTS ".in" file must be altered so that it's only looking for files with GPUBOX in the filename (see below).
After running the RTS (see below) and checking the solutions (using the python scripts supplied in VCSTools), I found that I also had to add 14 to the list of flagged channels, and 52, 53, 54, 59 and 61 to the list of flagged tiles.
...
| Code Block |
|---|
| language | bash |
|---|
| theme | Midnight |
|---|
| title | /astro/mwavcs/vcs/1318345216/cal/1318345216/rts/rts_1318345216.in |
|---|
| collapse | true |
|---|
|
ReadAllFromSingleFile=
BaseFilename=/astro/mwavcs/vcs/1318345216/cal/1318345216/vis/*_gpuboxGPUBOX
ReadGpuboxDirect=0
UseCorrelatorInput=1
ReadMetafitsFile=1
MetafitsFilename=/astro/mwavcs/vcs/1318345216/1318345216
DoCalibration=
doMWArxCorrections=1
doRawDataCorrections=1
doRFIflagging=0
useFastPrimaryBeamModels=0
generateDIjones=1
applyDIcalibration=1
UsePacketInput=0
UseThreadedVI=1
MaxFrequency=170.24
ObservationFrequencyBase=138.895
ObservationTimeBase=2455846.1252314816
ObservationPointCentreHA=-0.35952189154909175
ObservationPointCentreDec=-5.98757073636307
ChannelBandwidth=0.04
NumberOfChannels=32
CorrDumpsPerCadence=128
CorrDumpTime=1.0
NumberOfIntegrationBins=6
NumberOfIterations=1
StartProcessingAt=0
ArrayPositionLat=-26.7033194
ArrayPositionLong=116.6708152
ArrayNumberOfStations=128
ArrayFile=
SourceCatalogueFile=/astro/mwavcs/vcs/1318345216/cal/1318345216/srclist_pumav3_EoR0aegean_EoR1pietro+ForA_1318345216_patch1000.txt
NumberOfCalibrators=1
NumberOfSourcesToPeel=0
calBaselineMin=20.0
calShortBaselineTaper=40.0
FieldOfViewDegrees=1 |
...