CORBAM is a python package for Gaussian process based modeling of the archeomagnic field. It contains routines for hyperparameter optimization (model training), inversion and hyperparameter marginalization.
You can find a preprint of the related publication here.