A numerical software package for (2 plus 1)-dimensional simulation of optical limiters has been developed. The purpose of this effort is to provide the sensor protection community with a means to develop, understand, and optimize optical limiters. A graphics interface has been implemented to guide the user in both running the program and interpreting the results. Scientific visualization tools have been integrated with the package to facilitate a physical understanding of the numerical results. We have used it successfully to perform (2 plus 1)-dimensional nonlinear propagation of a focused optical beam with arbitrary intensity profile, such as Gaussian or tophat profile, in a f/5 system. The package can handle various nonlinear mechanisms, including intensity dependent refraction and absorption. Non-local nonlinearities such as thermal effects may be included as a modular subprogram. Temporal effects may also be included by simulating the integrated spatio-temporal evolution of the light and material properties, although this approach is memory- intensive and often requires the use of a super-computer.