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:       !! 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)
        
        
        if (( xd%tau_V(i,j) >= (T_VL + T_sh) ) .and. TESF)  then !.and. (TESF) RRD added 
           xd%tau_V(i,j)     = 0.0_ReKi
          ! LESF=.FALSE. !also added this
        end if
        
  !bjj: update sigma1 to value at t + dt      
        OtherState%sigma1(i,j) = 1.0_ReKi
  
negative.code:       else if (Kafactor < 0 ) then 
           OtherState%sigma3(i,j) = 4.0_ReKi
        end if
        
        if ((.not. TESF) .and. (Kq_f*dalpha0 < 0.0_ReKi)) then
           OtherState%sigma3(i,j) = 1.0_ReKi
  