anchor.diff:             Enorm = 0.0D0
              Enorm = DOT_PRODUCT(ui_temp,feqv)
              IF(Enorm .GT. Eref/TOLF) THEN
  -               WRITE(*,*) "Solution is diverging, exit N-R"
  +               call WrScr( "Solution is diverging, exit N-R" )
                  piter=niter
              ELSEIF(Enorm .LE. Eref) THEN
                  RETURN
  
positive.code:                               ,INode   = i            &
                                ,Pos     = TmpPos       &
                                ,ErrStat = ErrStat2      &
                                ,ErrMess = ErrMsg2       )
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
         
     ENDDO
  
negative.code:                                   ,INode   = temp_id      &
                                    ,Pos     = TmpPos       &
                                    ,ErrStat = ErrStat2      &
                                    ,ErrMess = ErrMsg2       &
                                    , Orient = ??? ) !bjj: add orient=DCM ! Orientation (direction cosine matrix) of node; identity by default
              CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
         ENDDO
     ENDDO
  