*note code coverage only applies to
python3
and C
code, CUDA
code is not currently supported by free code coverage software
WODEN
is C / GPU code designed to be able to simulate low-frequency radio interferometric data. It is written to be simplistic and fast to allow all-sky simulations. Although WODEN
was primarily written to simulate Murchinson Widefield Array (MWA, Tingay et al. 2013) visibilities, it is capable of bespoke array layouts, and has a number of primary beam options. WODEN
outputs uvfits
files. The WODEN
documentation lives here on readthedocs. If your internet has broken and you have already installed WODEN
, you can build a local copy by navigating into WODEN/docs/sphinx
and running make html
(you'll need to have doxygen
installed).
Note Jul 2024: Jack is back on a temporary contract to do some
WODEN
development. Keep your eyes peeled for new features!
New in version 2.3-alpha:
- Full polarisation is back! You can now simulate QUV visibilities, and they can be independent to Stokes I. In fact, you can set Stokes I to zero in the sky model, and still enter QUV information.
- Stokes V can either be a power-law, curved power-law, be a polarisation fraction of Stokes I, or be a list of flux densities.
- Linear polarisation (Q and U) can either be a power-law, curved power-law, be a polarisation fraction of Stokes I, or be a list of flux densities. The first three are always used in conjunction with a rotation measure. List-type fluxes can be given separately for Q and U, which makes Q and U independent of each other. Alternatively, a single list of linear polarisation fluxes can be given, which is used in conjunction with a rotation measure to determine Q and U.
New in version 2.2-alpha:
- Thanks to Marcin Sokolowski, and code from the PaCER Blink project, there is nominal support for AMD GPUs via ROCm. I've managed to compile and run on the Pawsey Setonix cluster, but this is the only AMD GPU I have access to. If you have an AMD GPU, please test and let me know how you go!
- There are some installation examples in
WODEN/templates/
for a CUDA and an AMD installation on superclusters. Hopefully you can adapt these to your needs if the Docker images don't work for you.- CUDA arch flag on compilation is changed to setting an environment variable
CUDAARCHS
to useCMakes
CMAKE_CUDA_ARCHITECTURES
feature.New in version 2.1:
- Dipole flagging for the MWA FEE beam is now enabled (
--use_MWA_dipflags
)- Dipole amplitudes for the MWA FEE beam is now enabled (
--use_MWA_dipamps
). Needs a modified metafits file- Things are still all Stokes I; plan is to implement QUV in version 2.3
New in version 2.0:
- Large swathes of the code have been transferred from
C
inpython
- The
C/CUDA
code that remains is called directly frompython
now, meaning the.json
and.dat
files are no longer needed- A new FITS sky model format is supported, meaning you can run simulations with greater than 25 million components via lazy loading
- I've done away with using Stokes QUV to speed things up. The plan is to implement some kind of rotation measure model to include Q/U, and reinstate Stokes V in a future release
Please be aware, before WODEN
version 1.4.0, in the output uvfits
files, the first polarisation (usually called XX) was derived from North-South dipoles, as is the labelling convention according to the IAU. However, most uvfits
users I've met, as well as the data out of the MWA telescope, define XX as East-West. So although the internal labelling and mathematics within the C/GPU code is to IAU spec, by default, run_woden.py
now writes out XX as East-West and YY as North-South. From version 1.4.0, a header value of IAUORDER = F
will appear, with F
meaning IAU ordering is False, so the polarisations go EW-EW, NS-NS, EW-NS, NS-EW. If IAUORDER = T
, the order is NS-NS, EW-EW, NS-EW, EW-NS. If there is no IAUORDER
at all, assume IAUORDER = T
.
If you have feature requests or want to contribute to WODEN
, have a read of the
guide to contributing to get started. I welcome the feedback and/or help!
Jack Line \ July 2024
Read the comprehensive installation guide on readthedocs.
The quickest way is to use a docker image (for a CUDA arch of 7.5 for example):
docker pull docker://jlbline/woden-2.3:cuda-75
and then run things through Docker or Singularity (more detail on that on ReadTheDocs). These versions come bundled with all the MWA FEE primary beam files which is nice.
Docker versions tested to run on various Australian superclusters are listed below:
jlbline/woden-2.3:cuda-60
- tested on Swinburne OzStarjlbline/woden-2.3:cuda-80
- tested on Swinburne Ngarrgu Tindebeekjlbline/woden-2.3:setonix
- tested on Pawsey Setonix (this is very experimental)To get the most control, install yourself. Again, go to ReadTheDocs for details, but in short, you will need the dependencies:
Once you have those, compilation is done via CMake
. Ideally, this will be enough:
$ cd WODEN
$ mkdir build && cd build
$ cmake ..
$ make -j 4
You then install the python3
code via something like
$ cd WODEN
$ pip install -r requirements.txt
$ pip install .
with a couple of post-compilation environment variables needed to use the MWA FEE primary beam model.
There are two routes to test WODEN:
ctest
, which will test your dependencies (see cmake testing on readthedocs)WODEN
is run via the script run_woden.py
, which will launch the woden
executable. All run_woden.py
arguments are explained here on readthedocs (or you can type run_woden.py --help
).
As WODEN
is written primarily for the MWA, the simplest way to feed in observational properties is via a metafits
file. You can obtain metafits files here. A very minimalistic example command using a metafits file looks like
run_woden.py \
--ra0=50.67 --dec0=-37.2 \
--cat_filename=srclist_msclean_fornaxA_phase1+2.fits \
--metafits_filename=1202815152_metafits_ppds.fits \
--primary_beam=MWA_FEE
where the array layout, observational settings, frequency and time resolution are all read in from the metafits file. All you have to specify is a phase centre (ra0
, dec0
), a sky model (cat_filename
), and a primary beam (in this case the MWA Fully Embedded Element beam - the delays are also read from the metafits). For full control, you can also specify all observational settings explicitly:
run_woden.py \
--ra0=50.67 --dec0=-37.2 \
--cat_filename=srclist_msclean_fornaxA_phase1+2.fits \
--primary_beam=MWA_FEE \
--MWA_FEE_delays=[6,4,2,0,8,6,4,2,10,8,6,4,12,10,8,6] \
--lowest_channel_freq=169.6e+6 \
--freq_res=10e+3 \
--num_time_steps=240 \
--time_res=0.5 \
--date=2018-02-28T08:47:06 \
--array_layout=MWA_phase2_extended.txt
WODEN
can read in an array layout as specified in local east, north, and height, and can simulate for any location on Earth (defaults to the MWA but you can specify Earth location with --longitude
/ --latitude
).
This sky model uses point sources, elliptical Gaussians, and shapelet models. All sources can either have power-law, curved power-law, or list of flux densities for Stokes I spectral behaviour. It can optionally also have a polarisation fraction, power-law, or curved power-law for linear (Q and U) and circular (V) polarisation.
A full breakdown of the sky model lives here on readthedocs, along with the possible formats to feed into WODEN
.