| Class | wa_base_module_sjpack |
| In: |
libsrc/wa_module_sjpack/wa_base_module_sjpack.f90
|
| Authors: | Shin-ichi Takehiro, Youhei SASAKI |
| Version: | $Id: wa_base_module_sjpack.f90 590 2013-08-19 08:48:21Z uwabami $ |
| Copyright&License: | See COPYRIGHT |
spml/wa_base_module_sjpack モジュールは球面上での流体運動を 球面調和函数を用いたスペクトル法によって数値計算するための モジュール wa_module_sjpack の下部モジュールであり, スペクトル計算の 基本的な Fortran90 関数を提供する.
球面上の 1 層モデル用 w_base_module_sjpack モジュールを多層モデル用に 拡張したものであり, 同時に複数個のスペクトルデータ, 格子点データに 対する変換が行える.
内部で ISPACK の SJPACK の Fortran77 サブルーチンを呼んでいる. スペクトルデータおよび格子点データの格納方法や変換の詳しい計算法に ついては ISPACK/SJPACK のマニュアルを参照されたい.
このモジュールを使うためには前もって w_base_initial を呼んで 切断波数, 格子点数の設定をしておく必要がある.
| Variable : | |||
| km =16 : | integer
|
| Subroutine : | |||
| wa_Psi(:,:) : | real(8), intent(in)
| ||
| wa_Chi(:,:) : | real(8), intent(in)
| ||
| xya_U(0:im-1,1:jm,size(wa_Psi,2)) : | real(8), intent(out)
| ||
| xya_V(0:im-1,1:jm,size(wa_Psi,2)) : | real(8), intent(out)
|
流線・ポテンシャル(スペクトルデータ)から速度場(格子データ)に (逆)変換する(多層用)
スペクトル変換を用いず微分を計算するために, 変換回数が 2 回ですむ.
u cosφ = ∂χ/∂λ - cosφ∂ψ/∂φ, v cosφ = cosφ∂χ/∂φ + ∂ψ/∂λ
subroutine wa_StreamPotential2Vector(wa_Psi, wa_Chi, xya_U, xya_V)
!
! 流線・ポテンシャル(スペクトルデータ)から速度場(格子データ)に
! (逆)変換する(多層用)
!
! スペクトル変換を用いず微分を計算するために, 変換回数が 2 回ですむ.
!
! u cosφ = ∂χ/∂λ - cosφ∂ψ/∂φ,
! v cosφ = cosφ∂χ/∂φ + ∂ψ/∂λ
!
real(8), intent(in) :: wa_Psi(:,:)
!(in) 流線関数((nm+1)*(nm+1),:)
real(8), intent(in) :: wa_Chi(:,:)
!(in) 速度ポテンシャル((nm+1)*(nm+1),:)
real(8), intent(out) :: xya_U(0:im-1,1:jm,size(wa_Psi,2))
!(out) 速度経度成分(0:im-1,1:jm,:)
real(8), intent(out) :: xya_V(0:im-1,1:jm,size(wa_Psi,2))
!(out) 速度緯度成分(0:im-1,1:jm,:)
integer :: k
do k=1,size(wa_Psi,2)
call w_StreamPotential2Vector ( wa_Psi(:,k), wa_Chi(:,k), xya_U(:,:,k), xya_V(:,:,k) )
enddo
end subroutine wa_StreamPotential2Vector
| Subroutine : | |||
| xya_U(0:,:,:) : | real(8), intent(in)
| ||
| xya_V(0:,:,:) : | real(8), intent(in)
| ||
| wa_Vor((nm+1)*(nm+1),size(xya_U,3)) : | real(8), intent(out)
| ||
| wa_Div((nm+1)*(nm+1),size(xya_U,3)) : | real(8), intent(out)
|
速度場(格子データ)から渦度・発散(スペクトルデータ)に (正)変換する(多層用)
スペクトル変換を用いず微分を計算するために, 変換回数が 2 回ですむ.
ζ = 1/cosφ∂v/∂λ - 1/cosφ ∂(u cosφ)/∂φ D = 1/cosφ∂u/∂λ + 1/cosφ ∂(v cosφ)/∂φ
subroutine wa_Vector2VorDiv(xya_U, xya_V, wa_Vor, wa_Div)
!
! 速度場(格子データ)から渦度・発散(スペクトルデータ)に
! (正)変換する(多層用)
!
! スペクトル変換を用いず微分を計算するために, 変換回数が 2 回ですむ.
!
! ζ = 1/cosφ∂v/∂λ - 1/cosφ ∂(u cosφ)/∂φ
! D = 1/cosφ∂u/∂λ + 1/cosφ ∂(v cosφ)/∂φ
!
real(8), intent(in) :: xya_U(0:,:,:)
!(in) 速度経度成分(0:im-1,1:jm,:)
real(8), intent(in) :: xya_V(0:,:,:)
!(in) 速度緯度成分(0:im-1,1:jm,:)
real(8), intent(out) :: wa_Vor((nm+1)*(nm+1),size(xya_U,3))
!(out) 流線関数
real(8), intent(out) :: wa_Div((nm+1)*(nm+1),size(xya_U,3))
!(out) 速度ポテンシャル
integer :: k
do k=1,size(xya_U,3)
call w_Vector2VorDiv (xya_U(:,:,k), xya_V(:,:,k), wa_Vor(:,k), wa_Div(:,k))
enddo
end subroutine wa_Vector2VorDiv
| Subroutine : | |||
| xya_UCosLat(0:,:,:) : | real(8), intent(in)
| ||
| xya_VCosLat(0:,:,:) : | real(8), intent(in)
| ||
| wa_Vor((nm+1)*(nm+1),size(xya_UCosLat,3)) : | real(8), intent(out)
| ||
| wa_Div((nm+1)*(nm+1),size(xya_UCosLat,3)) : | real(8), intent(out)
|
速度場(格子データ)から渦度・発散(スペクトルデータ)に (正)変換する(多層用)
スペクトル変換を用いず微分を計算するために, 変換回数が 2 回ですむ.
ζ = 1/cosφ∂v/∂λ - 1/cosφ ∂(u cosφ)/∂φ D = 1/cosφ∂u/∂λ + 1/cosφ ∂(v cosφ)/∂φ
subroutine wa_VectorCosLat2VorDiv(xya_UCosLat, xya_VCosLat, wa_Vor, wa_Div)
!
! 速度場(格子データ)から渦度・発散(スペクトルデータ)に
! (正)変換する(多層用)
!
! スペクトル変換を用いず微分を計算するために, 変換回数が 2 回ですむ.
!
! ζ = 1/cosφ∂v/∂λ - 1/cosφ ∂(u cosφ)/∂φ
! D = 1/cosφ∂u/∂λ + 1/cosφ ∂(v cosφ)/∂φ
!
real(8), intent(in) :: xya_UCosLat(0:,:,:)
!(in) 速度経度成分 * cos(lat) (0:im-1,1:jm,:)
real(8), intent(in) :: xya_VCosLat(0:,:,:)
!(in) 速度緯度成分 * cos(lat) (0:im-1,1:jm,:)
real(8), intent(out) :: wa_Vor((nm+1)*(nm+1),size(xya_UCosLat,3))
!(out) 流線関数
real(8), intent(out) :: wa_Div((nm+1)*(nm+1),size(xya_UCosLat,3))
!(out) 速度ポテンシャル
integer :: k
do k=1,size(xya_UCosLat,3)
call w_VectorCosLat2VorDiv (xya_UCosLat(:,:,k), xya_VCosLat(:,:,k), wa_Vor(:,k), wa_Div(:,k))
enddo
end subroutine wa_VectorCosLat2VorDiv
| Subroutine : |
このサブルーチンを単独で用いるのでなく, 上位サブルーチン wa_Finalize を使用すること.
subroutine wa_base_Finalize
!
!
! このサブルーチンを単独で用いるのでなく,
! 上位サブルーチン wa_Finalize を使用すること.
!
call MessageNotify('M','wa_base_finalize', 'wa_base_module_sjpack is finalized (dummy, 2013/02/23) ')
end subroutine wa_base_Finalize
| Subroutine : | |||
| k_in : | integer,intent(in)
|
SNPACK 用 wa_base_initial の互換のためのダミーサブルーチン
subroutine wa_base_initial(k_in)
!
! SNPACK 用 wa_base_initial の互換のためのダミーサブルーチン
!
integer,intent(in) :: k_in !(in) 最大データ数(層数)
km = k_in
call MessageNotify('M','wa_base_initial', 'No need to set maximum level number and in wa_base_module_sjpack (2013/02/23) ')
end subroutine wa_base_Initial
| Function : | |||
| wa_xya((nm+1)*(nm+1),size(xya_data,3)) : | real(8)
| ||
| xya_data(0:,:,:) : | real(8), intent(in)
| ||
| ipow : | integer, intent(in), optional
| ||
| iflag : | integer, intent(in), optional
|
格子点 -> 球面調和関数スペクトル
格子データからスペクトルデータへ(正)変換する(多層用).
function wa_xya(xya_data,ipow,iflag) ! 格子点 -> 球面調和関数スペクトル
!
! 格子データからスペクトルデータへ(正)変換する(多層用).
!
real(8), intent(in) :: xya_data(0:,:,:)
!(in) 格子点データ(0:im-1,1:jm,:)
real(8) :: wa_xya((nm+1)*(nm+1),size(xya_data,3))
!(out) スペクトルデータ((nm+1)*(nm+1),:)
integer, intent(in), optional :: ipow
!(in) 変換時に同時に作用させる 1/cosφ の次数. 省略時は 0.
integer, intent(in), optional :: iflag
! 変換の種類
! 0 : 通常の正変換
! -1 : 経度微分を作用させた正変換
! 1 : 緯度微分 1/cosφ・∂(f cos^2φ)/∂φ を作用させた正変換
! 2 : sinφを作用させた正変換
! 省略時は 0.
!
integer, parameter :: ipow_default = 0 ! スイッチデフォルト値
integer, parameter :: iflag_default = 0 ! スイッチデフォルト値
integer ipval, ifval
integer k
if (present(ipow)) then
ipval = ipow
else
ipval = ipow_default
endif
if (present(iflag)) then
ifval = iflag
else
ifval = iflag_default
endif
do k=1,size(xya_data,3)
wa_xya(:,k) = w_xy(xya_data(:,:,k),iflag=ifval,ipow=ipval)
enddo
end function wa_xya
| Function : | |||
| xya_wa(0:im-1,1:jm,size(wa_data,2)) : | real(8)
| ||
| wa_data(:,:) : | real(8), intent(in)
| ||
| ipow : | integer, intent(in), optional
| ||
| iflag : | integer, intent(in), optional
|
球面調和関数スペクトル -> 格子点
スペクトルデータから格子データへ変換する(多層用).
function xya_wa(wa_data,ipow,iflag) ! 球面調和関数スペクトル -> 格子点
!
! スペクトルデータから格子データへ変換する(多層用).
!
real(8), intent(in) :: wa_data(:,:)
!(in) スペクトルデータ((nm+1)*(nm+1),:)
!
real(8) :: xya_wa(0:im-1,1:jm,size(wa_data,2))
!(out) 格子点データ(0:im-1,1:jm,:)
!
integer, intent(in), optional :: ipow
!(in) 作用させる 1/cosφ の次数. 省略時は 0.
integer, intent(in), optional :: iflag
!(in) 変換の種類
! 0 : 通常の正変換
! -1 : 経度微分を作用させた逆変換
! 1 : 緯度微分 cosφ・∂/∂φ を作用させた逆変換
! 2 : sinφを作用させた逆変換
! 省略時は 0.
!
integer, parameter :: ipow_default = 0
integer, parameter :: iflag_default = 0
integer ipval, ifval
integer k
if (present(ipow)) then
ipval = ipow
else
ipval = ipow_default
endif
if (present(iflag)) then
ifval = iflag
else
ifval = iflag_default
endif
do k=1,size(wa_data,2)
xya_wa(:,:,k) = xy_w(wa_data(:,k),iflag=ifval,ipow=ipval)
enddo
end function xya_wa