anchor.diff:      character(len = 8)                                     :: number_of_cell
       character(len = 64)                                    :: fname
       integer                                                :: fid
  -    type(parameters)                                       :: param
       !------------------------------------------------------------------------
  -    call get_parameters(param)
  -    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)
       !-------------------------------------------------------------------------
  
positive.code:     real(kind=x_precision)                                 :: delta = 0.0d0
      !------------------------------------------------------------------------
      delta          = coeff_b**2 - 4_x_precision * coeff_a * coeff_c
      if (delta .lt. 0.) then
        write(*,*)'No solutions in the R field.'
     else
        sol_1         = -0.5_x_precision *(coeff_b + sign(sqrt(delta),coeff_b))/coeff_a
        if (sol_1 .lt. 0.d-12) then
           write(*,*)'There is a problem with the solution of the quadratic equation.'
           stop
        endif
  
negative.code:     write(*,*)'fmin = ',f_min
      write(*,*)'fmax = ',f_max
  
      if ( f_max * f_min > 0.) then
         write(*,*)'This function image does not switch its sign in this particular interval.'
      endif
  
      if( f_max * f_min > 0.) then
         iteration:do while ( dabs( Smax - Smin ) >= eps .and. j < 10000)
  
  
      call variables(T, Smin, Omega, H, rho, cs, nu, Q_plus, Q_minus, K_ff, K_e,&
  