PARFLOW AND MODFLOW6

19 views

Skip to first unread message

MABOSSANI AKARA

unread,

Sep 15, 2026, 7:15:37 PM (10 days ago) Sep 15

to ParFlow

Hello All,

I am trying to use PARFLOW (Richards mode) to replicate a MODFLOW6 test case. The model grid consists of 1 layer, 3 rows, and 10 columns (see figure).

Time-dependent specified head is applied on the eastern edge (blue cells).

model.Patch.right.BCPressure.Type = "DirEquilRefPatch"
model.Patch.right.BCPressure.Cycle = "pers"
model.Patch.right.BCPressure.RefGeom = "domain"
model.Patch.right.BCPressure.RefPatch = "bottom"
model.Patch.right.BCPressure.s1.Value = 200.130
model.Patch.right.BCPressure.s2.Value = 200.261
model.Patch.right.BCPressure.s3.Value = 200.391
model.Patch.right.BCPressure.s4.Value = 200.521

....

Time-dependent specified flux is applied in the western part of the model domain (red cell).

model.Patch.top21.BCPressure.Type = "FluxVolumetric"

model.Patch.top21.BCPressure.Cycle = "pers"
model.Patch.top21.BCPressure.s1.Value = -2521.173
model.Patch.top21.BCPressure.s2.Value = -3042.345
model.Patch.top21.BCPressure.s3.Value = -3563.518
model.Patch.top21.BCPressure.s4.Value = -4084.691
model.Patch.top21.BCPressure.s5.Value = -4605.863

....

The results show that simulated heads from PARFLOW are consistently higher than those from MODFLOW6.

I am wondering if the difference is related to the equations being solved (variably saturated vs. fully saturated).

I will appreciate any suggestions.

Thanks,

-Mahawa

Reply all

Reply to author

Forward

0 new messages