221   PetscScalar fwc, fthetac;  in BoundaryConditions()  local
228     fthetac = user->ws / (2 * user->H) * (user->PM_min - user->Pmax * PetscSinScalar(theta));  in BoundaryConditions()
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 …f[j][i] = fwc * (p[j][i + 1] - p[j][i - 1]) / (2 * user->dx) + fthetac * p[j][i] - user->disper_co…  in BoundaryConditions()
245 …f[j][i] = fwc * (p[j][i + 1] - p[j][i - 1]) / (2 * user->dx) + fthetac * p[j][i] - user->disper_co…  in BoundaryConditions()
336           PetscScalar fthetac, fwc;  in IJacobian()  local
337           fthetac = user->ws / (2 * user->H) * (user->PM_min - user->Pmax * PetscSinScalar(xi));  in IJacobian()
348             val[nc++] = -fwc / user->dx + fthetac + user->disper_coe / user->dy;  in IJacobian()
358             val[nc++] = -fwc / user->dx + fthetac - user->disper_coe / user->dy;  in IJacobian()
368             val[nc++] = fwc / user->dx + fthetac + user->disper_coe / user->dy;  in IJacobian()
378             val[nc++] = fwc / user->dx + fthetac - user->disper_coe / user->dy;  in IJacobian()
391             val[nc++] = -fwc / user->dx + fthetac;  in IJacobian()
404             val[nc++] = fwc / user->dx + fthetac;  in IJacobian()
417             val[nc++] = user->disper_coe / user->dy + fthetac;  in IJacobian()
430             val[nc++] = -user->disper_coe / user->dy + fthetac;  in IJacobian()