Before, the magnetic footpoint at the Sun had been calculated for a radial distance of zero, which is unphysical (though it doesn't make a big difference in the outcome). Now introduced is an (hidden) option to define to which radial distance (in solar radii) the backmapping should be carried out:
sm = SolarMACH(date, bodies, vsw, ref_long, ref_lat, coord_sys, target_solar_radius=2.5)
This is by default 1, so the radius of the Sun. Setting it to 0 provides the behaviour so far; 2.5 would correspond to the standard source surface in a PFSS model.
Before, the magnetic footpoint at the Sun had been calculated for a radial distance of zero, which is unphysical (though it doesn't make a big difference in the outcome). Now introduced is an (hidden) option to define to which radial distance (in solar radii) the backmapping should be carried out:
This is by default
1, so the radius of the Sun. Setting it to 0 provides the behaviour so far; 2.5 would correspond to the standard source surface in a PFSS model.