progress)

        - Program now issue a reasonable error message before exit

        - Implemented a proper trilinear interpolation of the field

        - operation of defining the maximum bin for fitting was done
        after the sampling, which uses this value leading to seg fault
        in some cases