f = netcdf.open('honshu2011ote_h.nc','nowrite');
lon=netcdf.getVar(f,netcdf.inqVarID(f,'lon'));
lat=netcdf.getVar(f,netcdf.inqVarID(f,'lat'));
h=netcdf.getVar(f,netcdf.inqVarID(f,'ha'));
t=netcdf.getVar(f,netcdf.inqVarID(f,'time'));
si=500 %Frame number
ha=double(h(:,:,si));
ha=ha';
loc=find(ha<-10000000);
ha(loc)=NaN;
surf(lon,lat,ha,'edgecolor','none')
