Comparison of model calculated >70 MeV proton fluxes with observations. Left panels: model calculations with g_b function. Right panels: model calculations with g_h function; Upper panels: fluxes along the 1-minute RPS-b trajectory in year 2013; Lower panels: fluxes at POES satellite altitude for each year between 1998–2013; y = x: the perfect fit line (red); Coefficients of determination R² are denoted.