#CORinfodef;fcidisp 1;fcioff .01;fciflip 1;hor;ver:start
#doc start "startup script (run w/ autoexecute)"
#
attach "/usr/local/phil/philib.anz" ulib2:atphil
atphil;scratch 5;setxterm;corinfodef;fciflip 1:autoexecute
#scratch 5;setxterm:autoexecute
doc autoexecute "autoexecute script"
#
0;sto 8;0;sto 28;corfrqcfr pr1 pr2 pr3 pr4;corfrq 0;sto 8:setfrq
doc setfrq "setfrq <fr1> <fr2> <fr3> <fr4>, stores it in reg 8"
#
s;1420.4058;rcl 8;sto 40;-;rcl 8;/;* 299792.458;sto 41:helvel
doc helvel "rcls frq from 8, stores velocities in reg. 41"
#
#p2=.31079595;p3=.00744506;p4=-.00174542;p5=.1498039e-3;p6=-.44974e-5:GZACOR%0
#gzacor%0;eval4 pr1:GZACOR
#
p2=0.278;p3=-2.473e-3;p4=6.11e-4;p5=-3.4e-5:tgzclw%0
tgzclw%0;eval4x:tgzclw
doc tgzclw "LW -compute T/G(za) for za in xreg. return in xreg"
#
s;p2=pr2*2 + 2;storec 100 pr1 p2:getdat%1
p1=100;p2=100+pr2;loop pr2;rcl p1;rcl p2;/;- 1;sto p1;p1=p1+1;p2=p2+1;end:gtdt%2
p7=100;loop pr2;rcl p7;rms;- pl(1);sto p7;p7=p7+1;end:getdat%2a
getdat%2a pr1 pr2;p8=p7-1;print "on/off -1 raw in 100 -",p8:getdat%2b
#
getdat%2b pr1 pr2;p2=102+2*pr2;corpwrl p2 pr1 pr2:getdat%3
rcl p1;@tmp2=x;gzacor;@tmp3=x;x=@tmp2;tzacor;@tmp4=x/@tmp3:getdat%3b
rcl p1;@tmp2=x;tgzclw;@tmp4=1 / x:gtdtlw%3b
i=1;p1=100;loop pr2;getdat%4b pr1 pr2;end;getdat%4c pr1 pr2:getdat%4
rcl p1;*@tmp4(i);p7=p1+pr2;sto p7;i=i+1;p1=p1+1:getdat%4b
p8=100+pr2;print "on/off -1 corr in",p8,"-",p7:getdat%4c
#
p1=100+pr2;x=0;loop pr2;rcl p1;+;p1=p1+1;end;/ pr2;p7=100+2*pr2;sto p7:gtdt%5
gtdt%5 pr1 pr2;print "on/off -1 corrected average in",p7:gtdt%5c
#
getdat%1 pr1 pr2;gtdt%2 pr1 pr2;getdat%3 pr1 pr2;getdat%3b pr1 pr2:getdat%6a
getdat%1 pr1 pr2;gtdt%2 pr1 pr2;getdat%3 pr1 pr2;gtdtlw%3b pr1 pr2:gtdtlw%6a
getdat%6a pr1 pr2;getdat%4 pr1 pr2;gtdt%5c pr1 pr2:getdat%6b
gtdtlw%6a pr1 pr2;getdat%4 pr1 pr2;gtdt%5c pr1 pr2:gtdtlw%6b
#
p6=p7;p7=p6+1;0;sto p7;p5=4*hdr(39);p4=2*p5:stodat%1a
rcl p6;shiftl p5;sto p7;rcl p6;rcl p7;+;/ 2;sto p7:stodat%1b
rcl p7;s p5 p4;0;sto p7;s;rcl p7;sto 42:stodat%1c
stodat%1a;stodat%1b;stodat%1c:stodat%1
print "on/off -1 corr, AB avg'd in",p7,42:stodat%3
#
#getdat%6b pr1 pr2;stodat%1 pr2;stodat%2;stodat%3:getdat
getdat%6b pr1 pr2;stodat%1 pr2;stodat%3:getdatln
doc getdatln "getdat(LN) <scan> <dumps>. uses registers 100 - 100+pr2*2 + p1"
#
gtdtlw%6b pr1 pr2;stodat%1 pr2;stodat%3:getdatlw
doc getdatlw "getdat(LW) <scan> <dumps>. uses registers 100 - 100+pr2*2 + p1"
#
rcl 42;rms;- pl(1);sto 43;rcl 43:mbase%1
mbase%1;print "rms subtracted data in 43":mbase
#
correw;getdatln pr1 pr2;s;mbase;rcl 43;p8=4*hdr(39);hor 1 p8;plot:redln
doc redln "data reduction (LN) - datred <scan> <dumps>"
#
correw;getdatlw pr1 pr2;p8=4*hdr(39);s 1 p8;mbase;rcl 43;hor 1 p8;plot:redlw
doc redlw "data reduction (LW) - datred <scan> <dumps>"
#
helvel;s pr1 pr2;sto 30;rcl 43;hor;ver;plot:setvel
doc setvel "puts vel. in 41, sections data from pr1 - pr2"
#
$baseh -5; -;sto 47;rcl 47;plot;plot;msrtest:basgauss
doc basgauss "Baseline red;Det. the central velocity, flux, and widths ofthe data"
#
pl(1,5) .5 .2 .5 .2 .2;$msrj 11223:msrtest
#
setvel pr1 pr2;basgauss:getres
doc getres "sections data, gets widths, etc - getres <sec1> <sec2>"
#
attach "/proj/a1366/corfile.18feb00.a1366.2" data:atdat
doc atdat "delete me"