CHIPS is available on Garrawarla
/astro/mwaeor/MWA/chips
Plotting scripts in Python available in /astro/mwaeor/MWA/chips:
plotchips_1d.py - 1D PS
plotchips.py - 2D PS
plotchipsratio.py - ratio in 2D of two different datasets
External cloud processing pared down version of CHIPS (v2.8)
CHIPS version for external computing
Full version CHIPS is a three-step data analysis suite:
Grid visibilities over time as a function of frequency (parallel over frequency with OpenMP)
Combine frequencies into a data cube
Perform spectral Fourier Transform to form power
This version of CHIPS includes Steps 1 and 2 only
CHIPS software:
ANSI C, compiled with gcc, using OpenMP parallelisation dependencies required: CFITSIO, LAPACK, BLAS, OpenMP, LIBSLA (astronomy) file input: UVFITS file outputs: complex floats, doubles (binary files, little endian)
USAGE:
Example usage: This example has 384 frequency channels split over 24 input uvfits files. Each channel is gridded separately, and then these frequencies are combined into a cube with prepare_diff:
#!/bin/bash
source env_variables.sh
ext='test_highband_eor0'
for i in 1 2 3 4 5 6 7 8 9
do
./gridvisdiff $OUTPUTDIR/uvdump_0${i}.uvfits 1160576480 ${ext} 1 -f 0 -e 1
echo ${i}
done
for i in 10 11 12 13 14 15 16 17 18 19 24
do
./gridvisdiff $OUTPUTDIR/uvdump_${i}.uvfits 1160576480 ${ext} 1 -f 0 -e 1
echo ${i}
done
./prepare_diff ${ext} 384 0 'yy' ${ext} 1 -p 2.0 -n 167.075e6 ./lssa_fg_thermal ${ext} 384 80 'yy' 300. ${ext} 0 1 0
IMPORTANT For data with Nchan neq 384, kriging cannot be used. Instead, run with kriging off. E.g., for 120 channels:
./prepare_diff ${ext} 384 0 'yy' ${ext} 1 -p 2.0 -n 167.075e6 ./lssa_fg_thermal ${ext} 120 80 'yy' 300. ${ext} 1 1 0
COMMAND LINE SWITCHES AVAILABLE FOR NON-STANDARD MODES:
fprintf(stderr,"Usage: %s <options> uvfits_filename obs_id extension band \n",argv[0]);
fprintf(stderr,"\t -p period (seconds)\n");
fprintf(stderr,"\t -ew flag 14m EW baselines\n");
fprintf(stderr,"\t -c chanwidth (Hz)\n");
fprintf(stderr,"\t -n lowest frequency (Hz)\n");
fprintf(stderr,"\t -field fieldnum (0=EoR0, 1=EoR1, >2 anything)\n");
fprintf(stderr,"\t -u umax\n");
exit(1);