if the model grid units, and diameter units, are km, then the depths should be calculated in km too. But imposing the depths on the grid appears to create (even small) craters that are in excess of 50, 100, 200km deep. This is not right - that would be cratering and exposing the mantle! But they are also not 50, 100, 200 m deep (that is too little!). So what is going wrong with the calculation of crater depth? is there a switch in units somewhere? are we using the wrong measures for crater diameter / radial distances (?)

Max Diameters should is a current of the model grid domain, and it should reflect x-length / 2 (or equally y-length /2).... but it currently does not reflect this and the biggest craters look small relative to the model domain.

This issue is probably within the code where the MinD and MaxD are set!

Need to fix "CrateringLesson.ipynb" (jupyter notebook) as it can't run: says the JSON formatting is corrupted. Not sure where the errors are.