A significant source of error in orbit modeling of Global Navigation Satellite System (GNSS) satellites is the effect of radiation pressure, both from the sun, and from the Earth. The new Empirical CODE Orbit Model (ECOM) of the Center for Orbit Determination in Europe (CODE), introduced in 2014, attempts to empirically account for the Solar Radiation Pressure (SRP) effects. No a priori effect, e.g., from a boxwing model, is applied because any error, uncertainty, or significant simplification has a comparably large effect on the SRP corrections for the orbits. For that reason we have used scaling factors to validate and adapt the parameters of the boxwing models. The scaling factors are computed from a long-term solution and applied when computing SRP and PRP effects. Reliable satellite characteristics have only been released for Galileo and QZSS, which can be used as a reference. For the other GNSS satellites, the true satellite models remain unknown, so the scaling factors have a bigger impact on the orbits. This study seeks to characterize box-wing models not only for GPS and GLONASS, but also for other GNSS systems. Due to their higher surface-to-mass ratio, it is expected that the adjustable boxwing model will also improve orbit modeling of Galileo satellites, even though corresponding models, including their surface properties, have been disclosed by their agencies. Additionally, we analyze improvements achieved by using the rescaled boxwing models, e.g., by evaluating the obtained ECOM parameters and with other usual orbit validation methods (satellite clock performance and SLR residuals as far as possible).