AaronH

Hi Paul, On further investigation I think you're right. I had actually tested truncating the result to N=10 and was getting different results, but I now realise that's due to the way that IGRF interpolates the Gauss coefficients (if you don't specify the IGRF initialization frequency, it defaults to the middle of the year rather than the start of the year). I get values much closer if I emulate this. I can try and modify the code to try and implement the full range of coefficients, although I don't...

I realise I forgot to specify that the above example was calculated for 2015/01/1, but the year used does not seem to matter -- the values are incorrect regardless.

IRBEM IGRF implementation producing incorrect values

For anyone else running into this error, the following will work around the error that Adam mentioned in the previous post without explicitly providing a myownfield_init function: In the file source/init_nouveau.f, comment out line 58, In the file source/onera_desp_lib.f, comment out lines 1267 and 1284 Remake the library. As long as you don't use the "own field" functionality, the library should work without issue now.

Error in conversion between SPH coordinates and RLL coordinates