资料同化能够为数值模式提供更高质量的初始场,改进预报效果。WRFDA 是应用较为广泛的资料同化系统,由于原始的观测资料往往格式各异,因此 WRFDA 采用 Little R 格式,解决原始资料到同化系统的过渡问题。Little R 也是一种 ASCII 的文本格式,从 MM5 模式时期就开始使用。将所需的观测资料转为Little R格式,是做资料同化的基础工作。
下面介绍 Little R 格式:
Little 格式包含三部分:Header record, Data record, Tail record。
Header record: 这行出现表示一个站点数据记录的开始。包含数据的经纬度、ID、name、FM-Code、海拔和日期等信息。
Data record: 记录数据本身包含气压、高度、温度、露点温度、风速、风向、U风、V风、相对湿度、厚度以及这些数据的质控码。
海洋浮标数据:
探空数据:
数据最后还包括一个结束的指示记录Ending recording,这行出现表明数据本身的结束,即气压和高度的变量值为-777777.00000。
PS:
多层数据,比如探空,可以有多条不同高度数据的记录。
GPSPW 和 GPSZTD 数据实际记录在 Header record 中, 而不是 Data record。Header record 的 "
Precipitable Water" 单位是厘米,"QC" 部分包含其观测误差。由于 QC 只能是整数,因此误差记录为毫米的十倍,比如误差是1.5毫米,那么 QC=15。
Tail record: 三个整型数(I7), 第一个整数记录此站点一共有多少条数据,第二第三个整数为0。Tail record的出现表示这个站点的数据结束。
更全面的信息可以参考 Little R 的官网介绍,有需要 python 语言编码 Little R 格式的demo可以关注公众号气海同途留言索取。
PS:GSI 也是一个应用比较广泛的资料同化系统,但是它只支持 Bufr 格式数据的同化,而 Bufr 格式的编码相对较难理解。WRFDA 的 obsproc.exe 处理观测资料时还可输出 prepbufr 格式的常规资料,可以将这个输出资料用于 GSI 同化系统。
参考:
https://www2.mmm.ucar.edu/wrf/users/wrfda/OnlineTutorial/Help/littler.html