matlab读取GFS的grib2文件的示例

又是那隻貓

关于用matlab一站式读取NetCDF, OPeNDAP, HDF5, GRIB, GRIB2, HDF4等等
http://bbs.06climate.com/forum.php?mod=viewthread&tid=47599
其实这里是补一下上次欠下的一个示例；

从型如GFS_Global_0p5deg_20140531_0000_anl.grib2的文件中读取温度等变量
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

grib=ncgeodataset(Data_File);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%time
utc_time=double(grib.data('time'));
Data.UTC=datestr(datenum([Data_File(end-22:end-15) Data_File(end-13:end-10)],'yyyymmddHHSS')+  utc_time/24  ,'dd-mmm-yyyy HH:MM:SS' );
%%%%%%%%%%%%%%%%%%%%%%%%%%%%grid isobaric
Data.isobaric=double(grib.data('isobaric1'))/100;%21 level and convert it into hPa
%%%%%%%%%%%%%%%%%%%%%%%%%%%%grid lat lon
Data.lat=double(grib.data('lat'));
Data.lon=double(grib.data('lon'));



Data.Temperature=squeeze(double(grib.data('Temperature_isobaric')));%K
Data.Temperature=Data.Temperature(1:21,:,:)-273;%convert it into degress Centigrade


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


上面是读取温度的示例，完整的我传成一个文件，收少许贡献。附件里包含了可降水率、相对湿度、风场、温度等
%float Precipitable_water_entire_atmosphere(time=1, lat=361, lon=720);
%float Relative_humidity_isobaric(time=1, isobaric3=25, lat=361, lon=720);
%float Temperature_isobaric(time=1, isobaric=26, lat=361, lon=720);
%float u-component_of_wind_isobaric(time=1, isobaric=26, lat=361, lon=720);
%float v-component_of_wind_isobaric(time=1, isobaric=26, lat=361, lon=720);
%float Vertical_velocity_pressure_isobaric(time=1, isobaric1=21, lat=361, lon=720);
%float Convective_available_potential_energy_surface(time=1, lat=361, lon=720);
%float Convective_inhibition_surface(time=1, lat=361, lon=720);
等变量的读取

Fun_GFS_Data.m (2.26 KB, 下载次数: 259, 售价: 3 贡献)

2014-10-22 21:46 上传

点击文件名下载附件

流行雨

请教楼主，我现在想解的数据名称是gfs.t12z.pgrb2.1p00.f006，
NCEP在年初才改的名称和格式，最近数据变更了，
我读了一下，提示错误：
Warning: Invalid escape sequence appears in format string. See help sprintf for valid escape sequences.
> In ncdataset>ncdataset.ncdataset at 89
  In cfdataset>cfdataset.cfdataset at 59
  In ncgeodataset>ncgeodataset.ncgeodataset at 74
??? Error using ==> ncdataset>ncdataset.ncdataset at 91
Failed to open D:
ead_grib

Error in ==> cfdataset>cfdataset.cfdataset at 59
            obj = obj@ncdataset(url);

Error in ==> ncgeodataset>ncgeodataset.ncgeodataset at 74
            obj = obj@cfdataset(url);

Caused by:
    Undefined variable "ucar" or class "ucar.nc2.dataset.NetcdfDataset.openDataset".

ningking1986

多谢分享啊

浪子

谢谢分享 贡献不足

afei790912

看来有钱就是好呀，

maccha

有用呀~~希望楼主多发

xull2010

雪中送炭啊！

diegu

谢谢分享 贡献不足

chinalys

感谢分享，谢谢楼主！

gszywfc

hao...........