【求助讨论】Grads里风向的表达式

红河岸

已知ugrdprs和vgrdprs，求风向的表达式，也就是算出来是度数。

有人是这么做的，就是把ugrdprs和vgrdprs分成正负值来讨论。
这里特别要提醒大家注意一下，grads里的atan2(a,b) 函数跟咱们平常用的不一样。

'define uv=mag(vgrdprs,ugrdprs)' *计算风速
'define pi=3.14159'
if (ugrdprs>=0.)
if (vgrdprs>=0.)
'define dir=180+atan2(ugrdprs,vgrdprs)*180/pi'
endif
endif
if (vgrdprs<=0.)
if (ugrdprs>=0.)
'define dir=270+atan2(ugrdprs,vgrdprs)*180/pi'
endif
endif
if (ugrdprs<=0.)
if (vgrdprs>=0.)
'define dir=180-atan2(ugrdprs,vgrdprs)*180/pi'
endif
endif
if (vgrdprs<0.)
if (ugrdprs<=0.)
'define dir=270-atan2(ugrdprs,vgrdprs)*180/pi'
endif
endif


不过我根据下面这个表达式，觉得不管什么情况，不管纬向、经向风是大于0或者是小于0，dir似乎都可以用一个表达式来用，就是

dir=180+atan2(ugrdprs,vgrdprs)*180/pi'


不知道思路对错与否，与大家探讨！

river

看了一下官方文档，是这么说的
atan2()

atan2 (expr1, expr2)

Applies the tan-1 function to the result of (expr1/expr2).

If expr1 and expr2 are both zero, the result is arbitrarily set to zero. The result of the atan2 function is in radians.

传说中的谁

好像是的，我最初写的时候想了各种情形，后来推敲了一下，一个算式就可以

红河岸

传说中的谁 发表于 2014-1-23 08:38

                               
登录/注册后可看大图

好像是的，我最初写的时候想了各种情形，后来推敲了一下，一个算式就可以

嘿嘿，如果确定一个算式就可以的话，那算起来简单许多了。

pawpaw丁

你好，我想问一下在表示风向和风速时风向的值是怎么定的？
图中的风由径向风分量和纬向风分量合成后指向了东北方向（就是图中的第一象限），风向为何是这种表示方式呢？

红河岸

pawpaw丁 发表于 2014-2-26 10:50
你好，我想问一下在表示风向和风速时风向的值是怎么定的？
图中的风由径向风分量和纬向风分量合成后指向了 ...

我这个风向是度数……

咖啡伴侣

u v都为0 的时候，公式dir=180+atan2(ugrdprs,vgrdprs)*180/pi 计算出来为180°
而u=0，v=0 为静风
我觉得需要特殊处理下

红河岸

咖啡伴侣 发表于 2014-9-10 15:43
u v都为0 的时候，公式dir=180+atan2(ugrdprs,vgrdprs)*180/pi 计算出来为180°
而u=0，v=0 为静风
我觉 ...

怎么处理？

咖啡伴侣

静风处理啊，想把静风的风向设为多少，或者取按前一时刻的风向值