anchor.diff:         &               n_inter_max_warn, &
          &               n_reabs_max,  &
          &               n_reabs_max_warn,  &
  +       &               forced_first_interaction, &
  +       &               forced_first_interaction_algorithm, &
          &               kill_on_scatter, &
          &               monochromatic_energy_threshold
   
  
positive.code:                 end if
  
                  ! Propagate until photon is absorbed again
                  call propagate(p, peeloff_scattering_only)
  
               end do
  
  
negative.code:                       end if
  
                        ! Propagate until photon is absorbed again
                        call propagate(p)
  
                     end if
  
  