Julia Simpson

Looking at the above, are you sure that you're calling different values to create ql_0, qh_0, tau_0, tauy_0, evap_0 as opposed to ql_1, qh_1, tau_1, tauy_1, evap_1? it looks like...

```python # Valid data range according to documentation # https://github.com/xgcm/aerobulk-python/blob/main/source/aerobulk/flux.py VALID_VALUE_RANGES = {'sst': [270, 320], 't_zt': [180, 330], 'hum_zt': [0, 0.08], 'u_zu': [-50, 50], 'v_zu': [-50, 50], 'slp': [80000, 110000],...

Or ideally, the below would run, with the xarray version ```python import numpy as np import pandas as pd import datetime as dt from datetime import datetime import xarray as...

The actual error message varies, but it usually is a problem with this line

Comparing the real observational dataset (which runs with both numpy and xarray wrappers) with the synthetic data created here:

> Comparing the real observational dataset (which runs with both numpy and xarray wrappers) with the synthetic data created here: > >

Inputs: ```python import numpy as np import pandas as pd import datetime as dt from datetime import datetime import xarray as xr import matplotlib.pyplot as plt from matplotlib import cm...

Could there be something telling in the way that data is created in the aerobulk-python test py file, here: (https://github.com/xgcm/aerobulk-python/blob/f1affeb275ae86ea21f3b5a59d83f9238d380cce/tests/create_test_data.py)?

[ocean_index_check.pdf](https://github.com/xgcm/aerobulk-python/files/14672846/ocean_index_check.pdf) when I run it like this, the ocean index/args_shrunk works for both the real and the synthetic data, but trying to call aerobulk with the synthetic data still crashes

Created this to outline the current state, where both np/xr are running the real data (ds_clean) but neither are running on the synthetic–though notably again, one timestep of the np...