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

qixiangren34

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

xiaolugui

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

chongzika

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

楼主的窃以为是对的吧

zhksoft

! 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

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

淡水河

验证了楼主是正确的

dqxiaoding

谢谢楼里各位的分享。

kelly

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

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

fgrace

正在计算风向，验证一下

天边的彩虹

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