Room for improvements

• Finer grid to resolve temperature gradients
• Greater extent of box
• Longer time sequences to cover thermal and growths time scale
• Faster (simpler and/or implicit) radiation transport
• Double precision
• Sphericity effects, 3D,...

Aim

• Comparison with 2C and other 2D models
• Step towards complete coverage of cool star surface convection zones