Infrared CO₂ cooling of the middle and upper atmosphere is increasing. We developed a new routine for very fast and accurate calculations of this cooling in general circulation models. The new algorithm accounts for non-local thermodynamic equilibrium and is about 1000 times faster than the standard matrix algorithms. It is based on advanced techniques for non-equilibrium emission calculations in stellar atmospheres, which so far have not been used in Earth’s and planetary atmospheres.