We construct an integrable colored six-vertex model whose partition function is a double Grothendieck polynomial. This gives an integrable systems interpretation of bumpless pipe dreams and recent results of Weigandt [arXiv:2003.07342] relating double Grothendieck polynomias with bumpless pipe dreams. For vexillary permutations, we then construct a new model that we call the semidual version model. We use our semidual model and the five-vertex model of Motegi and Sakai to given a new proof that double Grothendieck polynomials for vexillary permutations are equal to flagged factorial Grothendieck polynomials. Taking the stable limit of double Grothendieck polynomials, we obtain a new proof that the stable limit is a factorial Grothendieck polynomial as defined by McNamara. The states of our semidual model naturally correspond to families of nonintersecting lattice paths, where we can then use the Lindström-Gessel-Viennot lemma to give a determinant formula for double Schubert polynomials corresponding to vexillary permutations.