; routine to select the index identifying the polys to plot
; should be used in tandem with PlotPoly
; if anything changes below, change PlotPoly accordingly

pro SelectPoly, QPoly

while (QPoly LT 1 OR QPoly GT 399) do begin
  print, "  > Select the map required: "
  print, "  >   outlines: 1=UK, 2=Eire, 3=BritishIsles"
  print, "  >   counties: 4=GB, 5=Eng, 6=Scot, 7=Wales, 8=BritIs"  
  print, "  >   counties: 100=Eng, 200=Wales, 300=Scot"
  read, QPoly
  
  if (QPoly EQ 100) then begin
	print, '  > .avon_=101, .beds_=102, .berks=103'
	print, '  > .bucks=104, .cambs=105, .chesh=106'
	print, '  > .cleve=107, .cornw=108, .cumbr=109'
	print, '  > .derbs=110, .devon=111, .dorst=112'
	print, '  > .durhm=113, .esuss=114, .essex=115'
	print, '  > .glouc=116, .manch=117, .hants=118'
	print, '  > .herew=119, .herts=120, .humbr=121'
	print, '  > .kent_=122, .lancs=123, .leics=124'
	print, '  > .lincs=125, .londi=126, .londo=127'
	print, '  > .mersy=128, .norfk=129, .nyork=130'
	print, '  > .nhant=131, .nrhum=132, .notts=133'
	print, '  > .oxfrd=134, .shrop=135, .somer=136'
	print, '  > .syork=137, .staff=138, .suffk=139'
	print, '  > .surry=140, .tynew=141, .warwk=142'
	print, '  > .wmidl=143, .wsuss=144, .wyork=145'
	print, '  > .wight=146, .wilts=147'
	read, QPoly
  endif else if (QPoly EQ 200) then begin  
	print, '  > .clwyd=205, .dyfed=206, .gwent=204'
	print, '  > .gwynd=207, .mglam=202, .powys=208'
	print, '  > .sglam=201, .wglam=203'
	read, QPoly
  endif else if (QPoly EQ 300) then begin
	print, '  > .bordr=301, .centr=302, .dugal=303'
	print, '  > .fife_=304, .gramp=305, .highl=306'
	print, '  > .lothi=307, .orkny=308, .shetl=309'
	print, '  > .strat=310, .taysd=311, .westi=312'
	read, QPoly
  endif
endwhile

end