anchor.diff:    !-------------------------------------------------------------------------
     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(*,"(' Simga_0     =',1p,E12.4)") Sigma_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
  
positive.code:     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 )
  
    end subroutine initial_variables
  
  
negative.code:              call variables(T, Smin, Omega, H, rho, cs, nu, Q_plus, Q_minus, K_ff, K_e,&
                    tau_eff, P_rad, P_gaz,E_ff,Fz,f_min,Sigma_0, Omega_0,rs)
            endif
      
            dichotomy       = (Smin+Smax)/2.0d0
            call variables(T, dichotomy, Omega, H, rho, cs, nu, Q_plus, Q_minus, K_ff,&
                 K_e, tau_eff, P_rad, P_gaz,E_ff,Fz,f_dichotomy,Sigma_0, Omega_0,rs)
            j               = j + 1
      
         end do iteration
      endif
    
    end function dichotomy
    
    
  