Hi Mikael, why don't add the possibility to set the r and h parameter in the get_Bnec function as an array? Currently it is possible to use only a float. This could be useful when I have to calculate the Bnec from model at the r values taken from the Swarm MAG-L products rather than using a single r or h value per computation.