Max Columns of Data Issue

[tjfarnum]

Greetings,

I’m attempting to set the max columns of data that is written out by flopy. I have the following code blocks.

However, this still writes out the external files in ncol x nrows. What am I doing wrong?

[christianlangevin]

Hey @tjfarnum, does setting the max columns the way you've shown work if you don't write to external files? I wonder if this is not implemented for the case where arrays are written to external files.

[tjfarnum]

Looks like internally its still not writing it out as 1 column.

[tjfarnum]

To add a further wrinkle. I tried setting a few different options in that same code block and they worked. But the max columns does not. I'm using Flopy 3.9.1.

[christianlangevin]

I poked around a bit and got it to work. All seems fine if you create the simulation from scratch, but if you load an existing model, then it behaves a little differently. The following works for me:

Create a simulation and write it.

import numpy as np
import flopy
ncol = 10
nrow = 10
k = np.random.random((nrow, ncol))
ws = './mymodel'
name = 'mymodel'
sim = flopy.mf6.MFSimulation(sim_name=name, sim_ws=ws, exe_name='mf6')
tdis = flopy.mf6.ModflowTdis(sim)
ims = flopy.mf6.ModflowIms(sim)
gwf = flopy.mf6.ModflowGwf(sim, modelname=name, save_flows=True)
dis = flopy.mf6.ModflowGwfdis(gwf, nrow=10, ncol=10)
ic = flopy.mf6.ModflowGwfic(gwf)
npf = flopy.mf6.ModflowGwfnpf(gwf, k=k, save_specific_discharge=True)
chd = flopy.mf6.ModflowGwfchd(gwf, stress_period_data=[[(0, 0, 0), 1.],
                                                       [(0, 9, 9), 0.]])
budget_file = name + '.bud'
head_file = name + '.hds'
oc = flopy.mf6.ModflowGwfoc(gwf,
                            budget_filerecord=budget_file,
                            head_filerecord=head_file,
                            saverecord=[('HEAD', 'ALL'), ('BUDGET', 'ALL')])
sim.set_all_data_external()
sim.write_simulation()

Then load the simulation, make the changes, and write to a new folder. In this case, the arrays are written with one value per line. It looks like the max_columns_auto_set needs to be set to False for things to work. I can't explain the logic behind this, but hope this helps you move forward.

import numpy as np
import flopy
sim = flopy.mf6.MFSimulation.load(sim_ws=ws)
sim.simulation_data.max_columns_auto_set = False
sim.simulation_data.max_columns_of_data = 1
sim.set_sim_path("./mymodel2")
sim.set_all_data_external()
sim.write_simulation()

[tjfarnum]

Yep that worked. It seems the key is that sim.simulation_data.max_columns_auto_set = False.

I ran each of those in a separate block and once the sim.set_all_data_external() was run it saved out the .txt files even before the write_simulation() was called. Seems like a weird quirk to me, but it works and thats good enough for me.

Thanks Chris!