[原创]MATLAB矢量地理数据相关计算（涡度、散度、平流...）

hhh走走看看

牛人，谢谢分享

最后一个匈奴

感谢楼主分享

子达

请教一下楼主那个涡度函数第七行：“curlz=dv./dx-du./dy+u.*tan(latitude.*pi./180)./R;” 后面这个“+u.*tan(latitude.*pi./180)./R;”是什么意思？涡度函数不是只要前面“curlz=dv./dx-du./dy”就行了？

沙颖凯

子达 发表于 2016-4-6 08:15
请教一下楼主那个涡度函数第七行：“curlz=dv./dx-du./dy+u.*tan(latitude.*pi./180)./R;” 后面这个“+u.* ...

地球曲率的矫正，计算大尺度时不能忽略。

子达

沙颖凯 发表于 2016-4-10 11:59
地球曲率的矫正，计算大尺度时不能忽略。

原来如此，另外请教一下此校正的出处是哪里，有什么文献我可以参考一下？

沙颖凯

子达 发表于 2016-4-10 23:18
原来如此，另外请教一下此校正的出处是哪里，有什么文献我可以参考一下？

搜“球坐标系旋度公式”。
具体是把
x = Rsin(pi/2-lat)cos(lon)
y = Rsin(pi/2-lat)sin(lon)
z = Rcos(pi/2-lat)
带入 curl = delta \times V 同时对自变量和方向矢量求导。

子达

沙颖凯 发表于 2016-4-17 02:07
搜“球坐标系旋度公式”。
具体是把
x = Rsin(pi/2-lat)cos(lon)

谢谢

aas314

给力  大赞

就要唯一

可以，好用，图画得出来。

gtouc

楼主请问平流项计算中advh=-(u.*dTx./dx+v.*dTy./dy);为什么会有一个负号？平流量不都是(u.*dTx./dx+v.*dTy./dy)么？还有就是如图所示，用楼主的程序计算出的风应力旋度和网上下载的风应力旋度正负值以及高低值区对应比较好，只是数值上偏差不小，请问这是问什么？