获取海流数据并统计海流的变化规律

1. 海流源数据来自于：中国科学院海洋科学数据中心

http://msdc.qdio.ac.cn/data/metadata-special-detail?id=1462700665938763778&cnId=1462700665938763778&enId=1462700665968123905

2. 读取海流数据，并根据海流数据的格式，提取目标海域的表层海流数据。

close all;
clear all;
clc;
UE=load('u_eastward_monthly_1997_2016.mat');
VN=load('v_northward_monthly_1997_2016.mat');
sur_sh_E=squeeze(UE.data(30,:,:,:));
sur_sh_N=squeeze(VN.data(30,:,:,:));
% %目标位置：31°06' N,122°08' E
% % 117.01°E~131.66°E，分为288点
% % 29.04°N~42.09°N，分为324点
% % 由表到深30层，由远至近240月，
% % 文件中324×288表示模式在经向（主要为南北方向上）和纬向（东西方向上）的网格点数，与depth中的288×324意思类似，只是在维度存在转置的关系。
% % date中的6列数，应为年月日时分秒，可用MATLAB函数转为其他形式。希望有所帮助。
pa=round(((31+6/60)-29.04)/(42.09-29.04)*324);
pb=round(((122+8/60)-117.01)/(131.66-117.01)*288);
UE.depth(pb,pa,1);
UE.depth(pb,pa,30);
VN.depth(pb,pa,1);
VN.depth(pb,pa,30);
% 10/324*(42.09-29.04); bound
% 10/288*(131.66-117.01); bound
AE=sur_sh_E(pa-9:pa+10,pb-9:pb+10,:); %[-0.4°,0.4°]
BN=sur_sh_N(pa-9:pa+10,pb-9:pb+10,:); %[-0.51°,0.51°]bound ~ 1 knot
EN=sqrt(AE.^2+BN.^2);
DEN=atan(AE./BN);%与北向夹角

3.将预处理的数据进行可视化处理

[a1,b1,c1]=size(AE);%East