The FORTRAN code used to compute visco-elastic-plastic stresses in Io's crust resulting from temperature changes due to changing advective heat flow using finite difference methods. Details and supporting information are given in Kirchoff, M.R., 2006, Mountain building on Io: An unsteady relationship between volcanism and tectonism (Earth and Planetary Sci.). Washington University, Saint Louis, Missouri (Ph.D. Thesis). See also Kirchoff, M.R. & McKinnon, W.B., 2009, Formation of mountains on Io: Variable volcanism and thermal stresses. Icarus 201, 598–614.