anchor.diff:  
          call build_s_curve(sigma_real_1, sigma_real_0, temp_real_1,nb_it, index_scp, sigma_real, temp_real)
   
  -
  +        call display_critical_points(sigma_c_thin, temp_c_thin,sigma_c_thick, temp_c_thick, k)
  +
  +
  +       
             do l          = 1, nb_it
  -             write(fid_3,'(1p,E12.6,4x,1p,E12.6,4x,1p,E12.6)') sigma_real(l), temp_real(l)
  +             write(fid_3,'(1p,E12.6,4x,1p,E12.6,4x,1p,E12.6)')sigma_real(l),temp_real(l)
             enddo
          
          
  
positive.code:     P_gaz                = rho * T
      cs                   = Omega * H
      nu                   = params%alpha * cs * H
      K_ff                 = 6.13d22 * state_0%rho_0 * rho * (state_0%T_0 * T)**(-3.5_x_precision)
      K_e                  = params%kappa_e
      E_ff                 = 6.22d20 * (state_0%rho_0 * rho)**2 * sqrt(state_0%T_0 * T)
      tau_eff              = 0.5_x_precision * sqrt(K_e * K_ff) * Sigma * state_0%S_0
  
  
      !-------------------------------------------------------------------------
  
negative.code:     
      !------------------------------------------------------------------------
  
      call display_parameters()
      call initial_variables(rs, rmin, Mdot_0, Sigma_0, Omega_0, T_0, rho_0)
      call display_initial_variables(rs, rmin, Mdot_0, Sigma_0, Omega_0, T_0, rho_0)
      !-------------------------------------------------------------------------
      ! Test for 1 value of r
      !-------------------------------------------------------------------------
      do k              = 1 , n_cell
         !r = (rmax-rmin)/(n_cell-1)*(k-1) + rmin
         !r              = 10._x_precision*G*params%M/(c**2)
         r   = r_state%r(k)
         omega          = sqrt(G*params%M/r**3) / Omega_0
         write(number_of_cell,'(I5.5)') k
         fid = 20 + k
         fid_2 = 21 + k
         fid_3 = 22 + k
  
         fname = 's_curves/Temperature_Sigma_'//trim(number_of_cell)//'_1.dat'
         fname_2 = 's_curves/Temperature_Sigma_'//trim(number_of_cell)//'_0.dat'
         fname_3 = 's_curves/Temperature_Sigma_'//trim(number_of_cell)//'_tot.dat'
  
         open(fid,file  = fname, status='unknown',action='readwrite')
         open(fid_2,file  = fname_2, status='unknown',action='readwrite')
         open(fid_3,file  = fname_3, status='unknown',action='readwrite')
  
         do i          = 1, nb_it
            Smin        = 1d-2
            Smax        = 1d54
  
            temp        = (t_max-t_min)/(nb_it-1)*(i-1) + t_min
            optical_depth = 1
            sigma       = dichotomy(Smin, Smax, eps, temp, omega, sigma_0, Omega_0,rs, T_0, rho_0, optical_depth)
            
            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, optical_depth)
  
          !  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_real_1(i)   = log10(temp * T_0)
            sigma_real_1(i)  = log10(sigma * sigma_0)
  
            !do j          = 1, nb_it
            !   write(fid_2,'(1p,E12.6,4x,1p,E12.6,4x,1p,E12.6)')sigma_real_1(j),temp_real_1(j)
            !enddo
  
            optical_depth = 0
  
            sigma       = dichotomy(Smin, Smax, eps, temp, omega, sigma_0, Omega_0,rs, T_0, rho_0, optical_depth)
            
            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, optical_depth)
            
          !  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_real_0(i)   = log10(temp * T_0)
            sigma_real_0(i)  = log10(sigma * sigma_0)
  
            do l          = 1, nb_it
               write(fid,'(1p,E12.6,4x,1p,E12.6,4x,1p,E12.6)')sigma_real_0(l),temp_real_0(l)
            enddo
            
         enddo
  