This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
analysis:fourfit [2015/06/03 00:32] Jamie McCallum [Phase Cal estimation] |
analysis:fourfit [2017/06/06 01:26] (current) Jamie McCallum |
||
---|---|---|---|
Line 1: | Line 1: | ||
+ | Fourfit is the software that performs the multi-band delay fitting to the correlator visibilities. At present, magilla is set up with the latest version and has scripts installed which remove some of the donkey-work in preparing properly fringe-fitted data sets. The following recipe will work for data sets correlated at Curtin. | ||
+ | |||
+ | * Download the following files from the Pawsey server to ~/Downloads: AUG013.MARK4.tar.gz, aug013.vex and aug013.v2d | ||
+ | * Run "postcorr_prep.sh aug013" to copy the data across to /data/AUSTRAL/ extract them into the appropriate folders. It will also get you to check the ovex files station ID assignments and generate the apriori SNR file (delete all the leading text before the first two headers). | ||
+ | * After completing the initial fringe-fitting (using the default settings defined by /data/AUSTRAL/cf_aust), you will be prompted to enter a command like "cd /data/AUSTRAL/aug013/1344 ; export CORDATA=/data/AUSTRAL/aug013/ ; export DATADIR=/data/AUSTRAL/aug013/ ; aedit -f alist.ed.out". Do so, and you will be able to make a first inspection of the data. | ||
+ | |||
+ | ====== Aedit ====== | ||
+ | |||
+ | After running "edi dup snr", run "psplot" and check for either grey scans (uncorrelated data) or brown ones indicating non-detections. If found, make sure there's a sensible/known reason for them (module failure, antenna windstowed, etc). If there's no explanation, get in contact with the correlator ASAP. | ||
+ | |||
+ | ------------ | ||
+ | |||
+ | |||
+ | |||
+ | |||
There are two machines set up to run fourfit at present - hex2 (through the hex2:1 VNC session) or ACECOMPUTE (using acecompute.antcrc.utas.edu.au:1). The two systems give slightly different results, thought to be due to the difference in the 32/64-bit compilers. The 32-bit version on hex2 occasionally gives unreliable estimates of the SNR, while the 64-bit version on ACECOMPUTE does not. Bonn also uses the 64-bit version, so ACECOMPUTE in strongly recommended for use in processing the data. | There are two machines set up to run fourfit at present - hex2 (through the hex2:1 VNC session) or ACECOMPUTE (using acecompute.antcrc.utas.edu.au:1). The two systems give slightly different results, thought to be due to the difference in the 32/64-bit compilers. The 32-bit version on hex2 occasionally gives unreliable estimates of the SNR, while the 64-bit version on ACECOMPUTE does not. Bonn also uses the 64-bit version, so ACECOMPUTE in strongly recommended for use in processing the data. | ||
Line 15: | Line 30: | ||
A more concise version is included below: | A more concise version is included below: | ||
- | for i in $(cat /tmp/scans) ; do echo $i ; fourfit -c cf_1234 $i ; done | + | <code> |
+ | parallel fourfit -c cf_1234 ::: $(ls -d 12*) | ||
+ | </code> | ||
Once the fourfit processes are complete (15-30 minutes?), run alist with ''alist * '' and then edit the file with the following commands | Once the fourfit processes are complete (15-30 minutes?), run alist with ''alist * '' and then edit the file with the following commands | ||
sed -i 's/ 16383 / 1234 /g' alist.out | sed -i 's/ 16383 / 1234 /g' alist.out | ||
- | grep -v " LL " alist.out | grep -v " ii " | grep -v " xx " | grep -v " ee " | grep -v " HH " | grep -v " gg " | grep -v " WW " > alist.ed.out | + | grep -v '0.000000 0.000 -1 -1' alist.out > alist.ed.out |
Download the skd file and reprocess it into an ovex file. NB - ''nskd_ovex.pl'' needs to have internet access - you'll need to define the proxy settings for most machines. | Download the skd file and reprocess it into an ovex file. NB - ''nskd_ovex.pl'' needs to have internet access - you'll need to define the proxy settings for most machines. | ||
Line 77: | Line 94: | ||
alist * | alist * | ||
sed -i 's/ 16383 / 1234 /g' alist.out | sed -i 's/ 16383 / 1234 /g' alist.out | ||
- | grep -v " LL " alist.out | grep -v " ii " | grep -v " HH " | grep -v " ee " | grep -v " xx " | grep -v " gg " | grep -v " TT " | grep -v " WW " > alist.ed.out | + | grep -v '0.000000 0.000 -1 -1' alist.out > alist.ed.out |