Here is a simple tutorial for exporting NEWSTAR
FITS
Step-by-step procedure: from e.g. ASTE telescope
to CLASS
1. The ASTE telescope will produce output data
in the raw NEWSTAR format.
You will probably get
a DAT tape with all NEWSTAR fits files.
e.g. rawdata file name:
n1068
2. start NEWSTAR
3. in NEWSTAR: read in the DAT tape (use task
FROM_CMT with parameter
"1-99", set the tapedrive
name, and execute FROM_CMT).
4. in NEWSTAR: use task LINTEG/SPFits to write out
each raw NEWSTAR file for
one filterbank at a
time. i.e. specify the input rawdata file from the
previous step and an
output file name, check the box "indiv scan", check
one filterbank (i.e.
A1, A2, A3, A4) at a time, then execute task.
e.g. you will now have
n1068a1.fits n1068a2.fits n1068a3.fits n1068a4.fits
5. Use procedure NSCLASS.PRO to read in the
NEWSTAR fits files produced
in the above step,
and output CLASS compatible fits files. For help
on this step, see below.
6a. Use the CLASS procedure file generated by NSCLASS.PRO
to read in
and consolidate the fits
files produced in the previous step.
e.g. you will now have
one output file, "n1068.bur", which contains
all scans originally
in the raw NEWSTAR file "n1068".
6b. Start IRAF (change loginuser.cl to include the line
set imtype
=fits
so that IRAF will work
directly with FITS files.
7. Proceed with processing in CLASS or IRAF.
EXAMPLE USAGE OF THIS PROCEDURE:
IDL> .rnew NSCLASS.PRO
IDL> NSCLASS,'myinputlist','CO(3-2)','myoutputlist','classread.class'
'myinputlist' : contains a list
of NEWSTAR fits files written with SpFITS.
one file per line. E.g. contents of file <myinputlist>:
n1068a1.fits
n1068a2.fits
'CO(3-2)'
: this is the name of the transition observed.
'myoutputlist': is an optional
file which contains file names for the output
fits files. Each input file corresponds to n output files, where n
is the number of scans in the input fits file.
E.g. if each input fits file consists of 10 scans, then
the file 'myoutputlist' should contain 20 lines, as follows:
n1068a1_0001.fits
n1068a1_0002.fits
...
n1068a1_0010.fits
n1068a2_0001.fits
n1068a2_0002.fits
...
n1068a2_0010.fits
There is no restriction in the name of the files, as long as you
have one file per scan.
The program will then output fits files named as above.
If file <myoutputlist> does not exist, or is empty, then default
extensions are used for the output fits files; these consist of the
input name followed by "_sc0001" etc. E.g. output files
n1068a1_sc0001.fits
n1068a1_sc0002.fits
...
n1068a1_sc0010.fits
n1068a2_sc0001.fits
n1068a2_sc0002.fits
...
n1068a2_sc0010.fits
'classread.class': this file will
contain a class procedure to read in
and consolidate all fits files produced by NSCLASS.PRO.
To run it, enter CLASS (part
of GILDAS) and at the prompt type:
LAS> @classread.class
e.g. you will now have the files n1068a1.bur and n1068a2.bur.
n1068a1.bur will have all scans originally in the raw NEWSTAR
file n1068a1, and the same applies to n1068a2.
Any comments
and questions to Juan Cortés, crivasqu@das.uchile.cl, or
Neil Nagar, nagar@astro-udec.cl.
|