anchor.diff:        else if (Kafactor < 0 ) then 
            OtherState%sigma3(i,j) = 4.0_ReKi
         end if
  -      
  -      !   ! We are testing this heirarchical logic instead of the above block 5/29/2015
  -      !if ( (xd%tau_V(i,j) <= 2.0_ReKi*T_VL) .and. (xd%tau_V(i,j) >= T_VL) ) then
  -      !   OtherState%sigma3(i,j) =  3.0_ReKi
  -      !else if (.not. TESF) then
  -      !   OtherState%sigma3(i,j) =  4.0_ReKi
  -      !else if ( VRTX .and. (xd%tau_V(i,j) <= T_VL) ) then
  -      !   if (Kafactor < 0.0_ReKi) then
  -      !      OtherState%sigma3(i,j) =  2.0_ReKi
  -      !   else
  -      !      OtherState%sigma3(i,j) = 1.0_ReKi
  -      !    end if           
  -      !else if (Kafactor < 0 ) then 
  -      !   OtherState%sigma3(i,j) =  4.0_ReKi
  -      !end if
  -      
  -      !!!if ( (.not. VRTX) .OR. (.not. TESF) ) then !RRD: trying to emulate the old AD14 SEPAR.f90 with SHIFT=NOT(TESF), go back to original commented!!! out above when done
  -      !!!       OtherState%sigma3(i,j) =  2_ReKi
  -      !!!   else
  -      !!!      OtherState%sigma3(i,j) = 1.0_ReKi
  -      !!!endif
         
         if ((.not. TESF) .and. (Kq_f*dalpha0 < 0.0_ReKi)) then
            OtherState%sigma3(i,j) = 1.0_ReKi
  
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:             end if
              
           end if
        
  !bjj: update sigma3 to value at t + dt      
           
  