matlab计算水汽通量积分

zhchlue

AnneChen 发表于 2016-5-10 18:58
trapz（）函数是用梯形公式计算数值积分，对矩阵使用trapz（）函数时，相当于将行或列相邻元素相加除以2 ...

对的，中央差分法求面积

蒙面的羊驼

单位是: g/(s*hpa*cm)或者kg/(s*hpa*cm)
我是从：近10年夏季西北地区水汽空间分布和时间变化分析_李江林  这篇论文上看的

zddouc

caihuaxialove 发表于 2016-11-30 10:58
我最近在matlab 写了水汽通量散度的，但是总感觉怪怪的

您好 可否分享一下您的m文件？我最近在做这个，感觉自己做的不太对

qixianger

你好，我觉得你程序的积分不是从地面开始的，你每次的level除去最前面的地面气压其他的都还是从1000hPa开始的，所以有问题，我修改了一下，认为是这个样子的lev1=ncread('shum.mon.mean.nc','level');
ps=ncread('pres.mon.mean.nc','pres');

lps=size(ps);
for k=1:lps(1)
    for kk=1:lps(2)
        for kkk=1:lps(3)
            pss=ps(k,kk,kkk);
            lev = lev1;
            if pss>=0
                lev(lev>pss)=pss;
                lll=[pss;lev];
                dl=lll(1:end-1)-lll(2:end);
                dp(k,kk,:,kkk)=dl; %#ok<SAGROW>
            else
                lll=[1000;lev];
                dl=lll(1:end-1)-lll(2:end);
                dp(k,kk,:,kkk)=dl; %#ok<SAGROW>
            end
        end
    end
end

你可以算一下青藏高原的水汽通量结合文献，看一下你的。如果有什么不对的地方，欢迎指正。

dnh

qixianger 发表于 2021-1-8 21:42
你好，我觉得你程序的积分不是从地面开始的，你每次的level除去最前面的地面气压其他的都还是从1000hPa开始 ...

我把你们的代码都跑了一遍，你的和grads结果最像，主要区别就是青藏高原地区，楼主的结果是水汽通量偏大，自己也不是很懂，但我看文献中的水汽通量图，在青藏高原地区都是很低的

puck66

dnh 发表于 2021-4-13 21:55
我把你们的代码都跑了一遍，你的和grads结果最像，主要区别就是青藏高原地区，楼主的结果是水汽通量偏大 ...

大侠最后搞定Matlab水汽通量及其散度没？

微醉的阳光

qixianger 发表于 2021-1-8 21:42
你好，我觉得你程序的积分不是从地面开始的，你每次的level除去最前面的地面气压其他的都还是从1000hPa开始 ...

赞，优秀

微醉的阳光

楼主，你最后选择sum的方式结果可能不太好，由于层数太少，导致累加与积分结果可能有些误差，建议采用trapz

呆妹小霸王

除1000又是什么意思呀