function get_value, x, y1,  ipar1, icol1
@rdsumerext_block
@xraster_draw_block
y=ydim-y1-1
case ipar1 of 
         0: begin ; Intensity
             im = data_fit(1,y,x,icol1)*data_fit(3, y,x,icol1)/26.90865
             end
         1 : begin ; Width
              im = data_fit(3, y, x, icol1)*(1.414/100./0.6745)*deltas(0, icol1)
              end
         2 : begin ; Shift
   
              wl = float(strmid(data_names(icol1), 2, 6))
        im = (data_fit(2, y, x, icol1)-median(data_fit(2,y,*, icol1)))/100. $
                * deltas(0, icol1) /wl*3.e5
             
             end
endcase
return, im
end