PstFrom Assistance

[tjfarnum]

Greetings,

Apologies if this is somewhere in the GMDSI notebooks, I have looked but couldn't find what I'm looking for.

Here's a simplified example of what I'm trying to accomplish. If I have a 2 layer model, Layer 1 is a mix of Sand and Clay, Layer 2 is all sand, and where the sand of layer 1 and layer 2 meet they should have the same K properties. For traditional PEST what I would do is set up 2 different zones, so layer 1 would be zone 1 (clay) and zone 2 (sand), layer 2 would be just zone 2 (sand). Then create a k_by_node.txt file which is just a list of node numbers for all layers and what zone it is. That along with the pilot point information and which zone each pilot falls in, feeds into plproc which kicks out some factors.txt files. Then at model run time, within the forward run command, the plproc input file handles the kriging across different layers.

So, within the pyemu eco system, is there a way to replicate this behavior of having 1 set of pilot points for Layer 2, but where zone 2 is the same in layer 1 and 2, have fac2real krige the same K values across layers? In essence, is it possible to krige by zone regardless of layer and not just by layer?

Thanks for any assistance!