proc g3grid data=reflib.pond out=reflib.pondgrid;
   grid vdist*hdist=height / naxis1=100 naxis2=100;
run;