38   PetscScalar dx;         /* x 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()
140   user->dx = (user->xmax - user->xmin) / (M - 1);  in ini_bou()
183   *p1 = f * (p[j][i + 1] - p[j][i - 1]) / (2 * user->dx);  in adv1()
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()
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()
434         c1        = (yi - user->ws) / (2 * user->dx);  in IJacobian()