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

2023/10/2(月)
N88-BASICでコリオリ力 (1回目)
 
(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)の法線(垂直抗力の方向)
φ:地理緯度
Φ:地心緯度
β:パラメトリック緯度(更成緯度)
γ:天文経度(重力が指す方向との角γ≒φ)
 
▼ 地球半径
a = R_EE = 6378.1×10³(m) 地球赤道半径
b = R_PE = 6356.8×10³(m) 地球極　半径
 
β = Tan^-1{(b/a)tanφ}
R = √{(acosβ)² + (bsinβ)²}
 
または
tanβ = (b/a)tanφ
R = √{(a² + b²tan²β)/(1 + tan²β)}
 
■ 解説
φを入力して
βと地球半径Rを表示します

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

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

https://ulprojectmail.blogspot.com

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