绘图¶
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)