Importing Data from MIRIAD: Difference between revisions

From CASA Guides
Jump to navigationJump to search
No edit summary
No edit summary
Line 3: Line 3:
A filler routine is currently being developed to convert native MIRIAD format data into a CASA measurement set, but it has not yet been completed.  The following steps describe how to use FITS format files to transfer CARMA data into CASA, and should still be useful even after the filler is operational.  At that time, a description of how to use the CARMA filler routine will be added to this page.
A filler routine is currently being developed to convert native MIRIAD format data into a CASA measurement set, but it has not yet been completed.  The following steps describe how to use FITS format files to transfer CARMA data into CASA, and should still be useful even after the filler is operational.  At that time, a description of how to use the CARMA filler routine will be added to this page.


This page assumes you have MIRIAD already installed.  All MIRIAD commands are given in their shell-script (csh) form, assuming that the following have been defined at the top  
This page assumes you have MIRIAD already installed.  All MIRIAD commands are given in their shell-script (csh) form, assuming that a number of variables have been defined at the top of your script.
 
For this example, the following items would be defined at the beginning of the csh shell script:
 
<source lang="bash">
# in csh
 
#!/bin/csh
# FILE = name of original MIRIAD-format file
set FILE="../c0437.1E_115NGC694.1.miriad"
# PROJECT = base file name for processed data files
set PROJECT="c0437"
# ANTPOS = the location of the file which contains updates to the CARMA antenna positions
set ANTPOS="../antpos.090613"
</source>


== Preparing MIRIAD data for export to FITS ==
== Preparing MIRIAD data for export to FITS ==
To minimize the data you are transferring over, select only the astronomically useful data:
<source lang="bash">
# in csh
uvcat vis=$FILE  select='-source(NOISE),-auto' out=$PROJECT.vis
</source>


If you need to apply an update for the antenna positions, or linelength corrections, you should do these before exporting to FITS.
If you need to apply an update for the antenna positions, or linelength corrections, you should do these before exporting to FITS.
Line 13: Line 34:
<source lang="bash">
<source lang="bash">
# in csh
# in csh
uvedit vis=$PROJECT.rest out=$PROJECT.base apfile=$ANTPOS
uvedit vis=$PROJECT out=$PROJECT.ant apfile=$ANTPOS
</source>
</source>


Line 20: Line 41:
<source lang="bash">
<source lang="bash">
# in csh
# in csh
linecal vis=$PROJECT.base
linecal vis=$PROJECT.ant
uvcat vis=$PROJECT.base out=$PROJECT.lc
uvcat vis=$PROJECT.ant out=$PROJECT.lc
</source>
</source>
Currently, MIRIAD writes out the identity of the array as identified by Antenna 1, which is a formerly-OVRO antenna.  Thus, without correction, the dataset will identify the array as OVRO rather than CARMA.  To fix this problem, simply do a small header edit:
#####
# Write out windows of individual sidebands to export to CASA for
# testing--all sources in same file for common window
#####
#so CASA knows the correct telescope (not "OVRO")
puthd in=$PROJECT.lc/telescop value=CARMA
== Writing out the FITS files ==
foreach i (4 5 6)
    fits in=$PROJECT.lc op=uvout select="win($i)" out="./fits/$PROJECT.$i.fits"
end





Revision as of 15:33, 23 November 2009


A filler routine is currently being developed to convert native MIRIAD format data into a CASA measurement set, but it has not yet been completed. The following steps describe how to use FITS format files to transfer CARMA data into CASA, and should still be useful even after the filler is operational. At that time, a description of how to use the CARMA filler routine will be added to this page.

This page assumes you have MIRIAD already installed. All MIRIAD commands are given in their shell-script (csh) form, assuming that a number of variables have been defined at the top of your script.

For this example, the following items would be defined at the beginning of the csh shell script:

# in csh

#!/bin/csh
# FILE = name of original MIRIAD-format file
set FILE="../c0437.1E_115NGC694.1.miriad"
# PROJECT = base file name for processed data files
set PROJECT="c0437"
# ANTPOS = the location of the file which contains updates to the CARMA antenna positions
set ANTPOS="../antpos.090613"

Preparing MIRIAD data for export to FITS

To minimize the data you are transferring over, select only the astronomically useful data:

# in csh
uvcat vis=$FILE  select='-source(NOISE),-auto' out=$PROJECT.vis

If you need to apply an update for the antenna positions, or linelength corrections, you should do these before exporting to FITS.

To update antenna positions:

# in csh
uvedit vis=$PROJECT out=$PROJECT.ant apfile=$ANTPOS

To apply linelength corrections:

# in csh
linecal vis=$PROJECT.ant
uvcat vis=$PROJECT.ant out=$PROJECT.lc

Currently, MIRIAD writes out the identity of the array as identified by Antenna 1, which is a formerly-OVRO antenna. Thus, without correction, the dataset will identify the array as OVRO rather than CARMA. To fix this problem, simply do a small header edit:

  1. Write out windows of individual sidebands to export to CASA for
  2. testing--all sources in same file for common window
  1. so CASA knows the correct telescope (not "OVRO")

puthd in=$PROJECT.lc/telescop value=CARMA

Writing out the FITS files

foreach i (4 5 6)

   fits in=$PROJECT.lc op=uvout select="win($i)" out="./fits/$PROJECT.$i.fits"

end




--Mthornle 19:05, 23 November 2009 (UTC)