%Load a  data file
d=load('Jan_coads_lluvp.dat');

% get lon/lat's
lo=unique(d(:,1));
nlo=length(lo);
la=unique(d(:,2));
nla=length(la);

[LO,LA]=meshgrid(lo,la);
u_wind=reshape(d(:,3),size(LO));
v_wind=reshape(d(:,4),size(LO));
surf_pres=reshape(d(:,5),size(LO));


% MATLAB's interp2 commands can be used to map the
%1 deg grid data to arbitrary lon/lats
ui = interp2(LO,LA,u_wind,loi,lai);

% where loi,lai are used-defined locations to interpolate to.