trajectory

# SI units

3 # dimensions

# gravity: Earth
# GM=398,600.4418 => G=
newtonian 6.67384e-11 5.97219e24 6371e3 # G M R

# environment
constant 0.0 # rho(r) Function3D
constant 0.0 constant 0.0 constant 0.0 # u_i(r) Function3D[3]

# hard deck
sphere   0.0 0.0 0.0   6371e3

# rocket
constant 0.75 # Cd(u,t): typical model rocket vs x-sec area
constant 5.067e-4 # frontal area (mass,time): 1" dia
constant 0.0 # volume (mass, time)
constant 0.0 # F(t)
constant 1.0e10 # alpha: F = -alpha dm/dt cf F Dt=D(m u_x) => F = u_x dm/dt => alpha = u_x
constant 1.0e8 # max structural pressure (mass,time), ~2 x strength for uniform body
constant 0.0 # ablation rate (fdrag = pdrag * area, gas rho, rocket mass)

# launcher
none

# initial state
0.0   42164e3 0.0 0.0   0.0 1.6e3 0.0  0.10 # t r u m

# integration
5e6 100 # 259200.0 100.0 # dt dtsub

end