资料同化 | 将NCEP的prepbufr常规观测资料转little_r并合并国家站点数据

更多地学领域的文章，欢迎关注地学科学社公众号，加入交流群

WRF/WRFDA就是由NCAR中尺度和微尺度气象学 (MMM) 实验室提供的公开且供社区成员免费使用的中尺度气象模式和资料同化软件，WRF应用范围广泛非常广泛，从公里尺度的区域和中尺度建模到用于数千公里尺度的全球建模，WRFDA可用于为WRF模式预报提供初始条件，验证观测与模式之间的差异，进行观测系统的设计、监测和评估等。

NCEP使用全球数据同化系统（GDAS）将来自各种观测系统和仪器的数据内插到三维网格中, 用来初始化NCEP的全球预报系统模式（GFS），而这些数据是可供下载使用的，其中包括各种常规和非常规观测数据的bufr格式文件和处理过的常规观测数据prepbufr文件。

WRFDA可以直接同化这两种格式，非常方便，但由于这两类格式的解码和编码较为复杂，常看其数据的时空分布特征较为困难。本文旨在分享一个非常易于使用的常规观测prepbufr数据解码，以及与观测站点数据合并的方法。

WRFDA可以直接同化bufr和prepbufr数据，但仍存在一个缺陷，prepbufr中的观测数据为国际交换数据，可以非常方便的获取探空、飞机报、船舶等不容易获取的常规资料，虽然prepbufr中的地面站观测数据质量较高，但是相比于国家气象科学数据中心可获取的国家站数量来说要少很多，因此我们希望同时利用prepbufr和国家站数据作为WRFDA同化的常规观测数据文件。然而WRFDA只支持读取一个ob文件，即ob.bufr（prepbufr数据）或者ob.ascii，不像GSI可以分别指定，因此我们需要将两个数据合并。

笔者尝试了两个思路：一、将国家站数据处理为prepbufr格式，与下载的GDAS的prepbufr数据合并&