You need about 1 G of disk space. The data you need are
These can be downloaded and extracted using:
Download and install CASA; also see data and instructions for a quick test.
These VLA solar data are a small data set so they can be reduced quickly as a first introduction to calibration and imaging radio interferometry data in CASA. They comprise one scan on a primary flux scale calibrator B021+735, followed by scans alternating between the phase-reference B0945+664 and the Sun. The VLA was in D-array, with a maximum baseline of just over 1 km.
The most obvious observational differences with ALMA are the wavelength (20 cm) and the antenna diameter D=25 m, which means that the Sun can be observed in a single pointing. They are observed in circular polarization, whilst ALMA used linear polarization receivers (we only process total intensity here). These old VLA data also just have a single IF with a single channel, whilst ALMA data are multi-channel. As with ALMA, special observing modes are used and these data have been pre-cooked by applying applying flags reported by the on-line system and Tsys corrections (which will also be used, albeit measured differently, by ALMA).
Read the script sun-vla_new.py. Go to a terminal and subdirectory containing the script and the data DULK16081984.UVFITS and also sun_vla.flagcmds, and start CASA.
Run one step at a time by inserting the step number(s) in the script where it says mysteps = . You need to insert some missing values for parameters called SETTHIS1 - SETTHIS6, see below or in script. It is best to run one step at a time to start with.
To run the script, in CASA, type
CASA <>: execfile('sun-vla_new.py')
The steps in sun-vla_new.py (with main differences from ALMA) are:
CASA <>: !more listobs.last
CASA <>: default listobs
CASA <>: execfile('listobs.last')
================================================================================ MeasurementSet Name: /home/amsr/scratch/ALMA/CASANOTES/SOLAR/Script/sun_16081984.ms MS Version 2 ================================================================================ Observer: DULK Project: Observation: VLA Data records: 985344 Total elapsed time = 39590 seconds Observed from 16-Aug-1984/13:04:50.0 to 17-Aug-1984/00:04:40.0 (TAI) ObservationID = 0 ArrayID = 0 Date Timerange (TAI) Scan FldId FieldName nRows SpwIds Average Interval(s) ScanIntent 16-Aug-1984/13:04:50.0 - 13:08:30.0 1 0 0212+735 7374  [9.99] 13:23:40.0 - 13:44:50.0 2 1 QSUN16 33105  [9.99] 13:46:40.0 - 13:50:20.0 3 2 0945+664 7426  [9.99]
CASA <>: !more sun_vla.flagcmds mode='manual antenna='VA02' field='0' mode='manual' timerange='13:58:30~13:58:40' mode='manual' timerange='15:52:30~16:00:50' etc.
This is set up for the properties of an image made from the calibrated VLA data, qsun20cm_all.clean_ap.fits; you can use the downloaded version or make your own. Note the slightly smaller pixel size (to reduce artefacts due to resampling) and smaller image size (to reduce mosaicing time). You can change the clean parameters for the input image to select a time-range etc. if you want.
We will simulate observations at the lowest frequency which will eventually be available for ALMA, 36 GHz (lambda about 9 mm), in Band 1, to minimise the number of pointings needed. The usual continuum (TDM) mode of ALMA has four spectral windows (spw) of 2 GHz, each in 128 channels (for dual polarization). You can average these up to a limit determined by bandwidth smearing; for a mosaic one should be cautious in real life but here I use the full-sensitivity width of one spw, 1.875 GHz, which will be simulated as a single channel.
ALMA's field of view is given
FoV = 1.12*lambda/D
so you can estimate roughly how many pointings are needed to cover the Sun with 12-m or 7-m dishes (the task will do this for you exactly, at Nyquist spacing).
The VLA D-configuration is sensitive to scales up to at least 900 arcsec at 20 cm but the largest angular scale apparent in the VLA image is about 200 arcsec. The largest angular scale which the ALMA ACA is sensitive to is determined by the shortest antenna spacing, and is just over 200 arcsec but it is simplest in simulation just to use the main 12-m array (most compact configuration) and the ACA 7-m array.
The input image peak (if you use the suggested image) is 4.6 Jy/beam, where the primary beam correction has been applied. simalma interprets the input as Jy per pixel, so you can ask it to rescale the input if necessary. I used 10 Jy/pixel to get decent S/N in one second per pointing. I set the total time set to 21 min for each of the ACA and the main array. Although the main array is more sensitive, it needs 1003 pointings whilst the ACA only needs 351 pointings, so this gives approximately 1 visit per main array pointing and 3 visits per ACA pointing.
The simulated MS and diagnostic png's go in a directory here called qsun_sim. If you re-run the simulation, either set overwrite=T or use a different project name or delete the old one.
When it has finished, type 'ls qsun_sim'. You can inspect the .png's and the mosaic pointing files , *.ptg.txt. You will see two MS, and we will concatenate them using concat. The relative weights can be adjusted either more precisely (looking at the exact uv coverage) to optimise sensitivity, or to give finer/coarser resolution by giving more/less weight to the 12-m array. Here, we will roughly optimise sensitivity by weighting the 12-m v. 7-m arrays in proportion to the antenna areas. Enter suitable values at SETTHIS below. Enter the following in CASA:
default concat inp(concat) # to see the inputs, then enter: vis=['qsun_sim/qsun_sim.alma.cycle2.1.ms','qsun_sim/qsun_sim.aca.cycle2.i.ms'] concatvis='qsun_sim_concatvis.ms' # delete any old file of the same name first visweightscale=[X,Y] inp(concat) # check again and when happy: concat() # execute the task.
Set up clean yourself. Have a look at sun-vla_sim.py or at clean.last for initial guidance. Start by with
The main other things which you need to set are:
imagermode phasecenter cell imsize
If you have time, use help(clean) or the Cookbook and experiment with settings such as: