32f846e8b43c370e250a72c295afd1e7abf00837
diff --git a/source/physics/materialProperties/MagneticResistivity/MagneticResistivityMain/SpitzerHighZ/MagneticResistivity_fullState.F90 b/source/physics/materialProperties/MagneticResistivity/MagneticResistivityMain/SpitzerHighZ/MagneticResistivity_fullState.F90
index 7fe26d54c3..f85be6e37b 100644
--- a/source/physics/materialProperties/MagneticResistivity/MagneticResistivityMain/SpitzerHighZ/MagneticResistivity_fullState.F90
+++ b/source/physics/materialProperties/MagneticResistivity/MagneticResistivityMain/SpitzerHighZ/MagneticResistivity_fullState.F90
@@ -127,30 +127,4 @@ subroutine MagneticResistivity_fullState(solnVec,resPar, resPerp, resCross)
      end if
   end if
 
-contains
-
-  subroutine loglambda(tele, nele, zbar, ll)
-    implicit none
-
-    ! This subroutine computes the Coulomb logarithm. The formula used
-    ! comes from Atzeni.
-    
-    real, intent(in)  :: tele ! electron temperature [K]
-    real, intent(in)  :: nele ! electron number density [cm^-3]
-    real, intent(in)  :: zbar ! the average ionization [unitless]
-    real, intent(out) :: ll   ! the coulomb logarithm [unitless]
-
-    real :: bmax, bmin, bmin_classic, bmin_quantum
-    real, parameter :: ll_floor = 1.0
-    
-    bmax = sqrt(res_boltz * tele / (4*PI * res_qele**2 * nele))
-    
-    bmin_classic = zbar * res_qele**2 / (3*res_boltz*tele)
-    bmin_quantum = res_hbar / (2*sqrt(3*res_boltz*tele*res_mele))
-    bmin = max(bmin_classic, bmin_quantum)
-
-    ll = max(log(bmax/bmin), ll_floor)
-
-  end subroutine loglambda
-
 end subroutine MagneticResistivity_fullState