anchor.diff:            x     = x + (half*xp)*Red_S
             y     = y + (half*yp)*Red_S
   
  -          call calc_ion_loss_cry(is,pc,srefl,dest)
  +          call cry_calcIonLoss(is,pc,srefl,dest)
             pc = pc - dest*Srefl ! "added" energy loss before capture
   
  -          call calc_ion_loss_cry(is,pc,sdech,dest)
  +          call cry_calcIonLoss(is,pc,sdech,dest)
             pc = pc - (half*dest)*Sdech ! "added" energy loss while captured
   
  -          call calc_ion_loss_cry(is,pc,red_s,dest)
  -          call move_am(is,nam,red_s,dest,dlyi(is),dlri(is),xp,yp,pc)
  +          call cry_calcIonLoss(is,pc,red_s,dest)
  +          call cry_moveAM(is,nam,red_s,dest,dlyi(is),dlri(is),xp,yp,pc)
             x = x + (half*xp)*Red_S
             y = y + (half*yp)*Red_S
   
  
positive.code:         x     = x + (half*s_length)*xp
          y     = y + (half*s_length)*yp
          if(zn > zero) then
            call calc_ion_loss_cry(is,pc,s_length,dest)
            call move_am(is,nam,s_length,dest,dlyi(is),dlri(is),xp,yp,pc)
          end if
          x = x + (half*s_length)*xp
          y = y + (half*s_length)*yp
  
negative.code:       x  = x  + (half*s_length)*xp
        y  = y  + (half*s_length)*yp
  
        call cry_calcIonLoss(is,pc,s_length,dest)
        call cry_moveAM(is,nam,s_length,dest,dlyi(is),dlri(is),xp,yp,pc)
  
        x = x + (half*s_length)*xp
        y = y + (half*s_length)*yp
  