;
;
;
if (n_elements(sxc) eq 0) then begin
    input, 'Enter the start date to extract ', sttim, '1-sep-91'
    input, 'Enter the end date to extract   ', entim, gt_day(!stime, /str)
    rd_sxc, sttim, entim, index2, sxc
    ;rd_sxc, '1-sep-91', !stime, index2, sxc

    rr = get_rb0p(index2)
    rr = rr(0,*)/2.45

    ref_tim = index2(0)
    x = int2secarr(index2, ref_tim)		;do this once only
end
;
!p.multi = [0,1,4]
!p.charsize = 1.4
strver = 'PLOT_SXC.PRO  Ver 1.1'
;
d1 = sxc.suncenter(0)-sxc.sxt_center(0)
utplot, x, sxc.sxt_center(0), ref_tim, ytit='X-Ray E/W Center (pix)', /ynozero, timerange=timerange, psym=1
utplot, x, sxc.suncenter(0), ref_tim, ytit='GET_SUNCENTER E/W (pix)', /ynozero, timerange=timerange, psym=1
diff = sxc.suncenter(0)-sxc.sxt_center(0)
ss = where(abs(diff) lt 5) & dev = stdev(diff(ss), avg) & mtit='Avg = ' + string(avg) + ' Stdev = ' + string(dev)
utplot, x, d1, ref_tim, ytit='GET_SUNCENTER -XRay (E/W)', psym=1, tit=mtit, timerange=timerange
utplot, x, d1, ref_tim, ytit='GET_SUNCENTER -XRay (E/W)', psym=1, tit=mtit, yrange=[-10,10], timerange=timerange
plottime, 0, 0, strver & plottime & pause
;
d2 = sxc.suncenter(1)-sxc.sxt_center(1)
utplot, x, sxc.sxt_center(1), ref_tim, ytit='X-Ray N/S Center (pix)', /ynozero, timerange=timerange, psym=1
utplot, x, sxc.suncenter(1), ref_tim, ytit='GET_SUNCENTER N/S (pix)',/ynozero, timerange=timerange, psym=1
diff = sxc.suncenter(1)-sxc.sxt_center(1)
ss = where(abs(diff) lt 5) & dev = stdev(diff(ss), avg) & mtit='Avg = ' + string(avg) + ' Stdev = ' + string(dev)
utplot, x, d2, ref_tim, ytit='GET_SUNCENTER -XRay (N/S)', psym=1, tit=mtit, timerange=timerange
utplot, x, d2, ref_tim, ytit='GET_SUNCENTER -XRay (N/S)', psym=1, tit=mtit, yrange=[-10,10], timerange=timerange
plottime, 0, 0, strver & plottime & pause
;
d3 = rr - sxc.sxt_center(2)
utplot, x, sxc.sxt_center(2), ref_tim, ytit='X-Ray Radius (pix)', /ynozero, timerange=timerange, psym=1
utplot, x, rr,                ref_tim, ytit='GET_RB0P Radius (pix)',/ynozero, timerange=timerange, psym=1
utplot, x, d3,                ref_tim, ytit='GET_RB0P - XRay Radius', psym=1, timerange=timerange
utplot, x, d3,                ref_tim, ytit='GET_RB0P - XRay Radius', psym=1, yrange=[-10,10], timerange=timerange
plottime, 0, 0, strver & plottime & pause
;
stop
;
if (n_elements(sot) ne 0) then begin
    if (n_elements(xsot) eq 0) then begin
	ref_tim = index2(0)
	tout = int2secarr(index2, ref_tim)
	tin = int2secarr(index, ref_tim)
	xsot0 = sot.find_limb(0)
	ysot0 = sot.find_limb(1)
	xsot = interpol( xsot0, tin, tout )
	ysot = interpol( ysot0, tin, tout )
    end

    !p.multi = [0,1,3]
    utplot, index2, sxc.sxt_center(0), ytit='X-Ray E/W Center (pix)', /ynozero
    utplot, index2, xsot,              ytit='SOT E/W Center (pix)', /ynozero
    diff = xsot-sxc.sxt_center(0)
    ss = where(abs(diff) lt 5) & dev = stdev(diff(ss), avg) & mtit='Avg = ' + string(avg) + ' Stdev = ' + string(dev)
    utplot, index2, xsot-sxc.sxt_center(0), ytit='SOT - XRay (E/W)', psym=1, yrange=yrange_d2, tit=mtit
    plottime, 0, 0, strver & plottime & pause

    utplot, index2, sxc.sxt_center(1), ytit='X-Ray N/S Center (pix)', /ynozero
    utplot, index2, ysot,              ytit='SOT N/S Center (pix)', /ynozero
    diff = ysot-sxc.sxt_center(1)
    ss = where(abs(diff) lt 5) & dev = stdev(diff(ss), avg) & mtit='Avg = ' + string(avg) + ' Stdev = ' + string(dev)
    utplot, index2, ysot-sxc.sxt_center(1), ytit='SOT - XRay (N/S)', psym=1, yrange=yrange_d2, tit=mtit
    plottime, 0, 0, strver & plottime & pause
end

end