|
Документ взят из кэша поисковой машины. Адрес
оригинального документа
: http://www.naic.edu/~astro/spectral_line/handbook_old/app/old_analyz.html
Дата изменения: Wed Sep 20 02:21:10 2000 Дата индексирования: Sun Dec 23 05:24:20 2007 Кодировка: Поисковые слова: interstellar medium |
Note: There are other libraries available. However, all libraries function very closely to the below libraries. Additionally, these libraries are used in the HI analyz routines below. When information on the other libraries becomes available, it will be placed in the link at http://www.naic.edu/ ~ csalter/toppage.html
correlator library help summary.
this file is located (at AO) in /usr/local/src/Analyz/help/cat/other/corlibv1
expanded help available for SEE ALSO functions below the
expanded documentation is in /usr/local/src/Analyz/help/cat/other
last modified: 28dec99
FUNCTION
atoldata ..attach default online datafile
coravgp reg ..average two polarizations.
corazzal regOut scan numrecs... sto az regout,za regout+1
corclose ..close monitoring socket.
corcmbab rega regb .. move pola followed by polb to x-reg
cordisp .. display spectra in corget current display mode.
corfrqcfr cfrsbc1 cfrsbc2 cfrsbc3 cfrsbc4.. set center frequency each sbc
corfrq 0,1 .. stofreq array in x. 0-->a,1-->polb. use @frqX arrays
corget parm..input scan.grp to reg1-4. pola,b to reg 5,6.
corhdr .. print out cor header from last corget.
corhor sbcnum.. set horizontal scale for sbc n (0=all)
corinfodef ..set default values for corinfo
corld reg pol(0,1) {unp:1}..spectra->x,from reg.
corldf reg pol(0,1) {unp:1}..spectra->y,freq->x.from reg.
corlist .. list summary of current correlator file attached as data
corpl {reg}{sbc:1-4,0=all}..plt vs frq 2pol
corpl. {reg}{sbc:1-4,0=all}..plt. vs frq 2pol
corpla {reg}{sbc:1-4,0=all}..plt vs frq polA
corpla. {reg}{sbc:1-4,0=all}..plt. vs frq polA
corplb {reg}{sbc:1-4,0=all}..plt vs frq polB
corplb. {reg}{sbc:1-4,0=all}..plt. vs frq polB
corplit sbcNum(0=all) pol(=0a,1=b)(flag 0/1 donot reset cipl%).corpl uses.
corpos scan.rec ..position to scan.rec (rec is 5 digits)
corposnr scan.rec ..position to scan.rec (rec is 5 digits) no rewind.
corprint .. print hdr info to stdout
corpw regstart scan numrecs.. store polA,B,az,za to 8 registers.
corpwr ..move power info to x . 4 numbers each board
corpwrl regOut scan numrecs.. sto pwr info,az,za
correw ..reset pointer to start of current file
corsec sbcnum pol(0,1)..sec section for sbc 1-4, 0 all A or B,-1 all A+B
corseek recnum.. seek to record num..no input
corsmooth ..called before display.. data in x-reg
corsplab reg1 .. split xpolA->reg1, xpolB->reg1+1
corsum scan nrecs.. compute summary rec for scan. sto in 5,6
doit scan recs ..
domany regSrc regCal scanSt numSets recSrcOn..doset numSets times
doset regSrc regCal scanSt recSrcOn .on/off src then on/off cal
dosrc reg scan npairs rec/scan... load avg
gzacor ..compute gain(za) for za in xreg. return in xreg
gzacor1 za .. compute gain(za) for input za. return in pl(1)
medfonoff regStore scan numrecs.. medianfilter then on/off-1
medfscan regStore scan numrecs.. median filter recs of scan
ndoit scan numpairs numrecs ..
ndoitsto regSt scanSt numPairs numrecs..sto (on-off)/off
npair numpairs numrecs.. accumulate n pairs
pair numrecs.. process 1 pair
pairpl {sbc:1-4,0=all}..load,plt vs frq last onoff
pairpl. {sbc:1-4,0=all}..load,plt. vs frq last onoff
pairpla {sbc:1-4,0=all}..load,plt vs frq polA last onoff
pairpla. {sbc:1-4,0=all}..load,plt. vs frq polA last onoff
pairplb {sbc:1-4,0=all}..load,plt vs frq polB last onoff
pairplb. {sbc:1-4,0=all}..load,plt. vs frq polB last onoff
stokesc reg..load @i,@q,@u,@v,cmp delta,tau. plot delta. circdata
stokesl reg..load @i,@q,@u,@v,cmp delta,tau. plot delta. lin data
storec regSt scn nrecs .. sto bandpasses
storecnr regSt scn nrecs .. sto bandpasses
strippr reg1 reg2 offset step.. strip plt vs freq.
tzacor ..compute za Tsys correction. input za in x. return in x
tzacor1 za .. compute Tsys za correction for za ,return in p1
DESCRIPTION
----------------
register usage:
----------------
1 - band1 both pol raw data
2 - band2 both pol raw data
3 - band3 both pol raw data
4 - band4 both pol raw data
5 - flipped,and scaled (polA) for all brds concantenated
6 - flipped,and scaled (polB) for all brds concantenated
7 - default freq array for corget,cordisp displays
8 -
9 - scratch reg for corpwrl,corsum (9,10)
10-27 scratch for tsyscon
addon,addof for a single pair
20 - addon A
21 - addon B
22 - addoff A
23 - addoff B
24 - (on-off)/off A
25 - (on-off)/off B
28,29 - scratch for dosrcplx
--------------------------------------
variables used:
--------------------------------------
regS - start reg set by dosrc, used by dosrcplx
regE - end reg for storing pair on/off set by dosrc, used by dosrcplx
regA - reg for average all pairs. set by dosrc, used by dosrcplxx
stroff - offset for first strip. set by dosrcsclx, used by dosrcpl1,2a,2b
strstp - step for each strip. set by dosrcsclx, used by dosrcpl1,2a,2b
fcioff - function to change offset of polB when plotting. used by
corget,cordisp,dosrcpl3,pairpl,pairpl1
--------------------
basic functions
--------------------
atoldata .. attach online datafile default datafile
corinfo .. display current state
corget scan.rrrrr
.. read in scan.rrrrrr (default next rec) and load data:
reg 1-4 raw data, unflipped, unscaled
reg 5-6 flipped, scaled,and smoothed data.
5 1st sbc of each board (polA) concatenated
6 2nd sbc of each board (polB) concatenated (if used)
If fcidisp is not zero, cordisp will be called to display
the data.
ex: corget .. get next record
corget .00032 .. get record 32 of current scan
corget 123. .. get 1st rec scan 123
corget 123.00003 .. get record 3 of scan 123
note: use correw if you have positioned beyond the requested scan.
scannumber 815600012 can be abbreviated 12 until you
hit a day change (8157xxxxxx).
corsum scan.rrrrr
.. input and compute the summary record for the scan.
Store the summary record in registers 5 and 6.
The azimuth and zenith angle for each record is stored in
the scratch arrays @az, @za.
Use cordisp or corpl to disp the data.
ex: corsum .. compute summary at current position
ex: corsum 123 .. compute summary for scan 123
note: you need to do a correw if you are positioned beyond the
requested scan.
storec reg scan nrecs
.. input nrec records starting at scan. store each record
(bandpasses) starting at register reg. PolA and PolB are
concatenated into a single register. This will happily
cross scan boundaries. Use corpl1 reg to display the data.
^L
-----------------
general plotting
-----------------
With data stored in a register pair (say reg 5=polA, reg 6=polB) use
cordisp, corpl. When data is stored with all data in 1 register
(pol A followed by polB) use corpl1 to display the data.
cordisp .. display the data reg 1-4 (ficdisp 1) or reg 5-6
(fcidisp 2). Use plotxy with freq reg set by fcifreq reg.
polB is offset by the value from fcioff. You load whatever
you want in the freq reg.
corpl {reg} {sbcnum}
corpla {reg {sbcnum}optional reg}
corplb {reg} {sbcnum
display the spectra vs frequency. If sbcnum is given, then
just plot that sub correlator. If A/B is used then just plot
polA. reg can be:
none - default to 5,6 .. registers loaded by corget.
< 30 - use reg and reg+1 for pola,polB
>=30.. register has polA,polB packed in 1 register
(see storec, doset,domany, etc..)
see help corpl
ex: corpl .. plots 5 (polA), 6 (polB) vs corfrq
ex: corpl 24 .. plots 24 (polA), 25 (polB) vs corfrq
SEE ALSO help corpl
corpl. {optional reg}
same as corpl but uses plotxy. for plotting.
multiple calls within a function will overplot.
ex: corpl. .. plots 5 (polA), 6 (polB) vs corfrq
ex: corpl. 24 .. plots 24 (polA), 25 (polB) vs corfrq
corpos scan
.. position to start of scan. A rewind is automatically done
before the positioning.
ex: corpos 123 .. position to scan 123.
corposnr scan
.. position to start of scan with no rewind. Scan must be
after current position in file (if not use corpos). For
large files doing seqential positioning, this routine is
than corpos.
ex: corposnr 123 .. position to scan 123.
corseek recnum
.. seek to record recnum of current scan. The record is not input.
correw .. rewind to start of file
corlist .. list current file attached as data
corprint .. print out hdr info of last record input
corhdr .. print entire header of last record input
corpwr .. extract power info from header for current record. 16 numbers
in sbc order are stored in x reg:
sbc1pcA sbc1pcB sbc1lag0A sbc1lag0B ..
sbc2pcA sbc2pcB sbc2lag0A sbc2lag0B .. etc
corpwrl reg1 scan numRecs
.. read numRecs records starting at scan. extract the
power information,az,za and store in registers.
The nth entry in a register is from the nth record input.
register 9 is used as a scratch register.
reg1 : lag0A sbc1
reg1+ 1: lag0A sbc2
reg1+ 2: lag0A sbc3
reg1+ 3: lag0A sbc4
reg1+ 4: lag0B sbc1
reg1+ 5: lag0B sbc2
reg1+ 6: lag0B sbc3
reg1+ 7: lag0B sbc4
reg1+ 8: pwrCA sbc1
reg1+ 9: pwrCA sbc2
reg1+10: pwrCA sbc3
reg1+11: pwrCA sbc4
reg1+12: pwrCB sbc1
reg1+13: pwrCB sbc2
reg1+14: pwrCB sbc3
reg1+15: pwrCB sbc4
reg1+16: az position
reg1+17: za position dome
where reg(n) is the nth record from the start.
----------------------------------
utility routines
----------------------------------
coravgp regnum
average polA,B. If regnum <=30 then pola is in regnum and
polB is in regnum+1. If regnum > 30 then polA,polB are packed
in regnum. The averaged funtions are left in the x-reg.
SEE ALSO help coravgp
corfrqcfr cfrsbc1 cfrsbc2 cfrsbc3 cfrsbc4
Define the center frequency for the 4 subbands.
Used by corfrq, corpl...
ex. if cfr 1370 1390 1410 1420:
corfrqcfr 1370 1390 1410 1420:
if 1 pol/sub correlator.. 1370A 1370B 1390A 1390B then:
corfrqcfr 1370 1390
corfrqcfra/b cfrsb1 cfrsb2 cfrsb3 cfrsb4
define center frequency for polA or B sbc. Same as corfrqcfr
but allows polA and polB to not cover the same frequency range.
corfrq {0,1}
load the x register with the frequency valuse for polA or polB.
(0== polA, 1==polB). Used by corpl,pairpl, etc..
ex. corfrq 0 .. loads x-reg with polA freq.
corext sbcN {reg}
extract sbc N (N=1..4) from reg 5,6 or reg,reg+1r
x=polA,y=polB,z=freqA.
corhor {1..4}
set horizontal freqendy scale for sbc n =1 to 4.
SEE ALSO help corhor
corld {reg} {polA,B 0,1} {unpacked:1/0}
load spectra into x-reg. used mainly by other functions.
SEE ALSO help corld
corsbci sbcnum 0/1
return indices in reg 5,6 for sbcnum. 0==polA, 1==polB. return
info in pl array.
ex. corsbci 2 0
pl(1 2) index start,end in reg
pl(3 4) frequency start/end
pl(5) length this sbc.
corsec {sbc} {pola,b: 0,1}
set section size to display sbc.
SEE ALSO help corsec
corcmbab reg1 reg2
assume reg1,reg2 hold pola,polB. combine these in the x-reg
so you have polA followed by polB.
corsplab reg1
assume x-reg has polA followed by polb. Put polA in reg1
and polB in reg1+1.
medfscan regStore scan numrecs
median filter by channel the records in scan. There are
numrecs records in scan. Leave the median filter spectra in
x-reg. Use regStore to regStore+numrecs-1 as scratch.
medfscan 101 123 50 ..
medfonoff regStore scan numrecs
median filter by channel scan and scan+1 and then compute
on/off. Leave data in regs:
reg 20 - median filtered on position. polA followed by polB
reg 21 - median filtered on position. polB
reg 22 - median filtered off position. polA followed by polB
reg 23 - median filtered off position. polB
reg 24 - on/off-1 .
medfonoff2 regStore scan numrecs
compute on/off-1 by rec and then median filter.
used regStore to regStore +2*numrecs-1 as scratch
Leave data in regs:
reg 23 - median filtered off position. polB
reg 24 - on/off-1 .
^L
-----------------------------------
on/off position switching routines
-----------------------------------
general routines:
doit scan recs
.. position, do a single pair, then plot the results:
corpos scan; pair recs; pairpl
ex: doit 123 2 .. do a on/off pair starting at scan 123 with
2 recs in on and 2 in off.
data is left in registers 24 polA, 25 polB
pairpl {reg}
.. plot on/off-1 for the data in registers reg (polA) and
reg+1 (polB). offset polB by value from fcioff. Default
register is 24,25. This routine is the same as corpl (
except that the default register is different).
ex: pairpl .. plot on/off-1 from register 24,25 vs 8
pairpl 28 .. plot on/off-1 from register 28,29 vs 8
dosrc regSt scan numpairs recs/scan
.. process numpairs on/off pairs starting at scan. Store the
individual on/off-1 starting at regSt (polA followed by polB
in one register). Compute the average for all the on/off-1
and store in register regSt+numpairs. Set the variables
regS:start reg, regE:endReg regA:average reg. These
are used by the dosrcplXX routines.
ex: dosrc 31 123 10 4 .. do 10 pairs starting at scan 123.
store each on/off-1 starting at reg 31. There are 4 recs
in an on or off. The average will be stored in
reg 31+10=41. Plot the average with corpl1 41.
domany regSrc regCal scan numsets recsSrcon
.. call doset to process an on/off src starting at scan, then
process a cal on/off pair. repeat this for numSets times.
This routine should be used to process data taken with the
coronoffnewXX datataking routine.
register usage:
regSrc,regSrc+numSets-1.. the on/off-1 for the src scans
regSrc+numSets .. the avg for the src on/off-1
regCal,regCal+numSets-1.. the on/off-1 for the cal scans
regCal+numSets .. the avg for the Cal on/off-1
This routine assumes:
1. that there are recsSrcOn records for the source on or
off in each scan.
2. that the cal on /off each contain a single record.
3. that the scan numbers are contiguous.
4. the src average will be meaningfull if it was the
same src.
ex: domany 101 151 123 10 5
Process 10 sets of on/off, cal pairs.
The scan for the first on is 123
Store the src data in reg 101-110, avg in 111
Store the cal data in reg 151-160, avg in 161
there are 5 records in each src scan.
corpl1 111 will look at the src average.
^L
routines used by the above:
clear
.. clear on/off accum registers (20-23)
addon numrecs
.. add numrecs (default 1) to the on registers (20,21)
addoff numrecs
.. add numrecs (default 1) to the off registers (20,21)
pair numrecs
.. Do an on/off pair from the current disc position:
clear;addon numrecs;addoff numrecs; compute on/off -1
20,21 have the sum of ons, 22,23 have the sum of offs
24,25 have on/off-1 (polA,polB)
ex: pair 4 .. do on/off pairs with 4 records per scan
npair numpairs numrecs
.. accumulate numpairs pairs startin at the current position,
then compute on/off-1. The results are stored in reg 20-25
ex: npair 4 4 .. process 4 pairs. each pair has 4 recs in on,off.
ndoit scan numpairs numrecs
.. position to scan, accumulate numpairs pairs, then plot
the results. The data is stored in reg 20-25.The code is:
corpos;npair numpairs numrecs;pairpl
ex: ndoit 123 5 2 ... position to scan 123 then accumulate 5 pairs
each with 2 recs per on, off.
ndoitsto reg scan numpairs numrecs
.. process numpairs pairs. Store the on/off-1 starting at
register reg (polA followed by polB in 1 reg).
each on/off pair should have numrecs records in an on or off.
ex: ndoitsot 31 123 10 3 .. process 10 pairs starting at scan
123. Each pair has 3 recs in the on or off. Store the
individual on/off-1 starting at reg 31.
online:
corclose . close socket for monitoring
^L
-----------------------------------
misc functions
-----------------------------------
gzacor1 za ..compute the telescope gain at zenith angle=za. Return value
in pl(1). Units are Tsys(za=3deg).
gzacor ..compute the telescope gain fo zenth angles
passed in via the x-register. Return the gain the in x-reg.
Units are Tsys(za=3deg).
tzacor1 za ..compute the system temperature correction for za.
return value in pl(1). Normalized so Tsys(za=3) is unity.
tzacor ..compute the system temperature correction for zenth angles
passed in the x-register. Return the corrections in x-reg.
Normalized so Tsys(za=3) is unity.
********************************************************************************