code updated with Sigma interpolation.
Bug in the Weiss Field update solved by changing the
internal loop over time-indices and making use of pack_index_tri.
The sigma interpolation needs to be handles with some care, it can produces
problematic results when iterated.