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:       if (i .eq. 1) then
           state_out%v(i)  = 0._x_precision
        else if (i .eq. n_cell) then
           state_out%v(i)  = - 1._x_precision / state_out%S(i) / state_out%x(i)
        else
           state_out%v(i)  = - 1._x_precision / state_out%S(i) / state_out%x(i) * ( ((state_out%nu(i+1) * &
                             state_out%S(i+1)) - (state_out%nu(i-1) * state_out%S(i-1))) /  &
  