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

qixiangren34 · 发表于 2014-8-8 09:27:25

到底斜温图gs咋写？讨论结果呢？

xiaolugui · 发表于 2015-6-11 07:20:59

正好遇到要把u,v分量转化为角度，学习了！
强大的论坛,什么疑难问题都能找到答复.

chongzika · 发表于 2015-11-7 21:42:09

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

楼主的窃以为是对的吧

zhksoft · 发表于 2016-4-28 16:26:26

! Copyright (c) 2006 Yongjun ZHENG
! Author : Yongjun ZHENG
! Date : 5/9/2006

module wind_utils
  real(kind=8),parameter :: PI=3.14159265358979323846
contains
  subroutine uv2wind(u,v,spd,dir)
implicit none

! arguments
real,intent(in)  :: u,v
real,intent(out) :: spd,dir
real          :: tmp

spd=sqrt(u*u+v*v)
tmp=270.0-atan2(v,u)*180.0/PI
dir=mod(tmp,360.0)
  end subroutine uv2wind

  subroutine wind2uv(spd,dir,u,v)
implicit none

! arguments
real,intent(in) :: spd,dir
real,intent(out)  :: u,v

! local variables
real             :: tmp

tmp=(270.0-dir)*PI/180.0
u=spd*cos(tmp)
v=spd*sin(tmp)
  end subroutine wind2uv
end module wind_utils
来自：http://blog.chinaunix.net/uid-523040-id-2388503.html?/1799.html

zhksoft · 发表于 2016-4-28 16:26:41

本帖最后由 zhksoft 于 2016-4-28 16:28 编辑

从别的博客找到的：http://blog.chinaunix.net/uid-523040-id-2388503.html?/1799.html

淡水河 · 发表于 2016-7-16 15:24:20

验证了楼主是正确的

dqxiaoding · 发表于 2016-9-10 20:38:30

谢谢楼里各位的分享。

kelly · 发表于 2017-4-10 10:44:20

define dire=180+atan2(u,v)*180/pi

是define dire=180+atan2(v,u)*180/pi' 吧

fgrace · 发表于 2017-5-24 16:53:07

正在计算风向，验证一下

天边的彩虹 · 发表于 2017-6-13 12:44:30

向楼主学习，楼主真厉害。

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

[分享资料] 传说版主，利用u，v计算风向的代码似乎有bug