Introduction - If you have any usage issues, please Google them yourself
This is actually a generic varimax routine and knows nothing special about
% EOFs. It expects a matrix of "loadings". Typically (in state space
% rotation), these loadings are the expansion coefficients (aka Principal
% Component Time Series) for the truncated basis of eigenvectors (EOFs), but
% they could also be the EOFs*diag(L)^(1/2) (in the case of rotation in
% sample space).