Junya Okajima

I want to give different head values for each layer when inputting the head value for each cell as a side boundary condition. If using CHD package, can't input them for each layers. How can I input ?

I am in trouble because I can not create complex model including faults in the area. At hand, there are data of the positions of the 6 geological boundary surfaces represented by the distribution of points (x, y, z) and the positions of the 6 faults also represented by distribution of points(x, y, z). How can I make appropriate modeling? Which package should I use?