In radio occultation (RO) data processing and data assimilation, the forward model (FM) is used to calculate bending angle (BA) from refractivity (N). The accuracy and precision of forward modeled BA are affected by refractivity profiles and FM methods, including Abel integral algorithms (direct, exp, exp_T, linear) and methods of interpolating refractivity during integral (log-cubic spline and log-linear). Experiment 1 compares these forward model methods by comparing the difference and relative difference (RD) of the experimental value (forward modeled ECMWF analysis) and the true value (BA of FY3D RO data). Results suggested that the exp with log-cubic spline (log-cubic) interpolation is the most accurate FM because it has better integral accuracy (less than 2%) to inputs, especially when the input is lower than an order of magnitude of 1 × 10−2 (that is, above 60 km). By contrast, the direct induced a 10% error, and the improvement of exp T to exp is limited. Experiment 2 simulated the exact errors of an FM (exp) based on inputs on different vertical resolutions. The inputs are refractivity profiles on model levels of three widely used analyses, including ECMWF 4Dvar analysis, final operational global analysis data (FNL), and ERA5. Results demonstrated that based on exp and log-cubic interpolation, BA on model level of ECMWF 4Dvar has the highest accuracy, whose RD is 0.5% between 0–35 km, 4% between 35–58 km, and 1.8% between 58–80 km. By contrast, the other two analyses have low accuracy. This paper paves the way to better understanding the FM, and simulation errors on model levels of three analyses can be a helpful FM error reference.