使用GRADS读取转换nc文件的方法

下面给出使用Grads读取nc格式文件的一个范例脚本文件（.gs文件）：
以位势高度场为例
‘reinit’
‘sdfopen d:\xxxx\hgt.mon.mean.nc ‘ （打开需要的nc文件，注意路径）
‘set gxout fwrite’
‘set fwrite d:\xxxxxx\hgt.bin’ （指定要输出的二进制文件名和路径）
t=a （a为起始时次）
while（t<=b）（开始循环，b为结束时次）
‘set t’ t （设定时次）
‘set lev 1000’ （设定要输出的层次）
‘d hgt ‘ （输出）
‘set lev 925’ （设定要输出的层次）
‘d hgt ‘ （输出）
‘set lev 850’ （设定要输出的层次）
‘set lon 0 357.5’
‘set lat -90 90 ‘ （设定写出区域范围，这是全球范围，第一句是必需的，如果是用默认，则会在x方向输出为145个点，要注意）
‘set lev 1000’ （设定要输出的层次）
‘set t 1 640’ （设定时间范围）
‘d hgt’ （输出）
‘reinit’
这样就生成了需要的二进制文件，其对应的ctl文件需要另外给出。
以上是输出单层的二进制文件，把多层的数据放在同一个二进制文件中，就要使用循环，根据grads规定的数据存放次序（依次为经度、纬度、高度、变量、时次），如下：
‘reinit’
‘sdfopen d:\xxxx\hgt.mon.mean.nc’ （打开需要的nc文件，注意路径）
‘set gxout fwrite’
‘set fwrite d:\xxxxxx\hgt.bin’ （指定要输出的二进制文件名和路径）
‘set lon 0 357.5’
‘set lat -90 90’ （设定写出区域范围，这是全球范围，第一句是必需的，如果用默认set lon 0 360，则会在x方向输出为145个点，要注意）
‘d hgt ‘ （输出）
‘set lev 700’ （设定要输出的层次）
‘d hgt ‘ （输出）
‘set lev 600’ （设定要输出的层次）
‘d hgt ‘ （输出）
‘set lev 500’ （设定要输出的层次）
‘d hgt ‘ （输出）
（可以添加任意的层次）
（由于高度分布不连续而且数量有限，所以用这种方法）
t=t 1
endwhile （结束循环）
‘