登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 youngman 于 2020-5-20 18:17 编辑
GRIB1/2格式数据处理
2019/11/05
GRIB1/2主要用来存放数值天气分析和预报产品资料。wgrib/wgrib2命令行工具可高效处理GRIB1/GRIB2数据。以wgrib2安装为例:https://www.cpc.ncep.noaa.gov/products/wesley/wgrib2/
下载:ftp://ftp.cpc.ncep.noaa.gov/wd51we/wgrib2/wgrib2.tgz
安装:https://www.cpc.ncep.noaa.gov/products/wesley/wgrib2/compile_questions.html
注:官方推荐用gcc/gfortran编译器,icc和pgc都会出问题。在WSL上安装,亲测可用
tar -xzvf wgrib2.tgz
cd grib2
export CC=gcc
export FC=gfortran
make #注:默认只支持netcdf3因为netcdf4很大,修改makefile参数可选择需要的库,但不能同时支持netcdf3/4
wgrib2/wgrib2-config
帮助文档:https://www.ftp.cpc.ncep.noaa.gov/wd51we/wgrib2/tricks.ncep
结合正则表达式提取变量:-match, -not,-if
wgrib2 IN -match':TMP:' -grib ht.grb
wgrib2 IN | egrep“:(UGRD|VGRD):(2|3)00 mb” | sort -t: -k3,3 -k4,4 -k5,5 -k6,6 | wgrib2 IN -i-ieee sorted.bin
提取区域-small_grib:
wgrib2 in.grb2 -set_grib_type c2 -small_grib 109:110 35:37 region.grb2
插值到新的网格:
wgrib2 IN -new_grid_winds earth -new_grid latlon 0:360:1 90:181:-1 OUT
grib转nc文件:
wgrib2 in.grb2 -netcdf out.nc
参考:http://bbs.06climate.com/forum.php?mod=viewthread&tid=51744
---------------------------------------------------------------------------------------------------------------------------
大多常用气象数据处理软件均可处理GRIB数据:CDO,NCL,Grads,Python,Matlab等,简单介绍下:
CDO:https://code.mpimet.mpg.de/projects/cdo
安装:sudo apt install cdo
Grib转nc:cdo -f grb copy a.nc a.grb
nc转grib:cdo -f nc copy a.grb a.nc
按时间,区域,变量名提取:
cdo -selyear,2019 -selmonth,01 -selday,25 -selname,SSTK infile.nc outfile.nc
参考:https://www.cnblogs.com/jiangleads/p/11145774.html
---------------------------------------------------------------------------------------------------------------------------
NCL:http://www.ncl.ucar.edu/Applications/grib2netCDF.shtml
grib转nc:ncl_convert2nc file.grb 输出file.nc
查看grib文件信息:ncl_filedump file.grb
---------------------------------------------------------------------------------------------------------------------------
EC的Grib格式数据处理工具:https://confluence.ecmwf.int/display/ECC/GRIB+tools
安装ecodes库:conda install eccodes
参考:https://www.jianshu.com/p/683d3bb5b3cd
python处理GRIB数据常用包:pygrib
http://bbs.06climate.com/forum.php?mod=viewthread&tid=58258
|