当我将 Grib2 转换为 Netcdf4 时出现问题
Problems when i convert Grib2 to Netcdf4
当我尝试从 CFSR 转换一些 grib 文件时遇到了一些问题。
当我在 unix shell 中使用 ncl_convert2nc
时,nc 文件信息是这样的:
netcdf flxf01.gdas.197901 {
dimensions:
ygrid_0 = UNLIMITED ; // (0 currently)
lv_DBLL0 = 4 ;
variables:
// global attributes:
:creation_date = "Sun Nov 12 11:29:36 AEDT 2017" ;
:NCL_Version = "6.3.0" ;
:system = "Linux charcoalp-Lenovo-Erazer-Y50-70 4.10.0-38-generic #42-Ubuntu SMP Tue Oct 10 13:24:27 UTC 2017 x86_64 x86_64 x86_64 GNU/Linux" ;
:Conventions = "None" ;
:grib_source = "flxf01.gdas.197901.grb2" ;
:title = "NCL: convert-GRIB-to-netCDF" ;
}
这只会在我转换 CFSR 数据 grib2 文件时发生。
原始文件部分信息如下:
Nio file: flxf01.gdas.197901.grb2
global attributes:
dimensions:
lat_0 = 576
lon_0 = 1152
lv_DBLL0 = 4
variables:
float TMP_P8_L1_GGA0_avg744h [ lat_0, lon_0 ]
center : US National Weather Service - NCEP (WMC)
production_status : Operational products
long_name : Temperature
units : K
_FillValue : 1e+20
grid_type : Gaussian latitude/longitude
parameter_discipline_and_category : Meteorological products, Temperature
parameter_template_discipline_category_number : [8, 0, 0, 0]
level_type : Ground or water surface
level : 0
type_of_statistical_processing : Average
statistical_process_duration : 744 hours (ending at forecast time)
forecast_time : 745
forecast_time_units : hours
initial_time : 01/01/1979 (00:00)
float TMP_P8_L103_GGA0_avg744h [ lat_0, lon_0 ]
center : US National Weather Service - NCEP (WMC)
production_status : Operational products
long_name : Temperature
units : K
_FillValue : 1e+20
grid_type : Gaussian latitude/longitude
parameter_discipline_and_category : Meteorological products, Temperature
parameter_template_discipline_category_number : [8, 0, 0, 0]
level_type : Specified height level above ground (m)
level : 2
type_of_statistical_processing : Average
statistical_process_duration : 744 hours (ending at forecast time)
forecast_time : 745
forecast_time_units : hours
initial_time : 01/01/1979 (00:00)
如您所见,它能够从 grib 文件中找到变量。那么为什么nc文件是空的呢? ncl_convert2nc
为 JRA-55 再分析数据集的 Grib 文件工作。谁能帮我看看问题出在哪里?
出于兴趣,您是否尝试过使用 CDO?我认为命令是:
cdo -f nc4 copy in.grb out.nc
当我尝试从 CFSR 转换一些 grib 文件时遇到了一些问题。
当我在 unix shell 中使用 ncl_convert2nc
时,nc 文件信息是这样的:
netcdf flxf01.gdas.197901 {
dimensions:
ygrid_0 = UNLIMITED ; // (0 currently)
lv_DBLL0 = 4 ;
variables:
// global attributes:
:creation_date = "Sun Nov 12 11:29:36 AEDT 2017" ;
:NCL_Version = "6.3.0" ;
:system = "Linux charcoalp-Lenovo-Erazer-Y50-70 4.10.0-38-generic #42-Ubuntu SMP Tue Oct 10 13:24:27 UTC 2017 x86_64 x86_64 x86_64 GNU/Linux" ;
:Conventions = "None" ;
:grib_source = "flxf01.gdas.197901.grb2" ;
:title = "NCL: convert-GRIB-to-netCDF" ;
}
这只会在我转换 CFSR 数据 grib2 文件时发生。
原始文件部分信息如下:
Nio file: flxf01.gdas.197901.grb2
global attributes:
dimensions:
lat_0 = 576
lon_0 = 1152
lv_DBLL0 = 4
variables:
float TMP_P8_L1_GGA0_avg744h [ lat_0, lon_0 ]
center : US National Weather Service - NCEP (WMC)
production_status : Operational products
long_name : Temperature
units : K
_FillValue : 1e+20
grid_type : Gaussian latitude/longitude
parameter_discipline_and_category : Meteorological products, Temperature
parameter_template_discipline_category_number : [8, 0, 0, 0]
level_type : Ground or water surface
level : 0
type_of_statistical_processing : Average
statistical_process_duration : 744 hours (ending at forecast time)
forecast_time : 745
forecast_time_units : hours
initial_time : 01/01/1979 (00:00)
float TMP_P8_L103_GGA0_avg744h [ lat_0, lon_0 ]
center : US National Weather Service - NCEP (WMC)
production_status : Operational products
long_name : Temperature
units : K
_FillValue : 1e+20
grid_type : Gaussian latitude/longitude
parameter_discipline_and_category : Meteorological products, Temperature
parameter_template_discipline_category_number : [8, 0, 0, 0]
level_type : Specified height level above ground (m)
level : 2
type_of_statistical_processing : Average
statistical_process_duration : 744 hours (ending at forecast time)
forecast_time : 745
forecast_time_units : hours
initial_time : 01/01/1979 (00:00)
如您所见,它能够从 grib 文件中找到变量。那么为什么nc文件是空的呢? ncl_convert2nc
为 JRA-55 再分析数据集的 Grib 文件工作。谁能帮我看看问题出在哪里?
出于兴趣,您是否尝试过使用 CDO?我认为命令是:
cdo -f nc4 copy in.grb out.nc