gl.diagnostics.sim function

Comparing simulations against theoretical expectations