Thanks!
Maybe I’m just making assumptions about what model coords can do given my past experience with xarray. I was hoping to make use of automatic alignment/broadcasting/etc between model parameters and xarray datasets, but maybe that’s not possible?
For instance, in your example the dow_effect parameter has a weekday dim with named coordinates Monday, Tuesday, etc, but to actually use it in the likelihood it’s using integers to index into it. Is there any way to use the model coordinates to index into a parameter object? I can do dow_effect[0] but not dow_effect["Sunday"].