#
# 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


x1[1:8]   := [290, 290, 290, 290, 290, 290, 290, 290]
x1[9:12]  := [570, 570, 570, 570]
x1[13:14] := [1160, 1160]
x1[15]    := [2300]

x2[1:8]   := [310, 310, 310, 310, 310, 310, 310, 310]
x2[9:12]  := [630, 630, 630, 630]
x2[13:14] := [1240, 1240]
x2[15]    := [2500]

x3[1:8]   := [600, 600, 600, 600, 600, 600, 600, 600]
x3[9:12]  := [1200, 1200, 1200, 1200]
x3[13:14] := [2400, 2400]
x3[15]    := [4800]

y1 := [15.3,16.2,13.5,13.5,13.4,13.7,13.0,13.0,9.95,9.76,9.85,9.42,7.32,7.14,5.55]

y2 := [15.1,15.8,13.1,13.0,13.1,13.5,12.9,13.3,11.0,9.41,9.66,9.34,6.91,6.70,4.82] 

y3 := [9.79,12.1,9.27,9.37,9.37,9.30,9.35,9.51,7.64,6.90,6.80,6.73,5.09,4.91,3.77]

# convert to logs

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

xline := log([200,7200])
yline  := log( [(14.1 * sqrt(300/200) ), (14.1 * sqrt(300/7200) ) ] )

# make the plot

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

xlbl := "LOG10(Integration time in seconds)"
ylbl := "LOG10(Observed RMS in mK)"
toplbl := "TSPECTEST_03: RMS in 50MHz 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 = 19K')
pg.text(3.15,1.28,'dVel = 1.29 km/s')

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