anchor.diff:                           , INode = 1                &
                            , Pos = TmpPos             &
                            , ErrStat   = ErrStat2     &
  -                         , ErrMess   = ErrMsg2       )
  +                         , ErrMess   = ErrMsg2      &
  +                         , Orient = ??? ) !bjj: add orient=DCM ! Orientation (direction cosine matrix) of node; identity by default
         CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
   
   
  
positive.code: 
         CALL BD_StaticUpdateConfiguration(ui,uuNf,node_total,dof_node)
         IF(i==niter .OR. piter==niter) THEN
             WRITE(*,*) "Solution does not converge after the maximum number of iterations"
             EXIT !! USE EXIT INSTEAD OF STOP, NJ 3/18/2014
         ENDIF
     ENDDO
  
     END SUBROUTINE BD_StaticSolution
  
     SUBROUTINE BD_GenerateStaticElement(uuN0,uuNf,Mass0,Stif0,gravity,u,&
                                         elem_total,node_elem,dof_node,ngp,StifK,RHS)
  
negative.code:            u_temp%DistrLoad%Moment(:,:) = u_interp%DistrLoad%Moment(:,:)/i*j
             CALL BD_StaticSolution(p%uuN0,x%q,p%Mass0_GL,p%Stif0_GL,p%gravity,u_temp,&
                                    p%node_elem,p%dof_node,p%elem_total,&
                                    p%dof_total,p%node_total,p%ngp,niter,piter, ErrStat2, ErrMsg2)
                 call SetErrStat(ErrStat2,ErrMsg2,ErrStat, ErrMsg, RoutineName)
             IF(niter .EQ. piter) EXIT
         ENDDO
         IF(piter .LT. niter) THEN
             i=0
         ELSE
             i=i+1
             call WrScr( "Warning: Load may be too large, BeamDyn will attempt to solve with additional steps.")
             call WrScr( "  Load_Step="//trim(num2lstr(i)) )
             x%q(:) = 0.0D0
         ENDIF
     ENDDO
     
     call cleanup()
     return
     
  contains
     subroutine cleanup()
        CALL BD_DestroyInput(u_interp, ErrStat2, ErrMsg2 )
        CALL BD_DestroyInput(u_temp,   ErrStat2, ErrMsg2 )
     end subroutine cleanup
  END SUBROUTINE BD_Static
  
  SUBROUTINE BD_StaticSolution(uuN0,uuNf,Mass0,Stif0,gravity,u,&
                                  node_elem,dof_node,elem_total,&
                                  dof_total,node_total,ngp,niter,piter, ErrStat,ErrMsg)
  
     REAL(ReKi),        INTENT(IN   ):: uuN0(:,:)
     REAL(ReKi),        INTENT(IN   ):: Mass0(:,:,:)
  