Interpolate values for the contour plot. The interpolated data set is stored in REFLIB.CLAY2. 
proc g3grid data=reflib.clay out=reflib.clay2;
   grid y*x=pct_clay / naxis1=21
                       naxis2=21
                       join;
run;