Notes for reducting galfa data (Oct 27 06)

GENERAL:

- login to arecibo as user "madison" 
  using ssh -X remoto.naic.edu -l madison
  ssh to aolc1 (2,3, or 4), 
- Our working directories are in /share/galfa/madison (Bill and Samantha)
- Fits files are in /share/galfa
- Typical fits file name is:
  galfa.20050530.a2032.0002.fits 
  galfa.YearMonthDay.Project.Filenumber.fits

- Each fits file has its mh file.
  All mh files are in /share/galfa/galfamh

Each data (fits) file contains 10mins of data.
On each day we have several fits files -- this is one scan of data.
During the archiving process for each file we generate the so called
"mh" file which contains properly calculated header information,

PART I. 
Make a text file that lists all files per day and check that 
all fits files have appropriate mh files.

*** Write logs of everything you do! ***

Make a large table:

Date	 File #  MH		Region		Scan#	Start LST	End LST  	
20050529   0000  x		ms              13      22.607272
           0001  x 
.......
           0010  x              ms              13                      0.445346


Getting slst and elst:
- slst is start LST (Local Siderial Time) for that scan (=day)
In IDL open the first mh file for that scan, with:
restore, 'galfa.20050529.A2011.0000.mh.sav',/ver 
help, mh,/st
   and read
LST_MEANSTAMP   DOUBLE           23.605123
therefore slst=23.605123
or print,mh.lst_meanstamp --- take the first, smallest number
or print, minmax(mh.lst_meanstamp)
The smallest LST for that scan is your slst.

"elst" is the end LST for that scan. To find this value,
do the same by opening the last mh file for that scan and printing 
the latest LST. 

Also, by doing this you can get your "scan" 
by printing mh.object:
OBJECT          STRING    'ms_00_00   
This means your scan number is 00 (middle number). 


PART II. 
Run make_dirs.pro (see "The GALSPECT pipeline" document) or
generate appropriate directories manually.
This has to be done only once.


PART III. 
Run 1st stage of reduction: stg0.pro; need:

stg0, year, month, day, proj, region, root, startn, endn, slst,$
elst, scan, /nomh, mhdir=mdir, caldir=cldir,$
fitsdir=fdir

"The GALSPECT pipeline" document explains all these parameters.
See scripts: /share/galfa/madison/Samantha/ss_a2032_simple.idl 
and /share/galfa/madison/Bill/ss_a2034_simple.idl .
Use info from your log files, e.g. "slst" and "elst", scan.
 

IV. Check the data files after running stg0.pro
Output files are: *.sav and *.fits

- Some diagnostics:
Open a few fits files in IDL with:
data=readfits('galfa.20050516.a2032.0007.ms.fits')

DATA         FLOAT     = Array[8192, 2, 7, 600]
read a few data files and print minmax, note down if there are some huge 
numbers (like 1000 or more)
Think of other useful diagnostics as you go.

PART V: Write an IDL code to plot positions (RA and Dec) for all 
processed data files.

--------------------------------------------------------------------

PROJECT SPECIFIC INFORMATION:

PROJECT a2032:

Observing mode: basket-weave scans.
There is only one object "ms".
There were 17 scans in total, from 00 to 16. 

PROJECT a2034:

Observing mode: drift scans.
There are several objects, need to be careful in selecting them.
objects: G44p8, 36p2, G44p8, G30p5, G32, G35p6, G37p5.
Check Erik Muller's log info.
Each day's data will make one scan for a given object.
So, for each object will have a bunch of scans, starting with
0000 to 00XX.
For each scan must check objects and then select appropriate LST 
range per object.

Probably the best to start with: date 050504, worked!

To find lst where object changed can use:
IDL> help,mh[0].object
<Expression>    STRING    = 'G44p8_1a355     '
IDL> m=where(mh.object eq 'G44p8_1a355     ')
IDL> print,minmax(mh[m].lst_meanstamp) 


=========================================================

Note on getting display with VNC:

1st xterm: ssh -X remoto.naic.edu, then ssh fusion02, then vncserver -geometry
1280x1024
2nd xterm:
vncviewer -via sstanimi@remoto.naic.edu fusion02:1