/* project block group boundaries */
/* L. Hoyle */

libname bg '/homed/blockgroups';

data bgrad;
  set bg.bgbounds;
  x=x*arcos(-1)/180;
  y=y*arcos(-1)/180;


proc gproject data=bgrad out=bg.bgproj;
  id blkgrp;
run;