For each ocean regime, surface AT results are fitted to a
second-order polynomial in SSS and SST and are also forced
to satisfy equations for adjacent regime(s) which share
boundaries with the regime of interest. The use of this
functional form of equation is to ensure a smooth transition
in AT between adjacent regimes, and thus to minimize errors
in AT predictions at their mutual boundaries