clc;clear all;close all
addpath('/data3/yiwen/cesm_final_exp/matlab_toolbox/colorbrewer')

%%
px=linspace(1,800,800);
py=linspace(1,800,800);

set(0,'defaultfigurecolor','w')
figure(2)
pcolor(px(280:679),py(81:400),harbor(81:400,280:679));shading flat;colorbar
% pcolor(Xp,Yp,harbor);shading interp;colorbar
colormap(flipud(brewermap(30, 'RdYlBu')));
c=colorbar('location','eastoutside');
c.Title.String='meter'
% axis([500 700 50 350])
% axis([121.9100 121.9200 25.1220 25.1300])
title('Bitou Fishing Harbor water level','Fontsize',18)























