anchor.diff:     DEALLOCATE(temp_GLL)
      DEALLOCATE(SP_Coef)
   
  -! Compute sectional propertities ( 6 by 6 stiffness and mass matrices)
  -! at Gauss points
  +   ! Compute sectional propertities ( 6 by 6 stiffness and mass matrices)
  +   ! at Gauss points
      CALL AllocAry(temp_ratio,p%ngp,p%elem_total,'temp_ratio',ErrStat2,ErrMsg2)
         CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
         if (ErrStat >= AbortErrLev) then
  
positive.code:    p%GlbRot(1:3,1:3) = TRANSPOSE(InitInp%GlbRot(1:3,1:3))
     p%gravity(:) = MATMUL(TRANSPOSE(p%GlbRot),InitInp%gravity(:))
  
  ! Analysis type: 1 Static 2 Dynamic
     p%analysis_type  = InputFileData%analysis_type
  ! Numerical damping coefficient: [0,1].
  ! No numerical damping if rhoinf = 1; maximum numerical damping if rhoinf = 0.
     p%rhoinf = InputFileData%rhoinf
  ! Time step size
     p%dt = InputFileData%DTBeam
  ! Compute generalized-alpha time integrator coefficients given rhoinf
     CALL AllocAry(p%coef,9,'GA2 coefficient',ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
        IF( ErrStat >= AbortErrLev ) RETURN
     p%coef(:) = 0.0D0
     CALL BD_TiSchmComputeCoefficients(p%rhoinf,p%dt,p%coef, ErrStat2, ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  ! Maximum number of iterations in Newton-Ralphson algorithm
     p%niter = InputFileData%NRMax
  ! Tolerance used in stopping criterion
     p%tol = InputFileData%stop_tol
  ! Total number of elements
     p%elem_total = InputFileData%member_total
  ! Number of nodes per elelemt
     p%node_elem  = InputFileData%order_elem + 1   
  ! Number of Gauss points
     p%ngp        = p%node_elem - 1
  ! Degree-of-freedom (DoF) per node
     p%dof_node   = 6
  ! Total numbder of (finite element) nodes
     p%node_total  = p%elem_total*(p%node_elem-1) + 1         
  ! Total numbder of (finite element) dofs
     p%dof_total   = p%node_total*p%dof_node   
  ! Compute coefficients for cubic spline fit, clamped at two ends
     CALL AllocAry(SP_Coef,InputFileData%kp_total-1,4,4,'Spline coefficient matrix',ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
        if (ErrStat >= AbortErrLev) then
  
negative.code:          return
        end if
     
     ! Compute initial position vector uuN0 in blade frame
     temp_int = p%node_elem * p%dof_node
     CALL AllocAry(p%uuN0,temp_int,p%elem_total,'uuN0 (initial position) array',ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  