RAP 的绘图

RAP 的绘图由如下脚本控制:

${ROOT_DIR}/nwprod/rap.v4.0.12/jobs/JRAP_NCLPART1
${ROOT_DIR}/nwprod/rap.v4.0.12/jobs/JRAP_NCLPART2

调用:

${ROOT_DIR}/nwprod/rap.v4.0.12/scripts/exrap_nclpart1.sh.ecf
${ROOT_DIR}/nwprod/rap.v4.0.12/scripts/exrap_nclpart2.sh.ecf

这两个脚本调用具体的 *ncl* 脚本

具体的 ncl 脚本在:

cd ${ROOT_DIR}/nwprod/rap.v4.0.12/ncl/ncl/nclrr_rll

[~/nwprod/rap.v4.0.12/ncl/ncl/nclrr_rll]$ls *.ncl
WRFUserARW.ncl         rr_htxs_grmr.ncl     rr_sfc_ltng.ncl     rr_sfc_temp.ncl    yyrr_2m_temp.ncl    yyrr_sfc_pchg.ncl
griblook.ncl           rr_htxs_rwmr.ncl     rr_sfc_mcc.ncl      rr_sfc_totp.ncl    yyrr_500_temp.ncl   yyrr_sfc_ptyp.ncl
plot_hodo.ncl          rr_htxs_sh.ncl       rr_sfc_mnvv.ncl     rr_sfc_ulwrf.ncl   yyrr_700_temp.ncl   yyrr_sfc_pwtr.ncl
rr_10m_gust.ncl        rr_htxs_snmr.ncl     rr_sfc_mucp.ncl     rr_sfc_vis.ncl     yyrr_80m_wind.ncl   yyrr_sfc_rhpw.ncl
rr_10m_wind.ncl        rr_htxs_temp.ncl     rr_sfc_mulcp.ncl    rr_sfc_weasd.ncl   yyrr_850_temp.ncl   yyrr_sfc_shtfl.ncl
rr_10m_wind_track.ncl  rr_htxs_wind.ncl     rr_sfc_mxcp.ncl     rr_ua_ceil.ncl     yyrr_925_temp.ncl   yyrr_sfc_snod.ncl
rr_20cm_soilt.ncl      rr_nta_ulwrf.ncl     rr_sfc_pchg.ncl     rr_ua_ctop.ncl     yyrr_sfc_1hsnw.ncl  yyrr_sfc_solar.ncl
rr_20cm_soilw.ncl      rr_prxs_temp.ncl     rr_sfc_profile.ncl  rr_ua_rh.ncl       yyrr_sfc_acp.ncl    yyrr_sfc_tcc.ncl
rr_2ds_temp.ncl        rr_sfc_1hsnw.ncl     rr_sfc_ptyp.ncl     rr_ua_rh8.ncl      yyrr_sfc_acpcp.ncl  yyrr_sfc_temp.ncl
rr_2m_dewp.ncl         rr_sfc_acp.ncl       rr_sfc_pwtr.ncl     rr_ua_vort.ncl     yyrr_sfc_acsnw.ncl  yyrr_sfc_totp.ncl
rr_2m_ptemp.ncl        rr_sfc_acp_3h.ncl    rr_sfc_rhpw.ncl     rr_ua_vvel.ncl     yyrr_sfc_cape.ncl   yyrr_sfc_vis.ncl
rr_2m_rh.ncl           rr_sfc_acpcp.ncl     rr_sfc_shtfl.ncl    rr_ua_wind.ncl     yyrr_sfc_cin.ncl    yyrr_sfc_weasd.ncl
rr_2m_temp.ncl         rr_sfc_acpcp_3h.ncl  rr_sfc_skewt.ncl    rr_ua_wmag.ncl     yyrr_sfc_cref.ncl   yyrr_ua_ceil.ncl
rr_500_temp.ncl        rr_sfc_acsnw.ncl     rr_sfc_skewt1.ncl   skewt_func.ncl     yyrr_sfc_flru.ncl   yyrr_ua_ctop.ncl
rr_5cm_soilt.ncl       rr_sfc_cape.ncl      rr_sfc_skewt2.ncl   soilmois.ncl       yyrr_sfc_hcc.ncl    yyrr_ua_rh.ncl
rr_5cm_soilw.ncl       rr_sfc_cin.ncl       rr_sfc_skewt3.ncl   wrfxsecs.ncl       yyrr_sfc_hpbl.ncl   yyrr_ua_rh8.ncl
rr_700_temp.ncl        rr_sfc_cref.ncl      rr_sfc_skewt4.ncl   yyrr_10m_gust.ncl  yyrr_sfc_lcc.ncl    yyrr_ua_vort.ncl
rr_80m_wind.ncl        rr_sfc_flru.ncl      rr_sfc_snod.ncl     yyrr_10m_wind.ncl  yyrr_sfc_ltng.ncl   yyrr_ua_vvel.ncl
rr_850_temp.ncl        rr_sfc_hcc.ncl       rr_sfc_soilt.ncl    yyrr_2ds_temp.ncl  yyrr_sfc_mcc.ncl    yyrr_ua_wind.ncl
rr_925_temp.ncl        rr_sfc_hpbl.ncl      rr_sfc_soilw.ncl    yyrr_2m_dewp.ncl   yyrr_sfc_mucp.ncl   yyrr_ua_wmag.ncl
rr_bsfc_soilt.ncl      rr_sfc_lcc.ncl       rr_sfc_solar.ncl    yyrr_2m_ptemp.ncl  yyrr_sfc_mulcp.ncl
rr_bsfc_soilw.ncl      rr_sfc_lhtfl.ncl     rr_sfc_tcc.ncl      yyrr_2m_rh.ncl     yyrr_sfc_mxcp.ncl

