[LON,LAT]=meshgrid(lon,lat);
LON=LON';
LAT=LAT'; 
pcolor(LON,LAT,vardata(:,:,1))
shading interp

%this line add the coastline
world_coastp('k','linewidth',2)