為什麼這篇經緯度距離matlab鄉民發文收入到精華區：因為在經緯度距離matlab這個討論話題中，有許多相關的文章在討論，這篇最有參考價值！作者DKer ()看板TY_Research標題[問題] 渦度計算 (經緯度距離修正)時間Thu M...
小弟最近由於研究的需要
必須自己算出整個北半球的渦度變化(relative vorticity)
目前手上已經有NCEP所提供的u、v風速場及渦度場的網格資料
但我自己算出的資料與NCEP提供的渦度場不相符
我利用cos修正NCEP的經度座標距離 A點經度座標=cos(A點緯度)*A點經度
然後直接計算每個網格上u、v風速的curl
在預想中這樣算出的渦度與NCEP應該只有量級(單位)的差異
整體趨勢應該是相同的，不過結果似乎不光是這樣而已~"~
想請教一般在渦度場的計算中，除了距離修正之外，還會進行別的處理嗎？
謝謝 m(_ _)m
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.167.189.233
→ DKer :忘了說 我是用MATLAB算的 05/05 00:22
推 tytony :我在想是否為NCEP的渦度場是直接以經緯度網格curl 05/05 07:07
→ tytony :並沒有修正緯度不同造成的距離差異 05/05 07:08
→ tytony :用NCEP的網格直接curl u,v 比對NCEP渦度場看看 05/05 07:10
推 shenerica :邊界條件問題 05/05 09:14
→ blc :渦度應該沒有邊界條件問題，速度位和流函數才會有。 05/05 09:22
→ blc :ncep資料是從北排到南，看是不是這個問題。 05/05 09:23
→ blc :有些資料… 05/05 09:24
→ tytony :渦度會有邊界問題 ζ= dv/dx - du/dy 05/05 09:39
→ blc :哦哦，我看錯了…不過球座標的話就只有南北極有。 05/05 11:53
→ blc :而且也只會讓渦度場內縮一個格點(中差法的話)。 05/05 11:54
→ DKer :(汗)不做距離修正算出來的真的比較接近NCEP.. 05/05 15:36
→ DKer :不過還是有些差異 請問有沒有其他計算渦度的工具呢 05/05 15:37
→ DKer :我有google到一個GEMPAK 但是對我來說實在太複雜orz 05/05 15:37
→ blc :NCAR commandline language (ncl) 05/06 14:29
→ DKer :萬分感謝樓上 m(_ _)m 05/06 21:47
→ DKer :感謝各位 目前討論的結果是決定用MATLAB不修正距離 05/07 01:32
→ DKer :算出的渦度場 感恩~~ 05/07 01:33