121 PetscErrorCode WaterReadData(WATERDATA *water, char *filename)  in WaterReadData()  argument
139   water->nvertex = water->nedge = 0;  in WaterReadData()
150       water->nvertex += nlines;  in WaterReadData()
151       water->njunction = nlines;  in WaterReadData()
156       water->nvertex += nlines;  in WaterReadData()
157       water->nreservoir = nlines;  in WaterReadData()
162       water->nvertex += nlines;  in WaterReadData()
163       water->ntank = nlines;  in WaterReadData()
168       water->nedge += nlines;  in WaterReadData()
169       water->npipe = nlines;  in WaterReadData()
174       water->nedge += nlines;  in WaterReadData()
175       water->npump = nlines;  in WaterReadData()
182   PetscCall(PetscCalloc1(water->nvertex, &water->vertex));  in WaterReadData()
183   PetscCall(PetscCalloc1(water->nedge, &water->edge));  in WaterReadData()
184   vert = water->vertex;  in WaterReadData()
185   edge = water->edge;  in WaterReadData()
189   for (i = 0; i < water->njunction; i++) {  in WaterReadData()
207   for (i = 0; i < water->nreservoir; i++) {  in WaterReadData()
223   for (i = 0; i < water->ntank; i++) {  in WaterReadData()
244   for (i = 0; i < water->npipe; i++) {  in WaterReadData()
271   for (i = 0; i < water->npump; i++) {  in WaterReadData()
297     for (j = water->npipe; j < water->npipe + water->npump; j++) {  in WaterReadData()
298       if (water->edge[j].pump.paramid == curve_id) {  in WaterReadData()
300         pump                                      = &water->edge[j].pump;  in WaterReadData()
312   for (j = water->npipe; j < water->npipe + water->npump; j++) {  in WaterReadData()
313     pump = &water->edge[j].pump;  in WaterReadData()