anchor.diff:      rmin               = 3._x_precision * rs
       Mdot_0             = param%Mdot 
       Sigma_0            = Mdot_0 /(state_0%Omega_0 * rs**2 * 2 * pi)
  -    Omega_0            = sqrt( G*param%M/rs**3 )
  +    Omega_0            = state_0%Omega_0
   
     end subroutine initial_variables
   
  
positive.code:        !do j              = 1 , n_cell
         !  r = (param%rmax-rmin)/(n_cell-1)*(j-1) + rmin
         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-8
            Smax        = 1.0d9
            temp        = (t_max-t_min)/(nb_it-1)*(i-1) + t_min
            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,&
  
negative.code:   !-------------------------------------------------------------------------
    subroutine display_initial_variables(rs, rmin, Mdot_0, Sigma_0, Omega_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
      real(kind = x_precision), intent(in)                     :: Omega_0
  
      Type(parameters)                                         :: param
      call get_parameters(param)
      !-----------------------------------------------------------------------
  
      write(*,*)'           Initial Variables            '
      write(*,*)'****************************************'
      write(*,"(' Temp_0      =',1p,E12.4)") state_0%T_0
      write(*,"(' Sigma_0     =',1p,E12.4)") Sigma_0
      write(*,"(' Omega_0     =',1p,E12.4)") Omega_0
      write(*,"(' Omega_max   =',1p,E12.4)") sqrt(G*param%M / rmin**3)
      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
  