Introduction - If you have any usage issues, please Google them yourself
Code is for 2D incompressible Rayleigh-Bernard convection, using staggered grids, middle difference, box boundary and Boussinesq approximation. The upper boundary is cooled while the lower heated.