anchor.diff:     ! input variables
   
      INTEGER(IntKi), INTENT(IN   ):: N           ! Order of spectral element
  -   REAL(BDKi),     INTENT(  OUT):: x(:)        ! location of GLL nodes
  -   REAL(BDKi),     INTENT(  OUT):: w(:)        ! quadrature weights at GLL nodes
  +   REAL(BDKi),     INTENT(  OUT):: x(n+1)        ! location of GLL nodes
  +   REAL(BDKi),     INTENT(  OUT):: w(n+1)        ! quadrature weights at GLL nodes
      INTEGER(IntKi), INTENT(  OUT):: ErrStat     ! Error status of the operation
      CHARACTER(*),   INTENT(  OUT):: ErrMsg      ! Error message if ErrStat /= ErrID_None
   
  
positive.code: ! the rotation matrix
  !--------------------------------------------------
  
     REAL(BDKi),    INTENT(IN   ):: Rr(:,:)       ! Rotation Matrix
     REAL(BDKi),    INTENT(  OUT):: cc(:)         ! Crv paramteres
     INTEGER(IntKi),INTENT(  OUT):: ErrStat       ! Error status of the operation
     CHARACTER(*),  INTENT(  OUT):: ErrMsg        ! Error message if ErrStat /= ErrID_None
  
  
negative.code: ! This subroutine computes the rotation tensor (RT)
  ! given Wiener-Milenkovic rotation parameters
  !--------------------------------------------------
     REAL(BdKi),    INTENT(IN   ):: cc(3)
     REAL(BdKi),    INTENT(  OUT):: Rr(3,3)
     INTEGER(IntKi),INTENT(  OUT):: ErrStat       ! Error status of the operation
     CHARACTER(*),  INTENT(  OUT):: ErrMsg        ! Error message if ErrStat /= ErrID_None
  
  