anchor.diff:    ! Pass to T*_crit and S*_crit
     temperature_c  = temperature_c / state_0%T_0
     sigma_c = sigma_c / state_0%S_0 * x_state%x
  -  
  +
     ! Initiate the S_curve
     ! call s_curve(foo, bar)
   
  
positive.code:   real(kind = x_precision), dimension(n_cell) :: dt_nu, dt_T
    real(kind = x_precision)                    :: min_dt_nu, min_dt_T
    logical                                     :: T_converged
    real (kind = x_precision), dimension(n_cell):: dist, dist_crit
    
  
    ! FIXME
    S_crit = 1.e99_x_precision
  
negative.code:   !----------------------------------------------
    open(13, file="output.dat", status="replace", iostat=ios)
    if (ios /= 0) then
       stop "Error while opening output file."
    end if
  
    !----------------------------------------------
  