************************************************************************
*City plotter
* Galebonwe Ramaphane - Department of Meteorological Services - Botswana
* 30 Nov. 2010
* Adapted from the WRF_EMS graphics package by:
* Dan Leins - NWS WFO Cleveland, OH
**************************************************************************
*Open the list of cities
cities = cities.txt

*start the while loop which reads all the points from the cities.txt file
start=0
while(start<100)
say 'start is 'start
 ret = read(cities)
 rc = sublin(ret,1)
 if (rc>0)
   if (rc!=2)
     say 'File I/O Error'
     return
   endif
   break
 endif

rec = sublin(ret,2)
 lat.start=subwrd(rec,1)
 lon.start=subwrd(rec,2)
 site.start=subwrd(rec,3)
say 'lat is 'lat.start' lon is 'lon.start' and the city is 'site.start
start=start+1
endwhile
last_city=start-1
say 'last city is 'last_city

*This is the section where we actually plot the cities
start=0
while (start <=last_city)
'q w2xy 'lon.start' 'lat.start
say result
xx.start=subwrd(result,3)
yy.start=subwrd(result,6)

*Mark each point with a circle
say 'This means x is 'xx.start' and y is 'yy.start
*'set string 0 c 10 0'
*'set strsiz 0.3'
*'draw string 'xx.start' 'yy.start' '.
'draw mark '3' 'xx.start' 'yy.start' '0.08

*Now label each period with the city name - the name is placed above the circle
*'set string 0 tc 5 0'
*'set strsiz 0.10'
'draw string 'xx.start+0.1' 'yy.start+0.1' 'site.start
start=start+1
endwhile