Orbital debris soliton simulation code
Alexis Truitt
May 2020
Contact info: alexis.truitt@gmail.com, hartzell@umd.edu

Conditions of Use:
If you use this code for your academic work, please cite this DRUM doi and this JSR journal paper: https://doi.org/10.2514/1.A34652. 

References:
JSR Paper 1: https://doi.org/10.2514/1.A34652
(Note Paper 1 figure 5 and figure 6 have the wrong captions assigned, were reversed in proofing/editing.
Read the paragraph for the description. I submitted a correction).
JSR Paper 2: https://doi.org/10.2514/1.A34674  

Code:

plotgen_undamped - reproduces simulations from the first NIAC paper and JSR Paper 1
plotgen_damped - reproduces simulation of Figure 5 from the damping paper (JSR Paper #2)

orbit_results - identification of orbital locations where solitons can be created
dist_results - propagation distances including landau damping for orbit_results file

Functions called:

inmode - defines the domain: N,L,x,k,t,dt for 1d simulations, can modify if needed
inmode2d - for 2d simulations

CCMC.mat - file with saved data from the IRI model, as described in Paper #1
ccmc_params - interpolates the CCMC.mat file and outputs plasma parameters for a given orbital location
ccmc_output_dlat5_dalt50_maxalt2000.mat - saved plasma parameters from IRI for given plasmamodes (local times/seasons), delta altitude = 50 km, delta latitude = 5 degrees, for faster global simulations

jdera - converts from Julian date to earth rotation angle for the given local times/seasons
wkspc-ERAdeg.mat - saved file converting julian date to earth rotation angle for given conditions
seztoijk - converts from SEZ to IJK reference frame

sol_force - define forcing function for 1D undamped simulations
sol_force_d - damped
sol_force_d2d - 2D

sol_steps - propagation of debris location for 1D undamped simulations
sol_steps_d - 1D damped
sol_steps_d2d - 2D damped

scale_hydrogen - debris surface potential in hydrogen dominated plasma, cite Bibhas
scale_oxygen - debris surface potential in oxygen dominated plasma, cite Bibhas

trendfit.mat - soliton paramters (amplitude, width, speed, generation interval) based on simulations, 21 debris sizes and 16 velocities
trendinterp.mat - interpolation of trendfit for faster runs
trendgen - simulation of 2d damped results for a series of debris sizes and speeds, for input into a trend analysis of soliton properties, for JSR Paper 3

loop_alt - JSR Fig. 12: As a function of debris size, the percent of grid points where the
median and maximum change in Debye length is greater than 1%, assuming radial propagation. 
loop_lat - same as altloop but for latitudinal propagation.






