anchor.diff:  !  note: otherwise ipart1_in,ipart2_in have to set to one if not used       !
   !                                                                           !
   !----------------------------------------------------------------------------
  +  call getarg(7, arg)
  +  read (arg, *, iostat=stat) ipart1_in
  +  if (.not. stat == 0) then
  +    write (0,*) 'ipart1_in must be integer'
  +    stop 1
  +  endif
  +  call getarg(8, arg)
  +  read (arg, *, iostat=stat) ipart2_in
  +  if (.not. stat == 0) then
  +    write (0,*) 'ipart2_in must be integer'
  +    stop 1
  +  endif
   
   !----------------------------------------------------------------------------
  -  isquark1_in = 0                                                           !
  -  isquark2_in = 0                                                           !
  +!  isquark1_in = 0                                                          !
  +!  isquark2_in = 0                                                          !
   !                                                                           !
   !  for LO with light-squark flavor in the final state                       !
   !  isquark1_in     =  -5,-4,-3,-2,-1,+1,+2,+3,+4,+5                         !
  
positive.code:   real(kind=double)                    :: energy_in
    logical                              :: lfinal
    character(len=2)                     :: final_state_in
  
  !----------------------------------------------------------------------------
    inlo = 1       ! specify LO only[0] or complete NLO (slower)[1]           !
  !                ! results: LO     - LO, degenerate squarks, decoupling on  !
  !                !          NLO    - NLO, degenerate squarks, decoupling on !
  !                !          LO_ms  - LO, free squark masses, decoupling off !
  
negative.code: !                                                                           !
  !  squark and antisquark added, but taking into account different sb or ss  !
  !----------------------------------------------------------------------------
    call getarg(6, final_state_in)
  
  !----------------------------------------------------------------------------
  !  ipart1_in = 5                                                            !
  !  ipart2_in = 7                                                            !
  !                                                                           !
  !  final_state_in = ng,ns,nn                                                !
  !  ipart1_in   = 1,2,3,4  neutralinos                                       !
  