Thanks for the response, I found that the unit used for SO2 is DU how do I derive this in code?
This is my base SO2 code:
product = harp.import_product(r"/Users/jpoetzscher/Downloads/SO2 OCT 1/*.nc", operations=‘SO2_column_number_density_validity>50;keep(latitude_bounds,longitude_bounds,SO2_column_number_density);bin_spatial(1801,-90,0.1,3601,-180,0.1);derive(SO2_column_number_density [DU])’, post_operations=“bin();squash(time, (latitude_bounds,longitude_bounds));derive(latitude {latitude});derive(longitude {longitude});exclude(latitude_bounds,longitude_bounds,latitude_bounds_weight,longitude_bounds_weight,count,weight)”)
Also when I try to convert from the base unit of CO (mol/m2) to ppbv I get an error: CLibraryError: unit ‘mol/m^2’ cannot be converted to unit ‘ppbv’ (in unit conversion of variable ‘CO_column_number_density’)
The code I use for CO is:
product = harp.import_product(r"/Users/jpoetzscher/Downloads/CO OCT 1/*.nc", operations=‘CO_column_number_density_validity>50;keep(latitude_bounds,longitude_bounds,CO_column_number_density);bin_spatial(1801,-90,0.1,3601,-180,0.1);derive(CO_column_number_density [ppbv])’, post_operations=“bin();squash(time, (latitude_bounds,longitude_bounds));derive(latitude {latitude});derive(longitude {longitude});exclude(latitude_bounds,longitude_bounds,latitude_bounds_weight,longitude_bounds_weight,count,weight)”)
Also, I’ve completed my monthly averages for NO2 and CH4 but i wanted to export them to Panoply where there is more map customizability and I get an issue when I try to plot either one:
"There was an error preparing the plot. Failed creating the data handler: Axis var $s is > 1.
I can plot the daily averages just fine but, when in Visan I run the command:
average = harp.import_product("/Users/jpoetzscher/Downloads/October CH4/*.nc", post_operations=“bin()”)
and then export this file, which is now the monthly average, to downloads with the following command:
harp.export_product(average, “/Users/jpoetzscher/Downloads/October CH4 Average.nc”)
and then import this monthly average file into panoply I get this error.
