Good you spot that! Well, we can make it bigger in principle, but it
is good to resolve the smallest scales with many modes. Then I just
picked it to be a power of 2 times kmin, and just a bit more than the
required value, but lower than the maximum by a factor of 2-4. In a
second simulation, you can try making it bigger to see if it has any
effect.