Ionospheric error modeling is necessary to create reliable global navigation satellite system (GNSS) signals using a GNSS simulator. In this paper we developed algorithms to generate Klobuchar coefficients ${alpha}_n$, ${eta}_n$ (n = 1, 2, 3, 4) for a GNSS simulator and verified accuracy of the algorithm. The eight Klobuchar coefficients were extracted from three years of global positioning system broadcast (BRDC) messages provided by International GNSS service from 2006 through 2008 and were fitted with Fourier series. The generated coefficients from our developed algorithms are referred to as Fourier Klobuchar model (FOKM) coefficients, while those coefficients from BRDC massages are named as BRDC coefficients. The correlation coefficient values between FOKM and BRDC were higher than 0.97. We estimated total electron content using the Klobuchar model with FOKM coefficients and compared the result with that from the BRDC model. As a result, the maximum root mean square was 1.6 total electron content unit.