trilinear interpolation works with all 0D,1D and 2D data decomposition, interpolation has become a module for that purpose
(using ixmin instead of ra(1) which is way easier to read). domain decomposition is now stored in ra_table as well as
in a 2D integer array (which is saved to disk)