anchor.diff:         return
      end if
   
  -   Nuuu(:)            = 0.0D0
  -   Nrrr(:)            = 0.0D0
  -   EStif0_GL(:,:,:)   = 0.0D0
  -   EMass0_GL(:,:,:)   = 0.0D0
  -   DistrLoad_GL(:,:)  = 0.0D0
  -   elf(:)             = 0.0D0
  -   elk(:,:)           = 0.0D0
  +   Nuuu(:)            = 0.0_BDKi
  +   Nrrr(:)            = 0.0_BDKi
  +   EStif0_GL(:,:,:)   = 0.0_BDKi
  +   EMass0_GL(:,:,:)   = 0.0_BDKi
  +   DistrLoad_GL(:,:)  = 0.0_BDKi
  +   elf(:)             = 0.0_BDKi
  +   elk(:,:)           = 0.0_BDKi
   
      DO nelem=1,elem_total
          CALL BD_ElemNodalDisp(uuNf,node_elem,dof_node,nelem,Nuuu,ErrStat2,ErrMsg2)
  
positive.code:    ! Total number of elements
     p%elem_total = InputFileData%member_total      
     ! Total number of key points
     p%kp_total = InputFileData%kp_total  !bjj: not used
        
     ! Number of nodes per elelemt
     p%node_elem  = InputFileData%order_elem + 1   
  
negative.code:    end if   
  ! END Actuator
     
  
        
     call SetInitOut(p, InitOut, errStat, errMsg)
        call SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  