I16293 Water VLA Observation Summary

TIP Last Modified: AlWootten - 20 Nov 2018

Archive Extraction and Download

All data can be extracted from the EVLA archive available from the NRAO User Portal. Since these data sets are extremely large, one needs to use the CV lustre to both store and process (with CASA) these data. Steps to extraction and download are:
  1. Make the archive extraction request from NRAO User Portal, and wait for the data to be extracted into the archive directory.
  2. Login to username@cvpost-master.cv.nrao.edu to request a lustre (cvpost) node: ssh username@cvpost-master.cv.nrao.edu
  3. List all used and free lustre nodes: qstat -a.
  4. Request a node: nodescheduler --request 45 1 to request one node for 45 days (the current maximum), for example. You will get an email telling you which node you have been assigned. Let's assume in the following that you have been assigned "cvpost005".
  5. ssh into your assigned node: ssh cvpost005
  6. If you are going to use a VNC connection later, start the VNC server: vncserver
  7. cd to your data directory (you should have already been given a login on the lustre). NOTE: Data directories are something like /lustre/cv/projects/awootten.
  8. Use the following wget command to copy files from the archive:
wget -b -e robots=off --wait 1 -rc -nH --cut-dirs=1 --header="accept-encoding: gzip" --user=Chicken --password=Butt https://archive.nrao.edu/secured/Chicken/

  • If you are copying files from the public archive, use the following command (use wildcards on the MSNAME if necessary):
wget -b -e robots=off --wait 1 -rc -nH --cut-dirs=1 --header="accept-encoding: gzip" ftp://ftp.aoc.nrao.edu/e2earchive/MSNAME

Links to data reduction details provided in link in table.

  • 2018-11-20: Updated for 18A-507


Session Log
Project Date Array Band Targets SBID Duration (hrs)
18A-507A 2018-11-21 D K IRAS16293A1 4982259 1.0

Comments: Only 26 antennas available. First session. Copied from EVLA pipeline and located at /.lustre/cv/projects/awootten/I16293/2018Water/. EVLA pipeline flagged most water maser channels. Made cont.dat file to use rerunning pipeline, defining continuum channels so line channels will not be flagged. In addition to creating a 'cont.dat' file, we advise to modify the following tasks in 'casa_pipescript.py' for spectral line data:

  1. hifv_hanning: Hanning smoothing lessens the Gibbs ringing from strong spectral features, usually strong, narrow RFI, or very strong spectral lines such as masers. Hanning smoothing, however, reduces the spectral resolution. Therefore, depending on the data and the science case, one may or may not choose to apply Hanning smoothing. Disable the application of Hanning-smoothing in the pipeline if such smoothing is not needed or desired by commenting out hifv_hanning or simply removing the step from 'casa_pipescript.py'. Do not need to modify.
  2. hifv_flagata: The pipeline, by default, flags 5% of the data on each spw edge as well as the first and last 10 channels of each baseband. In some cases, for example spectral surveys, lines may fall right on such frequencies. The edgespw, fracspw, and baseband parameters in hifv_flagata can be adjusted to flag different percentages of the edges. Do not need to modify
  3. hifv_targetflag: Flagging prior to this step was only applied to the calibrator scans, which should be line-free. But hifv_targetflag attempts to auto-flag all fields including target fields. The rflag mode in CASA's flagdata is designed to remove outliers that deviate from a mean level. Strong spectral lines can fulfill this criterion and be flagged. The 'cont.dat' file will ensure that rflag will only be applied to the continuum frequency ranges specified in it. Alternatively, hifv_targetflag can be turned off completely, or, by specifying intents='*CALIBRATE*' (and thereby omitting '*TARGET*'), one can restrict the flagging to the calibrator data only, leaving the target data (scans with the TARGET intent) untouched. In either case, we recommend manual flagging for the spectral line frequency ranges after the pipeline has finished processing.
  4. hifv_statwt: A similar argument applies to the hifv_statwt step, where the visibilities are weighted by the square of the inverse of their RMS noise. Strong spectral lines will increase the RMS and will therefore be down-weighted. The cont.dat file will restrict statwt to only use the continuum frequency ranges for the rms and weight calculations and thus prevent the inclusion of spectral features. Alternatively, hifv_statwt can be excluded from the pipeline altogether and the CASA task statwt can be executed manually after the pipeline has finished, where statwt's parameter fitspw should be set to continuum channels only.
  5. casa -r 5.1.2-4 --pipeline run in /lustre/cv/projects/awootten/I16293/2018Water/rawdata/UJ3510 using execfile('casa_pipescript-20181110S194136.py') (Nov 20, 2018) One needs to avoid spectral line channels--the 'cont.dat' construct only works with casa 5.4.2 or later, Tobin tells me. I just eliminated flagging of the source altogether.

18A-507A 2018-11-21 D K IRAS16293A1 4982259 1.0
Comments: Second session. Extracted from archive and located as above.

AmmoniaTalks -- AlWootten - 2018-11-20
Topic attachments
I Attachment Action Size Date Who Comment
casa_pipescript-20181110S194136.py.txttxt casa_pipescript-20181110S194136.py.txt manage 2 K 2018-11-20 - 14:46 AlWootten EVLA 18A_507 pipeline script modified for spectral lines
Topic revision: r4 - 2018-11-30, AlWootten
This site is powered by FoswikiCopyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding NRAO Public Wiki? Send feedback