anchor.diff:  
       call get_parameters(param)
      
  -    coeff_a              = (Omega**2 * Omega_0**2 * Sigma * Sigma_0)/2.
  -    coeff_b              = (-1._x_precision/3._x_precision) * cst_rad*T**4 * T_0**4
  -    coeff_c              = (-1._x_precision * param%RTM * T * T_0* Sigma * Sigma_0)/2._x_precision
  +    coeff_a              = (Omega**2 * Omega_0**2 * Sigma * Sigma_0)/2._x_precision
  +    coeff_b              = (-1._x_precision/3._x_precision) * cst_rad*T**4 * T_0**4 /rs
  +    coeff_c              = (-1._x_precision * param%RTM * T  *  Sigma * Sigma_0)/(2._x_precision * rs**2)
   
       call quadratic(coeff_a , coeff_b , coeff_c , H)
   
  -    H                    = H /rs
  +  !  H                    = H / rs
       rho                  = Sigma / H
       P_rad                = T**4
       P_gaz                = rho * T
       cs                   = Omega * H
       nu                   = param%alpha * cs * H
       K_ff                 = 6.13d22 * rho_0 * rho * (T_0 * T)**(-3.5_x_precision)
  -    K_e                  = 0.2_x_precision * (1_x_precision + param%X)
  +    K_e                  = 0.2_x_precision * (1._x_precision + param%X)
       E_ff                 = 6.22d20 * (rho_0 * rho)**2 * sqrt(T_0 * T)
       tau_eff              = 0.5_x_precision * sqrt(K_e * K_ff) * Sigma * Sigma_0    
       
  
positive.code:     real(kind = x_precision)                               :: Omega_0
      real(kind = x_precision)                               :: rho_0
      real(kind = x_precision)                               :: T_0
      real(kind = x_precision)                               :: temp_real=0.0d0
      real(kind = x_precision)                               :: sigma_real=0.0d0
      character(len = 8)                                     :: number_of_cell
      character(len = 64)                                    :: fname
      integer                                                :: fid
  
negative.code:     !-------------------------------------------------------------------------
      dichotomy             = (Smin+Smax)/2.
      j = 0
  
      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,T_0,rho_0)
      
  