Run pattern scaling on the input temperature fields. This amounts to doing a linear regression on each grid cell, with the global temperature as the sole regressor.

pscl_analyze(tin, tgavin)

Arguments

tin

Matrix of temperature fields (ntime x ngrid)

tgavin

Matrix of global mean temperature over time (ntime x 1)

Value

List of r, w, and b matrices (see details)

Details

The return value is a list with the residuals and pattern coefficients as r (resid), w (slope), and b (intercept).

dim(r) == dim(tin) == ntime x nlat*nlon dim(w) == dim(b) == nlat*nlon