运行PB2NC

1. 将 METplus/ush 加入 PATH 环境变量

   export PATH=/home/wrf/METplus-3.0.2/ush:$PATH
   

2. 设置结果输出路径和MET系统路径

   $ cd /home/wrf/METplus-3.0.2
   $ vi parm/metplus_config/metplus_system.conf
   
   ...
   #  Location of METplus output files
   OUTPUT_BASE = /home/wrf/METout
   STAGING_DIR = {OUTPUT_BASE}/stage
   
   ##  Commonly used base MET variables
   
   ## Met install location
   MET_INSTALL_DIR = /home/wrf/METplus-3.0.2/build_components
   MET_BASE = {MET_INSTALL_DIR}/share/met
   
   
   ##  Output directories
   LOG_DIR = {OUTPUT_BASE}/logs
   TMP_DIR = {OUTPUT_BASE}/tmp
   ...
   

3. 设置输入数据路径

   $ cd /home/wrf/METplus-3.0.2
   $ vi parm/metplus_config/metplus_data.conf
   
   ....
   # Input data directories
   
   # This is the location of your input files for METplus
   INPUT_BASE = /home/wrf/com/prod
   ....
   

   输入数据的目录结构为：

   $ tree -L 2 -P "*prepbufr.tm00" /home/wrf/com/prod/
   
   ├── rap.20200620
   │   ├── nclprd
   │   ├── rap.t00z.prepbufr.tm00
   │   ├── rap.t01z.prepbufr.tm00
   │   ├── rap.t02z.prepbufr.tm00
   │   ├── rap.t03z.prepbufr.tm00
   │   ├── rap.t04z.prepbufr.tm00
   │   ├── rap.t05z.prepbufr.tm00
   │   ├── rap.t06z.prepbufr.tm00
   │   ├── rap.t07z.prepbufr.tm00
   │   ├── rap.t08z.prepbufr.tm00
   │   ├── rap.t09z.prepbufr.tm00
   │   ├── rap.t10z.prepbufr.tm00
   │   ├── rap.t11z.prepbufr.tm00
   │   ├── rap.t12z.prepbufr.tm00
   │   ├── rap.t13z.prepbufr.tm00
   │   ├── rap.t14z.prepbufr.tm00
   │   ├── rap.t15z.prepbufr.tm00
   │   ├── rap.t16z.prepbufr.tm00
   │   ├── rap.t17z.prepbufr.tm00
   │   ├── rap.t18z.prepbufr.tm00
   │   ├── rap.t19z.prepbufr.tm00
   │   ├── rap.t20z.prepbufr.tm00
   │   ├── rap.t21z.prepbufr.tm00
   │   ├── rap.t22z.prepbufr.tm00
   │   ├── rap.t23z.prepbufr.tm00
   │   └── wmo
   ├── rap.20200621
   │   ├── nclprd
   │   ├── rap.t00z.prepbufr.tm00
   │   ├── rap.t01z.prepbufr.tm00
   │   ├── rap.t02z.prepbufr.tm00
   │   ├── rap.t03z.prepbufr.tm00
   │   ├── rap.t04z.prepbufr.tm00
   │   ├── rap.t05z.prepbufr.tm00
   │   ├── rap.t06z.prepbufr.tm00
   │   ├── rap.t07z.prepbufr.tm00
   │   ├── rap.t08z.prepbufr.tm00
   │   ├── rap.t09z.prepbufr.tm00
   │   ├── rap.t10z.prepbufr.tm00
   │   ├── rap.t11z.prepbufr.tm00
   │   ├── rap.t12z.prepbufr.tm00
   │   ├── rap.t13z.prepbufr.tm00
   │   ├── rap.t14z.prepbufr.tm00
   │   ├── rap.t15z.prepbufr.tm00
   │   ├── rap.t16z.prepbufr.tm00
   │   ├── rap.t17z.prepbufr.tm00
   │   ├── rap.t18z.prepbufr.tm00
   │   ├── rap.t19z.prepbufr.tm00
   │   ├── rap.t20z.prepbufr.tm00
   │   ├── rap.t21z.prepbufr.tm00
   │   ├── rap.t22z.prepbufr.tm00
   │   ├── rap.t23z.prepbufr.tm00
   │   └── wmo
   ├── rap.20200622
   │   ├── rap.t00z.prepbufr.tm00
   │   ├── rap.t01z.prepbufr.tm00
   │   ├── rap.t02z.prepbufr.tm00
   │   └── wmo
   

