利用fortran推算以北為0度的風場角度

當手上有x方向風速值(u)和y方向風速值(v)時，
可拿來推算風速大小(vel)及風向(theta)。

風速大小(vel)沒什麼問題，
就是使用畢氏定裡便可求得。
vel=sqrt(u**2+v**2)

我本來用arc tan推算風向，
發現忘了考量theta可能會有負值。
但是！！  太久沒玩三角函數了。
天殺的，過了半小時後我沒有任何進度，
我想在別人眼中，我是在發呆吧...  囧

大概是弟弟學長看不下去了，
過來扶了我一把~~ (大感謝弟弟學長~~ QQ)
我才解出來了~
原來用arc cos就不用想那麼多了... >"<

theta=(acos (v / vel))* 180.0/ pi
if (u .lt. 0.0) then
   theta= 360.0 - theta
endif

我要拋棄arc tan，
奔向arc cos....
(明明就是中學數學要重修了...)

sin(x) 與 cos(x)
圖片來源：wiki

tan(x)
圖片來源：wiki