anchor.diff:      real (kind = x_precision),              dimension(n_cell) :: Delta
       real (kind = x_precision),              dimension(n_cell) :: kappa_ff, tau, epsil
       real (kind = x_precision)                                 :: kappa_e !cgs
  -    
  +
       type(state),               intent(inout)                    :: state_out
  -    
  +
       call get_parameters(para)
       kappa_e = 0.2_x_precision * (1._x_precision + para%X)
       !-------------Compute trinome coeff for H----------------
  
positive.code:   private
  
    public :: compute_variables, dim_adim, init_variable_0, state_0
    
  contains
    !--------------------------------------------------------------
    subroutine compute_variables(state_out)
  
negative.code:     type(parameters)                 :: para
  
      c2 = c**2
  
      call get_parameters(para)
      rs = 2._x_precision*G*para%M/c2
  
      state_0%omega_0 = sqrt( G*para%M/(rs)**3 )
      state_0%temps_0 = 2._x_precision / state_0%omega_0
      state_0%x_0     = sqrt(rs)
  