clc; clear;clf;
sta = 'Boluo2';
yr = 2018 ;
mn = 1 ;
%========= read data ========%
fname=['C:\Users\Administrator\Desktop\discharge',sta,'_',num2str(yr),'_filled.csv'];
A = importdata(fname,' ');
tt = datenum(A.textdata);
aa = A.data;
%========== plot figure =========%
figure(1)
subplot(4,1,1)
xmin = datenum(yr,mn,1);
xmax = datenum(yr,mn+1,1);
ik = find(and(tt>=xmin,tt<=xmax));
plot(tt(ik),aa(ik),'.-');
grid on
set(gca,'xlim',[xmin,xmax],'xtick',[xmin:1:xmax]);
datetick('x','keepticks');
RemoveSubplotWhiteArea(gca, 4, 1, 1, 1);
subplot(4,1,2)
xmin = datenum(yr,mn+1,1);
xmax = datenum(yr,mn+2,1);
ik = find(and(tt>=xmin,tt<=xmax));
plot(tt(ik),aa(ik),'.-');
grid on
set(gca,'xlim',[xmin,xmax],'xtick',[xmin:1:xmax]);
datetick('x','keepticks');
RemoveSubplotWhiteArea(gca, 4, 1, 2, 1);