Dear forum,
I am trying to use avl in a Jupyter notebook to plot L3 data on a map with a command like: avl.Geo(l3product, varname, colorrange=(0.,6.e16))
l3product is the result of an harp.import_product command.
However, I get an error message, saying UnboundLocalError: cannot access local variable 'data_type' where it is not associated with a value plus a bunch of diagnostics that boil down to: raise Exception('invalid lat/lon/data dimensions')
According to: https://github.com/stcorp/avl-demo-lps2022/blob/main/geo2d.md, plotting of L3 should be possible. The L3 data dimensions look okay to me.
Having a time dimension for grids is still something we need to support better in AVL (specifically the length=1 case), but for now you can explicitly remove the time dimension of the grid by using the squash operation (as is already done for the lat/lon variables):
You can add a post_operations argument to your import having: squash(time, (tropospheric_HCHO_column_number_density, weight)).