Grads求偏导数问题

栤唥de嶶笑

'define pi=3.1415926/180'      *pi即π，取值3.14
'define r=6.371e6'                 *r是地球半径
'define deqt=cdiff(eqt,y)'        *deqt是对eqt是指假相当位温，求cidff中央差分
'define dy=cdiff(lat,y)*pi*r'     *dy是对lat纬度，求cidff中央差分
'define eqt=-deqt/dy'             *求偏导数：-偏eqt/偏y


上述脚本中，计算结果，量级比实际大了10倍，有大神指导问题出在哪里吗？

river

我觉得好像没有问题，大了10倍怎么判断的？

栤唥de嶶笑

river 发表于 2021-1-13 11:53
我觉得好像没有问题，大了10倍怎么判断的？

跟别人发表的文献对比发现的

伽蓝鸟

额，以下两段复制自GrADS官网，
算相对涡度：
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
照着他们的写法，试试看？

栤唥de嶶笑

伽蓝鸟 发表于 2021-1-13 17:39
额，以下两段复制自GrADS官网，
算相对涡度：
define dv = cdiff(v,x)

是照着上面的写法来写的，反复确认了好多次，高值区分布都没问题，就是量级大了10倍左右

kaiofchina

伽蓝鸟 发表于 2021-1-13 17:39
额，以下两段复制自GrADS官网，
算相对涡度：
define dv = cdiff(v,x)

官网上的这2种算法，试过，下边那个没问题。上边算相对涡度的这个，define du = cdiff(u*cos(lat*3.1416/180),y)，这句换成define du = cdiff(u,y)*cos(lat*3.1416/180)就跟下边这种结果一样了，不知道为何？

栤唥de嶶笑

问题已经解决了，是最后乘数错了