anchor.diff:                  DIRCOS=transpose(p%elemprops(K)%DirCos)! global to local
                   
                   CALL CALC_NODE_FORCES( DIRCOS,p%MOutLst3(I)%Me(:,:,1,J),p%MOutLst3(I)%Ke(:,:,1,J),(/uddout( L : L+5  ),uddout( L2 : L2+5 )/), &
  -                                 (/yout( L : L+5 ),       yout( L2 : L2+5 )/),  p%MoutLst3(I)%Fg(:,1,J),   K2,FM_elm,FK_elm) 
  +                                 (/yout( L : L+5 ),       yout( L2 : L2+5 )/),  p%MOutLst3(I)%Fg(:,1,J),   K2,FM_elm,FK_elm) 
                   !transform back to global, need to do 3 at a time since cosine matrix is 3x3
                   DO L=1,2  
                      FM_elm2((L-1)*3+1:L*3) = FM_elm2((L-1)*3+1:L*3) + matmul(p%elemprops(K)%DirCos,FM_elm((L-1)*3+1:L*3))  !sum forces at joint in GLOBAL REF
  
positive.code:                 L2=p%IDY((K3(2)-1)*6+1)! starting index for node K3(2) within yout
                  DIRCOS=transpose(p%elemprops(K)%DirCos)! global to local
                  CALL CALC_NODE_FORCES( DIRCOS,p%MOutLst2(I)%Me2(:,:,J),p%MOutLst2(I)%Ke2(:,:,J),(/uddout( L : L+5  ),uddout( L2 : L2+5 )/), &
                                   (/yout( L : L+5 ),       yout( L2 : L2+5 )/), p%MoutLst2(I)%Fg2(:,J),  K2,FM_elm,FK_elm) 
                   ! Store in All Outs
                   L=MaxOutPts+(I-1)*24+(J-1)*12+1!start index
                   L2=L+11
  
negative.code: SUBROUTINE SDOut_Init( Init, y,  p, misc, InitOut, WtrDpth, ErrStat, ErrMsg )
     TYPE(SD_InitType),        INTENT( INOUT ) :: Init                 ! data needed to initialize the output module     
     TYPE(SD_OutputType),      INTENT( INOUT ) :: y                    ! SubDyn module's output data
     TYPE(SD_ParameterType),   INTENT( INOUT ) :: p                    ! SubDyn module paramters
     TYPE(SD_MiscVarType),     INTENT( INOUT ) :: misc                 ! SubDyn misc/optimization variables
     TYPE(SD_InitOutputType ), INTENT( INOUT ) :: InitOut              ! SubDyn module initialization output data
     REAL(ReKi),               INTENT( IN    ) :: WtrDpth              ! water depth from initialization routine  
  