anchor.diff:  36 format(F6.4, 4X, F8.4, 4X, F6.4, 4X, F6.4, 4X)
   38 format(F6.4, 4X, F6.4, 4X, I8)
   
  -open(unit=10, file='L_40-T_1_01.dat', status='unknown')
  +open(unit=10, file='L_120-T_1_01.dat', status='unknown')
   write(10, *) "#m    e    m2    e2    "
   
   ! Inicializo la red de spines.
  
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
  