1 simple_packing 2 华东区域中尺度数值模式 grib2 数据分地面和高空 2 类数据 3 地面要素逐小时提供一个文件, 文件名按中华人民共和国气象行业标准 ( 详见附件 ) 如下 : Z_NAFP_C_BCSH_ _P_surface-warms-f09

Size: px

Start display at page:

Download "1 simple_packing 2 华东区域中尺度数值模式 grib2 数据分地面和高空 2 类数据 3 地面要素逐小时提供一个文件, 文件名按中华人民共和国气象行业标准 ( 详见附件 ) 如下 : Z_NAFP_C_BCSH_ _P_surface-warms-f09"

谐苍
5 years ago
Views:

1 附件 2 华东区域中尺度数值模式 grib2 数据说明系统名称产品种类格式华东区域中尺度数值天气预报模式 (SMS-WARMS) 2 米温度 2 米相对湿度 2 米露点温度地面温度东西方向 10 米风速南北方向 10 米风速地面气压海平面气压对流有效位能对流抑制系数能见度综合雷达回波累计降水累计降雪总云量高空温度高空相对湿度高空露点温度高空位势高度高空东西方向风速高空南北方向风速经度纬度 grib2 分辨率 9 公里范围涵盖全国区域时次 08 时 20 时两个时次时效 72 小时传输方式 FTP 上传数据服务区域上海数据生成时间 08 时起报的 grib2 数据一般在 14:30 分可以上传到服务器 20 时起报的 grib2 数据一般在 02:30 分可以上传达服务器投影方式及参数兰伯特投影 (Lambert Conformal) CEN_LAT=35.0; CEN_LON=105.0;TRUELAT1=30.0; TRUELAT2=60.0; STAND_LON=105.0; XLAT(0,0)= , XLONG(0,0)=

2 1 simple_packing 2 华东区域中尺度数值模式 grib2 数据分地面和高空 2 类数据 3 地面要素逐小时提供一个文件, 文件名按中华人民共和国气象行业标准 ( 详见附件 ) 如下 : Z_NAFP_C_BCSH_ _P_surface-warms-f09.BIN, 其中,Z: 为数据类型编码方式, 为不符合 WMO 编码格式的气象传输标识 ; NAFP: 为数据类型识别, 数值预报模式获得的预报产品 ; C: 为数据生成中心编码方式, 按编报中心进行编码 ; BCSH: 为数据生成中心标识, 上海 ; 为文件的生成时间 yyyymmddhhmmss, 使用国际协调时 (UTC); P: 为文件属性, 预报产品 ; surface-warms-f09: 为自由字段,surface 表示地面预报文件,warms 为华东区域中尺度数值模式简称,f09 为预报时间 ; BIN: 为文件类型标识, 按 WMO 规定的 GRIB 码格式编码的文件 ; 每个文件正常大小约 17M 地面要素及单位说明表要素变量名单位 2 米温度 TMP_P0_L103_GLC0 K 2 米相对湿度 SPFH_P0_L103_GLC0 % 2 米露点温度 DPT_P0_L103_GLC0 K 地面温度 TMP_P0_L1_GLC0 K 东西方向 10 米风速 UGRD_P0_L103_GLC0 m s-1 南北方向 10 米风速 VGRD_P0_L103_GLC0 m s-1 地面气压 PRES_P0_L1_GLC0 Pa 海平面气压 PRES_P0_L101_GLC0 Pa 对流有效位能 CAPE_P0_L1_GLC0 J kg-1 对流抑制系数 CIN_P0_L1_GLC0 J kg-1 能见度 VIS_P0_L1_GLC0 m 综合雷达回波 REFC_P0_L10_GLC0 db 逐小时累计降水 APCP_P8_L1_GLC0_acc kg m-2(mm)

