matlab计算水汽通量积分

huzengyun

高手，学习中

yuanyuan815

AnneChen 发表于 2016-11-7 22:04
代码里没有权重计算，你指的是什么意思，我不大懂。建议你先看看计算水汽通量的公式原理。

楼主，那dp(dp<0)=0;这个是什么意思呢？执行玩这个语句之后dp值只有0,1两个值了。

做积分的时候是从地面到300hpa的积分，但是数据是离散的，所以想问问楼主您，怎么处理的呢？麻烦了。

AnneChen

yuanyuan815 发表于 2016-11-8 22:21
楼主，那dp(dp

求水汽积分是先计算每一层层高积分计算的，层高是用下层减上层之间的差值。由于高海拔地区如青藏高原地表气压可能小于地面气压数据，因此在用地表气压订正的时候，该地区可能出现负值层高。这行代码用来置零。

yuanyuan815

AnneChen 发表于 2016-11-9 05:04
求水汽积分是先计算每一层层高积分计算的，层高是用下层减上层之间的差值。由于高海拔地区如青藏高原地表 ...

dp（dp<0)=0在MATLAB会判断是否小于0，小于0赋值为0，大于等于0赋值为1，再者微积分的思想是进行微分再积分，直接计算层高作为积分的dp，每个气压之间的差值在100~200多hpa，跟微积分不像，更像是离散的去的计算累积的和。

yuanyuan815

yuanyuan815 发表于 2016-11-9 11:39
dp（dp

我一直纠结与这个问题，不知道楼主是怎么处理的，盼回复，谢谢。

AnneChen

yuanyuan815 发表于 2016-11-9 11:39
dp（dp

你说的对，这应该算是离散求和。但是dp这个判断只会给小于0的赋值为0，大于0的都是原值。

yuanyuan815

AnneChen 发表于 2016-11-9 20:33
你说的对，这应该算是离散求和。但是dp这个判断只会给小于0的赋值为0，大于0的都是原值。

楼主，那你是不是就是用这个代码计算的水汽通量呢？

yuanyuan815

yuanyuan815 发表于 2016-11-11 21:28
楼主，那你是不是就是用这个代码计算的水汽通量呢？

之后有没有什么改动呢，因为离散的这种不知道是不是正确的，所以想请教请教楼主是怎么做的

海小宝

{:eb502:}{:eb502:}{:eb502:}

AnneChen

caihuaxialove 发表于 2016-11-30 10:58
这位大神，请问有没有写过散度的代码？可否分享

没有呢。抱歉