function range_x,x1,x2,dx
   
;- make a scale based on the input
   
   nx   = abs( fix( (x2 - x1) / dx + 1) )
   x    = x2 - findgen(nx) * dx 
   
   return,x
   end