164   PetscInt i, j, l, Mx, My, xs, ys, xm, ym;  in FormInitialGuess()  local
174       for (l = 0; l < N_SPECIES; l++) {  in FormInitialGuess()
176           if (l == 0) x[j][i].sp[l] = (ctx->x_inflow.sp[l] * ((PetscScalar)j) / (My - 1));  in FormInitialGuess()
177 …       else if (l == 1) x[j][i].sp[l] = (ctx->x_inflow.sp[l] * (1. - ((PetscScalar)j) / (My - 1)));  in FormInitialGuess()
178           else x[j][i].sp[l] = ctx->x_0.sp[l];  in FormInitialGuess()
189   PetscInt    i, j, l, m;  in FormIFunctionLocal()  local
216       for (l = 0; l < N_SPECIES; l++) {  in FormIFunctionLocal()
219           sxp = (x[j][i + 1].sp[l] - x[j][i].sp[l]) * dhx;  in FormIFunctionLocal()
222           sxp = (x[j][i].sp[l] - x[j][i - 1].sp[l]) * dhx;  in FormIFunctionLocal()
225           sxm = (x[j][i + 1].sp[l] - x[j][i].sp[l]) * dhx;  in FormIFunctionLocal()
226           sxp = (x[j][i].sp[l] - x[j][i - 1].sp[l]) * dhx;  in FormIFunctionLocal()
229           syp = (x[j + 1][i].sp[l] - x[j][i].sp[l]) * dhy;  in FormIFunctionLocal()
232           syp = (x[j][i].sp[l] - x[j - 1][i].sp[l]) * dhy;  in FormIFunctionLocal()
235           sym = (x[j + 1][i].sp[l] - x[j][i].sp[l]) * dhy;  in FormIFunctionLocal()
236           syp = (x[j][i].sp[l] - x[j - 1][i].sp[l]) * dhy;  in FormIFunctionLocal()
240 …if (l == 0) f[j][i].sp[l] = (x[j][i].sp[l] - ctx->x_inflow.sp[l] * ((PetscScalar)j) / (info->my - …  in FormIFunctionLocal()
241 …else if (l == 1) f[j][i].sp[l] = (x[j][i].sp[l] - ctx->x_inflow.sp[l] * (1. - ((PetscScalar)j) / (…  in FormIFunctionLocal()
242           else f[j][i].sp[l] = x[j][i].sp[l];  in FormIFunctionLocal()
245           f[j][i].sp[l] = xt[j][i].sp[l] * scale;  in FormIFunctionLocal()
246           u             = x[j][i].sp[l];  in FormIFunctionLocal()
247           if (j == 0) uyy = u - x[j + 1][i].sp[l];  in FormIFunctionLocal()
248           else if (j == info->my - 1) uyy = u - x[j - 1][i].sp[l];  in FormIFunctionLocal()
249           else uyy = (2.0 * u - x[j - 1][i].sp[l] - x[j + 1][i].sp[l]) * hxdhy;  in FormIFunctionLocal()
251           if (i != info->mx - 1) uxx = (2.0 * u - x[j][i - 1].sp[l] - x[j][i + 1].sp[l]) * hydhx;  in FormIFunctionLocal()
252           else uxx = u - x[j][i - 1].sp[l];  in FormIFunctionLocal()
258           f[j][i].sp[l] += f_advect + ctx->dispersivity * (uxx + uyy);  in FormIFunctionLocal()
266           for (l = 0; l < N_SPECIES; l++) {  in FormIFunctionLocal()
267             if (stoich(m, l) < 0) {  in FormIFunctionLocal()
269               rate *= PetscPowScalar(x[j][i].sp[l], PetscAbsScalar(stoich(m, l)));  in FormIFunctionLocal()
273           for (l = 0; l < N_SPECIES; l++) {  in FormIFunctionLocal()
274             f[j][i].sp[l] += -scale * stoich(m, l) * rate; /* Reaction term */  in FormIFunctionLocal()
287   PetscInt    i, j, l, Mx, My, xs, ys, xm, ym;  in ReactingFlowPostCheck()  local
306       for (l = 0; l < N_SPECIES; l++) {  in ReactingFlowPostCheck()
307         if (x[j][i].sp[l] < 0.) x[j][i].sp[l] = 0.;  in ReactingFlowPostCheck()