subroutine rezero_spin_diff() character(*), parameter :: this_routine = "rezero_spin_diff" inst_spin_diff = 0.0_dp end subroutine rezero_spin_diff