anchor.diff:           N(3, jsize6) = 0.0D0
            
           END DO
  -        
  -        !--------------------------------------------------
  -        
  -        w_w_w_det = w_w_lx*2.0D0*det
  -        
  -        !--------------------------------------------------
  -        
  +                
           DO isize = 1, ndof*nn
            
  -         !vect(isize)                               &
  -         != vect(isize)                             &
  -         ! +w_w_w_det*( N(1, isize)*e3_hat(1)       &
  -         !             +N(2, isize)*e3_hat(2)       &
  -         !             +N(3, isize)*e3_hat(3) )*val 
            vect(isize)                                    &
            = vect(isize)                                  &
  -          +w_w_w_det*( N(1, isize)*g1_cross_g2(1)       &
  +          +w_w_lx*( N(1, isize)*g1_cross_g2(1)       &
                         +N(2, isize)*g1_cross_g2(2)       &
                         +N(3, isize)*g1_cross_g2(3) )*val 
            
  
positive.code:        S(3, 1) = Sv(5)
         S(1, 3) = Sv(5)
  	  
         !--------------------------------------------------
         
         stress(lx,1)                  &
          = ( S(1, 1)*g1(1)*g1(1)   &
  
negative.code:        !--------------------------------------------------------
         
        END DO
             
  !--------------------------------------------------------------------
        
        ! MITC4
  