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

[xyzemc]

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: Fix scanline calculation for FFT used by ATMS remapping  bufr-query#24

[xyzemc]

ATMS_N20 data Offline Domain Check

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

[xyzemc]

Demo the result from channel 10:
OmB before QC:

OmB after QC:

Bias

[guoqing-noaa]

@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]

@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]

Analysis Potential Temperature increment at Level 20

Analysis Potential Temperature increment at Level 25