diff --git a/source/KKRnano/source/bfield/torque.f90 b/source/KKRnano/source/bfield/torque.f90 index 485515e17495bb7bc56351ce7fa9fe563d365828..b80eb32f70ac47838544d2dd0d7d438c5ecb10ab 100644 --- a/source/KKRnano/source/bfield/torque.f90 +++ b/source/KKRnano/source/bfield/torque.f90 @@ -108,7 +108,7 @@ contains call intcheb_cell(integrand,integrate_result,rpan_intervall,ipan_intervall,npan_tot,ncheb,irmd) mag(ilm,i) = real(integrate_result * rfpi) ! Same for only the muffin tin: Set to zero outside, integrate again - integrand(imt:) = (0., 0.) + integrand(imt+1:) = (0., 0.) call intcheb_cell(integrand(:),integrate_result,rpan_intervall,ipan_intervall,npan_tot,ncheb,irmd) mag_mt(ilm,i) = real(integrate_result * rfpi) end do @@ -135,7 +135,7 @@ contains call intcheb_cell(integrand(:),integrate_result,rpan_intervall,ipan_intervall,npan_tot,ncheb,irmd) torque(i) = torque(i) + real(integrate_result * rfpi) ! rfpi only contained in convoluted quantities ! Same for mt - integrand(imt:) = 0. + integrand(imt+1:) = 0. call intcheb_cell(integrand(:),integrate_result,rpan_intervall,ipan_intervall,npan_tot,ncheb,irmd) torque_mt(i) = torque_mt(i) + real(integrate_result * rfpi) end do