anchor.diff:             ui(j+6) = ui_temp(j)
          ENDDO
   
  -       temp = BD_Norm(feqv)
  +       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)
  
positive.code:    ENDDO
  
     e3 = 0.0D0
     e3 = BD_CrossProduct(e1,e2)
     DO i=1,3
         Rr(i,3) = e3(i)
     ENDDO
  
negative.code:        e1(i) = Coef(i,2) + 2.0D0*Coef(i,3)*temp + 3.0D0*Coef(i,4)*temp*temp
     ENDDO
  !   WRITE(*,*) "e1",e1(:)
     e1(:) = e1(:)/TwoNorm(e1)
     Twist_Angle = 0.0D0
     Twist_Angle = Coef(4,1) + Coef(4,2)*temp + Coef(4,3)*temp*temp + Coef(4,4)*temp*temp*temp
  !   WRITE(*,*) "e1",e1(:)
  
  
     END SUBROUTINE BD_ComputeIniNodalPosition
  
     SUBROUTINE BD_CrvExtractCrv(Rr,cc,ErrStat,ErrMsg)
     !--------------------------------------------------
  