anchor.diff:              end if
               
            end if
  -      !!!if (.not. LESF ) then  !RRD: trying to emulate the old AD14 SEPAR.f90 with SHIFT=NOT(LESF), go back to original commented!!! out above when done
  -      !!!        OtherState%sigma1(i,j) = 0.667_ReKi
  -      !!!end if
         
   !bjj: update sigma3 to value at t + dt      
            
  
positive.code:          xd%tau_V(i,j)          = xd%tau_V(i,j) + 2.0_ReKi*p%dt*u%U / p%c(i,j)  
        end if
     
        
  #ifdef UA_OUTS
     m%TESF(i,j) = TESF  
  
negative.code:       !! NEW CODE 2/19/2015
        VRTX = (xd%tau_V(i,j) <= 2.0_ReKi*T_VL) .and. (xd%tau_V(i,j) > 0.0_ReKi)
        
  
        
  !bjj: update sigma1 to value at t + dt      
        OtherState%sigma1(i,j) = 1.0_ReKi
  