anchor.diff:         return
      end if
   
  -   Eref = 0.0D0
  +   Eref = 0.0_BDKi
      DO i=1,niter
          piter=i 
          CALL BD_GenerateStaticElement(uuNf,Mass0,Stif0,gravity,u,&
  
positive.code:                           ErrStat2, ErrMsg2) 
            CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  
         ui(:) = 0.0D0
         DO j=1,dof_total-6
             ui(j+6) = RHS_LU(j)
         ENDDO
  
negative.code:    ENDIF
     DO i=1,p%node_total
         temp_id = (i-1)*p%dof_node
         temp6(:) = temp_Force(temp_id+1:temp_id+6)
         temp_cc(:) = MATMUL(p%GlbRot,temp6(1:3))
         y%BldForce%Force(1,i) = temp_cc(2)
  