Bug #16031

pycama does not generate so2 total column data for MPC QC report and website

Added by Mark ter Linden 8 months ago. Updated 8 months ago.

In Progress
Target version:
Start date:
Due date:
% Done:



The sulfurdioxide_total_vertical_column data is not extracted from NRTI files, the field is not present in the output file.

Attached are joborder, log file and ncdump.

so2.log View (2.04 MB) Mark ter Linden, 02/13/2019 03:23 PM

so2.ncdump (2.72 MB) Mark ter Linden, 02/13/2019 03:23 PM

so2.xml View (40.3 KB) Mark ter Linden, 02/13/2019 03:23 PM


#1 Updated by Maarten Sneep 8 months ago

Same symptoms observed in HCHO.

#2 Updated by Maarten Sneep 8 months ago

  • Status changed from New to In Progress
  • Assignee set to Maarten Sneep
  • Target version set to PyCAMA 0.10.0
  • % Done changed from 0 to 80

Cause found: The limit for the solar zenith angle where no retrievals are expected was set to 88 degrees. This is too high for both SO2 and HCHO, as they retrieve data for solar zenith angles up to 85 (?) degrees. A NRT granule may fall fully beyond 85 degrees, and will therefore only contain fill values in SO2 and HCHO. If the granule is fully beyond 88 degrees the granule will be removed, but if part of the granule has a SZA between 85 and 88 degrees, then the variable was removed. By setting the limit to 80 degrees, we still have the desired effect, but the condition is not triggered as frequently.

This will require a code update. Note that only NRT is affected by this.

#3 Updated by Maarten Sneep 8 months ago

  • % Done changed from 80 to 50

Need to rethink some assumptions, as the start of a day is over the pacific, and this may result in a full empty orbit, or at least empty for a significant part of the orbit.

Also available in: Atom PDF