每个脚本负责画出一个图形。绘制的图形统一存放在:

${ROOT_DIR}/com/rap/prod/rap.ccyymmdd/nclprd

[/com/rap/prod/rap.20200209/nclprd]$ls -la
total 448
drwxr-xr-x 4 hailang work    512 Feb  9 20:54 .
drwxrwxr-x 4 hailang work 131072 Feb 10 00:12 ..
drwxr-xr-x 2 hailang work  32768 Feb 10 00:02 dfi
drwxr-xr-x 2 hailang work  65536 Feb 10 00:14 full
[hailang@sya01n01 ~/chinawise/com/rap/prod/rap.20200209/nclprd]$ls dfi
cape_sfc.t17z.f00.png  mucp_sfc.t17z.f00.png   pwtr_sfc.t17z.f00.png  temp_sfc.t17z.f00.png   wind_850.t17z.f00.png
cape_sfc.t23z.f00.png  mucp_sfc.t23z.f00.png   pwtr_sfc.t23z.f00.png  temp_sfc.t23z.f00.png   wind_850.t23z.f00.png
cin_sfc.t17z.f00.png   mulcp_sfc.t17z.f00.png  rh_2m.t17z.f00.png     weasd_sfc.t17z.f00.png  wmag_250.t17z.f00.png
cin_sfc.t23z.f00.png   mulcp_sfc.t23z.f00.png  rh_2m.t23z.f00.png     weasd_sfc.t23z.f00.png  wmag_250.t23z.f00.png
cref_sfc.t17z.f00.png  mxcp_sfc.t17z.f00.png   snod_sfc.t17z.f00.png  wind_10m.t17z.f00.png   wmag_850.t17z.f00.png
cref_sfc.t23z.f00.png  mxcp_sfc.t23z.f00.png   snod_sfc.t23z.f00.png  wind_10m.t23z.f00.png   wmag_850.t23z.f00.png
dewp_2m.t17z.f00.png   ptemp_2m.t17z.f00.png   temp_2ds.t17z.f00.png  wind_250.t17z.f00.png
dewp_2m.t23z.f00.png   ptemp_2m.t23z.f00.png   temp_2ds.t23z.f00.png  wind_250.t23z.f00.png
gust_10m.t17z.f00.png  ptyp_sfc.t17z.f00.png   temp_2m.t17z.f00.png   wind_80m.t17z.f00.png
gust_10m.t23z.f00.png  ptyp_sfc.t23z.f00.png   temp_2m.t23z.f00.png   wind_80m.t23z.f00.png
[hailang@sya01n01 ~/chinawise/com/rap/prod/rap.20200209/nclprd]$ls full
acpcp_sfc.t17z.f02.png  hcc_sfc.t17z.f03.png    ptemp_2m.t23z.f01.png   temp_2ds.t23z.f00.png  vis_sfc.t23z.f03.png
acpcp_sfc.t17z.f03.png  hcc_sfc.t17z.f04.png    ptemp_2m.t23z.f02.png   temp_2ds.t23z.f01.png  vis_sfc.t23z.f04.png
acpcp_sfc.t17z.f04.png  hcc_sfc.t17z.f05.png    ptemp_2m.t23z.f03.png   temp_2ds.t23z.f02.png  vis_sfc.t23z.f05.png
acpcp_sfc.t17z.f05.png  hcc_sfc.t17z.f06.png    ptemp_2m.t23z.f04.png   temp_2ds.t23z.f03.png  vis_sfc.t23z.f06.png
acpcp_sfc.t17z.f06.png  hcc_sfc.t23z.f00.png    ptemp_2m.t23z.f05.png   temp_2ds.t23z.f04.png  vvel_700.t17z.f00.png
acpcp_sfc.t23z.f02.png  hcc_sfc.t23z.f01.png    ptemp_2m.t23z.f06.png   temp_2ds.t23z.f05.png  vvel_700.t17z.f01.png
acpcp_sfc.t23z.f03.png  hcc_sfc.t23z.f02.png    pwtr_sfc.t17z.f00.png   temp_2ds.t23z.f06.png  vvel_700.t17z.f02.png
acpcp_sfc.t23z.f04.png  hcc_sfc.t23z.f03.png    pwtr_sfc.t17z.f01.png   temp_2m.t17z.f00.png   vvel_700.t17z.f03.png
acpcp_sfc.t23z.f05.png  hcc_sfc.t23z.f04.png    pwtr_sfc.t17z.f02.png   temp_2m.t17z.f01.png   vvel_700.t17z.f04.png
acpcp_sfc.t23z.f06.png  hcc_sfc.t23z.f05.png    pwtr_sfc.t17z.f03.png   temp_2m.t17z.f02.png   vvel_700.t17z.f05.png
cape_sfc.t17z.f00.png   hcc_sfc.t23z.f06.png    pwtr_sfc.t17z.f04.png   temp_2m.t17z.f03.png   vvel_700.t17z.f06.png
cape_sfc.t17z.f01.png   lcc_sfc.t17z.f00.png    pwtr_sfc.t17z.f05.png   temp_2m.t17z.f04.png   vvel_700.t23z.f00.png
cape_sfc.t17z.f02.png   lcc_sfc.t17z.f01.png    pwtr_sfc.t17z.f06.png   temp_2m.t17z.f05.png   vvel_700.t23z.f01.png
cape_sfc.t17z.f03.png   lcc_sfc.t17z.f02.png    pwtr_sfc.t23z.f00.png   temp_2m.t17z.f06.png   vvel_700.t23z.f02.png
cape_sfc.t17z.f04.png   lcc_sfc.t17z.f03.png    pwtr_sfc.t23z.f01.png   temp_2m.t23z.f00.png   vvel_700.t23z.f03.png

