座標変換に関するサブルーチンパッケージ. 変換する座標系はルーチン名の末尾2文字に略号であらわされる. ここで扱う座標系の略号は以下の通り.
放物線座標との変換は直角双曲線座標の逆変換である.次元 英語名 略号 直角座標 2 cartesian coordinates C 極座標 2 polar coordinates P 楕円座標 2 elliptic coordinates E 直角双曲線座標 2 rectangular hyperbolic H 双極座標 2 bipolar coordinates B 直角座標 3 cartesian coordinates C 球座標 3 spherical coordinates S 楕円面座標 3 ellipsoidal E coordinates
角度に関する引数の単位はすべてラジアンである.
CT2PC(R,THETA,X,Y) 2次元極座標を直角座標に変換する. CT2CP(X,Y,R,THETA) 直角座標を2次元極座標に変換する. CT2EC(U,V,X,Y) 楕円座標を直角座標に変換する. CT2BC(U,V,X,Y) 双極座標を直角座標に変換する. CT2HC(U,V,X,Y) 直角双曲線座標を直角座標に変換する. CT2CH(X,Y,U,V) 直角座標を直角双曲線座標に変換する. CT3SC(R,THETA,PHI,X,Y,Z) 3次元球面座標を直角座標に変換する. CT3CS(X,Y,Z,R,THETA,PHI) 直角座標を3次元球面座標に変換する.
CR2C(THETA,X0,Y0,X1,Y1) 2次元直角座標を回転する. CR3C(THETA,PHI,PSI,X0,Y0,Z0,X1,Y1,Z1) 3次元直角座標を回転する. CR3S 球面座標を回転する. (THETA,PHI,PSI,THETA0,PHI0,THETA1,PHI1)
2次元極座標と直角座標の変換をする.
CALL CT2PC(R,THETA,X,Y)
CALL CT2CP(X,Y,R,THETA)
X,Y (R) 直角座標. R, (R) 極座標. THETA
x = r (θ) y = r (θ)
楕円座標と直角座標の変換をする.
CALL CT2EC(U,V,X,Y)
X, (R) 直角座標. Y U, (R) 楕円座標. V
x = (u)(v) y = (u)(v)
双極座標と直角座標の変換をする.
CALL CT2BC(U,V,X,Y)
X, (R) 直角座標. Y U, (R) 双極座標. V
x = (v)(v) + (u) y = (u)(v) + (u)
直角双曲線座標と直角座標の変換をする.
CALL CT2HC(U,V,X,Y)
CALL CT2CH(X,Y,U,V)
X, (R) 直角座標. Y U, (R) 直角双曲線座標. V
u = x^2 - y^2 v = 2xy
3次元球面座標と直角座標の変換をする.
CALL CT3SC(R,THETA,PHI,X,Y,Z)
CALL CT3CS(X,Y,Z,R,THETA,PHI)
X,Y,Z (R) 直角座標. R,THETA, (R) 極座標. PHI
x = r (θ)(φ) y = r (θ)(φ) z = r (θ)
2次元直角座標を回転する.
CALL CR2C(THETA,X0,Y0,X1,Y1)
THETA (R) 回転角. X0, (R) 回転前の座標値. Y0 X1, (R) 回転後の座標値. Y1
x_1 = (θ)x_0 + (θ)y_0 y_1 = -(θ)x_0 + (θ)y_0
3次元直角座標を回転する.
CALL CR3C(THETA,PHI,PSI,X0,Y0,Z0,X1,Y1,Z1)
THETA,PHI, (R) Euler の回転角 (θ, φ, ψ) PSI . X0,Y0,Z0 (R) 回転前の座標値. X1,Y1,Z1 (R) 回転後の座標値.
球面座標を回転する.
CALL CR3S(THETA,PHI,PSI,THETA0,PHI0,THETA1,PHI1)
THETA,PHI, (R) Euler の回転角 (θ, φ, ψ) PSI . THETA0,PHI0 (R) 回転前の座標値. THETA1,PHI1 (R) 回転後の座標値.
NUMAGUTI Atusi <a1n@gfdl.gov> Last Modified: Thu Aug 31 13:00:27 EDT 1995