FNL_1x1资料认识与应用(GrADS)- 兰溪整合版(grib1 grib2)

FNL 1X1 资料认识与应用

兰溪之水整合版2011-12-16 14:34:05（第一版）

2012-05-01 22:42:49（第二版）FNL 1.0X1.0数据下载地址：/datasets/ds083.2/需要注册通过方可下载。

FNL 1.0X1.0 资料包含的物理量及其缩写

缩写参数名称

No4LFTXsfc 近地表四层等压面的抬升指数

No5WAVAprs 500 hPa等压面位势高度距平

No5WAVHprs 500 hPa等压面位势高度

ABSVprs 绝对涡度

CAPE 对流有效位能

CIN 对流抑制能

CLWMRprs 云水

CWATclm 气柱云水

GPAprs 位势高度距平

HGT 位势高度

HPBLsfc 地表行星边界层高度

ICECsfc 海冰密集度

LANDsfc 陆地覆盖

LFTXsfc 地表抬升指数

O3MRprs 臭氧层混合比

POTsig995 位温

PRE 气压

PWATclm 可降水量

RH 相对湿度

SOILW 土壤体积含水量

SPFH 比湿

TCDCcvl 对流云总云量

TM 温度

TOZNEclm 臭氧含量

UGRD u分量

VGRD v 分量

VVEL 垂直速度

VWSH 垂直风切变

WEASDsfc 累积雪量

GrADS处理FNL 1.0X1.0（grib1）数据

处理前需要先将grib2ctl.exe放到GrADS安装文件夹GrADS19\win32（1.9版本）或OpenGrADS\Contents\Cygwin\Versions\2.0.a9.oga.1\i686（2.0版本）下，方便操作。gribmap.exe（这个GrADS自带的）

第一步：先生成一个描述文件ctl

打开命令提示符，

C:\Documents and Settings\Administrator>grib2ctl grib_file > grib_file.ctl

如：

或者进入GrADS：

gs>!grib2ctl grib_file > grib_file.ctl （！表示调用外部的程序，注意路径用“/”）

如：

这样fnl_20101013_00_00_c.ctl描述文件就生成了！

第二步：利用GrADS自带的gribmap.exe生成索引文件：

C:\Documents and Settings\Administrator>gribmap –v –i grib_file.ctl > b.txt

如：

或进入GrADS：

这样就会在H:\test路径下生成名为fnl_20101013_00_00_c.idx的索引文件，并会把整个映射过程写入到b.txt文件中，当然这里也可省略查看映射过程，即

C:\Documents and Settings\Administrator>gribmap –i grib_file.ctl

这样就可以开始用GrADS读取FNL文件画图了！

不过我们处理FNL资料的时候一般都是不止一个文件，所以我们就需要用到批处理了！

第一步：利用grib2ctl.exe生成初始时刻ncep数据的ctl文件；示例中生成的是fnl_20101013_00_00_c数据的文件。

第二步：将ctl中相应的信息进行修改：

a、将dset的数据集文件名写成代换模块格式；如：dset H:/test/fnl_201010%d2_%h2_00_c

％y2 代表两位数年

％y4 代表四位数年

％m1 代表一位或者两位数的月

％m2 代表两位数月（用0补齐1位数）

％mc 3个字符月份的缩写

％d1 1或2位天

％d2 两位天

％h1 1或者2位时

％h2 2位时

b、在第一行后添加一行“options template”定义时间范围和增量；

c、tdef number LINEAR start increment，将其中的number和increment（ncep时间间隔为6小时，可写成6hr）设成与所拥有的ncep资料一致；（修改后的信息如图所示,里面的路径需自己设置）；

第三步：利用gribmap.exe生成已改过后的ctl相对应的idx文件；

这样后，我们就可以用GrADS打开多个fnl文件进行数据处理。

############################################################################### ### 一种更简单的方法### ############################################################################## 第一步：

第二步：

FNL 1*1 资料认识与应用兰溪之水整合版

############################################################################### ############################################################################### ### 还有一种方法 ### ############################################################################## 就是把多个文件合并成一个，然后再使用grib2ctl.exe 和gribmap.exe生成对应的ctl文件和idx文件：

Windows下一般使用copy命令，linux下一般使用cat命令,

copy /b file1+file2 newfile, 文件比较多可以使用通配符：copy /b fnl* newfile

或

cat file1 file2 > newfile, 使用通配符：cat fnl* > newfile

然后对合并后的newfile用grib2ctl.exe 和gribmap.exe生成对应的ctl文件和idx文件，这个上面的重复，就不详细说明了！

附：如果我们只需要fnl文件中的某些变量，我们可以编写gs文件提取所需要时刻的ncep 数据了（注意：提取数据要很注意数据的存储顺序，时间-变量-层次，水平维设置建议使用’set x …’, ‘set y …’以便识别格点数和准确经纬度）。

gs如下：

'reinit'

'open H:/test/fnl_201010Megi.ctl' 'set gxout fwrite'

'set fwrite H:/test/ Megi.dat'

'set lat 15 40'

'set lon 110 145'

tt=1

while(tt<=60)

'set t ' tt

zz=1

while(zz<=26)

'set z ' zz

'd UGRDprs'

zz=zz+1 endwhile

zz=1

while(zz<=26)

'set z ' zz

'd VGRDprs'

zz=zz+1 endwhile

tt=tt+1 endwhile

'disable fwrite' 'reinit'

;