PRO ryr, x, y, y1, y2

; Restrict the y-range of a 1-d function

if y1 GE y2 then return

ok = where (y GE y1 AND y LE y2, n_ok)
IF n_ok GT 0 THEN begin
    x = x(ok)
    y = y(ok)
ENDIF

END