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

海蓝蓝

也不知道传说版主看到此贴没有

pqman3

海蓝蓝 发表于 2012-1-24 13:07

                               
登录/注册后可看大图

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

佩服得五体投地！

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轴开始逆时针算不一样。

river

@传说中的谁你怎么看？

Forcast

是要好好验证一下再用

传说中的谁

很可能是我算错了，当时只是初略的跟MICAPS的探空对比了一下。

善良猪是我

这个贴很有意思啊

海蓝蓝

楼上，你的回复好有意思啊

红河岸

请教下LZ，如果我有模式输出的u10和v10风场，要插值到站点上（也就是站点要插出一个u，v值），但是我用define dire=180+atan2(u,v)*180/pi' 是算不出来的。

不知道能不能跟LZ讨论一下代码？

海蓝蓝

红河岸 发表于 2013-3-2 12:38

                               
登录/注册后可看大图

请教下LZ，如果我有模式输出的u10和v10风场，要插值到站点上（也就是站点要插出一个u，v值），但是我用defi ...

你这个u和v是插值之后的？
难道是插值除了问题