41   PetscScalar disper_coe; /* Dispersion coefficient */  member
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 …f[j][i] = fwc * (p[j][i + 1] - p[j][i]) / user->dx + fthetac * p[j][i] - user->disper_coe * (p[j +…  in BoundaryConditions()
233 …f[j][i] = fwc * (p[j][i + 1] - p[j][i]) / user->dx + fthetac * p[j][i] - user->disper_coe * (p[j][…  in BoundaryConditions()
235 …f[j][i] = fwc * (p[j][i] - p[j][i - 1]) / user->dx + fthetac * p[j][i] - user->disper_coe * (p[j +…  in BoundaryConditions()
237 …f[j][i] = fwc * (p[j][i] - p[j][i - 1]) / user->dx + fthetac * p[j][i] - user->disper_coe * (p[j][…  in BoundaryConditions()
239 …f[j][i] = fwc * (p[j][i + 1] - p[j][i]) / (user->dx) + fthetac * p[j][i] - user->disper_coe * (p[j…  in BoundaryConditions()
241 …f[j][i] = fwc * (p[j][i] - p[j][i - 1]) / (user->dx) + fthetac * p[j][i] - user->disper_coe * (p[j…  in BoundaryConditions()
243 …[j][i + 1] - p[j][i - 1]) / (2 * user->dx) + fthetac * p[j][i] - user->disper_coe * (p[j + 1][i] -…  in BoundaryConditions()
245 …[j][i + 1] - p[j][i - 1]) / (2 * user->dx) + fthetac * p[j][i] - user->disper_coe * (p[j][i] - p[j…  in BoundaryConditions()
282 …user->disper_coe = PetscPowScalar((user->lambda * user->ws) / (2 * user->H), 2) * user->q * (1.0 -…  in IFunction()
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()