worldmap([20 27.5],[115 125])
load coast
geoshow(lat,long)