Without self-interaction corrections or the use of hybrid functionals, approximations to the density-functional theory (DFT) often favor intermediate spin systems over high-spin systems. In this paper, we apply the recently proposed Fermi-Löwdin-orbital self-interaction corrected density functional formalism to a simple tetra-coordinated Fe(II)-porphyrin molecule and show that the energetic orderings of the S = 1 and S = 2 spin states are changed qualitatively relative to the results of Generalized Gradient Approximation (developed by Perdew, Burke, and Ernzerhof, PBE-GGA) and Local Density Approximation (developed by Perdew and Wang, PW92-LDA). Because the energetics, associated with changes in total spin, are small, we have also calculated the second-order spin-orbit energies and the zero-point vibrational energies to determine whether such corrections could be important in metal-substituted porphins. Our results find that the size of the spin-orbit and vibrational corrections to the energy orderings are small compared to the changes due to the self-interaction correction. Spin dependencies in the Infrared (IR)/Raman spectra and the zero-field splittings are provided as a possible means for identifying the spin in porphyrins containing Fe(II).