Dear all, we would like to simulate natural convection using the Boussinesq approximation. In the example rayleigh the buoyancy force is implemented in userf, however there are the parameters BETAG and GTHETA, which appear not to be used. Is there any reason for that? Regards, clio