anchor.diff:    !----------------------------------------------
     open(13, file="output.dat", status="replace", iostat=ios)
     if (ios /= 0) then
  -    stop "Error while opening output file."
  +     stop "Error while opening output file."
     end if
   
     !----------------------------------------------
  
positive.code:   ! 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)
  
  
negative.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
  
  
    ! FIXME
    S_crit = 1.e99_x_precision
  