How to modulate the resolution of Julia DifferentialEquations solver

I ran it without tstops but in this case I get the second graph even with t=1000…
But according to this post just using tstops=t_in, works just fine. (The present post was a special case of that post…)