README: description of files in zenodo archive for Lalich et al. 2024. Results file includes results of random generation used in the paper. Included Matlab codes will allow the user to generate and analyze their own random stratigraphies following the methods described in the paper.

archive_results.mat:
 
	All results of simulations/random generation described in the paper. All figures and analysis should be reproducible using these values. Variables are defined as follows:

	PssPs3/4/5: Subsurface to surface power ratios at 3/4/5 MHz for all initial random stratigraphies
	
	PssPs3/4/5_match: each distribution of Pss/Ps ratios at 3/4/5 MHz that matches expected behavior INSIDE bright region

	PssPs3/4/5dimvar: Pss/Ps distributions generated from above stratigraphies corresponding to behavior OUTSIDE bright region

	R3/4/5: effective reflectivity for each initial random stratigraphy

	dim_eps: individual permittivity profiles used to seed distributions outside bright region

	dim_strats: individual layer thickness profiles used to seed distributions outside bright region

	eps_strat: permittivity profiles defining initial random stratigraphies

	epsvar_dim: full permittivity profiles defining distributions outside bright region

	epsvar_match: full permittivity profiles defining distributions inside bright region

	ind_common_match: indices of initial random scenarios that produced distributions matching observed behavior inside bright region

	normpow3/4/5: Observed Pss/Ps ratios INSIDE bright region for 3/4/5 MHz

	normpow3/4/5_dim: Observed Pss/Ps ratios OUTSIDE bright region for 3/4/5 MHz

	strat: layer thickness profiles defining initial random stratigraphies

	stratvar_dim: full layer thickness profiles defining distributions outside bright region
	
	stratvar_match: full layer thickness profiles defining distributions inside bright region

rand_reflectivity.m: creates random stratigraphies, calculates their effective reflectivity, and converts to Pss/Ps for each frequency

generate_strat.m: generates a single random stratigraphy, used by rand_reflectivity.m

marsis_reflectivity.m: calculates effective reflectivity for a given stratigraphy at each MARSIS frequency, used by rand_reflectivity.m

strat_vary_loop.m: creates variant stratigraphies for each profile and determines if resulting distributions match observed behavior inside bright region

strat_vary_analysis.m: creates variants for a single stratigraphy, then calculates Pss/Ps for each frequency, used by strat_vary_loop.m and get_dim_strats.m

get_dim_strats.m: determines which stratigraphies from INSIDE bright region distributions most closely match OUTSIDE bright region, then uses those as seeds to produce new variant distribution representing behavior outside bright region.

match_metric: calculates metric for determining which profiles create reflections most similar to observations at all frequencies, used by get_dim_strats.m

vary_layers.m: creates a single variant on a single stratigraphy, used by strat_vary_analysis.m

Data folder:
	
	bright_area_analysis.m: re-analyzes MARSIS data from Orosei et al. 2018. See the archive associated with that paper for individual data files. Flip the inequality in line 30 to select inside/outside bright region

	datlists.mat: lists of data file names from Orosei et al. 2018 archive

	summary_3/4/5MHz.mat: MARSIS observations from inside bright region used as comparison in paper