Galerkin finite element discretizations for atmospheric modeling often require the solution of ill-conditioned, saddle point equations which can be efficiently solved using a hybridized method. By extending Firedrake's domain-specific abstraction, we provide a mechanism for the rapid implementation of hybridization methods for a wide class of methods. In this paper, we show that hybridization is an effective alternative to traditional block solvers for simulating geophysical flows.