anchor.diff:         rho(1:3,1:3) = EMass0_GL(4:6,4:6,igp)
   
          CALL BD_GaussPointDataMass(hhx,hpx,Nvvv,temp_Naaa,RR0,node_elem,dof_node,&
  -                                  vvv,temp_aaa,vvp,mmm,mEta,rho)
  -       CALL BD_ElasticForce(E1,RR0,kapa,Stif,cet,Fc,Fd,Oe,Pe,Qe)
  -       CALL BD_GravityForce(mmm,mEta,gravity,Fg)
  +                                  vvv,temp_aaa,vvp,mmm,mEta,rho,ErrStat2,ErrMsg2)
  +       CALL BD_ElasticForce(E1,RR0,kapa,Stif,cet,Fc,Fd,Oe,Pe,Qe,ErrStat2,ErrMsg2)
  +       CALL BD_GravityForce(mmm,mEta,gravity,Fg,ErrStat2,ErrMsg2)
          Fd(:) = Fd(:) - Fg(:) - DistrLoad_GL(:,igp)
   
          DO i=1,node_elem
  
positive.code:    DO i=1,p%ngp*p%elem_total+2
         TmpPos(1:3) = p%GlbPos(1:3) + MATMUL(p%GlbRot,temp_L2(1:3,i))
         temp_CRV(:) = MATMUL(p%GlbRot,temp_L2(4:6,i))
         CALL BD_CrvMatrixR(temp_CRV,TmpDCM)
         TmpDCM(:,:) = TRANSPOSE(TmpDCM)
         CALL MeshPositionNode ( Mesh    = u%DistrLoad  &
                                ,INode   = i            &
  
negative.code:    CHARACTER(ErrMsgLen)        :: ErrMsg2                      ! Temporary Error message
     CHARACTER(*), PARAMETER     :: RoutineName = 'BD_GaussPointData'
  
     ErrStat = ErrID_None
     ErrMsg  = ""
  
     uuu = 0.0D0
     uup = 0.0D0
  