anchor.diff:                   s%Mdot(n_cell) = s%Mdot(n_cell) * 2._x_precision
                    print*, 'Converged in mean!', maxval(abs(mean)), S_rel_diff_max
                    print*, 'Mdot kick! YOLOOOOO', s%Mdot(n_cell)
  +                 mean = 0._x_precision
                 else
                    if (mod(iteration, output_freq) < 10) then
                       print*, maxval(abs(mean)), S_rel_diff_max
  
positive.code:               s%Mdot(n_cell) = s%Mdot(n_cell) * 2._x_precision
                print*, 'Converged totally!'
                print*, 'Mdot kick! YOLOOOOO', s%Mdot(n_cell), S_steps
             else
                mean = S_rel_diff / mean_size + &
                     (mean_size - 1._x_precision) / mean_size * mean
  
negative.code:   real(x_precision), dimension(n_cell) :: dist_crit
  
    call getarg(1, arg)
    mean        = 0._x_precision
    !----------------------------------------------
    ! Convergence criteria for S and T
    !----------------------------------------------
  