clc;clear %先对分量进行调和分析,然后利用elli_para2程序计算椭圆要素,ap2ep程序画潮流椭圆
%%
xlsend = [3477 3897 3477 3478];
lat_point = [39+39.922/60 39+42.356/60 39+37.024/60 39+25.933/60];
tidename = ['O1 ';'K1 ';'M2 ';'S2 ';'M4 ';'MS4 '];
color={'r';'g';'b';'k';'c';'m';'r--';'g--';'b--';'k--';'c--';'m--'};
ef_l=[11 17 31 14];%有效水层
depth=[15.5 22 19.8 18.6] %局地水深
%% 流速分析
%层循环
for i_site = 1:4 %站位循环
% subplot(3,4,(layer-1)*4+i_site);
time = xlsread('.\QHD_ADCP.xlsx', i_site, 'B2:G2');
if i_site==2
N=xlsread('.\QHD_ADCP.xlsx',i_site,['J2:BZ',num2str(xlsend(i_site))]