anchor.diff:  
   ! dar vuelta la red por completo una vez
   ! de forma ordenada
  -do tmc =1, 100000
  +do tmc =1, 1000000
       do k = 0, (L * L) -1
           i = mod(k, L)
           j = int(k/L)
  
positive.code: 
  implicit none
  integer                     :: i, j, k
  integer, parameter          :: L=40, u_max=100
  integer                     :: s(0:(L*L)-1), tmc
  integer                     :: n_down, n_up, n_left, n_right, spinadd, loc
  ! real(dp), parameter       :: k_b=1.3806488e-23
  
negative.code:         end if
      end do
  
      if (tmc > 10000) then  ! descarto transitorio
          m = sum(s)/real(L*L, pr)
          e_w = e/real(L*L, pr)
          ! actualizo e_c, m_c y calculo correlaciones corr_m,e
  