In this study, a three-dimensional numerical model (PersianGulfOceanicModel (ZSF974)) based on the primitive equations, in the Earth’s Spherical Coordinates System with a sigma vertical coordinate, has been developed with the aim of predicting and calculating oceanographic parameters in the marine medium of the Persian Gulf. The finite difference method has been used for the numerical solution of the model equations. For discretization, the Lax–Wendroff scheme was applied to the advection terms, the DuFort–Frankel scheme to the diffusion terms, and the Matsuno scheme was used to eliminate instabilities arising from certain calculations within the program. The mesh used is the modified Arakawa C grid. This model, in addition to accommodating any type of non-level bottom, has the capability to vary resolution in both horizontal and vertical directions. The superiority of the process implemented in its design, the proper application of Nihoul's theory (1977) regarding the effect of surface stress induced by wind in the lower layers has enhanced the accuracy of the model's outputs. In this study is presented the process of designing the base model and the dimensions of its validation through models in laboratorial mediums. In all laboratorial oceanic mediums, the stationary state, the effect of neglecting the Coriolis force and earth's curvature, as well as the effect of applying a constant force have been investigated. Hypothetical (and approximate real) data and forced tide, have been used as the program's input. The results of the model’s execution are consistent with ocean physics principles and the findings of previous researchers and can represent the overall behavior of an oceanic medium. This model can be used as a base model to examine the overall behavior of an oceanic medium. It also has the capability to be developed for implementation and conclusion in a real-world medium.