%let name=surf2;
filename odsout '.';

data nums;
   keep X Y Z;
   do i=1 to 30;
      X=10*ranuni(33)-5;
      Y=10*ranuni(35)-5;
      Z=sin(sqrt(x*x+y*y));
      output;
   end;
run;

proc g3grid data=nums out=a;
   grid y*x=z / spline
                axis1=-5 to 5 by .02
                axis2=-5 to 5 by .02;
run;

data a; set a;
z=z+2;
run;


 GOPTIONS DEVICE=gif;
 
 ODS LISTING CLOSE;
 ODS HTML path=odsout body="&name..htm" (title="Surface Chart (top view)") style=minimal;
 goptions noborder;

goptions gunit=pct htitle=6 ftitle="albany amt/bold" htext=4.25 ftext="albany amt/bold";

goptions xpixels=720 ypixels=600;

pattern1 v=solid color=cx9999ff;  /* light blue */
pattern2 v=solid color=cx993366;  /* purplish */

axis1 label=(a=90 'Y Axis') offset=(0,0) order=(-5 to 5 by 2.5) minor=none;
axis2 label=('X Axis') offset=(0,0) order=(-5 to 5 by 2.5) minor=none;

title link="area1.htm" "Surface Chart (top view)";
proc gcontour data=a; 
plot y*x=z / 
 pattern
 join
 levels= 1.1 2.3
 caxis=black
 vaxis=axis1 
 haxis=axis2
 nolegend
 des="" name="&name" ;  
run;


 quit;
 ODS HTML CLOSE;
 ODS LISTING;