;;RESTORE,'/data/2023sep07/aligned_results/vis-aligned-xu-correct-level-0.5.sav'


img=RF(hadata[5,0,*,*])
dopp=RF(hadata[5,1,*,*])
img=img[250:1800,250:1800]
dopp=dopp[250:1800,250:1800]
nx=(SIZE(img,/DIM))[0]
ny=(SIZE(img,/DIM))[1]

SET_PLOT,'X'

tmap=(img-0.533)/(6.12*1e-5)

LOADCT, 62
WINDOW, xs=nx, YS=ny, /PIXMAP
tv, bytscl(-tmap,min=-12000,max=-7000)

LOADCT, 0
SQUARE, /FILL

xb=793 & yb=597
PLOTS, [2*xb-250,2*xb-250],[2*yb-250,2*yb-250], PSYM=8, SYMSIZE=1.2, COLOR=0, /DEVICE

xr=805 & yr=571
PLOTS, [2*xr-250,2*xr-250],[2*yr-250,2*yr-250], PSYM=8, SYMSIZE=1.2, COLOR=0, /DEVICE

xb=656 & yb=576
PLOTS, [2*xb-250,2*xb-250],[2*yb-250,2*yb-250], PSYM=8, SYMSIZE=1.2, COLOR=0, /DEVICE

;xr=845 & yr=529
xr=636 & yr=585
PLOTS, [2*xr-250,2*xr-250],[2*yr-250,2*yr-250], PSYM=8, SYMSIZE=1.2, COLOR=0, /DEVICE

x0=670 & y0=546
PLOTS, [2*x0-250,2*x0-250],[2*y0-250,2*y0-250], PSYM=8, SYMSIZE=1.2, COLOR=0, /DEVICE

img=TVRD(/TRUE)


LOADCT, 62
scale=BYTARR(256,20)
for i=0,19 do scale[*,i]=-findgen(256)
WINDOW, xs=256, YS=20, /PIXMAP
TV, scale
scale=TVRD(/TRUE)

LOADCT, 0

SET_PLOT,'PS'

!P.CHARSIZE=1.5
aspect=FLOAT(ny)/nx
psx=18
psy=psx*aspect

!P.FONT=0
!P.THICK=3
DEVICE, XS=psx, YS=psy, /ENCAPSULATED,BITS_PER_PIXEL=8, /COLOR, filename='T-map-fwhm.eps', /avantgard,/book

x1=0.12 & x2=0.98 & y1=0.12 & y2=0.98

TV, img, /TRUE, x1, y1, XS=x2-x1, YS=y2-y1, /NORMAL

;LOADCT, 16
;CONTOUR,bz, POS=[x1,y1,x2,y2],XS=5,YS=5,XR=[0,720],YR=[0,720],LEVELS=[75],COLOR=250, $
;        XTICKFORMAT='(a1)', YTICKFORMAT='(a1)', /NORMAL, /NOERASE
;LOADCT, 0

xr1=0 & yr1=0
xr2=nx*0.029 & yr2=ny*0.029

PLOT, [xr1,xr2],[yr1,yr2],XR=[xr1,xr2],YR=[yr1,yr2],XS=1, YS=1,XTIT='arcsec', YTIT='arcsec', /NODATA, $
POSITION=[x1,y1,x2,y2], /NORM, XTICKLEN=-0.02, YTICKLEN=-0.02*aspect, /NOERASE, /ISO, XMINOR=10, YMINOR=10

TV, scale, /TRUE, 0.5, 0.943, XS=0.3, YS=0.02, /NORM
LOADCT, 0

PLOT, [0,255],[0,1], POSITION=[0.5,0.943,0.8,0.963], YTICKFORMAT='(a1)', yTICKLEN=0.00001, $
XR=[7,12], XS=1,/NODATA, XTICkLEN=0.2, XTICKS=5, XMINOR=2, CHARSIZE=1.0, /NORM, /NOERASE, THICK=0, COLOR=255
XYOUTS, 12.4,-1.4,'x 10!U3!N, K', COLOR=255, CHARSIZE=1.0, /DATA

;XYOUTS, xr1+4+0.1, yr1+4+0.1, 'GST/VIS 0.08 nm Doppler 2017 June 19 18:45-18:58 UT', FONT=1, COLOR=50, CHARSIZE=1.75, /DATA
;XYOUTS, xr1+4, yr1+4,         'GST/VIS 0.08 nm Doppler 2017 June 19 18:45-18:58 UT', FONT=1, COLOR=255, CHARSIZE=1.75, /DATA

XYOUTS, 0.87, 0.65, '1', COLOR=0, CHARSIZE=1.6, /NORM
XYOUTS, 0.88, 0.62, '2', COLOR=0, CHARSIZE=1.6, /NORM
XYOUTS, 0.72, 0.62, '3', COLOR=0, CHARSIZE=1.6, /NORM
XYOUTS, 0.68, 0.64, '4', COLOR=0, CHARSIZE=1.6, /NORM
XYOUTS, 0.71, 0.56, '5', COLOR=0, CHARSIZE=1.6, /NORM

