% lookat at 4x4 netcdf output

f=netcdf('ECCO_MaskAreaBathy.nc');
for i=1:11
  f{i}
  minmax(f{i}(:))
  pause
end

f=netcdf('ECCO_Stationary_1_phys.nc')
for i=1:10
  f{i}
  minmax(f{i}(:))
  pause
end

f=netcdf('ECCO_Stationary_3001_phys.nc')
for i=1:10
  f{i}
  minmax(f{i}(:))
  pause
end

f=netcdf('ECCO_Timedep_phys.nc')
for i=1:10
  f{i}
  minmax(f{i}(:))
  pause
end

f=netcdf('ECCO_Stationary_BasisFNCTNS_3001.nc')
tmp=f{1}(:);
for i=1:30
  clf
  subplot(221),mypcolor(squeeze(tmp(i,1,:,:)));colorbar,title(i)
  subplot(222),mypcolor(squeeze(tmp(i,3,:,:)));colorbar,title(i)
  subplot(223),mypcolor(squeeze(tmp(i,6,:,:)));colorbar,title(i)
  subplot(224),mypcolor(squeeze(tmp(i,9,:,:)));colorbar,title(i)
  pause(1)
end

f=netcdf('ECCO_Timedep_BasisFNCTNS_1775.nc')
tmp=f{1}(:);
for i=1:30
  clf
  subplot(221),mypcolor(squeeze(tmp(i,1,:,:)));colorbar,title(i)
  subplot(222),mypcolor(squeeze(tmp(i,3,:,:)));colorbar,title(i)
  subplot(223),mypcolor(squeeze(tmp(i,6,:,:)));colorbar,title(i)
  subplot(224),mypcolor(squeeze(tmp(i,9,:,:)));colorbar,title(i)
  pause(1)
end

f=netcdf('ECCO_Timedep_BasisFNCTNS_2005.nc')
tmp=f{1}(:);
for i=1:30
  clf
  subplot(221),mypcolor(squeeze(tmp(i,1,:,:)));colorbar,title(i)
  subplot(222),mypcolor(squeeze(tmp(i,3,:,:)));colorbar,title(i)
  subplot(223),mypcolor(squeeze(tmp(i,6,:,:)));colorbar,title(i)
  subplot(224),mypcolor(squeeze(tmp(i,9,:,:)));colorbar,title(i)
  pause(1)
end

f=netcdf('ECCO_Stationary_diag_2D.nc')
tmp1=f{2}(:);
tmp2=f{3}(:);
for i=1:30
  clf
  subplot(211),mypcolor(squeeze(tmp1(i,:,:)));colorbar,title(i)
  subplot(212),mypcolor(squeeze(tmp2(i,:,:)));colorbar,title(i)
  pause(.1)
end

f=netcdf('ECCO_Timedep_diag_0D.nc')
tmp1=f{1}(:);
tmp2=f{2}(:);
tmp3=f{3}(:);
tmp4=f{4}(:);
for i=1:30
  clf
  subplot(311),plot(tmp1,tmp2(:,i))
  subplot(312),plot(tmp1,tmp3(:,i))
  subplot(313),plot(tmp1,tmp4(:,i))
  title(i), pause(1)
end

f=netcdf('ECCO_Timedep_diag_2D.nc')
tmp1=f{2}(:);
tmp2=f{3}(:);
for i=1:30
  clf
  subplot(221),mypcolor(squeeze(tmp1(1,i,:,:)));colorbar,title(i)
  subplot(222),mypcolor(squeeze(tmp1(55,i,:,:)));colorbar
  subplot(223),mypcolor(squeeze(tmp2(1,i,:,:)));colorbar,title(i)
  subplot(224),mypcolor(squeeze(tmp2(55,i,:,:)));colorbar
  pause(.1)
end

f=netcdf('ECCO_Stationary_diag_0D.nc')
tmp1=f{1}(:);
tmp2=f{2}(:);
tmp3=f{3}(:);
tmp4=f{4}(:);
for i=1:30
  clf
  subplot(311),plot(tmp1,tmp2(:,i))
  subplot(312),plot(tmp1,tmp3(:,i))
  subplot(313),plot(tmp1,tmp4(:,i))
  title(i), pause(1)
end