READ ME

The following documents are as follows

RALF_SupplementaryInfo_v3.docx - This document covers the more in depth derivations of the solutions we present in the main text. This document also provides a more exhaustive set of results covering explicitly all of the focus scenarios. We also provide a direct closed form solution for nonlinear binding consider only advective transport of pesticide compounds. Lastly, we also include an analytic solution for the  time-dependent the linear binding scenario.

SYMS.m - This is the symbolic code that is used to get the analytic solution for the nonlinear binding. In particular, the code inverts the matrix required to get the coefficients for the polynomial expression in the main coefficient of the solution.

Variable_movin_averae_NonLin_cleaner - Code runs ensemble simulations for a given focus scenario for different factors for Se coefficient.

ConcNLTempLaer_ss.m - Function that uses the solution for nonlinear binding considering multiple soil layers

f_Temp.m - Temperature dependency on biodegredation

ConcPestNL_ss.m - Core nonlinear binding solution

region_params.m - Sets up parameters for a given focus scenario and calls weather files

SW_top_box.m - Soil water model for estimating effective water flux

get_met_data_input.m - Pulls weather data from weather data file (rain and evaporation)

front_speed.m - Estimates  the amount fo time for compounds to reach the leaching depth based on linear binding

front_speed_nl.m - Estimates  the amount fo time for compounds to reach the leaching depth based on non-linear binding (didn't work well)

MET_Files.zip (Very important) - These are the .met files for the 100 year (repeated) runs of the focus scenarios. Importantly, these are where we take temperature, rainfall, and evaporation.
Note that you will have to unzip and place these in a directory. It should run if you keep the folder together with the other codes. Howver, if you place the directory elsewhere, you will have to change line 18 of region_params.m to correspond to their location. For example,
if you place the met files in C:\home\met_files, you will have to define 
fname = strcat('C:\home\met_files',region,'-M.met');