Abstract
Depending on the configuration (geometry, boundary conditions, and physical properties), a solution of the equation set describing thermoconvection within a rigid porous structure may become unstable. The onset of this instability can be determined using linear stability theory. The program presented in this article has been designed to perform this task for complex geological structures whose geometry permits 2D simulation. Two different class of problems are addressed: the onset of free convection and the stability of 2D convective solutions to 3D perturbations. The instability is caused by critical perturbations that can also be computed.