181   PetscInt comp;  in f0_zero_u()  local
182   for (comp = 0; comp < dim; ++comp) f0[comp] = u[comp];  in f0_zero_u()
188   PetscInt    comp;  in f0_constant_u()  local
191   for (comp = 0; comp < dim && comp < 3; ++comp) f0[comp] = u[comp] - wind[comp];  in f0_constant_u()
196   PetscInt comp;  in f1_constant_u()  local
197   for (comp = 0; comp < dim * dim; ++comp) f1[comp] = 0.0;  in f1_constant_u()
213   PetscInt comp;  in f0_lap_u()  local
214   for (comp = 0; comp < dim; ++comp) f0[comp] = 4.0;  in f0_lap_u()
219   PetscInt comp, d;  in f1_lap_u()  local
220   for (comp = 0; comp < dim; ++comp) {  in f1_lap_u()
221     for (d = 0; d < dim; ++d) f1[comp * dim + d] = u_x[comp * dim + d];  in f1_lap_u()