trajectory

# SI units

3 # dimensions

# gravity
constant 0.0 0.0 -10.0

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

# hard deck
nowhere

# 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 1.77345e-4 # 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   0.0 0.0 100.0   0.0 0.0 0.0  0.10 # t r u m

# integration
100.0 0.1 # dt dtsub

end