注解

所有的脚本使用 ${ROOT_DIR}/nwprod/rap.v4.0.12/ncl/ncl/nclrr_rll/names.txt 中的变量名字典列表来读取所需的变量。可以利用 ncl 工具 ncl_filedump, 来列出 grib 文件中所包含的列表名,例如:

$ncl_filedump /sya/u/hailang/chinawise/com/rap/prod/rap.20200209/rap.t17z.wrfprsf06.grib2
 Copyright (C) 1995-2013 - All Rights Reserved
 University Corporation for Atmospheric Research
 NCAR Command Language Version 6.1.2
 The use of this software is governed by a License Agreement.
 See http://www.ncl.ucar.edu/ for more details.

Variable: f
Type: file
filename:       rap.t17z.wrfprsf06
path:   /sya/u/hailang/chinawise/com/rap/prod/rap.20200209/rap.t17z.wrfprsf06.grib2
   file global attributes:
   dimensions:
      ygrid_0 = 462
      xgrid_0 = 525
      lv_ISBL0 = 40
      lv_HTGL1 = 2
      lv_SPDL2 = 6
      lv_HTGL3 = 2
      lv_HTGL4 = 2
      lv_TMPL5 = 2
      lv_SPDL6 = 3
      lv_HTGL7 = 2
      lv_HTGL8 = 2
      lv_DBLL9 = 9
   variables:
      float TMP_P0_L1_GLC0 ( ygrid_0, xgrid_0 )
         center :       US National Weather Service - NCEP (WMC)
         production_status :    Operational products
         long_name :    Temperature
         units :        K
         _FillValue :   1e+20
         coordinates :  gridlat_0 gridlon_0
         grid_type :    Lambert Conformal can be secant or tangent, conical or bipolar
         parameter_discipline_and_category :    Meteorological products, Temperature
         parameter_template_discipline_category_number :        ( 0, 0, 0, 0 )
         level_type :   Ground or water surface
         level :         0
         forecast_time :        6
         forecast_time_units :  hours
         initial_time : 02/09/2020 (17:00)

         float TMP_P0_L105_GLC0 ( ygrid_0, xgrid_0 )
            center :       US National Weather Service - NCEP (WMC)
            production_status :    Operational products
            long_name :    Temperature
            units :        K
            _FillValue :   1e+20
            coordinates :  gridlat_0 gridlon_0
            grid_type :    Lambert Conformal can be secant or tangent, conical or bipolar
            parameter_discipline_and_category :    Meteorological products, Temperature
            parameter_template_discipline_category_number :        ( 0, 0, 0, 0 )
            level_type :   Hybrid level
            level :         1
            forecast_time :        6
            forecast_time_units :  hours
            initial_time : 02/09/2020 (17:00)

         float TMP_P0_2L108_GLC0 ( lv_SPDL2, ygrid_0, xgrid_0 )
            center :       US National Weather Service - NCEP (WMC)
            production_status :    Operational products
            long_name :    Temperature
            units :        K
            _FillValue :   1e+20
            coordinates :  gridlat_0 gridlon_0
            grid_type :    Lambert Conformal can be secant or tangent, conical or bipolar
            parameter_discipline_and_category :    Meteorological products, Temperature
            parameter_template_discipline_category_number :        ( 0, 0, 0, 0 )
            level_type :   Level at specified pressure difference from ground to level (Pa)
            forecast_time :        6
            forecast_time_units :  hours
            initial_time : 02/09/2020 (17:00)

