Joint Astronomy Centre
Show document only
JAC Home
JCMT
UKIRT
Contact info
JAC Divisions
OMP
Outreach
Seminars
Staff-only Wiki
Weather
Web Cameras
____________________

JCMT home
Observing at JCMT
OMP Observation Manager
Telescope
Spectral Line Observing
Continuum Observing
Schedule
Data Archive
Future Developments
Legacy Surveys
Newsletter & Publications
How to reduce DAS data with CLASS

Reducing DAS data with CLASS

When it seems (at the JCMT) that the GAG software (cfits, class, greg, graphic, gfits, astro) is not there,
you should type first:

setenv PATH /jac_sw/gag/linux-intel/bin:$PATH

The latest installed version (December 2005) can be reached via:

source /local/progs/gildas/gildas-exe/pc-redhatrelease-g95/etc/login

If you prefer to use CLASS and other Gildas software to reduce your spectra these are the steps to take:
Unfortunately you have to use Specx to write the spectra to disk as fits files.

For single spectrum #131:

>> o-fil (open output file)
>> r-g-d 131
>> das-merge
(other commands if needed, e.g. fsw 8.3)
>> wr-sp 1 (to output file 1)

For raster/grid maps in scan #132:

>> o-fil (open output file)
>> r-g-r 132 (you will be asked for the output file and if you want to do das-merge, answer Y)

>> tofits (writes the data in a specified output file to disk as fits files).

Alternatively you can use g2f (gsd2fits). It writes a range of scan numbers to fits files doing a das-merge with default parameters (see the warning shown when typing the command). It usually works well for single channel spectra, and RxB dual channel spectra at the same frequency (samples and patterns).

For more complicated cases (frequency switching, 2 channels) some information is given here

The fits files will have these names:

jcmt_0131_001_a.fits
jcmt_0132_001_a.fits
jcmt_0132_002_a.fits
...
jcmt_0132_015_a.fits
(there will also be files jcmt_0131_001_b.fits if you observe with two mixers)

These files can be converted into class format to file jcmtdata with Cfits (in the newest version use class):

%cfits
cfits>file out jcmtdata new (add new if the file does not yet exist)
cfits>read jcmt_0131_001a.fits
(or alternatively)
cfits>read "/dir/subdir/jcmt_0131_001a.fits"
cfits>las\write

In the example above all spectra 132_001a, 132_002a, etc will get the number 132, so it might be of use to modify the numbers later in class.
If you have many large raster maps it is useful to do the conversion to class format with a commandfile, e.g.
%ls *.fits > files.cfits
and then edit files.cfits with your favourite editor from
jcmt_0131_001_a.fits
to
read jcmt_0131_001_a.fits
las\write
etc
then
cfits>@files.cfits

%class
to start class

For help within cfits, class, etc type help, or look at the GILDAS homepage.

In the newest version of Greg/Graphic the Explain command is not working. Here is a list of available tasks (run task).

Jan Wouterloot

Contact: Per Friberg. Updated: Tue Apr 10 14:39:32 HST 2007

Return to top ^