[xi2,yi2]=meshgrid(-40:.05:-36,-54:.05:-50);%GASEX III site
ZI = griddata(lon,lat,tsnk,xi2,yi2);
figure;[c,h] = contour(xi2,yi2,ZI); clabel(c,h), colorbar;
 hold on;plot(lon,lat,'-')
 title('Sea T C)');
 

[xi2,yi2]=meshgrid(-40:.05:-36,-52:.05:-50);%GASEX III site
ZI = griddata(lon,lat,tsnk,xi2,yi2);
figure;[c,h] = contour(xi2,yi2,ZI); clabel(c,h), colorbar;
 hold on;plot(lon,lat,'-')
 title('Sea T C)');
 
[xi2,yi2]=meshgrid(-40:.05:-36,-52:.05:-50);%GASEX III site 
ZI = griddata(lon,lat,U,xi2,yi2);
figure;[c,h] = contour(xi2,yi2,ZI); clabel(c,h), colorbar;
hold on;plot(lon,lat,'-')
title('Wind Speed (m/s)');