anchor.diff:  
   
      END SUBROUTINE BD_ComputeIniCoef
  -
  -   FUNCTION BD_OuterProduct(vec1,vec2)
  -
  -   REAL(ReKi),INTENT(IN):: vec1(:),vec2(:)
  -   REAL(ReKi)::BD_OuterProduct(SIZE(vec1),SIZE(vec2))
  -
  -   INTEGER(IntKi)::i,j,n1,n2
  -
  -   n1=SIZE(vec1)
  -   n2=SIZE(vec2)
  -
  -   DO i=1,n1
  -       DO j=1,n2
  -           BD_OuterProduct(i,j) = vec1(i) * vec2(j)
  -       ENDDO
  -   ENDDO
  -
  -   END FUNCTION BD_OuterProduct
   
      SUBROUTINE BD_Static(t,n,u,utimes,p,x,xd,z,OtherState,ErrStat,ErrMsg)
   
  
positive.code:    ENDDO
  
     e3 = 0.0D0
     e3 = BD_CrossProduct(e1,e2)
     DO i=1,3
         Rr(i,3) = e3(i)
     ENDDO
  
negative.code:            ui(j+6) = ui_temp(j)
         ENDDO
  
         temp = TwoNorm(feqv)
         WRITE(13,*) i,temp
         piter=i !! ADDED BY NJ 3/18/14
  !       IF(i==1) Eref = SQRT(DOT_PRODUCT(ui_temp,feqv)*TOLF)
  