anchor.diff:  !-----------------------------------------------------------------------------------------------------------------------------------
   SUBROUTINE BD_CrvMatrixH(cc,Hh)
   
  -   REAL(BDKi),INTENT(IN) ::cc(:)
  -   REAL(BDKi),INTENT(OUT)::Hh(:,:)
  +   REAL(BDKi),INTENT(IN) ::cc(3)
  +   REAL(BDKi),INTENT(OUT)::Hh(3,3)
   
      REAL(BDKi):: cf1,cf2,cf3,cq,ocq,aa,cb0,cb1,cb2,cb3
   
  
positive.code:    ! 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
     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
  
  