Hello,
I have two problems :

1. Is Modflow solver is only for saturated/parially saturated zones.Means it
will not support simulation if initial head is less than bottom elavation of a
layer.

2. What is the significance of maximum stage in Lake package . My lake stage
during simulation crosses above this maximum stage. Why does it happen???

I am working with MODFLOW using GMS 6.0 interface.