cdiff的用法

278803532

麦田_smile 发表于 2014-5-8 15:23
请教下 关于4楼给的例子 为什么前后不一致？有6楼那样的疑惑

因为每个纬圈上的经向距离是不一样的，所以需要乘以后面的那个东西，不知道这样说明白不。

麦田_smile

278803532 发表于 2014-5-8 19:54
因为每个纬圈上的经向距离是不一样的，所以需要乘以后面的那个东西，不知道这样说明白不。

谢谢 这一点我也想到了。但为什么两个例子乘的位置不一样。第一个
define dv = cdiff(v,x)
define dx = cdiff(lon,x)*3.1416/180
define du = cdiff(u*cos(lat*3.1416/180),y)
define dy = cdiff(lat,y)*3.1416/180
display (dv/dx-du/dy)/(6.37e6*cos(lat*3.1416/180))
第二个
define dtx = cdiff(t,x)
define dty = cdiff(t,y)
define dx = cdiff(lon,x)*3.1416/180
define dy = cdiff(lat,y)*3.1416/180
display -1*( (u*dtx)/(cos(lat*3.1416/180)*dx) + v*dty/dy )/6.37e6

麦田_smile

278803532 发表于 2014-5-8 19:54
因为每个纬圈上的经向距离是不一样的，所以需要乘以后面的那个东西，不知道这样说明白不。

我现在要算这两项

                               
登录/注册后可看大图

我是这样计算的：
第一项：
'define dv=cdiff(vor,x)'
'define dx=cdiff(lon,x)*3.1416/180'
'd (dv/dx)/6.37e6'

第二项：
'define dv=cdiff(vor,y)'
'define dy=cdiff(lat,y)*3.1416/180'
'd (dv/dy)/6.37e6'

按照例子的说法，似乎应该在哪里加上cos(lat*3.1416/180)项，但具体怎么加，我自己比较乱，两个例子好像提供了矛盾的信息。

278803532

这个我还真的没有深入研究过，
我看过看个cidffi的，也有同样的疑问，一个是涡度的，一个是平流的，
你这个是算涡度的，应该是按你第一个算法那样算的吧
这个;display (dv/dx-du/dy)/(6.37e6*cos(lat*3.1416/180))

278803532

看错了，不好意思。，我想到一个办法，就是用涡度和散度那个两个函数算，你这个应该是按平流按个算法算的。
hdivg()+hcurl就是对X偏导，两个函数相减就是对Y偏导

一片落叶惊秋风

楼主的高度那个算错了，位势米和米相差很小，基本上可以用位势米代替米。
9.8*H=g*z
H=g*z/9.8

星火

刚好遇到和楼主差不多的问题，很有启发，在此感谢楼主。

另外官网的例子，其实把涡度公式换算成球坐标下就明白了，例子中为了整洁做了一个通分，所以会多出来一项（感谢我的舍友好兄弟）

dcx

请问楼主垂直方向生如何求偏导，这个问题解决了么

karenlk

siljechen 发表于 2013-7-12 14:31
继续说下我对cdiff的理解，和师妹讨论后，认为要用中央差分函数cdiff来写，应该是写成cdiff(thetase,z)/cdi ...

请问，cdiff对z方向的使用不能够定义的问题最后怎么解决的？

shenqiankai

正在研究这个用法，过来学习