DEVICE, /CLOSE
SET_PLOT,'X'


;------------------------;

LOADCT, 70
WINDOW, xs=nx, YS=ny, /PIXMAP
tv, bytscl(-dopp,min=-0.1,max=0.1)

LOADCT, 0, /SILENT
SQUARE, /FILL

xb=793 & yb=597
PLOTS, [2*xb-250,2*xb-250],[2*yb-250,2*yb-250], PSYM=8, SYMSIZE=1.2, COLOR=255, /DEVICE

xr=805 & yr=571
PLOTS, [2*xr-250,2*xr-250],[2*yr-250,2*yr-250], PSYM=8, SYMSIZE=1.2, COLOR=255, /DEVICE

xb=656 & yb=576
PLOTS, [2*xb-250,2*xb-250],[2*yb-250,2*yb-250], PSYM=8, SYMSIZE=1.2, COLOR=0, /DEVICE

;xr=845 & yr=529
xr=636 & yr=585
PLOTS, [2*xr-250,2*xr-250],[2*yr-250,2*yr-250], PSYM=8, SYMSIZE=1.2, COLOR=0, /DEVICE

x0=670 & y0=546
PLOTS, [2*x0-250,2*x0-250],[2*y0-250,2*y0-250], PSYM=8, SYMSIZE=1.2, COLOR=0, /DEVICE

LOADCT, 70, /SILENT

img=TVRD(/TRUE)

scale=BYTARR(256,20)
for i=0,19 do scale[*,i]=-findgen(256)
WINDOW, xs=256, YS=20, /PIXMAP
TV, scale
scale=TVRD(/TRUE)

LOADCT, 0

SET_PLOT,'PS'

!P.CHARSIZE=1.5
aspect=FLOAT(ny)/nx
psx=18
psy=psx*aspect

!P.FONT=0
!P.THICK=3
DEVICE, XS=psx, YS=psy, /ENCAPSULATED,BITS_PER_PIXEL=8, /COLOR, filename='Dopp-fwhm.eps', /avantgard,/book

x1=0.12 & x2=0.98 & y1=0.12 & y2=0.98

TV, img, /TRUE, x1, y1, XS=x2-x1, YS=y2-y1, /NORMAL

;LOADCT, 16
;CONTOUR,bz, POS=[x1,y1,x2,y2],XS=5,YS=5,XR=[0,720],YR=[0,720],LEVELS=[75],COLOR=250, $
;        XTICKFORMAT='(a1)', YTICKFORMAT='(a1)', /NORMAL, /NOERASE
;LOADCT, 0

xr1=0 & yr1=0
xr2=nx*0.029 & yr2=ny*0.029

PLOT, [xr1,xr2],[yr1,yr2],XR=[xr1,xr2],YR=[yr1,yr2],XS=1, YS=1,XTIT='arcsec', YTIT='arcsec', /NODATA, $
POSITION=[x1,y1,x2,y2], /NORM, XTICKLEN=-0.02, YTICKLEN=-0.02*aspect, /NOERASE, /ISO, XMINOR=10, YMINOR=10

TV, scale, /TRUE, 0.5, 0.943, XS=0.3, YS=0.02, /NORM
LOADCT, 0

;ps_grid

PLOT, [0,255],[0,1], POSITION=[0.5,0.943,0.8,0.963], YTICKFORMAT='(a1)', yTICKLEN=0.00001, $
XR=[-0.01,0.01], XS=1,/NODATA, XTICkLEN=0.2, XTICKS=4, XTICKFORMAT='(F6.3)',CHARSIZE=1.,/NORM, /NOERASE, THICK=0, COLOR=0
XYOUTS, 0.012,-1.4,' nm', FONT=1, COLOR=0, CHARSIZE=1.5, /DATA


;XYOUTS, xr1+4+0.1, yr1+4+0.1, 'GST/VIS 0.08 nm Doppler 2017 June 19 18:45-18:58 UT', FONT=1, COLOR=50, CHARSIZE=1.75, /DATA
;XYOUTS, xr1+4, yr1+4,         'GST/VIS 0.08 nm Doppler 2017 June 19 18:45-18:58 UT', FONT=1, COLOR=255, CHARSIZE=1.75, /DATA

XYOUTS, 0.87, 0.65, '1', COLOR=0, CHARSIZE=1.6, /NORM
XYOUTS, 0.88, 0.62, '2', COLOR=0, CHARSIZE=1.6, /NORM
XYOUTS, 0.72, 0.62, '3', COLOR=0, CHARSIZE=1.6, /NORM
XYOUTS, 0.68, 0.64, '4', COLOR=0, CHARSIZE=1.6, /NORM
XYOUTS, 0.71, 0.56, '5', COLOR=0, CHARSIZE=1.6, /NORM

DEVICE, /CLOSE
SET_PLOT,'X'

END