matlab计算大洋经圈翻转流函数怎么都算不对

xueyang-20

用的一套再分析资料里的经向海流速度，在垂直方向上积分得到后再在纬向上做平均，得到moc的流函数图。不知道为什么画出来的结果就是不对，怎么找都找不出错误，求求有没有大佬能帮忙看一眼呢

以下是代码和图

vc = squeeze(mean(vomecrty,4,'omitnan'));   % vc is the climatological meridional velocity

% integrate from the surface

mocs = zeros(360,180,75);  % lon x lat x depth

for ii = 1:360

%     disp(['ii=',num2str(ii)])

    for jj = 1:180

        for kk = 1:75

            if(kk==1)

                v_mid = squeeze(0.5*vc(ii,jj,kk));

                dz = deptht(kk);

                if(~isnan(v_mid))

                    mocs(ii,jj,kk) = v_mid*dz;

                end

            else

                v_mid = squeeze(0.5*(vc(ii,jj,kk)+vc(ii,jj,kkk-1)));

                dz = deptht(kk)-deptht(kk-1);

                if(~isnan(v_mid))

                    mocs(ii,jj,kk) = mocs(ii,jj,kk-1) + v_mid*dz;

                end

            end

        end

        mocs(ii,jj,:) = mocs(ii,jj,:)*cos(lat(jj)*pi/180);

    end

end

mocs = mocs*radius*10^(-6);    % unit: sv

moc_zm = squeeze(mean(mocs,1,'omitnan'));   % zonal mean

zhangn288

楼主解决了吗  该怎么做啊？