RAP 的绘图样本

../_images/1hsnw_sfc.t21z.f30.png ../_images/3hap_sfc.t21z.f30.png ../_images/acpcp_sfc.t21z.f30.png ../_images/acsnw_sfc.t21z.f30.png ../_images/cape_sfc.t21z.f30.png ../_images/ceil.t21z.f30.png ../_images/cin_sfc.t21z.f30.png ../_images/cref_sfc.t21z.f30.png ../_images/ctop.t21z.f30.png ../_images/dewp_2m.t21z.f30.png ../_images/flru_sfc.t21z.f30.png ../_images/gust_10m.t21z.f30.png ../_images/hcc_sfc.t21z.f30.png ../_images/hpbl_sfc.t21z.f30.png ../_images/lcc_sfc.t21z.f30.png ../_images/ltng_sfc.t21z.f30.png ../_images/mcc_sfc.t21z.f30.png ../_images/mucp_sfc.t21z.f30.png ../_images/mulcp_sfc.t21z.f30.png ../_images/mxcp_sfc.t21z.f30.png ../_images/pchg_sfc.t21z.f30.png ../_images/ptemp_2m.t21z.f30.png ../_images/ptyp_sfc.t21z.f30.png ../_images/pwtr_sfc.t21z.f30.png ../_images/rh_2m.t21z.f30.png ../_images/rh_500.t21z.f30.png ../_images/rh_850.t21z.f30.png ../_images/rhpw_sfc.t21z.f30.png ../_images/snod_sfc.t21z.f30.png ../_images/solar_sfc.t21z.f30.png ../_images/tcc_sfc.t21z.f30.png ../_images/temp_2ds.t21z.f30.png ../_images/temp_2m.t21z.f30.png ../_images/temp_500.t21z.f30.png ../_images/temp_700.t21z.f30.png ../_images/temp_850.t21z.f30.png ../_images/temp_925.t21z.f30.png ../_images/temp_sfc.t21z.f30.png ../_images/totp_sfc.t21z.f30.png ../_images/vis_sfc.t21z.f30.png ../_images/vort_500.t21z.f30.png ../_images/vvel_700.t21z.f30.png ../_images/weasd_sfc.t21z.f30.png ../_images/wind_10m.t21z.f30.png ../_images/wind_250.t21z.f30.png ../_images/wind_80m.t21z.f30.png ../_images/wind_850.t21z.f30.png ../_images/wmag_250.t21z.f30.png ../_images/wmag_850.t21z.f30.png