Generate the surface plot.

proc g3d data=iris;
    scatter PetalLength*PetalWidth=SepalLength/
      color=color
      shape=shape
      noneedle
      grid;
run;
quit;