anchor.diff:    if (ios /= 0) then
        stop "Error while opening output file."
     end if
  -
  + 
     ! Save initial conditions
  -  write(15, '(7(A16))')'r', 'x', 'T*', 'S*', 'T', 'Sigma', 'H*'
  +  write(15, '(7(A16))')'r', 'x', 'T*', 'S*', 'T', 'Sigma', 'H'
     do i= 1, n_cell
        write(15, '(7(e16.6e2))')r_state%r(i), x_state%x(i), s%T(i), s%S(i),&
  -                                            s%T(i)*state_0%T_0, s%S(i) * state_0%s_0 / x_state%x(i), CI%H_over_r(i)
  +          s%T(i)*state_0%T_0, s%S(i) * state_0%s_0 / x_state%x(i), CI%H_over_r(i)
     enddo
     close(15)
   
  
positive.code:   !----------------------------------------------
    ! Set the initial conditions for S, T
    !----------------------------------------------
    s%T = 8._x_precision*CI%T_ci / state_0%T_0
    s%S = CI%Sig_ci / state_0%S_0 * x_state%x
  
    ! H_over_r becoming H*
  
negative.code:   implicit none
  
    !--------------------------- Parameters for s_curve-----------------------
    real(x_precision), parameter         :: eps_in = 1.e-6_x_precision ! Precision required for the dichotomy
    real(x_precision), parameter         :: Tmin   = 2.5e-2_x_precision / 2.
    real(x_precision), parameter         :: Tmax   = 4.49e0_x_precision * 2.
    real(x_precision), parameter         :: Smin   = 2.36e1_x_precision / 5.
  