Layer array input packages (READASARRAYS) not supported

[mjreno]

Grid input (e.g. RCHA) packages are currently not supported:

• test(rcha): support writing layer array input package #238
• Current failure is in structure.py structure_array(), which assumes list style formatted input, e.g. https://github.com/modflowpy/pyphoenix-project/actions/runs/18855284829/job/53801578235?pr=238
• Current test is formatting layer array based input as:

boundaries = np.full((20, 30), DNODATA, dtype=float)
for row in range(5, 15):
    boundaries[row, 0] = 100.0
for row in range(8, 12):
    boundaries[row, 29] = 95.0
for col in range(10, 20):
    boundaries[19, col] = 98.0

rch = Rcha(
    parent=gwf, recharge={0: boundaries}, print_input=True, save_flows=True, dims={"nper": 1}
)

[wpbonelli]

There are separate issues here, right? What syntax the initializer accepts for arrays, vs whether we can successfully load a grid-array or layer-array package from and/or write it to file. Likewise with #234.

The former is a duplicate of #235 and #189 so maybe you can move your initialization code examples to #189 @mjreno, and update this and #234 to focus on writing? I created #255 regarding the reading side and will work on that.