anchor.diff:    real(kind=x_precision),intent(in)                        :: eps         
     real(kind=x_precision),intent(in)                        :: T           
     real(kind=x_precision),intent(in)                        :: omega
  +  real(kind=x_precision)                                   :: H         
  +  real(kind=x_precision)                                   :: rho           
  +  real(kind=x_precision)                                   :: cs
  +  real(kind=x_precision)                                   :: nu        
  +  real(kind=x_precision)                                   :: Q_plus           
  +  real(kind=x_precision)                                   :: Q_minus
  +  real(kind=x_precision)                                   :: K_ff        
  +  real(kind=x_precision)                                   :: K_e           
  +  real(kind=x_precision)                                   :: tau_eff
  +  real(kind=x_precision)                                   :: P_rad        
  +  real(kind=x_precision)                                   :: P_gaz
  +  real(kind=x_precision)                                   :: f_min          
  +  real(kind=x_precision)                                   :: f_max
  +  real(kind=x_precision)                                   :: f_dichotomy
  +  real(kind=x_precision)                                   :: E_ff
  +  real(kind=x_precision)                                   :: Fz
  +  real(kind=x_precision)                                   :: Sigma_0
  +  
     !-------------------------------------------------------------------------
     ! N-> Number of iterations for the dichotomy
     ! Smin, Smax -> Starting range
  
positive.code:   write(*,"(' Temp        =',1p,E12.4)") temp
    write(*,"(' Sigma       =',1p,E12.4)") sigma
    write(*,"(' Omega       =',1p,E12.4)") Omega
    !write(*,*)'H         = ', sol
    write(*,*)'****************************************'
    read(*,*)
  
  
negative.code: 
  end subroutine display_parameters
  
  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
    write(*,"(' rmin        =',1p,E12.4)") rmin
    write(*,"(' rs          =',1p,E12.4)") rs
    write(*,*)'****************************************'
    read(*,*)
  end subroutine display_initial_variables
  
  
  
  subroutine 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)
    implicit none
    real(kind = x_precision), intent(in)                     :: temp
    real(kind = x_precision), intent(in)                     :: Omega
    real(kind = x_precision), intent(in)                     :: r
    real(kind = x_precision), intent(in)                     :: sigma
    real(kind = x_precision), intent(in)                     :: H     
    real(kind = x_precision), intent(in)                     :: rho        
    real(kind = x_precision), intent(in)                     :: cs  
    real(kind = x_precision), intent(in)                     :: nu    
    real(kind = x_precision), intent(in)                     :: Q_plus          
    real(kind = x_precision), intent(in)                     :: Q_minus  
    real(kind = x_precision), intent(in)                     :: K_ff         
    real(kind = x_precision), intent(in)                     :: K_e            
    real(kind = x_precision), intent(in)                     :: tau_eff  
    real(kind = x_precision), intent(in)                     :: P_rad        
    real(kind = x_precision), intent(in)                     :: P_gaz    
    real(kind = x_precision), intent(in)                     :: f        
    real(kind = x_precision), intent(in)                     :: E_ff     
    real(kind = x_precision), intent(in)                     :: Fz  
  
    write(*,*)'               Variables                '
    write(*,*)'****************************************'
  