传说版主，利用u，v计算风向的代码似乎有bug

海蓝蓝

窃以为，不管什么象限，算风向，下面这一句就够了。
define dire=180+atan2(u,v)*180/pi'

一下是版主在画斜温图的gs里面的一段，应该是计算风向吧。似乎有点问题？
比如，u=10m/s,v= -1m/s,这时候，风向应该是276度，西北风。可是下面代码算出来是366度。
----------------------------------------
'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

bighug

atan2函数的值域是多少？我从网上找到一个fortran函数的日志，说此值域是-π~π，但正常反正切函数的值域应该是-π/2~π/2。对atan2函数不够了解，所以不知道你的答案对不对，我个人认为不对。我是用正常的反正切函数atan(v/u)来算的：
FORTRAN：
if (u>0..and.v>0.) dir=270-atan(v/u)*180/pi
if (u<0..and.v>0.) dir=90-atan(v/u)*180/pi
if (u<0..and.v<0.) dir=90-atan(v/u)*180/pi
if (u>0..and.v<0.) dir=270-atan(v/u)*180/pi
if (u==0..and.v>0.) dir=180
if (u==0..and.v<0.) dir=0
if (u>0..and.v==0.) dir=270
if (u<0..and.v==0.) dir=90
if (u==0..and.v==0.) dir=999
其中uv等于零的五种情况要单独挑出来，不然程序会有瑕疵。atan函数换成atand函数的话直接是度数，不用*180/pi
我四个象限和轴都试了，应该没错。
最需要注意的问题，一个是函数值域，另一个是uv矢量方向和风向是反着的，并且风向角度数是从正Y轴开始顺时针算，和三角函数里度数从正X轴开始逆时针算不一样。

hm_others

我觉得用matlab的话可以这样：
wd=90-atan2(v,u)*180/pi
wd(wd<0)=wd(wd<0)+360

孤蓝et

uv 转风向风速的确需要用4个象限分别计算。
楼主，你可以自己用坐标系画示意图，然后用三角函数计算下。其实上述函数就是实现这个步骤，你可以查查u＞0，v＜0的条件，看看是不是三角函数用错了。

孤蓝et

我以前用VB写过的，忘记了这个是不是最终的版本，楼主可以参考下。当时我就是把每个象限的画出来，然后用三角函数。呵呵。

1. 
   
2.                 temp2(st_nn,3)=mid(read_line(i),20,6)  'U
   
3.                 temp2(st_nn,4)=mid(read_line(i),26,6)  'V
   
4. 
   
5. 
   
6. if temp2(st_nn,3)<0 and temp2(st_nn,4)<0 then '第1象限
   
7. deg=cint(atn(temp2(st_nn,3)/temp2(st_nn,4))*180/3.1415926)
   
8. end if
   
9. if temp2(st_nn,3)<=0 and temp2(st_nn,4)>0 then '第2象限
   
10. deg= 180-cint(atn(abs(temp2(st_nn,3))/abs(temp2(st_nn,4)) )*180/3.1415926)
    
11. end if
    
12. if temp2(st_nn,3)>0 and temp2(st_nn,4)>=0 then '第3象限
    
13. deg=270-cint(atn(abs(temp2(st_nn,4))/abs(temp2(st_nn,3)) )*180/3.1415926)
    
14. end if
    
15. if temp2(st_nn,3)>0 and temp2(st_nn,4)<0 then '第4象限
    
16. deg=cint(atn(abs(temp2(st_nn,4))/abs(temp2(st_nn,3)) )*180/3.1415926)+270
    
17. end if

复制代码

孤蓝et

上面的代码，是2009年写的了。。。现在没有数据了，没办法测试到底对不对。
你就用这个参考下UV转风向风速的思路吧。希望别误导你。

海蓝蓝

已经验证了啊。
比如，u=10m/s, v= - 1m/s,  这时候，风向应该是276度，西北风。可是那个代码算出来是366度!

海蓝蓝

孤蓝et 发表于 2012-1-24 12:30

                               
登录/注册后可看大图

上面的代码，是2009年写的了。。。现在没有数据了，没办法测试到底对不对。
你就用这个参考下UV转风向风速 ...

已经验证了啊。
比如，u=10m/s, v= - 1m/s,  这时候，风向应该是276度，西北风。可是那个代码算出来是366度!

海蓝蓝

孤蓝et 发表于 2012-1-24 12:20

                               
登录/注册后可看大图

uv 转风向风速的确需要用4个象限分别计算。
楼主，你可以自己用坐标系画示意图，然后用三角函数计算下。其 ...

窃以为，不管什么象限，算风向，下面这一句就够了。
define dire=180+atan2(u,v)*180/pi'

孤蓝et

海蓝蓝 发表于 2012-1-24 13:07
窃以为，不管什么象限，算风向，下面这一句就够了。
define dire=180+atan2(u,v)*180/pi'.{:eb302:}真不知道啊

海蓝蓝

孤蓝et 发表于 2012-1-24 21:19

                               
登录/注册后可看大图

海蓝蓝 发表于 2012-1-24 13:07
窃以为，不管什么象限，算风向，下面这一句就够了。
define dire=180+at ...

哎，你试试就知道了。随便举几个例子，一个象限一个。

孤蓝et

海蓝蓝 发表于 2012-1-24 22:05

                               
登录/注册后可看大图

哎，你试试就知道了。随便举几个例子，一个象限一个。

早点发觉这个就好了。。。。。
看来当年数学、几何都没有学好。