individual result timeseries for cells

[Sonic]

Is there a way to get the timeseries of each state variable for each cell?

I see getting the 'final' state as:

cells = region_model.get_cells()
c0 = cells[0]
elv = np.array([c.geo.mid_point().z for c in cells])
lwc = np.array([c.state.gs.lwc for c in cells])

but how would you get the timeseries of the state variable?

[sigbjorn.helset]

Hi,

You would first enable state-collection prior to calling the run-method.

ref. set_state_collection :

from shyft.api.pt_gs_k import PTGSKCellModel

help(PTGSKCellModel)

After the run, you can just pull the state-timeseries (they have n+1 values), by looking into the cell properties,

or use the region-model cell statistics method for each stack level to access e.g. the gs -state parameters.

Br Sigbjørn