It is known that the moments of the maximum value of a one-dimensional conditional Brownian motion, the three-dimensional Bessel bridge with duration 1 started from the origin, are expressed using the Riemann zeta function. We consider a system of two Bessel bridges, in which noncolliding condition is imposed. We show that the moments of the maximum value is then expressed using the double Dirichlet series, or using the integrals of products of the Jacobi theta functions and its derivatives. Since the present system will be provided as a diffusion scaling limit of a version of vicious walker model, the ensemble of 2-watermelons with a wall, the dominant terms in longtime asymptotics of moments of height of 2-watermelons are completely determined. For the height of 2-watermelons with a wall, the average value was recently studied by Fulmek by a method of enumerative combinatorics.