anchor.diff:      end if
       etaplus = 0.0d0
       etaminus= 0.0d0
  -   !   xnplus  = xbplus
  -    !  xnminus = xbminus
  -  !    dxnplus = dxbplus
  -!	dxnminus= dxbminus
  +		xnplus  = xbplus   ! 0.0d0  ! xbplus
  +    xnminus = xbminus  ! xbminus
  +    dxnplus = dxbplus  ! dxbplus
  +    dxnminus= dxbminus ! dxbminus
  +
     	do inbeam=1,nbeam
     		if (kluch_beam.eq.2 .and. k.eq.1) then 	   ! режим заданного тока
       		call beam_calc (tau)
  
positive.code: 
    integer inr,is,k,kluch1
  	real power,powerenter,powerexit,powersum,powerim,powerplus,powerminus, dva_d_na_lambda
  !	double complex xrab1(20)
    if (kluch1.eq.1) then
      do is=1,sk
  
negative.code: 
  
  subroutine manager (stepwrite,iw)
  	use array_work
  	use beam_nes
  	use field_nes
    implicit none
  
    integer k,i11,is,in,stepwrite,iw,i,im,istwr
    real*8 tau,rel_coord
    complex fffw(256),fffb(256)
  
    dxnplus = 0.0d0
    dxnminus= 0.0d0
  