x cd /pub/data/sf32000/Tools;  %let pgm=usgeos_reports;  filename pgm "&pgm..sas";
  *---Genrate reports based on the sf3.usgeos data set.----;
 libname sf32000 '/pub/data/sf32000';
  libname sf3 (sf32000);

 proc sort data=sf3.usgeos(keep=sumlev geoid state county geocode areaname where=(state ne '72')) out=usgeos; by geoid;  run;
 data sf3.usgeos_sumlev_summary;
  set usgeos;
  by sumlev;
  if first.sumlev then n=1; else n+1;
  if last.sumlev then do;
    geoidlen=length(geoid);
    sumlev_label=put(sumlev,$sumlev.);
    output;
    end;
   run;
   ods html file="&pgm..html";
   footnote .h=1 "From setup in usgeos_reports.sas";

 proc print data=sf3.usgeos_sumlev_summary;
   id sumlev;
   format sumlev_label $40.;
   title 'sf3.usgeos:  SumLev Values with Typical geoid and geocode Values';
    title2 '(sumlev_label has been truncated to 40 chars)';
   var sumlev_label  N  geoid  geocode areaname state geoidlen;
   format n comma7.;
   run;

 ods html close;
 %include sascode(notify);