怎样用Fortran或者NCL读取JRA55的资料

bhy

数据文件如上
数据描述文件（.ctl）如下：
dset ^fcst_phy2m125.%y4%m2
index ^fcst_phy2m125.monthly.idx
undef 9.999E+20
title fcst_phy2m125
*  produced by grib2ctl v0.9.12.5p41
dtype grib 255
options yrev template
ydef 145 linear -90.000000 1.25
xdef 288 linear 0.000000 1.250000
tdef 678 linear 00Z01Jan1958 1mo
zdef 1 linear 1 1
vars 26
CPRATsfc  0 63,1,0  ** surface Convective precipitation [mm/day]
CSDLFsfc  0 163,1,0  ** surface Clear sky downward long wave radiation flux [W/m^2]
CSDSFsfc  0 161,1,0  ** surface Clear sky downward solar radiation flux [W/m^2]
CSULFtoa  0 162,8,0 ** top of atmos Clear sky upward long wave radiation flux [W/m^2]
CSUSFsfc  0 160,1,0  ** surface Clear sky upward solar radiation flux [W/m^2]
CSUSFtoa  0 160,8,0 ** top of atmos Clear sky upward solar radiation flux [W/m^2]
DLWRFsfc  0 205,1,0  ** surface Downward long wave radiation flux [W/m^2]
DSWRFsfc  0 204,1,0  ** surface Downward solar radiation flux [W/m^2]
DSWRFtoa  0 204,8,0 ** top of atmos Downward solar radiation flux [W/m^2]
EVPsfc  0 57,1,0  ** surface Evaporation [mm/day]
FGLUsfc  0 147,1,0  ** surface Zonal momentum flux by long gravity wave [N/m^2]
FGLVsfc  0 148,1,0  ** surface Meridional momentum flux by long gravity wave [N/m^2]
FGSUsfc  0 159,1,0  ** surface Zonal momentum flux by short gravity wave [N/m^2]
FGSVsfc  0 154,1,0  ** surface Meridional momentum flux by short gravity wave [N/m^2]
LHTFLsfc  0 121,1,0  ** surface Latent heat flux [W/m^2]
LPRATsfc  0 62,1,0  ** surface Large scale precipitation [mm/day]
PRESsfc  0 1,1,0  ** surface Pressure [Pa]
SHTFLsfc  0 122,1,0  ** surface Sensible heat flux [W/m^2]
SRWEQsfc  0 64,1,0  ** surface Snowfall rate water equivalent [mm/day]
TPRATsfc  0 61,1,0  ** surface Total precipitation [mm/day]
UFLXsfc  0 124,1,0  ** surface Momentum flux, u-component [N/m^2]
ULWRFsfc  0 212,1,0  ** surface Upward long wave radiation flux [W/m^2]
ULWRFtoa  0 212,8,0 ** top of atmos Upward long wave radiation flux [W/m^2]
USWRFsfc  0 211,1,0  ** surface Upward solar radiation flux [W/m^2]
USWRFtoa  0 211,8,0 ** top of atmos Upward solar radiation flux [W/m^2]
VFLXsfc  0 125,1,0  ** surface Momentum flux, v-component [N/m^2]
ENDVARS
求高人指点

lysx

这是grib数据啊，Fortran不能直接读取的（除非安装点东西，我没有，也不懂）。你可以用wgrib将其解压为一般的二进制数据就 可以用fortran读了

kongfeng0824

yongMATLAB处理会更好

bhy

lysx 发表于 2014-8-30 20:54
这是grib数据啊，Fortran不能直接读取的（除非安装点东西，我没有，也不懂）。你可以用wgrib将其解压为一般 ...

谢谢，能具体说一下怎么解压吗？

lysx

bhy 发表于 2014-9-1 17:13
谢谢，能具体说一下怎么解压吗？

grib2ctl，论坛里有很多介绍！新手上路吧，多浏览下论坛中的技术帖，对你只有好处！

bhy

lysx 发表于 2014-9-2 07:45
grib2ctl，论坛里有很多介绍！新手上路吧，多浏览下论坛中的技术帖，对你只有好处！

找到了，谢了