3 逐小时累计降雪 WEASD_P0_L1_GLC0 kg m-2(mm) 总云量 TCDC_P0_L10_GLC0 % 经度 ELON_P0_L1_GLC0 deg 纬度 NLAT_P0_L1_GLC0 deg 4 高空数据逐 3 小时提供 1 个文件, 文件名如下 : 按中华人民共和国气象行业标准 ( 详见附件 ) 格式如下 : Z_NAFP_C_BCSH_ _P_high-warms-f09.BIN, 其中,Z: 为数据类型编码方式, 为不符合 WMO 编码格式的气象传输标识 ; (UTC); NAFP: 为数据类型识别, 数值预报模式获得的预报产品 ; C: 为数据生成中心编码方式, 按编报中心进行编码 ; BCSH: 为数据生成中心标识, 上海 ; 为文件的生成时间 yyyymmddhhmmss, 使用国际协调时 P: 为文件属性, 预报产品 ; high.warms.f09: 为自由字段,high 表示高空预报文件,warms 为华东区域中尺度数值模式简称,f09 为预报时间 ; BIN: 为文件类型标识, 按 WMO 规定的 GRIB 码格式编码的文件 ; 每个文件正常大小约在 68 兆垂直层次有 10 层, 分别是 1000hPa 925hPa 850hPa 700hPa 600hPa 500hPa 400hPa 300hPa 200hPa 100hPa 高空要素及单位说明表要素变量名单位温度 TMP_P0_L100_GLC0 K 露点温度 DPT_P0_L100_GLC0 K 相对湿度 SPFH_P0_L100_GLC0 % 东西方向风速 UGRD_P0_L100_GLC0 m s-1 南北方向风速 VGRD_P0_L100_GLC0 m s-1 位势高度 HGT_P0_L100_GLC0 gpm 经度 ELON_P0_L1_GLC0 deg 纬度 NLAT_P0_L1_GLC0 deg 5 数据中的经纬度变量为 ELON_P0_L1_GLC0 和 NLAT_P0_L1_GLC0, 请勿使用变量

4 griblon_0 和 griblat_0 6 提供模式静态数据 : 模式经纬度地形高度和地表分类数据, 存放于下载网址中的 /warms9km/readme/ 目录, 二进制文件为 latlonluhgt.grd,ascii 文件为 latlonluhgt.dat, 生成上述两个文件的程序为 readlatlon2.f( 供读取数据时参考 ) 7 在 /warms9km/readme 提供了 NCL 脚本和 FORTRAN 程序, 可将 grib2 数据提取为二进制数据样本地面数据请参考 mk_sur.ncl readsurface.f, 高空数据请参考 mk_levdata.ncl readlevdata.f 以下提供简单的数据读取或绘图说明 ( 上海为例 ) 8 提供利用 wgrib2 方式直接提取数据转为二进制码和使用 GRADS 软件直接绘图的说明 ( 附件 3) 9 在 /warms9km/readme 提供 NCL 绘图范例 draw_case.ncl( 附件 4) 附件 3 应用 WGRIB2 解码和 GRADS 画图的方法 : 由于地图投影的问题, 不能直接采用 GRADS 提供的处理 grib2 格式资料的方式来查看资料, 这里提供一种间接的方法地面资料 : 如 Z_NAFP_C_BCSH_ _P_surface-warms-f00.BIN 首先使用 wgrib2 命令将其转换成 2 进制格式的文件 test.dat wgrib2 -s -header Z_NAFP_C_BCSH_ _P_surface-warmsf00.BIN wgrib2 -i Z_NAFP_C_BCSH_ _P_surface-warmsf00.BIN -bin test.dat 然后使用以下的 CTL 文件 open 即可 dset test.dat undef 9.999E+20 options sequential pdef lcc xdef 1101 linear

5 ydef 551 linear tdef 1 linear 12Z06apr2016 1hr zdef 1 levels 1 vars 15 SLP 0 99 Sea level Pressure [Pa] REFC 0 99 Composite reflectivity [db] VIS 0 99 Visibility [m] PRES 0 99 Surface Pressure [Pa] TMPS 0 99 Surface Temperature [K] WEASD 0 99 Water Equivalent of Accumulated Snow Depth [kg/m^2] TMP2m m Temperature [K] DPT2m m Dew Point Temperature [K] RH2m m Relative Humidity [%] U10m m U-Component of Wind [m/s] V10m m V-Component of Wind [m/s] APCP 0 99 Accumulated Precipitation [kg/m^2] CAPE 0 99 Convective Available Potential Energy [J/kg] CIN 0 99 Convective Inhibition [J/kg] CLD 0 99 Total Cloud Cover [%] ENDVARS 对于高空资料由于 grib2 格式的层次资料的顺序为 do k = 1, kk write(10) var1(:,:,k) write(10) var2(:,:,k).. end do 将多个变量输出一个文件的脚本如下 : #!/bin/csh wgrib2 -s -header Z_NAFP_C_BCSH_ _P_high-warms-f00.BIN grep ":HGT" wgrib2 -i Z_NAFP_C_BCSH_ _P_high-warmsf00.BIN -bin test.dat match ':TMP' -append -bin test.dat match ':RH' -append -bin test.dat match ':DPT' -append -bin test.dat match ':UGRD' -append -bin test.dat match ':VGRD' -append -bin test.dat 然后使用以下 CTL 文件 open 即可

