The presence of strong magnetic fields in neutron stars, such as in magnetars, may significantly affect their crust-core transition properties and the crust size. This knowledge is crucial in the correct interpretation of astrophysical phenomena involving magnetars, such as glitches in observed rotation frequencies, cooling, bursts and possibly tidal polarizabilities. A recently developed metamodelling technique allows exploring the model dependence of density functional theory equation of state calculations. In this work, we extend this metamodel to investigate the effect of strong magnetic fields on spinodal instabilities of neutron star matter and the associated crust-core properties. Both Tolman-Oppenheimer-Volkov and a full self-consistent numerical calculations are performed for the neutron star structure, the results being quantitatively different for strong magnetic fields. * dchatterjee@lpccaen.in2p3.fr