39   PetscScalar dy;         /* y step size */  member
119 …WORLD, "sum(p)*dw*dtheta at t = %3.2f = %3.6f\n", (double)t, (double)(sum * user->dx * user->dy)));  in PostStep()
141   user->dy = (user->ymax - user->ymin) / (N - 1);  in ini_bou()
207   *p2 = f * (p[j + 1][i] - p[j - 1][i]) / (2 * user->dy);  in adv2()
215   *p_diff = user->disper_coe * ((p[j - 1][i] - 2 * p[j][i] + p[j + 1][i]) / (user->dy * user->dy));  in diffuse()
231 … - p[j][i]) / user->dx + fthetac * p[j][i] - user->disper_coe * (p[j + 1][i] - p[j][i]) / user->dy;  in BoundaryConditions()
233 … - p[j][i]) / user->dx + fthetac * p[j][i] - user->disper_coe * (p[j][i] - p[j - 1][i]) / user->dy;  in BoundaryConditions()
235 …[j][i - 1]) / user->dx + fthetac * p[j][i] - user->disper_coe * (p[j + 1][i] - p[j][i]) / user->dy;  in BoundaryConditions()
237 …[j][i - 1]) / user->dx + fthetac * p[j][i] - user->disper_coe * (p[j][i] - p[j - 1][i]) / user->dy;  in BoundaryConditions()
239 …/ (user->dx) + fthetac * p[j][i] - user->disper_coe * (p[j + 1][i] - p[j - 1][i]) / (2 * user->dy);  in BoundaryConditions()
241 …/ (user->dx) + fthetac * p[j][i] - user->disper_coe * (p[j + 1][i] - p[j - 1][i]) / (2 * user->dy);  in BoundaryConditions()
243 …1]) / (2 * user->dx) + fthetac * p[j][i] - user->disper_coe * (p[j + 1][i] - p[j][i]) / (user->dy);  in BoundaryConditions()
245 …1]) / (2 * user->dx) + fthetac * p[j][i] - user->disper_coe * (p[j][i] - p[j - 1][i]) / (user->dy);  in BoundaryConditions()
345             val[nc++] = -user->disper_coe / user->dy;  in IJacobian()
348             val[nc++] = -fwc / user->dx + fthetac + user->disper_coe / user->dy;  in IJacobian()
355             val[nc++] = user->disper_coe / user->dy;  in IJacobian()
358             val[nc++] = -fwc / user->dx + fthetac - user->disper_coe / user->dy;  in IJacobian()
365             val[nc++] = -user->disper_coe / user->dy;  in IJacobian()
368             val[nc++] = fwc / user->dx + fthetac + user->disper_coe / user->dy;  in IJacobian()
375             val[nc++] = user->disper_coe / user->dy;  in IJacobian()
378             val[nc++] = fwc / user->dx + fthetac - user->disper_coe / user->dy;  in IJacobian()
385             val[nc++] = -user->disper_coe / (2 * user->dy);  in IJacobian()
388             val[nc++] = user->disper_coe / (2 * user->dy);  in IJacobian()
398             val[nc++] = -user->disper_coe / (2 * user->dy);  in IJacobian()
401             val[nc++] = user->disper_coe / (2 * user->dy);  in IJacobian()
414             val[nc++] = -user->disper_coe / user->dy;  in IJacobian()
417             val[nc++] = user->disper_coe / user->dy + fthetac;  in IJacobian()
427             val[nc++] = user->disper_coe / user->dy;  in IJacobian()
430             val[nc++] = -user->disper_coe / user->dy + fthetac;  in IJacobian()
435 …= (user->ws / (2.0 * user->H)) * (user->PM_min - user->Pmax * PetscSinScalar(xi)) / (2 * user->dy);  in IJacobian()
436 …/ (2 * user->H), 2) * user->q * (1.0 - PetscExpScalar(-t / user->lambda))) / (user->dy * user->dy);  in IJacobian()