anchor.diff:  
           T_converged = .false.
   
  -        t_steps = 0
  +        call compute_variables(s)
  +        
  +        T_steps = 0
           do while (.not. T_converged)
  -           t_steps = t_steps + 1
  -           if (t_steps >= 20 .and. mod(t_steps, 20) == 0) then
  -              print*, 'W: T is taking time to converge, already', t_steps, 'iterations'
  +           T_steps = T_steps + 1
  +           if (T_steps >= 20 .and. mod(T_steps, 20) == 0) then
  +              print*, 'W: T is taking time to converge, already', T_steps, 'iterations'
              end if
   
              call compute_variables(s)
  
positive.code:         iteration = iteration + 1
  
          if (mod(iteration, output_freq) == 0 ) then
             if (arg == 'load' .and. iteration > stp_value) then
                stop
             endif
             call snapshot(s, iteration, t, 13)
  
negative.code: 
             ! Do a snapshot
             if (mod(iteration, output_freq) == 0) then
                if ((arg == 'load' .or. arg == 'restart') &
                     .and. iteration-start_iteration > stp_value) then
                   stop
                endif
                call snapshot(s, iteration, t, 13)
  