6 dset test.dat options sequential undef 9.999E+20 pdef lcc xdef 1101 linear ydef 551 linear tdef 1 linear 12Z06apr2016 1h zdef 10 levels vars 6 HGT Height [m] TMP Temperature [K] RH Relative Humidity [%] DPT Dew Point Temperature [K] U U-Component of Wind [m/s] V V-Component of Wind [m/s] ENDVARS 附件 4 应用 NCL 绘图脚本范例 (draw_case.ncl): ******* ; draw_case.ncl ******* load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" begin ; open file and read in data fili = "Z_NAFP_C_BCSH_ _P_surface-warms-f00.BIN.grb" f = addfile (fili, "r") vpt = f->pres_p0_l1_glc0 lat2d = f->nlat_p0_l1_glc0 lon2d = f->elon_p0_l1_glc0 nlat = dimsizes(lat2d(:,0)) nlon = dimsizes(lat2d(0,:)) ; create plot wks = gsn_open_wks("png","lcnative") ; send graphics to PNG file res = True ; plot mods desired res@timainstring = "Native Lambert Conformal" res@mplimitmode = "Corners" ; choose range of map

7 = lat2d(0,0) res@mpleftcornerlonf = lon2d(0,0) res@mprightcornerlatf = lat2d(nlat-1,nlon-1) res@mprightcornerlonf = lon2d(nlat-1,nlon-1) res@mpprojection = "LambertConformal" res@mplambertparallel1f = 30.0 res@mplambertparallel2f = 60.0 res@mplambertmeridianf = ; ncl adds from grib file res@pmtickmarkdisplaymode = "Always" res@mpfillon = False ; turn off map fill res@mpgeophysicallinecolor = "black" res@mpgeophysicallinethicknessf = 1. res@mpdatasetname = "Earth..4" res@mpoutlinedraworder = "PostDraw" ; draw continental outline last res@mpgridspacingf = 20 res@mpprovinciallinecolor = "Black" res@mpnationallinecolor = "Black" res@mpoutlineboundarysets = "AllBoundaries" res@tfdondcoverlay = True res@cnfillon = True ; color plot desired res@cnlineson = False ; turn off contour lines res@cnfillpalette = "gui_default" ; set color map res@gsnaddcyclic = False ; regional data plot = gsn_csm_contour_map(wks,vpt(:,:),res) end

<4D F736F F F696E74202D20B1D7B7C5B9CFB0F2A5BBB1D0A7F7>

<4D F736F F F696E74202D20B1D7B7C5B9CFB0F2A5BBB1D0A7F7> 斜溫圖簡介怎麼畫斜溫圖如何求取大氣物理參數逆溫層現象大氣穩定度探空熱力圖的種類探空熱力圖種類繁多, 主要有下面幾種能量圖 (Emagram)((-lnp)-T) 溫熵圖 (Tephigram)(lnθ-T) Refsdal 圖 ((-Tlnp)-lnT ) 司徒夫圖 (Stuve

1 simple_packing 2 华东区域中尺度数值模式 grib2 数据分地面和高空 2 类数据 3 地面要素逐小时提供一个文件, 文件名按 中华人民共和国气象行业标准 ( 详见附件 ) 如下 : Z_NAFP_C_BCSH_ _P_surface-warms-f09

1 simple_packing 2 华东区域中尺度数值模式 grib2 数据分地面和高空 2 类数据 3 地面要素逐小时提供一个文件, 文件名按中华人民共和国气象行业标准 ( 详见附件 ) 如下 : Z_NAFP_C_BCSH_ _P_surface-warms-f09