return


% Process Andrew's data.

load KOE_GRD.mat
grd.T=200;
grd.L=399;
grd.P=249;

files=rnt_getfilenames('/nas/adavis44/KOE_3day', 'uv');


U=zeros(399,249, grd.T);
V=zeros(399,249, grd.T);

k=0;
for i=1:length(files)
    file=files{i};
    k=k+1;
    month(k)=str2num(file(27:28));
    day(k)=str2num(file(30:31));
    year(k)=str2num(file(33:36));
    load(file);
    U(:,:,k)=uf;
    V(:,:,k)=vf;
        
    if k==grd.T      
      time=datenum(year, month, day);
      filename=['uv_',num2str(i),'.nc'];
      disp(filename);
	Convert2NC(grd,filename);
      nc=netcdf(filename, 'w');
      nc{'time'}(:)=time;
      nc{'u'}(:,:,:)=perm(U);
      nc{'v'}(:,:,:)=perm(V);
      close(nc)
      U=zeros(399,249, grd.T);
      V=zeros(399,249, grd.T);
	clear day month year
	k=0;
    end 
    
    
end    

grd.T=k;
      time=datenum(year, month, day);
      filename=['uv_',num2str(i),'.nc'];
      disp(filename);
	Convert2NC(grd,filename);
      nc=netcdf(filename, 'w');
      nc{'time'}(:)=time;
      nc{'u'}(:,:,:)=perm(U(:,:,1:grd.T));
      nc{'v'}(:,:,:)=perm(V(:,:,1:grd.T));
      close(nc)
 


return
%% Load GRID:
% Process IPRC data.
user='theme2';
pass='apdrc_best';
ufile='http://theme2:apdrc_best@apdrc.soest.hawaii.edu:80/dods/iprc_esc/OfES/ncep_0.1_global_3day/uvel';
vfile='http://theme2:apdrc_best@apdrc.soest.hawaii.edu:80/dods/iprc_esc/OfES/ncep_0.1_global_3day/vvel';

nc=netcdf(ufile);
lon=nc{'lon'}(:)-360;
lat=nc{'lat'}(:);
[lat,lon]=meshgrid(lat, lon);
xlim=[125 165]-360;
ylim=[25 50];
[I,J]=rgrd_FindIJ(lon,lat, xlim, ylim);

% Define GRD
grd.lon=lon(I,J);
grd.lat=lat(I,J);
mask=perm(nc{'uvel'}(1,1,J,I));
in=find(mask < -10000);
mask(:)=1; mask(in)=nan;
grd.mask=mask;

% load time and adjust timescale
time=nc{'time'}(:) + datenum(1,-1,-1);
[year,month, day]=dates_datenum(time);
time=datenum(year,month,day);
grd.time=time;
grd.month=month;
grd.year=year;
grd.day=day;
grd.z=nc{'lev'}(:);
grd.cstfile='/dods/matlib/rgrd/rgrd_CoastlineWorld.mat';
close(nc);

[L,P]=size(grd.lon);
grd.T=10;
grd.L=L;
grd.P=P;

for it=920:10:2820
  str=['../uv_',num2str(it),'.mat'];
  load (str);
  
  grd.time=time(it-9:it);
  filename=['../nc/uv_',num2str(it),'.nc'];
  disp(filename);
  Convert2NC(grd,filename);
  nc=netcdf(filename, 'w');
  nc{'time'}(:)=grd.time;
  nc{'u'}(:,:,:)=perm(U);
  nc{'v'}(:,:,:)=perm(V);
  close(nc)
end