Generate the surface plot.

proc g3d data=sashelp.humid;
   scatter airtemp*bulbtemp=humidity/
      shape=pillar
       color=colorvar
       caxis=blue
       rotate=-15
       yticknum=5
       xticknum=2
       zticknum=4
       zmin=0
       zmax=100
run;
quit;