戻る|進む

DCL:MATH２:SHTLIB : 球面調和関数:サブルーチンの説明

4.3.19 SHTSWZ

1.

機能

スペクトル逆変換の前半部分であるルジャンドル逆変換を 指定された一つの緯度円上で指定された波数区間のみについて行う.

2.

定義

スペクトルデータS^m_nから, ウエーブデータ W^m_(φ） 指定された一つの緯度円φj上で  M1≤|m| ≤M2 の波数範囲のみについて行う. 変換式はSHTS2Wを参照.

3.

呼び出し方法

SHTSWJ(MM,JM,ISW,J,M1,M2,S,WJ,WORK)

4.

パラメーターの説明

MM	(I)	入力. 切断波数(M).
JM	(I)	入力. 南北分割数の1/2(J)
ISW	(I)	入力. 変換の種類の指定( SHTS2Wを参照).
J	(I)	入力. 変換を行う緯度円の指定(備考を参照).
M1	(I)	入力. 変換する波数区間の最小値(上記定義を参照)
M2	(I)	入力. 変換する波数区間の最大値(上記定義を参照)
S	(R)	入力. スペクトルデータ (長さおよび並び方は SHTNMLを参照).
WJ	(R)	出力. Wm(φj）  が格納される. 長さ 2*MM+1の配列(並び方は備考を参照).
WORK	(R)	SHTINTで初期化された作業領域.

5.

備考 

(a)

WJ(-MM:MM)と宣言しておけば, WJ(M) (M>0)には Re(W^m_(φj)) が, WJ(-M) (M>0)には Im(W^m_(φj)) がそれぞれ格納され, WJ(0)にはW⁰_(φj） が格納される (j=J).

(b)

指定された波数区間の外すなわちWJ(M) (|M|<M1, |M|>M2) には0が代入される.