4. 修改 PB2NC 参数文件

   $ vi /home/wrf/METplus-3.0.2/parm/use_cases/met_tool_wrapper/PB2NC/PB2NC.conf
   

   • 起始时间和间隔

     ...
     # Start time for METplus run - must match VALID_TIME_FMT
     VALID_BEG = 2020062000
     
     # End time for METplus run - must match VALID_TIME_FMT
     VALID_END = 2020062123
     
     # Increment between METplus runs (in seconds if no units are specified)
     #  Must be >= 60 seconds
     VALID_INCREMENT = 1H
     ...
     

   • 时间窗口

     ...
     # Time relative to each input file's valid time (in seconds if no units are specified) for data within the file to be
     # considered valid. Values are set in the 'obs_window' dictionary in the PB2NC config file
     PB2NC_WINDOW_BEGIN = -1800
     PB2NC_WINDOW_END = 1800
     ...
     

   • 区域和类型选择

     ...
     # Values to pass to pb2nc config file using environment variables of the same name.
     # See MET User's Guide for more information
     PB2NC_GRID =
     PB2NC_POLY =
     PB2NC_STATION_ID =
     PB2NC_MESSAGE_TYPE =
     ...
     

   • Offsets

     ...
     # list of offsets in the prepBUFR input filenames to allow. List is in order of preference
     # i.e. if 12, 6 is listed, it will try to use a 12 offset file and then try to use a 6 offset
     # if the 12 does not exist
     PB2NC_OFFSETS = 0
     ...
     

   • 输入和输出

     ...
     # directory containing input to PB2NC
     PB2NC_INPUT_DIR = {INPUT_BASE}
     
     # directory to write output from PB2NC
     PB2NC_OUTPUT_DIR = {OUTPUT_BASE}/pb2nc
     
     # End of [dir] section and start of [filename_templates] section
     [filename_templates]
     # Template to look for forecast input to PB2NC relative to PB2NC_INPUT_DIR
     PB2NC_INPUT_TEMPLATE = rap.{da_init?fmt=%Y%m%d}/rap.t{da_init?fmt=%2H}z.prepbufr.tm{offset?fmt=%2H}
     
     # Template to use to write output from PB2NC
     PB2NC_OUTPUT_TEMPLATE = rap.{da_init?fmt=%Y%m%d}/rap.t{da_init?fmt=%2H}z.prepbufr.tm{offset?fmt=%2H}.nc
     

5. 运行

   $ master_metplus.py -c /home/wrf/METplus-3.0.2/parm/use_cases/met_tool_wrapper/PB2NC/PB2NC.conf
   

