pc_init Subroutine

  subroutine pc_init(self, fieldmat)
    implicit none
    class(pc_type), intent(inout) :: self
    class(fieldmat_type), intent(in) :: fieldmat

    !First allocate arrays
    call self%allocate

    !Now find cell locality
    call self%find_locality

    !Count avail data
    call self%count_avail(fieldmat)
  end subroutine pc_init