anchor.diff:             ! Output things here
           end if
        else
  -        print*, 'Mdot kick'
  +       ! print*, 'Mdot kick'
           ! give a Mdot kick
        end if
        
  
positive.code:   ! Copy the value of state_0 into state vector s
       s%T = CI%T_ci / state_0%T_0
       s%S = CI%Sig_ci / state_0%S_0
  
    call compute_variables(s)
  
  
negative.code: 
    implicit none
  
    integer                                     :: iteration, ios, j, i
    type(state)                                 :: s
    real(kind = x_precision)                    :: delta_S_max, delta_T_max, t
    real(kind = x_precision), dimension(n_cell) :: prev_S, prev_T, S_crit
  