anchor.diff:      x_F = (-B_F-sqrt(B_F**2-four*A_F*C_F))/(two*A_F)
       s_F = alpha_F*x_F + beta_F
       write(*,*) B_F**2, A_F*C_F, four*A_F*C_F, B_F**2-four*A_F*C_F
  -    write(*,*) x_F, s_F
  +    write(*,*) 'new F', x_F, s_F
       
     endif
   
     ! MISCUT fourth step: deflection and length
     a = sqrt(s_F**2+(x-x_F)**2)
  +  tP = acos_mb((2*r**2-a**2)/(2*r**2))
     tdefl = asin_mb((s_f-s_P)/r)
  -  L_chan = r*tdefl
  +  L_chan = r*tP
     write(*,*) "Results", tdefl, L_chan
   
     xp_rel = xp - c_miscut
  
positive.code:   real(kind=fPrec) alpha                ! Par for new chann prob
    real(kind=fPrec) Pvr                  ! Prob for VR->AM transition
  
    real(kind=fPrec) const_dech,xpin,ypin,tchan,tdefl,L_chan,mep
    real(kind=fPrec) s_K,x_K,s_M,x_M,s_F,x_F,r,a
    real(kind=fPrec) A_F,B_F,C_F,alpha_F,beta_F
  
  
negative.code:           x = x + (half*xp)*(s_length - Srefl)
            y = y + (half*yp)*(s_length - Srefl)
          else
            write(*,*) 'TRAM'
            iProc = proc_TRAM
            x = x + xp*Srefl
            y = y + yp*Srefl
            Dxp = ((((-one*(13.6_fPrec/pc))*sqrt(s_length/dlri(is)))*c1m3)*xp_rel)/(two*xpcrit) + &
              (((13.6_fPrec/pc)*sqrt(s_length/DLRi(is)))*c1m3)*(one+(tdefl-c_miscut)/(two*xpcrit))
            xp = xp+Dxp
            x  = x + (half*xp)*(s_length-Srefl)
            y  = y + (half*yp)*(s_length-Srefl)
  