[原创]MATLAB流函数与势函数的数值计算程序

yeah...

请问这个可以计算区域的吗

yeah...

我看很多都是计算全球的 但是计算区域的话会出现很大的问题 清问楼主您的这个程序计算边界可行吗

一滴雨

用附件中的程序计算得到的辐散风全为nan是什么原因呀？程序如下，拜托大神帮忙解决下吧
clc;clear;clc;
% 88888888888888888888888888读取u8888888888888888888888888888888888
fileu='F:\omega\uwnd.mon.mean.194801to201612.nc';
% ncdisp(fileu);
lonn=ncread(fileu,'lon');
latt=ncread(fileu,'lat');
lonsta=find(lonn==0);
lonend=find(lonn==357.5);
latsta=find(latt==90);
latend=find(latt==-90);
nnx=lonend-lonsta+1;
nny=latend-latsta+1;
lon=lonn(lonsta:lonend);
lat=latt(latend:-1:latsta);
clear lonn;
clear latt;
t=ncread(fileu,'time',541,264);
level=ncread(fileu,'level',1,10);
level1=level;
uzong1=ncread(fileu,'uwnd',[lonsta latsta 1 541],[nnx nny 10 264]);
for j=1:nny
    uzong(:,j,:,:)=uzong1(:,nny+1-j,:,:);
end
clear fileu;
clear uzong1;
% 88888888888888888888888888读取v8888888888888888888888888888888888
filev='F:\omega\vwnd.mon.mean.194801to201612.nc';
vzong1=ncread(filev,'vwnd',[lonsta latsta 1 541],[nnx nny 10 264]);
for j=1:nny
    vzong(:,j,:,:)=vzong1(:,nny+1-j,:,:);
end
clear filev;
clear vzong1;
[lat,lon]=meshgrid(lat,lon);
for k=1:264
    for j=1:10
        u(:,:)=uzong(:,:,j,k);
        v(:,:)=vzong(:,:,j,k);
        [chi,Uchi,Vchi]=chi_potential(lon, lat, u, v);
        udiv(:,:,j,k)=Uchi(:,:);
        vdiv(:,:,j,k)=Vchi(:,:);
        shi(:,:,j,k)=chi(:,:);
    end
end

lqe2004

谢谢分享！正好需要呢

雪泉

还是牛人多啊

lqe2004

多谢楼主！下下来好好研究一下。提醒一下想用这个程序的海洋专业小朋友记得在准备数据的时候把u和v的岛屿和陆地都设为零，然后得到的psi要把岛屿和陆地设为nan。

kc121

计算部分是有了，画图部分和数据读取部分楼主都省了。让人还是看不懂。

路痴草莓

感谢楼主分享

liujimin025

感谢分享，楼主好人啊

轻拢烟云

大佬谢谢分享