anchor.diff:             ENDDO
          ENDDO
      ENDDO
  -
      DEALLOCATE(temp_GL)
      DEALLOCATE(temp_ratio)
   ! Physical damping flag and 6 damping coefficients
      p%damp_flag  = InputFileData%InpBl%damp_flag
      CALL AllocAry(p%beta,6,'Damping coefficient',ErrStat2,ErrMsg2)
         CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  +   p%beta(:)  = 0.0D0
      p%beta(:)  = InputFileData%InpBl%beta(:)
   
      IF( ErrStat >= AbortErrLev ) RETURN
  
positive.code: ! Compute blade/member/segment lengths and the ratios between member/segment and blade lengths
     CALL AllocAry(p%member_length,InputFileData%member_total,2,'member length array',ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
     CALL AllocAry(p%segment_length,InputFileData%kp_total-1,3,'segment length array',ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
     CALL BD_ComputeMemberLength(InputFileData%member_total,InputFileData%kp_member,SP_Coef,&
                                 p%segment_length,p%member_length,p%blade_length,&
                                 ErrStat2,ErrMsg2)
  
negative.code: ! 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 )
     SP_Coef(:,:,:) = 0.0D0
     temp_id = 0
     temp_id2 = 0
     DO i=1,InputFileData%member_total
  