6. 检查结果

   $ tree -L 2 /home/wrf/METout/pb2nc/ -h
   
   /home/wrf/METout/pb2nc/
   ├── [4.0K]  rap.20200620
   │   ├── [ 26K]  rap.t00z.prepbufr.tm00.nc
   │   ├── [111K]  rap.t01z.prepbufr.tm00.nc
   │   ├── [140K]  rap.t02z.prepbufr.tm00.nc
   │   ├── [ 92K]  rap.t03z.prepbufr.tm00.nc
   │   ├── [126K]  rap.t04z.prepbufr.tm00.nc
   │   ├── [111K]  rap.t05z.prepbufr.tm00.nc
   │   ├── [129K]  rap.t06z.prepbufr.tm00.nc
   │   ├── [134K]  rap.t07z.prepbufr.tm00.nc
   │   ├── [125K]  rap.t08z.prepbufr.tm00.nc
   │   ├── [116K]  rap.t09z.prepbufr.tm00.nc
   │   ├── [112K]  rap.t10z.prepbufr.tm00.nc
   │   ├── [110K]  rap.t11z.prepbufr.tm00.nc
   │   ├── [113K]  rap.t12z.prepbufr.tm00.nc
   │   ├── [ 95K]  rap.t13z.prepbufr.tm00.nc
   │   ├── [ 72K]  rap.t14z.prepbufr.tm00.nc
   │   ├── [ 76K]  rap.t15z.prepbufr.tm00.nc
   │   ├── [ 45K]  rap.t16z.prepbufr.tm00.nc
   │   ├── [ 32K]  rap.t17z.prepbufr.tm00.nc
   │   ├── [ 40K]  rap.t18z.prepbufr.tm00.nc
   │   ├── [ 31K]  rap.t19z.prepbufr.tm00.nc
   │   ├── [ 31K]  rap.t20z.prepbufr.tm00.nc
   │   ├── [ 32K]  rap.t21z.prepbufr.tm00.nc
   │   ├── [ 38K]  rap.t22z.prepbufr.tm00.nc
   │   └── [ 79K]  rap.t23z.prepbufr.tm00.nc
   ├── [4.0K]  rap.20200621
   │   ├── [112K]  rap.t00z.prepbufr.tm00.nc
   │   ├── [122K]  rap.t01z.prepbufr.tm00.nc
   │   ├── [139K]  rap.t02z.prepbufr.tm00.nc
   │   ├── [109K]  rap.t03z.prepbufr.tm00.nc
   │   ├── [112K]  rap.t04z.prepbufr.tm00.nc
   │   ├── [128K]  rap.t05z.prepbufr.tm00.nc
   │   ├── [112K]  rap.t06z.prepbufr.tm00.nc
   │   ├── [124K]  rap.t07z.prepbufr.tm00.nc
   │   ├── [131K]  rap.t08z.prepbufr.tm00.nc
   │   ├── [112K]  rap.t09z.prepbufr.tm00.nc
   │   ├── [120K]  rap.t10z.prepbufr.tm00.nc
   │   ├── [ 97K]  rap.t11z.prepbufr.tm00.nc
   │   ├── [123K]  rap.t12z.prepbufr.tm00.nc
   │   ├── [ 93K]  rap.t13z.prepbufr.tm00.nc
   │   ├── [ 90K]  rap.t14z.prepbufr.tm00.nc
   │   ├── [ 90K]  rap.t15z.prepbufr.tm00.nc
   │   ├── [ 35K]  rap.t16z.prepbufr.tm00.nc
   │   ├── [ 28K]  rap.t17z.prepbufr.tm00.nc
   │   ├── [ 28K]  rap.t18z.prepbufr.tm00.nc
   │   ├── [ 26K]  rap.t19z.prepbufr.tm00.nc
   │   ├── [ 27K]  rap.t20z.prepbufr.tm00.nc
   │   ├── [ 27K]  rap.t21z.prepbufr.tm00.nc
   │   ├── [ 39K]  rap.t22z.prepbufr.tm00.nc
   │   └── [ 81K]  rap.t23z.prepbufr.tm00.nc
   

7. 绘图

   MET提供了 plot_point_obs 绘制观测的分布情况

   $ ~/METplus-3.0.2/build_components/bin/plot_point_obs \
     ./rap.20200622/rap.t00z.prepbufr.tm00.nc \
     rap.t00z.prepbufr.tm00.nc.ps \
     -data_file rap.t09z.wrfprsf00.grib2 \
     -v 2
   

   rap.t00z.prepbufr.tm00.nc.ps 为观测分布的PS文件。