Motivated by an application to LDPC (low density parity check) algebraic geometry codes described by Voloch and Zarzar, we describe a computational procedure for establishing an upper bound on the arithmetic or geometric Picard number of a smooth projective surface over a finite field, by computing the Frobenius action on p-adic cohomology to a small degree of p-adic accuracy. We have implement...