load MeanVpark3x3_060227.dat
%load MeanVsurf3x3_060227.dat

lon = MeanVpark3x3_060227(:,1);
lat = MeanVpark3x3_060227(:,2);
zonal = MeanVpark3x3_060227(:,3);
merid = MeanVpark3x3_060227(:,4);
num = MeanVpark3x3_060227(:,5);
%%
LON = unique(lon);
LAT = unique(lat);

ZON = NaN*zeros(length(LON),length(LAT));
MER = ZON;
NUM = ZON;
for ndo=1:length(zonal)
    ilon = find(lon(ndo)==LON);
    ilat = find(lat(ndo)==LAT);
    ZON(ilon,ilat) = zonal(ndo);
    MER(ilon,ilat) = merid(ndo);
    NUM(ilon,ilat) = num(ndo);
end
%%
pcolor(LON,LAT,(ZON.^2+MER.^2)')
shading flat
caxis([0 100])
colorbar
%% 
%quiver(lon,lat,zonal,merid)
ilat = find(lat>=50&lat<=70);
ilon = find(lon>=-70&lon<=-10);
myrange = intersect(ilon,ilat);
quiver(lon(myrange),lat(myrange),zonal(myrange),merid(myrange));
%%
if ~exist('yomaha070612','var')
load yomaha070612.dat
end
data = yomaha070612;
lon = data(:,1);
lat = data(:,2);
lons = data(:,9);
lats = data(:,10);
%% Reduce range
latrange = [52 67];
lonrange = [-70 -42];
ilat = find(lat>=min(latrange)&lat<=max(latrange));
ilat2 = find(lats>=min(latrange)&lats<=max(latrange));
ilon = find(lon>=min(lonrange)&lon<=max(lonrange));
ilon2 = find(lons>=min(lonrange)&lons<=max(lonrange));
myrange1 = intersect(ilon,ilat);
myrange2 = intersect(ilon2,ilat2);
myrange = union(myrange1,myrange2);

%%
lon = lon(myrange);
lat = lat(myrange);
lons = lons(myrange);
lats = lats(myrange);

park = data(myrange,3);
time = data(myrange,4);
zonal = data(myrange,5);
merid = data(myrange,6);
errx = data(myrange,7);
erry = data(myrange,8);

lons = data(myrange,9);
lats = data(myrange,10);
times = data(myrange,11);
zonals = data(myrange,12);
merids = data(myrange,13);
errxs = data(myrange,14);
errys = data(myrange,15);
lastlat = data(myrange,17);
lastlon = data(myrange,16);
lasttime = data(myrange,18);
nextlat = data(myrange,20);
nextlon = data(myrange,19);
nexttime = data(myrange,21);
nextlat2 = data(myrange,23);
nextlon2 = data(myrange,22);
nexttime2 = data(myrange,24);
numsurf = data(myrange,25);
float = data(myrange,26);
cycle = data(myrange,27);
dup  = data(myrange,28);
%% 
zrange = find(lon<-999);
lon(zrange) = NaN;
lat(zrange)=NaN;
zonal(zrange) =NaN;
merid(zrange) = NaN;
zrange = find(zonal<-999);
zonal(zrange)=NaN;
merid(zrange)=NaN;
errx(zrange)=NaN;
zrange = find(park>1800);
lon(zrange)=NaN;
lat(zrange)=NaN;
zrange=find(park<1200);
lon(zrange)=NaN;
lat(zrange)=NaN;
%%
zrange = find(lons<-999);
lons(zrange)=NaN;
lats(zrange)=NaN;
zonals(zrange)=NaN;
merids(zrange)=NaN;
errxs(zrange)=NaN;
%% Do a gridding
dlat=.25;
dlon=.25;
latgrid = min(latrange):dlat:max(latrange);
longrid = min(lonrange):dlon:max(lonrange);
zonalgrid = NaN*zeros(length(longrid),length(latgrid));
meridgrid = zonalgrid;
zonalgrids=zonalgrid;
meridgrids=zonalgrid;
for odo=1:length(longrid)
    lon1=longrid(odo);
    ilon = find(abs(lon-lon1)<dlon);
    ilons = find(abs(lons-lon1)<dlon);
    for ado=1:length(latgrid)
        lat1 = latgrid(ado);
        ilat = find(abs(lat-lat1)<dlat);
        ilats = find(abs(lats-lat1)<dlat);
        myrange = intersect(ilon,ilat);
        myranges = intersect(ilons,ilats);
        if length(myrange)
            zonalgrid(odo,ado) = meanNaN(zonal(myrange));
            meridgrid(odo,ado) = meanNaN(merid(myrange));
            numgrid(odo,ado) = sum(~isnan(zonal(myrange)));
            maxerr(odo,ado) = max(errx(myrange));
        end
        if length(myranges)
            zonalgrids(odo,ado) = meanNaN(zonals(myranges));
            meridgrids(odo,ado) = meanNaN(merids(myranges));
            numgrids(odo,ado) = sum(~isnan(zonals(myranges)));
            maxerrs(odo,ado) = max(errxs(myranges));
    end
end    
end
%% 
[LAT,LON]=meshgrid(latgrid,longrid);
[OO AA]=size(LAT);
lonlong = reshape(LON,OO*AA,1);
latlong = reshape(LAT,OO*AA,1);
zonallong = reshape(zonalgrid,OO*AA,1);
meridlong = reshape(meridgrid,OO*AA,1);
zonallongs = reshape(zonalgrids,OO*AA,1);
meridlongs = reshape(meridgrids,OO*AA,1);
%%
figure(1);
quiver(lonlong,latlong,zonallong,meridlong,5);
figure(2)
subplot(121)
pcolor(numgrid');
shading flat
colorbar
subplot(122)
pcolor(maxerr');
shading flat
caxis([0 1.35]);
colorbar
%%
create_xyzz(lonlong,latlong,zonallong,meridlong,'deep1500_uv.xyz');
%%
figure(3)
quiver(lonlong,latlong,zonallongs,meridlongs,5);
figure(4)
subplot(121)
pcolor(numgrids');
shading flat
colorbar
subplot(122)
pcolor(maxerrs')
shading flat
caxis([0 20]);
colorbar
%%
create_xyzz(lonlong,latlong,zonallongs,meridlongs,'surf_uv.xyz');