search

NOAA-EMC / RDASApp

Regional DAS
GNU Lesser General Public License v2.1
1 stars 11 forks source link

Assimilate ATMS DBNet and ERAs data with the Normal feed data together in MPAS_JEDI #203

Open xyzemc opened 1 week ago

xyzemc commented 1 week ago

After confirmed the bufr2ioda converts the combined DBNet, ERAs and Normal feed ATMS bufr data correctly, I am starting to test and compare the result of assimilating the combined IODA data in MPAS JEDI with only assimilate Normal feed data. For this testing, we will separately test the atms_n20.yaml and atms_npp.yaml. To assimilate three ATMS bufr data together, the data conversion needs 2 steps:

  1. Read and write out bufr. In the python code , the Normal feed, DBNet, and ERAs bufr data will be read in with their own bufr table. Then write out to a combined bufr file with a new modified bufr table.
  2. Convert the new bufr file to IODA with bufr2ioda.x with the bug fixed bufr_query in PR: https://github.com/NOAA-EMC/bufr-query/pull/24
xyzemc commented 3 days ago

ATMS_N20 data Offline Domain Check image

Total data count read in JEDI is 11479, but only channel 7, 8, 9, 10, 11, 12, 13,14 have data passed the QC : image

xyzemc commented 3 days ago

Demo the result from channel 10: OmB before QC: image

OmB after QC: image

Bias image

guoqing-noaa commented 3 days ago

@xyzemc FYI, for the Offline Domain Check figure, it looks like you used an old version of the plotting script. Donnies updated that script so that the bottom box line takes the concave shape.

xyzemc commented 3 days ago

@xyzemc FYI, for the Offline Domain Check figure, it looks like you used an old version of the plotting script. Donnies updated that script so that the bottom box line takes the concave shape.

Thanks. I will check with Donny.

xyzemc commented 3 days ago

Analysis Potential Temperature increment at Level 20 image Analysis Potential Temperature increment at Level 25 image