function get_image, ipar1, icol1
@rdsumerext_block
@xraster_draw_block
case ipar1 of 
         0: begin ; Intensity
             im = data_fit(1,*,*,icol1)*data_fit(3, *,*,icol1)/26.90865
             end
         1 : begin ; Width
              im = data_fit(3, *, *, icol1)*(1.414/100./0.6745)*deltas(0, icol1)
              end
         2 : begin ; Shift
              wl = float(strmid(data_names(icol1), 2, 6))
               im=fltarr(ydim, xdim)
    
     for k=0,  ydim-1 do   $
     im(k, *) = (data_fit(2, k, *, icol1)-median(data_fit(2,k,*, icol1)))/100. $
                * deltas(0, icol1) /wl*3.e5
             end
endcase
 im1 = reform(data_fit(1,*,*,icol1)*data_fit(3, *,*,icol1)/26.90865)
 im=reform(im)
s=where(im1 le 0., count) 
;s1=where(im1 gt 100.)
if count ge 1 then im(s) = median(im)
return, rotate(reform(im),5)
end