anchor.diff:  
      ErrStat = ErrID_None
      ErrMsg  = ""
  -   RHS(:)  = 0.0D0
  +   RHS(:)  = 0.0_BDKi
   
      dof_elem = dof_node * node_elem
      rot_elem = (dof_node/2) * node_elem
  
positive.code:        call Cleanup()
         return
     end if
     Nuuu(:)  = 0.0D0
     Nrrr(:)  = 0.0D0
     Nvvv(:)  = 0.0D0
     EStif0_GL(:,:,:)   = 0.0D0
     EMass0_GL(:,:,:)   = 0.0D0
     DistrLoad_GL(:,:)  = 0.0D0
     elf(:)   = 0.0D0
  
     DO nelem=1,elem_total
         CALL BD_ElemNodalDisp(uuN,node_elem,dof_node,nelem,Nuuu,ErrStat2,ErrMsg2)
  
negative.code:           CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
     ENDIF
  
     IF(p%analysis_type .EQ. 2) THEN
         temp6(1:6) = -OS_tmp%Acc(1:6)
         temp_cc(:) = MATMUL(p%GlbRot,temp6(1:3))
  