anchor.diff: -module mod_S_curve
  +module mod_s_curve
     use mod_constants
     use mod_read_parameters
     use mod_variables
  
positive.code:     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)
      !-------------------------------------------------------------------------
  
negative.code:        K_ff, K_e, tau_eff, P_rad, P_gaz,E_ff,Fz,f,Sigma_0, Omega_0,rs, T_0, rho_0)
      implicit none
  
      real(kind = x_precision),intent(in)                      :: T,Sigma,Omega
      real(kind = x_precision)                                 :: coeff_a=0.,coeff_b=0.,coeff_c=0.
      real(kind = x_precision),intent(out)                     :: H
  