#
# glish script to plot rms results for standard spectra checkout
#
# R. Prestage 9th September 2002
#
#
# set up data arrays. These are for multiples of 300 seconds, plus 
# additional multiples of two for the averages of the two receivers
# stagger the X values slightly so points don't overlap
#

# y1 = receiver X, y2 = receiver Y, y3 = average

#        1    2    3    4    5    6    7    8    9   10   11   12   13   14
x1 :=[ 290, 290, 290, 290, 290, 290, 290, 570, 570, 570, 870, 870,1160,2000]
x2 :=[ 310, 310, 310, 310, 310, 310, 310, 630, 630, 630, 930, 930,1240,2200]
x3 :=[ 600, 600, 600, 600, 600, 600, 600,1200,1200,1200,1800,1800,2400,4200]

y1 :=[27.0,30.0,26.7,28.4,26.8,27.3,30.0,20.5,19.6,19.2,16.0,16.2,14.1,10.7]
y2 :=[27.4,29.6,26.4,28.3,26.4,26.8,33.0,20.4,19.8,19.0,16.1,15.9,13.8,11.2]
y3 :=[19.2,22.6,19.3,20.5,19.0,19.2,22.1,15.1,14.3,13.7,11.6,11.7,10.2,7.94] 

# convert to logs

x1 := log(x1); x2 := log(x2); x3 := log(x3); 
y1 := log(y1); y2 := log(y2); y3 := log(y3); 

xline := log([200,5500])
yline  := log( [(28.8 * sqrt(300/200) ), (28.8 * sqrt(300/5500) ) ] )

# make the plot

pg.clear()
pg.pap(7.5, 0.8)
pg.sch(1.5)
pg.env(2.3,3.7,0.8,1.6,0,0)

xlbl := "LOG10(Integration time in seconds)"
ylbl := "LOG10(Observed RMS in mK)"
toplbl := "TSPECTEST_02: RMS in 12.5MHz scans"
pg.lab(xlbl,ylbl,toplbl)

pg.sch(0.8)
pg.pt(x1,y1,4)
pg.pt(x2,y2,-9)
pg.sch(1.0)
pg.pt(x3,y3,3)
pg.line(xline, yline)

pg.pt(3.1,1.49,4)
pg.pt(3.1,1.44,-9)
pg.sch(1.0)
pg.pt(3.1,1.39,3)

pg.sch(1.5)
pg.text(3.15,1.48,'Receiver X')
pg.text(3.15,1.43,'Receiver Y')
pg.text(3.15,1.38,'Both')
pg.text(3.15,1.33,'Tsys = 19.5K')
pg.text(3.15,1.28,'dVel = 0.32 km/s')

pg.postscript('figure6.ps',color=F,landscape=F)