Step1:利用grib2ctl.exe生成初始时刻ncep数据的ctl文件；示例中生成的是fnl_20090801_00_00_c数据的文件。 Step2:将ctl中相应的信息进行修改：
a、将dset的数据集文件名写成代换模块格式；如：dset D:\\fog\\fnl_200212ò_%h2_00_c %y2 ===> 两位数年；%y4 ===> 四位数年；%m1 ===> 1或2位数月；
%m2 ===> 2位数月（用0补齐一位数）；%mc ===> 3字符月份缩写；
? ===> 1或2位天；ò ===> 2位天；%h1 ===> 1或2位时；%h2 ===> 2位时； b、在第一行后添加一行“options template”定义时间范围和增量；
c、tdef number LINEAR start increment，将其中的number和increment（ncep时间间隔为6小时，可写成6hr）设成与所拥有的ncep资料一致；（修改后的信息如图所示,里面的路径需自己设置）；
Step3：利用gribmap.exe生成与改过后的ctl相对应的idx文件； 这样就可以编写gs文件提取所需要时刻的ncep数据了。 gs如下：
'reinit'
'open E:\\Morakot_fnl\\fnl200908Morakot.ctl' 'set lat 15 40' 'set lon 110 145'
'set fwrite E:\\program\\data.dat' 'set gxout fwrite' t=1
while(t<=61) 'set t ' t z=1
while(z<=21) 'set z ' z 'd HGTprs' 'd RHprs' 'd TMPprs' 'd UGRDprs'
'd VGRDprs' 'd VVELprs' z=z+1 endwhile t=t+1 endwhile
'disable fwrite' 'reinit'
(注意：此过程为基于GrADs 2.0.a3.oga.1版本的试验，其它版本是否一致不敢保证)