Ocean surface waves have been utilized as fundamental information in various fields of oceanic research. In this paper, we suggest a simulation and modelling technique for generating an ocean surface wave using an inverse Fast Fourier Transform (iFFT), and we subsequently verify the wave’s accuracy. The conventional method, linear superposition, requires recursive calculation because of the double summation and the time variable; to circumvent this issue, the new algorithm is presented. The Joint North Sea Wave Project (JONSWAP) spectrum is utilized for the ocean surface wave simulation example, and the parameters are the significant wave height (HS) and the zero-crossing wave period (TZ). A coordinate transform for the wavenumber domain was used to apply the inverse FFT algorithm. To verify the accuracy of the simulation result, the relative error between the input condition and the analysis result was calculated. The result for TZ is below 4% relative error, and the maximum relative error for HS is 7%. To avoid the Nyquist frequency for wave-field analysis and simulation, the minimum grid size was calculated by twice applying the maximum wavenumber.