wgrib提取出的fnl资料中的变量数值的存储顺序是什么样的

thukd35

       GrADS的小白一枚，最近做的一些事情需要将fnl资料中的位势高度、风速等变量的数值提取出来，用wgrib试着将HGT在1000mb的数据提取了出来：wgrib -s F:/2006/fnl_20060101_00_00.grib1 | find ":HGT:1000mb:" | wgrib -i -text F:/2006/fnl_20060101_00_00.grib1 -o out.txt
然后用matlab读了一下out.txt，发现是一个65160×1的向量（即有360×181个位势高度数据），应该能说明提取方法是对的吧。但问题是这个数据是如何按经纬度排列的呢？我用grib2ctl给出了fnl_20060101_00_00.grib1的数据描述文件ct：
dset F:/2006/fnl_200601%d2_%h2_00.grib1
options template
index F:/2006/fnl_20060101_00_00.grib1.idx
undef 9.999E+20
title F:/2006/fnl_20060101_00_00.grib1
*  produced by grib2ctl v0.9.12.5p33k
dtype grib 3
options yrev
ydef 181 linear -90.000000 1
xdef 360 linear 0.000000 1.000000
tdef 5 linear 00Z01jan2006 6hr
zdef 26 levels
1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 70 50 30 20 10
vars 85

        能根据ctl文件判断我提取出的HGT数据的排序方式吗？还是用wgrib提出变量的数值的排序方式是固定的，如果是固定的，是什么顺序呢？求各位大神指教哈

schliezer

你好，你解决了这个问题吗？我也碰上了

一枝梅

http://blog.sina.com.cn/s/blog_86cf6ab60101qujf.html
2．建立描述文件
解码后的数据文件要配上 ctl 文件才能使用，这要手工建立。注意 ctl 中要有 OPTIONS
yrev，表示 y轴数据是反向的，因为 Grads默认 y轴（纬度）从南到北，而解码出的数据（纬
度）排列是从北到南。前面生成的数据文件对应的 ctl 内容应为：