The thermodynamic processes undergone by olivine, the chief mineral of subducting slabs, are important in studying the mechanisms of deep-focus earthquakes and the interaction between slabs and the surrounding mantle.We propose a modified coupling scheme to calculate the thermal and phase structures of slabs, which presents a high computational accuracy at an accessible CPU expense. Using the new code of computation, we have calculated the thermal and phase structures of the Kurile subduction zone.