221   PetscScalar fwc, fthetac;  in BoundaryConditions()  local
229     fwc     = (w * w / 2.0 - user->ws * w);  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
338           fwc     = (yi * yi / 2.0 - user->ws * yi);  in IJacobian()
342             val[nc++] = fwc / user->dx;  in IJacobian()
348             val[nc++] = -fwc / user->dx + fthetac + user->disper_coe / user->dy;  in IJacobian()
352             val[nc++] = fwc / user->dx;  in IJacobian()
358             val[nc++] = -fwc / user->dx + fthetac - user->disper_coe / user->dy;  in IJacobian()
362             val[nc++] = -fwc / user->dx;  in IJacobian()
368             val[nc++] = fwc / user->dx + fthetac + user->disper_coe / user->dy;  in IJacobian()
372             val[nc++] = -fwc / user->dx;  in IJacobian()
378             val[nc++] = fwc / user->dx + fthetac - user->disper_coe / user->dy;  in IJacobian()
382             val[nc++] = fwc / user->dx;  in IJacobian()
391             val[nc++] = -fwc / user->dx + fthetac;  in IJacobian()
395             val[nc++] = -fwc / user->dx;  in IJacobian()
404             val[nc++] = fwc / user->dx + fthetac;  in IJacobian()
408             val[nc++] = fwc / (2 * user->dx);  in IJacobian()
411             val[nc++] = -fwc / (2 * user->dx);  in IJacobian()
421             val[nc++] = fwc / (2 * user->dx);  in IJacobian()
424             val[nc++] = -fwc / (2 * user->dx);  in IJacobian()