% R6

constants

ve = 3000;

CD = 2.2*10; % 10X at solar maximum.
rho = 4e-10; % at 200 km
rho = 1e-13 % at 600 km
A = 1;
z = 600e3;
v = sqrt(G*M_E/(R_E+z))

% Force of atmos. drag.
FD = 1/2*CD*A*rho*v^2

dmdt = FD/ve % mass flow rate

% 1 year
yr = 60*60*24*365

% Mass flow per year.
m_per_yr = dmdt*yr