anchor.diff:      !------------------------------------------------------------------------
       call get_parameters(param)
       call display_parameters()
  -    call initial_variables(rs, rmin, Mdot_0, Sigma_0)
  -    call display_initial_variables(rs, rmin, Mdot_0, Sigma_0)
  +    call initial_variables(rs, rmin, Mdot_0, Sigma_0, Omega_0)
  +    call display_initial_variables(rs, rmin, Mdot_0, Sigma_0, Omega_0)
       !-------------------------------------------------------------------------
       ! Test for 1 value of r
       !-------------------------------------------------------------------------
       do j              = 1 , 1
          !do j              = 1 , n_cell
          !  r = (param%rmax-rmin)/(n_cell-1)*(j-1) + rmin
  -       r              = 6*G*param%M/(c**2)     
  +       r              = 10*G*param%M/(c**2)     
          omega          = G*param%M/r**3
          write(number_of_cell,'(I5.5)') j
          fid = 20 + j
          fname = 'Temperature_Sigma_'//trim(number_of_cell)//'.dat'
  -       
          open(fid,file   = fname, status='unknown',action='readwrite') 
  -       write(fid,*)"T              Sigma"
  -       do i           = 1 , nb_it
  -          Smin        = 1.0d-5
  -          Smax        = 1.0d-0
  +       !  write(fid,*)"T              Sigma"
  +        do i           = 1 , nb_it
  +          Smin        = 1.0d-8
  +          Smax        = 1.0d9
             temp        = (t_max-t_min)/(nb_it-1)*(i-1) + t_min
  -          sigma       = dichotomy(Smin, Smax, eps, temp, omega)
  +          sigma       = dichotomy(Smin, Smax, eps, temp, omega, sigma_0, Omega_0,rs)
             call variables(temp, sigma, omega, H, rho, cs, nu, Q_plus, Q_minus, K_ff,&
  -               K_e, tau_eff, P_rad, P_gaz,E_ff,Fz,f,Sigma_0)
  +               K_e, tau_eff, P_rad, P_gaz,E_ff,Fz,f,Sigma_0, Omega_0,rs)
             call display_variables(temp,Omega,r, sigma, H, rho, cs, nu, Q_plus, Q_minus,&
                  K_ff, K_e, tau_eff, P_rad, P_gaz,E_ff,Fz,f)
  -          
             write(fid,'(1p,E12.6,4x,1p,E12.6)')temp,sigma
          enddo
          close(fid)
  
positive.code:   !-------------------------------------------------------------------------
    !Subroutine in order to display initial variables
    !-------------------------------------------------------------------------
    subroutine display_initial_variables(rs, rmin, Mdot_0, Sigma_0)
      implicit none
      real(kind = x_precision), intent(in)                     :: rs
      real(kind = x_precision), intent(in)                     :: rmin
      real(kind = x_precision), intent(in)                     :: Mdot_0
      real(kind = x_precision), intent(in)                     :: Sigma_0
      !-----------------------------------------------------------------------
  
      write(*,*)'           Initial Variables            '
      write(*,*)'****************************************'
      write(*,"(' Temp_0      =',1p,E12.4)") state_0%T_0
      write(*,"(' Simga_0     =',1p,E12.4)") Sigma_0
      write(*,"(' Omega_0     =',1p,E12.4)") state_0%Omega_0
      write(*,"(' H_0         =',1p,E12.4)") state_0%H_0
      write(*,"(' Mdot_0      =',1p,E12.4)") Mdot_0
      write(*,"(' rho_0       =',1p,E12.4)") state_0%rho_0
  
negative.code:     real(kind=x_precision)                                   :: f_dichotomy
      real(kind=x_precision)                                   :: E_ff
      real(kind=x_precision)                                   :: Fz
      real(kind=x_precision),intent(in)                        :: Sigma_0
      real(kind=x_precision),intent(in)                        :: Omega_0
      real(kind=x_precision),intent(in)                        :: rs
  
      !-------------------------------------------------------------------------
      ! N-> Number of iterations for the dichotomy
      ! Smin, Smax -> Starting range
  