An ionospheric error simulation is needed for creating precise Global Positioning System (GPS) signal using GPS simulator. In this paper we developed Klobuchar coefficients n ${alpha}_n$ and ${eta}_n$ (n = 1, 2, 3, 4) generation algorithms for simulator and verified accuracy of the algorithm. The algorithm extract those Klobuchar coefficients from broadcast (BRDC) messages provided by International GNSS Service during three years from 2006 through 2008 and curve-fit them with sinusoidal and linear functions or constant. The generated coefficients from our developed algorithms are referred to as MODL coefficients, while those coefficients from BRDC messages are named as BRDC coefficients. The maximum correlation coefficient between MODL and BRDC coefficients was found for ${alpha}_2$ and the value was 0.94. On the other hand, the minimum correlation was 0.64 for the case of ${alpha}_1$. We estimated vertical total electron content using the Klobuchar model with MODL coefficients, and compared the result with those from the BRDC model and global ionosphere maps. As a result, the maximum RMS was 3.92 and 7.90 TECU, respectively.