anchor.diff:      real(kind=x_precision),intent(out)                      :: Omega_0
       real(kind=x_precision),intent(out)                      :: T_0
       real(kind=x_precision),intent(out)                      :: rho_0
  -    real(kind=x_precision)                                  :: nu_0
  -
  -    type(parameters)                                        :: param
  +    !real(kind=x_precision)                                  :: nu_0
  +
       !------------------------------------------------------------------------
   
  -    call get_parameters(param)
  -    rs                 = 2._x_precision * G * param%M/(c**2)
  +    rs                 = 2._x_precision * G * params%M/(c**2)
       rmin               = 3._x_precision * rs
  -    Mdot_0             = param%Mdot
  -    Omega_0            = sqrt(G * param%M / rs**3 )
  +    Mdot_0             = params%Mdot
  +    Omega_0            = sqrt(G * params%M / rs**3 )
       Sigma_0            = Mdot_0 /(Omega_0 * rs**2 * 2_x_precision * pi)
       !nu_0                 = 2._x_precision * rs**2 * Omega_0/3._x_precision
   
  
positive.code:   end subroutine display_variables
  
  
  end module mod_S_curve
  
negative.code:     character(len = 8)                                     :: number_of_cell
      character(len = 64)                                    :: fname
      integer                                                :: fid
      !------------------------------------------------------------------------
      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)
      !-------------------------------------------------------------------------
  