In order to understand the influence of magnetic fields on the propagation properties of waves, as derived from different local helioseismology techniques, forward modeling of waves is required. Such calculations need a model in magnetohydrostatic equilibrium as an initial atmosphere through which to propagate oscillations. We provide a set of the IDL routines to create a magneto-static sunspot model in equilibrium from the sub-photospheric to chromospheric layers for a wide range of parameters. The method combines the advantages of self-similar solutions and current-distributed models. A set of models already calculated by this method are also provided. The steps are described in The parameters of the magneto-static equilibrium model can be adjusted directly in this routine. More information in Khomenko, E., & Collados, M. 2008, ApJ, 689, 1379