ANALYZ instruction sequences HI in galaxies

This appendix contains example analyz sequences used to look at HI spectra.

This routine takes the average of the (gain/za corrected) dumps in each scan, and analysis is done on the individual correlator sections. Thus if you have the correlators overapping in frequency, and the galaxy therefore shows up in more than one correlator this routine will only look at one `detection' at a time. The commands to by typed are in boldface, followed by a description of the command. For help with these routines contact Karen O'Neil (koneil@naic.edu). All this assume you have run once analyz before. If you have not, please go to http://www.naic.edu/~astro/spectral_line/analyz and load hi_now.cmd, mjsr.o, and baseh.o into you analyz directories. Additionally, if you have not run analyz before, instruction (2) should read:

DESCRIPTION of the above commands:
 

• setvel: Uses the "optical" definition to convert frequencies to velocity, and then stores the velocities in register 41. The data is then "sectoned" for the rest of the analysis.
• baseh -5: This runs the baseline reduction routine, allowing you to choose which order fit you would like (from 1 to 5). The baseline is then subtracted from the data, and the results are stored in register 47.
• msrtest: This runs the msrj program, and outputs the folloqing results:

Line 1 - The area, central velocity, and width at 50% of the mean flux.
Line 2 - The area, central velocity, and width at 20% of the mean flux.
Line 3 - The area, central velocity, and width at 50% of the peak flux.
Line 4 - The area, central velocity, and width at 20% of the peak flux.
Line 5 - The area, central velocity, and width at 20% of each of the two horns,
searched from the inside out.
NOTE: Ignore the minus signs on the widths and areas - they come about because the data is flipped from what the msrj program expects (so the high velocity is subtracted from the low, instead of vice-versa.)

EXAMPLE using NEW routines (this data will be available for public use):

FIRST time:

        % cd analyz
        % analyz
        Enter your initials  [array size and number of STO/RCL registers optional]
        > temp 16284 800
        File ANZ_BASE/site/syslib.anz attached as SYSLIB.

        Enter <CR> or select an FLIB from the following:
         corv1    CORLIB for data prior to:07mar00
         corv2    CORLIB for data after   :06mar00
        FLIB> corv2
        File ANZ_BASE/site/corv2lib.anz attached as FLIB.

        >>> Creating a file for the stack and sto/rcl registers <<<
        Attaching temppanz.anz as STO/RCL file.
        Your arrays are 16384 elements long, and you have 800 sto/rcl registers.
        None of your sto/rcl registers are kept in main memory.
        Loading... PL X Y Z T HDR 
        File templib.anz attached as ULIB1
 
        Welcome to Analyz, version 5.16 patchlevel 2
 
 
        >~hi_n0w.cmd
        Error 144: ANALYZ fault: unrecognized command structure

        >exit

        STOP: ANALYZ exit
         Note: Following IEEE floating-point traps enabled; see ieee_handler(3M): 
         Overflow;  Division by Zero;  Invalid Operand; 
         Sun's implementation of IEEE arithmetic is discussed in 
         the Numerical Computation Guide.

        %analyz

        Enter your initials  [array size and number of STO/RCL registers optional]
        >temp

        Enter <CR> or select an FLIB from the following:
         corv1    CORLIB for data prior to:07mar00
         corv2    CORLIB for data after   :06mar00
        FLIB>corv2

        File ANZ_BASE/site/corv2lib.anz attached as FLIB.
 
        Attaching tempanz.anz as STO/RCL file.
        Your arrays are 16384 elements long, and you have 800 sto/rcl registers.
        None of your sto/rcl registers are kept in main memory.
        Loading... PL X Y Z T HDR 
        File templib.anz attached as ULIB1.
 
        Welcome to Analyz, version 5.16 patchlevel 2
 
 
        File /usr/local/phil/philib.anz attached as ULIB2.
        Created new array named @calv

        > attach "/proj/a1366/corfile.21mar00.a1366.1" data
        File /proj/a1366/corfile.21mar00.a1366.1 attached as DATA.
        Status = OLD; access = SEQUENTIAL; form = UNFORMATTED

        > redata 13 50 1.88
        loading via: ANZ_BASE/tools/ld -N -x -A /usr/local/bin/analyz.exe -T 037bc00 -o a.out
        ANZ_BASE/external/extern_main.o  ANZ_BASE/external/libextern.a -LANZ_BASE/external -lF77 -lm -lc
        New program text loaded successfully.


                .....

        Created new array named @ci
        on/off -1 raw in 100 -   149.0
        Created new array named @tmp2
        Created new array named @tmp3
        Created new array named @tmp4
        on/off -1 corr in   150-   199.0
        on/off -1 corrected average in   202.0
        CAL (in K) in reg nbsp; 200.0
        Calibrated data in register    203.0
        Temperature is    31.08
        on/off -1 corr, AB avg'd in   204.0        42.00
        rms subtracted data in 43

        > rcl 43;ver -0.02 0.01;plot 
 
        > hor 0 600;plot

        > getres 0 600

        loading via: ANZ_BASE/tools/ld -N -x -A /usr/local/bin/analyz.exe -T 0579400 -o a.out baseh.o
          -LANZ_BASE/external -lF77 -lm -lc
        New program text loaded successfully.

         Set base regions, with crosshair:
         Hit space key after each positioning,
         DO NOT touch mouse button!
         Hit Q after last (rightmost) limit.
         Limits are:   10  583  186  389  0  0  0  0  0  0
        Are limits OK? (y/n) [y] y
 
         Order     rms    Ftest
           1      0.001    27.985
           2      0.001     1.258
           3      0.001    11.416
           4      0.001    10.127
           5      0.001     2.408
        What order do you want?  2
        loading via: ANZ_BASE/tools/ld -N -x -A /usr/local/bin/analyz.exe -T 059e000 -o a.out msrj.o
          -LANZ_BASE/external -lF77 -lm -lc
        New program text loaded successfully.
        ow flag LEFT side chan for area bounds.
         Hit any key.
         Flag RIGHT side channel.
         Hit any key.
        Happy with the choice of bounds?  (y/n) [y] y
         area(230, 335) =-.3657E+01 Center = 9438.65884    Width(239, 335) = -532.81340
         area(230, 335) =-.3657E+01 Center = 9440.86550    Width(232, 335) = -569.33955
         area(230, 335) =-.3657E+01 Center = 9453.49903    Width(241, 329) = -487.90375
         area(230, 335) =-.3657E+01 Center = 9440.68466    Width(238, 335) = -537.20839
         area(230, 335) =-.3558E+01 Center = 9458.75213    Width(235, 333) = -528.72213
         area = -.3657E+01 Ok? (y/n or e (for exit))e
        MSR Done...
        > 

________________________________________________

> rcl 43;sto 500
>rcl 8;sto 501

>rcl 500;sto 43;rcl 501;sto 8

If you want to smooth your data and work with the smoothed version, just recall your data from 43, smooth it, ad store it back in 43. I.e.

>rcl 43;smo 19;sto 43

**************************************************************************************

---

Next:Trouble Shooting at the Telescope Up:Spectral Line Appendices Previous:Basic Analyz Routines -

Karen O'Neil