anchor.diff:      real(kind = x_precision)                               :: Omega_0
       real(kind = x_precision)                               :: rho_0
       real(kind = x_precision)                               :: T_0
  -    real(kind = x_precision)                               :: temp_reel=0.0d0
  -    real(kind = x_precision)                               :: sigma_reel=0.0d0
  +    real(kind = x_precision)                               :: temp_real=0.0d0
  +    real(kind = x_precision)                               :: sigma_real=0.0d0
       character(len = 8)                                     :: number_of_cell
       character(len = 64)                                    :: fname
       integer                                                :: fid
  
positive.code:     do j              = 1 , 1
         !do j              = 1 , n_cell
         !  r = (param%rmax-rmin)/(n_cell-1)*(j-1) + rmin
         r              = 10._x_precision*G*param%M/(c**2)     
         omega          = sqrt(G*param%M/r**3) / Omega_0 
         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        = 1d-3
            Smax        = 1d51
            
            temp        = (t_max-t_min)/(nb_it-1)*(i-1) + t_min
            sigma       = dichotomy(Smin, Smax, eps, temp, omega, sigma_0, Omega_0,rs, T_0, rho_0)
            
            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, Omega_0,rs, T_0, rho_0)
  
            
            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)
            
            temp_reel   = temp * T_0
            sigma_reel  = sigma * sigma_0
            write(fid,'(1p,E12.6,4x,1p,E12.6,4x,1p,E12.6)')sigma_reel,temp_reel, tau_eff
         enddo
         close(fid)
      enddo
      
    end subroutine curve
  
  
  
negative.code:   subroutine variables(T, 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, Omega_0,rs, T_0, rho_0)
      implicit none
  
      type(parameters)                                         :: param
      real(kind = x_precision),intent(in)                      :: T,Sigma,Omega
      real(kind = x_precision)                                 :: coeff_a=0.,coeff_b=0.,coeff_c=0.
      real(kind = x_precision),intent(out)                     :: H
      real(kind = x_precision),intent(out)                     :: rho
      real(kind = x_precision),intent(out)                     :: cs
      real(kind = x_precision),intent(out)                     :: nu
      real(kind = x_precision),intent(out)                     :: Q_plus
      real(kind = x_precision),intent(out)                     :: K_ff
      real(kind = x_precision),intent(out)                     :: K_e
      real(kind = x_precision),intent(out)                     :: E_ff
      real(kind = x_precision),intent(out)                     :: tau_eff
      real(kind = x_precision),intent(out)                     :: Fz
      real(kind = x_precision),intent(out)                     :: Q_minus
      real(kind = x_precision),intent(out)                     :: P_rad
      real(kind = x_precision),intent(out)                     :: P_gaz
      real(kind = x_precision),intent(in)                      :: Sigma_0
      real(kind = x_precision),intent(in)                      :: Omega_0
      real(kind = x_precision),intent(in)                      :: T_0
  