Dear Sander,
We are trying to ingest the S5P_PAL__L2__AER_OT product using harp 1.20, running via
Python 3.11.6 | packaged by conda-forge | (main, Oct 3 2023, 10:40:35) [GCC 12.3.0]
IPython 8.16.1 – An enhanced Interactive Python.
Using the command:
operations = “;”.join([
“aerosol_optical_depth>50”, “keep(latitude_bounds,longitude_bounds,datetime_start,datetime_length,aerosol_optical_depth)”,
“bin_spatial(161,34,0.05,301,15,0.05)”,
“derive(longitude {longitude});derive(latitude {latitude})”])
reduce_operations = “squash(time, (latitude, longitude, latitude_bounds, longitude_bounds));bin()”
mean_no2=harp.import_product(files, operations, reduce_operations=reduce_operations)
harp.export_product(mean_no2, output_dir+‘s5p-NO2_L3_averaged_20Aug8Sep2023.nc’)
We saw in the other post the trouble for harp 1.19, which is solved by excluding the qa filter. Even though we are using the newest harp, we also tried that and are still getting the error:
create_daily_period_files_AoT.py’, wdir=‘/mnt/lapsat/groups/lapsat/s5p_no2_from_fires/s5p_aot/Programs’)
Traceback (most recent call last):
File ~/.conda/envs/newHarp/lib/python3.11/site-packages/spyder_kernels/py3compat.py:356 in compat_exec
exec(code, globals, locals)
File /mnt/lapsat/groups/lapsat/s5p_no2_from_fires/s5p_aot/Programs/create_daily_period_files_AoT.py:84
mean_no2=harp.import_product(files, operations, reduce_operations=reduce_operations)
File ~/.conda/envs/newHarp/lib/python3.11/site-packages/harp/_harppy.py:1257 in import_product
raise CLibraryError()
CLibraryError: /mnt/lapsat/groups/lapsat/s5p_no2_from_fires/s5p_aot/S5P_PAL__L2__AER_OT_20230821T093131_20230821T111300_30333_03_010100_20230831T000636.nc: unsupported product
Any ideas?
Thanks
MariLiza