proc gcontour data=reflib.swirl;
   plot y*x=z;
run;