anchor.diff:        ! must initialize these because BD_AssembleStiffK and BD_AssembleRHS are INOUT
      m%RHS    =  0.0_BDKi
      m%StifK  =  0.0_BDKi
  -   
  +
         ! These values have not been set yet for the QP
      CALL BD_QPData_mEta_rho( p,m )            ! Calculate the \f$ m \eta \f$ and \f$ \rho \f$ terms
   
  
positive.code:       y%WriteOutput(i) = p%OutParam(i)%SignM * AllOuts( p%OutParam(i)%Indx )
     end do             ! i - All selected output channels
  
     
     call cleanup()
     return
     
  CONTAINS
        SUBROUTINE Cleanup()
           CALL BD_DestroyContState(x_tmp, ErrStat2, ErrMsg2 )
  
negative.code: !-----------------------------------------------------------------------------------------------------------------------------------
  !> This subroutine is just to clean up the code a bit.  This is called between the BD_InputGlobalLocal and BD_BoundaryGA2 routines.
  !! It could probably live in the BD_InputGlobablLocal except for the call just before the BD_CalcIC call (though it might not matter there).
  !! FIXME: explore moving this up to InputGlobalLocal: Not easy because BD_InputGlobalLocal is used in Init_ContinousStates, which is
  !!          before Init_MiscVars.
  ! NOTE: This routine could be entirely removed if the u%DistrLoad arrays are used directly, but that would require some messy indexing.
  SUBROUTINE BD_DistrLoadCopy( p, u, m )
  