anchor.diff:  
      ErrStat    = ErrID_None
      ErrMsg     = ""
  -   StifK(:,:) = 0.0D0
  -   RHS(:)     = 0.0D0
  +   StifK(:,:) = 0.0_BDKi
  +   RHS(:)     = 0.0_BDKi
   
      dof_elem = dof_node * node_elem
      rot_elem = (dof_node/2) * node_elem
  
positive.code: 
     ErrStat  = ErrID_None
     ErrMsg   = ""
     elf(:)   = 0.0D0
  
     fact = .FALSE.
  
  
negative.code:       end if   
        
     SP_Coef(:,:,:) = 0.0_BDKi
     temp_id = 1
     temp_id2 = 0
     DO i=1,InputFileData%member_total
         temp_id2= temp_id + InputFileData%kp_member(i) - 1
         CALL BD_ComputeIniCoef(InputFileData%kp_member(i),InputFileData%kp_coordinate(temp_id:temp_id2,1:4),&
                                SP_Coef(temp_id:temp_id2-1,1:4,1:4), ErrStat2, ErrMsg2)
  