MeteoInfo 读集合预报文件问题

gszywfc · 发表于 2016-4-11 16:37:02

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

看似只能读1个成员数据，请问王老师在文件信息中维数排序是对的吧：
File Name: D:\grib\pl_20160328_0.grib1
File type: GRIB1 Collection (GRIB-1)
Dimensions: 5
lon = 181;
lat = 111;
time = 1;
isobaric = 11;
ens = 51;
X Dimension: Xmin = 0.0; Xmax = 180.0; Xsize = 181; Xdelta = 1.0
Y Dimension: Ymin = -20.0; Ymax = 90.0; Ysize = 111; Ydelta = 1.0

MeteoInfo · 发表于 2016-4-11 16:39:45

没看明白你想问什么？

gszywfc · 发表于 2016-4-11 16:39:59

是这样吗？data(ens ，isobaric，time，lat，lon)

MeteoInfo · 发表于 2016-4-11 16:42:51

gszywfc 发表于 2016-4-11 16:39
是这样吗？data(ens ，isobaric，time，lat，lon)

你1楼的信息没给全吧，看看下面的变量中维的设置，不同的变量可能会不同。

gszywfc · 发表于 2016-4-11 16:59:40

File Name: D:\grib\pl_20160328_0.grib1
File type: GRIB1 Collection (GRIB-1)
Dimensions: 5
lon = 181;
lat = 111;
time = 1;
isobaric = 11;
ens = 51;
X Dimension: Xmin = 0.0; Xmax = 180.0; Xsize = 181; Xdelta = 1.0
Y Dimension: Ymin = -20.0; Ymax = 90.0; Ysize = 111; Ydelta = 1.0
Global Attributes:
: Originating_or_generating_Center = "European Centre for Medium Range Weather Forecasts (ECMWF) (RSMC)"
: Originating_or_generating_Subcenter = "0"
: GRIB_table_version = "0,128"
: file_format = "GRIB-1"
: Conventions = "CF-1.6"
: history = "Read using CDM IOSP GribCollection v3"
: featureType = "GRID"
Variations: 16
int LatLon_Projection);
LatLon_Projection: grid_mapping_name = "latitude_longitude"
LatLon_Projection: earth_radius = 6367470.0
float lat(lat);
lat: units = "degrees_north"
float lon(lon);
lon: units = "degrees_east"
double reftime);
reftime: units = "Hour since 2016-03-28T12:00:00Z"
reftime: standard_name = "forecast_reference_time"
reftime: long_name = "GRIB reference time"
reftime: calendar = "proleptic_gregorian"
double time(time);
time: units = "Hour since 2016-03-28T12:00:00Z"
time: standard_name = "time"
time: long_name = "GRIB forecast or observation time"
time: calendar = "proleptic_gregorian"
float isobaric(isobaric);
isobaric: units = "hPa"
isobaric: long_name = "Isobaric surface"
isobaric: positive = "down"
isobaric: Grib_level_type = 100
int ens(ens);
ens: _CoordinateAxisType = "Ensemble"
float Temperature_isobaric(time,ens,isobaric,lat,lon);
Temperature_isobaric: long_name = "Temperature @ Isobaric surface"
Temperature_isobaric: units = "K"
Temperature_isobaric: missing_value = NaNf
Temperature_isobaric: grid_mapping = "LatLon_Projection"
Temperature_isobaric: coordinates = "reftime time ens isobaric lat lon "
Temperature_isobaric: Grib_Variable_Id = "VAR_98-0-128-130_L100"
Temperature_isobaric: Grib1_Center = 98
Temperature_isobaric: Grib1_Subcenter = 0
Temperature_isobaric: Grib1_TableVersion = 128
Temperature_isobaric: Grib1_Parameter = 130
Temperature_isobaric: Grib1_Parameter_Name = "t"
Temperature_isobaric: Grib1_Level_Type = 100
Temperature_isobaric: Grib1_Level_Desc = "Isobaric surface"
float U_component_of_wind_isobaric(time,ens,isobaric,lat,lon);
U_component_of_wind_isobaric: long_name = "U component of wind @ Isobaric surface"
U_component_of_wind_isobaric: units = "m s**-1"
U_component_of_wind_isobaric: missing_value = NaNf
U_component_of_wind_isobaric: grid_mapping = "LatLon_Projection"
U_component_of_wind_isobaric: coordinates = "reftime time ens isobaric lat lon "
U_component_of_wind_isobaric: Grib_Variable_Id = "VAR_98-0-128-131_L100"
U_component_of_wind_isobaric: Grib1_Center = 98
U_component_of_wind_isobaric: Grib1_Subcenter = 0
U_component_of_wind_isobaric: Grib1_TableVersion = 128
U_component_of_wind_isobaric: Grib1_Parameter = 131
U_component_of_wind_isobaric: Grib1_Parameter_Name = "u"
U_component_of_wind_isobaric: Grib1_Level_Type = 100
U_component_of_wind_isobaric: Grib1_Level_Desc = "Isobaric surface"
float V_component_of_wind_isobaric(time,ens,isobaric,lat,lon);
V_component_of_wind_isobaric: long_name = "V component of wind @ Isobaric surface"
V_component_of_wind_isobaric: units = "m s**-1"
V_component_of_wind_isobaric: missing_value = NaNf
V_component_of_wind_isobaric: grid_mapping = "LatLon_Projection"
V_component_of_wind_isobaric: coordinates = "reftime time ens isobaric lat lon "
V_component_of_wind_isobaric: Grib_Variable_Id = "VAR_98-0-128-132_L100"
V_component_of_wind_isobaric: Grib1_Center = 98
V_component_of_wind_isobaric: Grib1_Subcenter = 0
V_component_of_wind_isobaric: Grib1_TableVersion = 128
V_component_of_wind_isobaric: Grib1_Parameter = 132
V_component_of_wind_isobaric: Grib1_Parameter_Name = "v"
V_component_of_wind_isobaric: Grib1_Level_Type = 100
V_component_of_wind_isobaric: Grib1_Level_Desc = "Isobaric surface"
float Specific_humidity_isobaric(time,ens,isobaric,lat,lon);
Specific_humidity_isobaric: long_name = "Specific humidity @ Isobaric surface"
Specific_humidity_isobaric: units = "kg kg**-1"
Specific_humidity_isobaric: missing_value = NaNf
Specific_humidity_isobaric: grid_mapping = "LatLon_Projection"
Specific_humidity_isobaric: coordinates = "reftime time ens isobaric lat lon "
Specific_humidity_isobaric: Grib_Variable_Id = "VAR_98-0-128-133_L100"
Specific_humidity_isobaric: Grib1_Center = 98
Specific_humidity_isobaric: Grib1_Subcenter = 0
Specific_humidity_isobaric: Grib1_TableVersion = 128
Specific_humidity_isobaric: Grib1_Parameter = 133
Specific_humidity_isobaric: Grib1_Parameter_Name = "q"
Specific_humidity_isobaric: Grib1_Level_Type = 100
Specific_humidity_isobaric: Grib1_Level_Desc = "Isobaric surface"
float Vertical_velocity_isobaric(time,ens,isobaric,lat,lon);
Vertical_velocity_isobaric: long_name = "Vertical velocity @ Isobaric surface"
Vertical_velocity_isobaric: units = "Pa s**-1"
Vertical_velocity_isobaric: missing_value = NaNf
Vertical_velocity_isobaric: grid_mapping = "LatLon_Projection"
Vertical_velocity_isobaric: coordinates = "reftime time ens isobaric lat lon "
Vertical_velocity_isobaric: Grib_Variable_Id = "VAR_98-0-128-135_L100"
Vertical_velocity_isobaric: Grib1_Center = 98
Vertical_velocity_isobaric: Grib1_Subcenter = 0
Vertical_velocity_isobaric: Grib1_TableVersion = 128
Vertical_velocity_isobaric: Grib1_Parameter = 135
Vertical_velocity_isobaric: Grib1_Parameter_Name = "w"
Vertical_velocity_isobaric: Grib1_Level_Type = 100
Vertical_velocity_isobaric: Grib1_Level_Desc = "Isobaric surface"
float Vorticity_relative_isobaric(time,ens,isobaric,lat,lon);
Vorticity_relative_isobaric: long_name = "Vorticity relative @ Isobaric surface"
Vorticity_relative_isobaric: units = "s**-1"
Vorticity_relative_isobaric: missing_value = NaNf
Vorticity_relative_isobaric: grid_mapping = "LatLon_Projection"
Vorticity_relative_isobaric: coordinates = "reftime time ens isobaric lat lon "
Vorticity_relative_isobaric: Grib_Variable_Id = "VAR_98-0-128-138_L100"
Vorticity_relative_isobaric: Grib1_Center = 98
Vorticity_relative_isobaric: Grib1_Subcenter = 0
Vorticity_relative_isobaric: Grib1_TableVersion = 128
Vorticity_relative_isobaric: Grib1_Parameter = 138
Vorticity_relative_isobaric: Grib1_Parameter_Name = "vo"
Vorticity_relative_isobaric: Grib1_Level_Type = 100
Vorticity_relative_isobaric: Grib1_Level_Desc = "Isobaric surface"
float Divergence_isobaric(time,ens,isobaric,lat,lon);
Divergence_isobaric: long_name = "Divergence @ Isobaric surface"
Divergence_isobaric: units = "s**-1"
Divergence_isobaric: missing_value = NaNf
Divergence_isobaric: grid_mapping = "LatLon_Projection"
Divergence_isobaric: coordinates = "reftime time ens isobaric lat lon "
Divergence_isobaric: Grib_Variable_Id = "VAR_98-0-128-155_L100"
Divergence_isobaric: Grib1_Center = 98
Divergence_isobaric: Grib1_Subcenter = 0
Divergence_isobaric: Grib1_TableVersion = 128
Divergence_isobaric: Grib1_Parameter = 155
Divergence_isobaric: Grib1_Parameter_Name = "d"
Divergence_isobaric: Grib1_Level_Type = 100
Divergence_isobaric: Grib1_Level_Desc = "Isobaric surface"
float Geopotential_Height_isobaric(time,ens,isobaric,lat,lon);
Geopotential_Height_isobaric: long_name = "Geopotential Height @ Isobaric surface"
Geopotential_Height_isobaric: units = "gpm"
Geopotential_Height_isobaric: missing_value = NaNf
Geopotential_Height_isobaric: grid_mapping = "LatLon_Projection"
Geopotential_Height_isobaric: coordinates = "reftime time ens isobaric lat lon "
Geopotential_Height_isobaric: Grib_Variable_Id = "VAR_98-0-128-156_L100"
Geopotential_Height_isobaric: Grib1_Center = 98
Geopotential_Height_isobaric: Grib1_Subcenter = 0
Geopotential_Height_isobaric: Grib1_TableVersion = 128
Geopotential_Height_isobaric: Grib1_Parameter = 156
Geopotential_Height_isobaric: Grib1_Parameter_Name = "gh"
Geopotential_Height_isobaric: Grib1_Level_Type = 100
Geopotential_Height_isobaric: Grib1_Level_Desc = "Isobaric surface"
float Relative_humidity_isobaric(time,ens,isobaric,lat,lon);
Relative_humidity_isobaric: long_name = "Relative humidity @ Isobaric surface"
Relative_humidity_isobaric: units = "%"
Relative_humidity_isobaric: missing_value = NaNf
Relative_humidity_isobaric: grid_mapping = "LatLon_Projection"
Relative_humidity_isobaric: coordinates = "reftime time ens isobaric lat lon "
Relative_humidity_isobaric: Grib_Variable_Id = "VAR_98-0-128-157_L100"
Relative_humidity_isobaric: Grib1_Center = 98
Relative_humidity_isobaric: Grib1_Subcenter = 0
Relative_humidity_isobaric: Grib1_TableVersion = 128
Relative_humidity_isobaric: Grib1_Parameter = 157
Relative_humidity_isobaric: Grib1_Parameter_Name = "r"
Relative_humidity_isobaric: Grib1_Level_Type = 100
Relative_humidity_isobaric: Grib1_Level_Desc = "Isobaric surface"

gszywfc · 发表于 2016-4-11 17:00:46

比如说变量温度的维数顺序

MeteoInfo · 发表于 2016-4-11 17:09:39

gszywfc 发表于 2016-4-11 17:00
比如说变量温度的维数顺序

float Temperature_isobaric(time,ens,isobaric,lat,lon);
不是很清楚吗？

gszywfc · 发表于 2016-4-11 17:09:45

float Temperature_isobaric(time,ens,isobaric,lat,lon);
这个应该是正确的吧？

gszywfc · 发表于 2016-4-11 17:14:23

明白了，多谢了

MeteoInfo · 发表于 2016-4-11 17:14:36

gszywfc 发表于 2016-4-11 17:09
float Temperature_isobaric(time,ens,isobaric,lat,lon);
这个应该是正确的吧？

为什么会怀疑其正确性呢？

		自动登录	找回密码
密码			立即注册