search

patternizer / ENSEMBLE_SST

Code to generate the FCDR ensemble used as input to surface temp CDR
Other
0 stars 0 forks source link

Create a 10-member ensemble distribution by sampling the covariance matrix #1

Closed patternizer closed 5 years ago

patternizer commented 5 years ago

While calc_ensemble.py produces the eigenvalues and eigenvectors of the harmonisation covariance matrix for a single AVHRR infrared channel, it is necessary to finalise the mathematical framework that uses the eigen-decomposition to generate 10 perturbed values of each harmonisation parameter in such a way that they ensemble is equi-probable.

patternizer commented 5 years ago

Now using numpy.random.multivariate_normal to sample from covariance matrix using EIV-calculated harmonisation parameters as "best" case mean.

patternizer commented 5 years ago

Method calc_ensemble_draws() calculates the ensemble from the minimum of the Frobenius norm between a draw and the decile values of the distribution of each harmonisation parameter.