anchor.diff:  
       Lrefl = xp_rel*r ! Distance of refl. point [m]
       Srefl = sin_mb(xp_rel/two + c_miscut)*Lrefl
  -
  -    if(Lrefl > zero .and. Lrefl < Length) then ! VR point inside
  +    write(*,*) 'Lrefl', Lrefl
  +    write(*,*) 'L_chan', L_chan
  +
  +    if(Lrefl > zero .and. Lrefl < L_chan) then ! VR point inside
   
         ! 2 options: volume capture and volume reflection
   
  
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)
  