real*8 function MAGNE(S, L)
  implicit none
  integer*4, intent(in) :: L
  integer*2, intent(out) :: S(1:L,1:L)
  integer*4 :: i, j
  real*8 :: mag
  mag = 0d0

  do i = 1, L
    do j = 1, L
      mag = mag + S(i, j)
    enddo
  enddo

  magne = mag
  return
end function MAGNE