Dear Sander,
Please don’t laugh, but I am wondering if one can call coda from the command line, to fetch some parameter that is not already included in the ingestion def of the GOME2_L2_O3MOTO — HARP 1.24 documentation product. Namely the FittingRMS DETAILED_RESULTS/FittingRMS.
I have bash scripts that read in the GOME2 SO2-related info, using for e.g.
harpmerge -o ‘detailed_results=SO2;corrected=true;so2_column=2.5km’ -ap ‘bin_spatial(721,-90.0,0.25,1441,-180.,0.25);derive(latitude {latitude});derive(longitude {longitude});derive(SO2_column_numbe
r_density [DU]);derive(SO2_slant_column_number_density [DU]);exclude(weight,solar_zenith_angle_weight,SO2_column_number_density_count)’ -a ‘SO2_column_number_density_validity!& 255;keep(latitude, longitude,
latitude_bounds,longitude_bounds,SO2_column_number_density,solar_zenith_angle, cloud_fraction,SO2_column_number_density_validity, scan_direction_type, SO2_slant_column_number_density);cloud_fraction <= 0.2
;solar_zenith_angle <= 75; scan_direction_type==0’ -l $filename /mnt/lapaero_b/groups/lap_aero/ACSAF/Output/daily_gridded/metopab_anthro/$year/$date.nc
Which I think run faster than via the python interface.
If not possible, I’ll translate into python and call coda from there, no worries.
Thank you as always,
MariLiza