;
/*
Countries near the Mediteranean sea

This SAS program selects countries
around the Mediteranean sea for a map.

It sorts IDs form SAS/GRAPH MAP data set
WORLDMAP, version 6.10, for use in proc
GMAP.

 Arnold Schick
*/

*libname maps v610 'c:\sas\maps';

data one;
  set maps.worldmap;
  where cont in (93 94 95);
  x=-x;
run;

data two;
  set one;
  where x < 1.0 and x > -0.3 and
        y < 0.8 and y >  0.5;
run;

data three;
  set two;
  id_found = id;
  if id ^= lag(id) then output;
  keep id id_found;
run;

proc sort data=three out=three;
  by id;
run;
proc sort data=maps.worldmap out=four;
  by id;
run;

data five;
  merge four three;
  if id=id_found and id ^= 825;
  drop id_found;
  by id;
  x=-x;
run;

title1 h=5 pct f=zapf 'Countries';
title2 h=4 pct f=zapf 'near the';
title3 h=5 pct f=zapf 'Mediteranean Sea';
pattern v=e r=123;
proc gmap data=five map=five;
  id id;
  choro id / discrete nolegend;
run;
title1; title2; title3; pattern;

*
;