anchor.diff:      real(kind = x_precision)                               :: sigma = 0.0d0
       real(kind = x_precision)                               :: Smin  = 0.0d0
       real(kind = x_precision)                               :: Smax  = 0.0d0
  +    real(kind = x_precision)                               :: omega = 0.0d0
       real(kind = x_precision)                               :: r     = 0.0d0
  -    real(kind = x_precision)                               :: omega = 0.0d0
       real(kind = x_precision)                               :: H     = 0.0d0
       real(kind = x_precision)                               :: rho   = 0.0d0
       real(kind = x_precision)                               :: cs    = 0.0d0
  
positive.code:   end function dichotomy
  
  
    !Subroutine in order to display variables
    !-------------------------------------------------------------------------
    subroutine display_variables(temp,Omega,r,sigma, H, rho, cs, nu, Q_plus, &
  
negative.code:     character(len = 64)                                    :: fname
      character(len = 64)                                    :: fname_2
      character(len = 64)                                    :: fname_3
      integer                                                :: fid
      integer                                                :: fid_2
      integer                                                :: fid_3
      !integer                                               :: n_cell = 1
      integer                                                :: index_fcp
      real(kind = x_precision)                               :: sigma_c_thick
      real(kind = x_precision)                               :: temp_c_thick
  