FNL 1X1 资料认识与应用2

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

'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'
;

然后编写描述文件：Megi.ctl

 

 

 

 

还可以写成批处理.bat文件，如： @echo off rem 此程序用于处理NECP FNL资料，并生成.ctl和.idx文件，然后使用grads绘图 rem 此程序使用的为Grads2.0 echo 正在生成.ctl文件,Please wait ... rem ---------------------------------------------------------------------------

pushd C:\

cd C:\OpenGrADS\Contents\Cygwin\Versions\2.0.a9.oga.1\i686\

grib2ctl H:\test\fnl_20101013_00_00_c > H:\test\fnl_20101013_00_00_c.ctl

rem ---------------------------------------------------------------------------

echo 正在生成.idx索引文件，请稍后...

pushd C:\

cd C:\OpenGrADS\Contents\Cygwin\Versions\2.0.a9.oga.1\i686\ gribmap -v -i H:\test\fnl_20101013_00_00_c.ctl

rem ---------------------------------------------------------------------------

pushd H:\test\

if exist H:\test\fnl_500hpa.gs (del H:\test\fnl_500hpa.gs) else (echo sorry!there is no file)

rem ---------------------------------------------------------------------------

rem 下面开始编写GS文件...

set f=H:\test\

set n=fnl_20101013_00_00_c

set l=.ctl

echo 'reinit' >> fnl_500hpa.gs

echo 'open %f%%n%%l% ' >>H:\test\fnl_500hpa.gs

echo 'set grads off' >>H:\test\fnl_500hpa.gs

echo 'set grid on' >>H:\test\fnl_500hpa.gs

echo 'set gxout contour' >>H:\test\fnl_500hpa.gs

echo *自定义颜色>> fnl_500hpa.gs

echo 'set rgb 16 51 255 0'>> fnl_500hpa.gs

echo 'set rgb 17 255 153 0'>> fnl_500hpa.gs

echo 'set rgb 18 255 0 0'>> fnl_500hpa.gs

echo 'set rgb 19 255 102 153'>> fnl_500hpa.gs

echo 'set rgb 20 102 0 102'>> fnl_500hpa.gs

echo 'set rgb 21 51 0 51'>> fnl_500hpa.gs

echo **设置绘图区域>> fnl_500hpa.gs

echo 'set lat 10 70'>> fnl_500hpa.gs

echo 'set lon 20 150'>> fnl_500hpa.gs

echo 'set lev 500'>> fnl_500hpa.gs

echo *设置等值线>> fnl_500hpa.gs

echo *echo 'setcterp on'>> fnl_500hpa.gs

echo *echo 'setcsmooth on'>> fnl_500hpa.gs

echo 'set clevs 520 524 528 532 536 540 544 548 552 556 560 564 568 572 576 580 584 588 592'>> fnl_500hpa.gs

echo *'set ccolor 18'>> fnl_500hpa.gs

echo 'set cmin 4'>> fnl_500hpa.gs

echo 'd HGTprs/10'>> fnl_500hpa.gs

echo 'set gxout barb'>> fnl_500hpa.gs

echo 'd skip(UGRdprs,3);skip(VGRdprs,3)'>> fnl_500hpa.gs

echo *设置地图>> fnl_500hpa.gs

echo 'set mpdset cnriver' >> fnl_500hpa.gs

echo 'set map 4 1 15' >> fnl_500hpa.gs

echo 'draw map' >> fnl_500hpa.gs

echo 'set mpdset cnworld' >> fnl_500hpa.gs

echo 'set map 15 1 5' >> fnl_500hpa.gs

echo 'draw map' >> fnl_500hpa.gs

echo 'set mpdset xinjiang' >> fnl_500hpa.gs

echo 'set map 14 1 7' >> fnl_500hpa.gs

echo 'draw map' >> fnl_500hpa.gs
echo 'set mproj scaled' >> fnl_500hpa.gs

echo 'set mpdset cnworld' >> fnl_500hpa.gs

echo *'draw title 2007/07/18 02:00(BJT) 500hPa HGT' >> fnl_500hpa.gs

echo 'printim %f%%n%UTC_hgt_wind.png white x1000 y800' >> fnl_500hpa.gs

echo ; >> fnl_500hpa.gs

rem ---------------------------------------------------------------------------

echo 正在绘图，马上完成，请稍等....

pushd C:\

cd C:\OpenGrADS\Contents\Cygwin\Versions\2.0.a9.oga.1\i686\

rem cd C:\OpenGrADS\Contents\Cygwin\Versions\2.0.a9.oga.1\i686\grads -cl

H:\test\fnl_500hpa_hgt_wind.gs

C:\OpenGrADS\Contents\Cygwin\Versions\2.0.a9.oga.1\i686\grads -cl

H:\test\fnl_500hpa.gs

将以上复制到记事本，保存为FNL2Graph4GrADS.bat，双击运行即可！
GrADS