Wednesday, November 11, 2009

以歐中預報溫度

時間關係, 預報方程只以厚度(500-1000)為單一參數. 把各資料點劃了出來:


由於歐中沒有 1000 高度, 故要以 slp 代替. 好在 slp 和 h1000 基本上是一樣的(相關系數非常接近 1). 預報方程為:
接著下載今早歐中 00Z 的資料(約在香港時間每天下午四時左右有). 以線性插值(可用更高 order)求出京士柏(~ 22.3 N; 114.2 E)的預測 slp 和 h500, 再代入上面方程.

3-7 天早上八時的預報溫度為 18.6, 18.7, 19.8, 17.5, 16.7, 跟官方數字很大出入. 誤差可能有三, 包括系統性, 有限資料, 或方程不適合極端值. 而天文台預報的最低溫是否就指 00Z 也是疑問.
另外, 雖然相關系數高, 但按照上圖, 誤差仍然可以達兩度. 下回再以京士柏 925-1000 hPa 厚度配合 GFS 看看能否改善.

No comments: