function [LON,LAT]=make_meshgrid(outfile,S_LAT,W_LON,N_LAT,E_LON,nlon,nlat)
  
  lon=linspace(W_LON,E_LON,nlon);
  lat=linspace(S_LAT,N_LAT,nlat);
  
  [LAT,LON]=meshgrid(lat,lon);
  mg=[LON(:) LAT(:)];
  dlmwrite(outfile,mg,'delimiter','\t','precision',6)