Eliminate logic for particles not at surface
Preload data from netCDF into arrays to make the model steps faster
Add test for traits model code to compare with run.py output