半拉格朗日法平流内插的计算过程研究.pdfVIP

半拉格朗日法平流內插的計算過程研究 張博雄 1 曾忠一 2 1 中央氣象局 2 中央研究院物理研究所 發點必須運用到 O(p 3 ) 次計算。如果是 摘要 三次內插函數，那麼p就等於 4，對於 半拉格朗日平流格式近來普遍受 更高次的內插函數而言 ，計算次數將 到採用因為這個方法有穩定， 、準確 、 呈現倍數成長 。Purser and Leslie 且有效率的優點。但是由於出發點通 (1991)提出可以改變數值計算的概 常不在固定的網格點 ，平流過程必須 念，透過循序的一維內插過程 ，進行 使用大量的內插計算 ，這個過程佔去 半拉格朗日平流內插的計算。他們將 大量的計算時間，而且內插格式越高 複雜的高次內插過程簡化到只需 O(p ) 階則計算內插的時間也越多 。透過串 次計算。 聯內插方法 (cascadeinterpolation) Purser and Leslie (1991)所提 的計算 ，可以縮短內插計算的時間 。 出的概念，首先定義一組固定網格與 一般的無輻散流場平流過程顯示串聯 拉格朗日網格 ，其中拉格朗日網格式 內插格式是省時且結果接近於傳統方 由出發點所組成 。這兩個網格系統的 法的格式。雲模式模擬顯示 ，串聯內 交點為中繼點 （intermediate 插結果與其他結果相近 。在這實驗中 point），通常位於固定網格上 ，利用 我們也提出三度空間的串聯內插流 一般的內插方法可以很快的得到中繼 程 ，結果顯示 ，串聯內插可以大幅減 點的值 。而通常中繼點又位於出發點 少三度空間的內插計算，且計算結果 所組成的曲線上 ，透過延曲線的內插 甚至較單調格式好。 就可以得到出發點的值 ，如此 ，一個 關鍵字 ：半拉格朗日法 ，串聯內插 ， 出發點僅需兩次的計算即可。當然 ， 雲模式 這些內插格式可以應用線性 、二次 、 三次或更高次格式來進行。Nair 等人 一、前言 (1999)應用單調格式與串聯內插方 半拉格朗日法近年來成為受注目 法 ，並且提出二維座標轉換的計算過 的數值方法 ，在許多模式中 ，半拉格 程 ，顯示守恆格式可以輕易的加入 。 朗日法表現出極佳的計算速率 ，但是 在串聯內插的方法中 ，求得出發點的 內插的過程仍然是限制半拉格朗日法 值必須利用到固定網格與拉格朗日網 快速計算的因素 。以三度空間計算過 格的交點 。若交點過於集中 ，或時步 程的例子來說 ，每個預報變數每個出 過大造成出發點曲線相交 ，則串聯內 1 插可能無法計算 ，而必須採用傳統方 過程利用 INTEL Fortran 編譯器 法求解。更有可能是由於 overshoot (LINUX版 )進行編譯 ，故計算時間以 或 undershoot造成的流場扭曲 ，使得 Pentium4為主 。 串聯內插無法計算 。 1. 二維無輻散風場及變形流場平流 本研究將採用無輻散流場測試串 測試及單調格式 聯內插格式的穩定性 、保守性及效 我們利用線性平流方程來平流一 率 ，對於串聯內插