;RESTORE,'vis-aligned-xu-products-level-0.5.sav' ;dmap=reform(hadata[5,1,*,*]) ;RESTORE,'vis-aligned-xu-products-level-0.1.sav' ;tmap=reform(hadata[5,0,*,*]) dim=SIZE(haData,/DIM) WINDOW, 0, XS=dim[2]/2., YS=dim[3]/2. ; --- plot line width LOADCT, 9 TV, bytscl(myresize(smooth(tmap,3)),min=0.2,max=0.7) LOADCT, 1 CONTOUR, MYRESIZE(dmap), POS=[0,0,dim[2]/2.,dim[3]/2.],/DEVICE,XR=[0,dim[2]/2.],XS=1,YR=[0,dim[3]/2.],YS=1,/NOERASE, $ LEVELS=[-0.03], C_COLOR=100 LOADCT, 3 CONTOUR, MYRESIZE(dmap), POS=[0,0,dim[2]/2.,dim[3]/2.],/DEVICE,XR=[0,dim[2]/2.],XS=1,YR=[0,dim[3]/2.],YS=1,/NOERASE, $ LEVELS=[0.03], C_COLOR=180 X2PNG, filename='Temp-Dopp.png', /COLOR LOADCT, 33 TV, bytscl(myresize(smooth(dmap,3)),min=-0.1,max=0.1) X2PNG, filename='Dopp.png', /COLOR ;WINDOW, 1, XS=dim[2]/2., YS=dim[3]/2. ; --- plot doppler shift at the line core ;imgm=tmap*dmap ;imgd=tmap/FLOAT(dmap) ;tv, bytscl(myresize(img),min=-1,max=1) ; --- plot intensity at the core ;WINDOW, 2, XS=dim[2]/2., YS=dim[3]/2. ;tv, bytscl(myresize(reform(hadata[5,2,*,*])),min=-5,max=5) ; --- to conver the width into Kelvins (from Leenaarts et al. 2012, ApJ 749, 136) dx=0.8-0.48 ; A dy=8.5-4 ; K k=dy/dx y0=k*0.8-8.5 ;T=y0+k*line width (i.e., hadata[i,0,*,*]) END