anchor.diff:           return
         end if
         
  -
  -   
  +      
      ! 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
            call cleanup()
            return
  -      end if
  -      
  -   SP_Coef(:,:,:) = 0.0D0
  +      end if   
  +      
  +   SP_Coef(:,:,:) = 0.0_BDKi
      temp_id = 0
      temp_id2 = 0
      DO i=1,InputFileData%member_total
  
positive.code: ! nodes along beam axis
  ! Nodal forces = K u
  !***************************************************************************************
     REAL(ReKi),        INTENT(IN   ):: Stif0(:,:,:) ! Element stiffness matrix
     REAL(ReKi),        INTENT(IN   ):: Mass0(:,:,:) ! Element stiffness matrix
     TYPE(BD_InputType),INTENT(IN   ):: u           ! Inputs at t
     REAL(ReKi),        INTENT(IN   ):: uuN(:) ! Displacement of Mass 1: m
     REAL(ReKi),        INTENT(IN   ):: vvN(:) ! Displacement of Mass 1: m
     INTEGER(IntKi),    INTENT(IN   ):: node_elem ! Node per element
     INTEGER(IntKi),    INTENT(IN   ):: dof_node ! Degrees of freedom per element
     INTEGER(IntKi),    INTENT(IN   ):: elem_total ! Total number of elements
     INTEGER(IntKi),    INTENT(IN   ):: ngp
     INTEGER(IntKi),    INTENT(IN   ):: quadrature
     REAL(ReKi),        INTENT(IN   ):: gw(:)
     REAL(ReKi),        INTENT(IN   ):: hhx(:,:)
     REAL(ReKi),        INTENT(IN   ):: hpx(:,:)
     REAL(ReKi),        INTENT(IN   ):: Jaco(:,:)
     REAL(ReKi),        INTENT(IN   ):: uu0(:,:)
     REAL(ReKi),        INTENT(IN   ):: E10(:,:)
     REAL(ReKi),        INTENT(  OUT):: Force(:)
     INTEGER(IntKi),    INTENT(  OUT):: ErrStat       ! Error status of the operation
     CHARACTER(*),      INTENT(  OUT):: ErrMsg        ! Error message if ErrStat /= ErrID_None
  
  
negative.code:                            p%Stif0_GL(1:6,1:6,temp_id) = InputFileData%InpBl%stiff0(1:6,1:6,k)
                             p%Mass0_GL(1:6,1:6,temp_id) = InputFileData%InpBl%mass0(1:6,1:6,k)
                         ELSE
                             temp66(:,:) = 0.0_BDKi
                             temp66(1:6,1:6) = (InputFileData%InpBl%stiff0(1:6,1:6,k)-InputFileData%InpBl%stiff0(1:6,1:6,k-1)) / &
                                               (InputFileData%InpBl%station_eta(k) - InputFileData%InpBl%station_eta(k-1))
                             p%Stif0_GL(1:6,1:6,temp_id) = temp66(1:6,1:6) * temp_ratio(j,i) + &
                                                           InputFileData%InpBl%stiff0(1:6,1:6,k-1) - &
                                                           temp66(1:6,1:6) * InputFileData%InpBl%station_eta(k-1)
                             temp66(:,:) = 0.0_BDKi
                             temp66(1:6,1:6) = (InputFileData%InpBl%mass0(1:6,1:6,k)-InputFileData%InpBl%mass0(1:6,1:6,k-1)) / &
                                               (InputFileData%InpBl%station_eta(k) - InputFileData%InpBl%station_eta(k-1))
                             p%Mass0_GL(1:6,1:6,temp_id) = temp66(1:6,1:6) * temp_ratio(j,i) + &
  