求助：用NCL根据U风V风如何计算风向

chennizi · 发表于 2016-8-11 12:16:12

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 chennizi 于 2016-8-18 08:15 编辑

求助：用NCL根据U风V风如何计算风向

已有格点的U V风，求该点的风向风速数据，不是要画图，是要求出数据

6.1.0以后的版本提供wind_speed和wind_direction函数计算风速和风向，但是我运行程序提示不存在这两个函数，于是只能自己计算，
计算风速应该就是用sqrt(u^2+v^2)
那么请问计算风向应该用什么公式呢？

1649518749 · 发表于 2016-8-11 15:34:28

wind=(u*u+v*v)^(0.5)

松下风 · 发表于 2016-8-11 15:58:44

风向应该是通过arcsin（u/wind）计算

chennizi · 发表于 2016-8-18 08:13:49

松下风发表于 2016-8-11 15:58
风向应该是通过arcsin（u/wind）计算

感谢，已经测试过可以通过以下两种方式计算
1\ wdir_p = wind_direction(uu10_p,vv10_p,0)

2\ PI = 3.14159265358979323846
tmp = 270.0-atan2(v,u)*180.0/PI
wdir_p = mod(tmp,360.0)
你这个方法我没有测试

zly4814624 · 发表于 2016-8-18 08:45:10

u10 = wrf_user_getvar(fl,"u10",-1)
v10 = wrf_user_getvar(fl,"v10",-1)
wd10m = atan2(u10,v10)/0.01745329 +180

用这个。

chennizi · 发表于 2016-8-18 10:54:47

zly4814624 发表于 2016-8-18 08:45
u10 = wrf_user_getvar(fl,"u10",-1)
v10 = wrf_user_getvar(fl,"v10",-1)
wd10m = atan2(u10,v10)/0.017 ...

看起来公式差不多，但是我没有测试你这个，看着觉得可行的恩恩。那个函数wrf_user_getvar是针对WRF输出数据的函数吗

zly4814624 · 发表于 2016-8-18 10:57:00

chennizi 发表于 2016-8-18 10:54
看起来公式差不多，但是我没有测试你这个，看着觉得可行的恩恩。那个函数wrf_user_getvar是针对WRF输出数 ...

是的。我是用这个做站点时序的。没问题的应该。

beesfly · 发表于 2017-5-10 16:23:40

谢谢各位大神，学习了

梁苏洁 · 发表于 2017-5-11 10:08:56

zly4814624 发表于 2016-8-18 08:45
u10 = wrf_user_getvar(fl,"u10",-1)
v10 = wrf_user_getvar(fl,"v10",-1)
wd10m = atan2(u10,v10)/0.017 ...

学习了

小肥妞 · 发表于 2017-11-19 16:07:04

终于找到学习了

		自动登录	找回密码
密码			立即注册