N88-BASICでコリオリ力 (2回目)

2023/10/5(木)
N88-BASICでコリオリ力 (2回目)
 
(Coriolis force)
 
緯度と重力加速度の方向
 
■ 前提
▼ 参照
https://ulprojectmail.blogspot.com/2023/09/coriolis-2.html
コリオリ力 (2回目)
 
▼ 定義
T:地球の自転周期[86164.098903691(s)]
ω:地球の角速度[2π/T (rad/s)]
R_EE:地球赤道半径[6378.1×10³(m)]
R_PE:地球極　半径[6356.8×10³(m)]
G:重力定数[6.67430×10^-11(Nm²/kg²)]
M:地球質量[5.972×10²⁴(kg)]
R:地球の半径[6367.5×10³(m)北緯45°]
g:重力加速度[g = GM/R² ≒ 9.83077(m/s²)]
n :P(x,y)の法線(垂直抗力の方向)
φ:地理緯度
Φ:地心緯度
β:パラメトリック緯度(更成緯度)
γ:天文経度(重力が指す方向との角γ≒φ)
 
原点(x, y, z) = (0, 0, 0)
 
▼ φ,β,Φ
tanφ = (a/b)tanβ = (a²/b²)tanΦ
(b²/a²)tanφ = (b/a)tanβ = tanΦ
 
▼ φ,β,γ
r = √{(a² + b²tan²β)/(1 + tan²β)}
tanβ = (b/a)tanφ
tanγ = {g(b/a) / (g - ω²r)}tanβ
 
■ 解説
緯度と重力加速度の方向の表示

VL,NL,XLとblg～.zip(cori002.bas)は

以下のリンクからダウンロードできます

https://ulprojectmail.blogspot.com

Readme.txtを読んで遊んで下さい