search

dtcenter / MET

Model Evaluation Tools
https://dtcenter.org/community-code/model-evaluation-tools-met
Apache License 2.0
77 stars 24 forks source link

PB2NC will only output HGT observations for the ADPSFC message type in GDAS PrepBUFR files #125

Closed dwfncar closed 12 years ago

dwfncar commented 12 years ago

During testing, pb2nc would not write ADPSFC message type observations in GDAS files for any forecast field other that HGT (grib code 7) to the output NetCDF file. This bug was reproduced by both John and Paul in previous versions (v3.1 and v3.0.1). However, correspondence with Jonathan Case on and around 3/1 appeared to show pb2nc working correctly. The problem does not occur when converting NDAS PrepBUFR observations, nor does it appear to happen with other message types. Several tests were carried out including all/limited message types, all/limited station ids, text dumps and several GDAS files and all tests reproduced the problem. [MET-125] created by pgoldenb

dwfncar commented 12 years ago

Upon further investigation, GDAS PrepBUFR files appear to have only HGT (P) and Z observations with a quality mark of less than 9 on the top of the observation stack for each station. This was determined by running pb2nc with the -dump option and studying the observation stacks. When a quality_mark_threshold of 9 is used, the other variable observations are exported, as expected. by pgoldenb