Refactor if statements
Add diffusivity to output
Move reset of species grid to 0 to fitSpeciesMassOnGrid
Fix, pass simulator to species fit
TODO FIX remove call to calculateComputationDomainGrids before agentgrid is defined