fwrite提取出来fnl的高度场画的图跟原始图不一样- -，求救！

Mr.Yes

大神们，本人大三，在做北京721特大降水的实习，资料idx格式的，想提取700hpa的高度场资料做barnes滤波，但是提取出来的图像跟直接用fnl画的图就是不一样- -一个多星期了，就是不知道问题出现在那里，左侧为直接用idx画的图，右侧为提取出来的数据画的。
附程序：
1.提取资料的gs文件
'reinit'
'open D:\Barnes\fnl_20120719_00_00.ctl'
'set mpdset cnworld cnriver'
'set lon  0 360.000000'
'set lat -90.00000 90.000000'
'set fwrite d:\barnes\out\hgtprs.grd'
'set gxout fwrite'
lev.1=1000;lev.2=975;lev.3=950;lev.4=925;lev.5=900;lev.6=850;lev.7=800;lev.8=750;lev.9=700;lev.10=650;lev.11=600;
lev.12=550;lev.13=500;lev.14=450;lev.15=400;lev.16=350;lev.17=300;lev.18=250;lev.19=200;lev.20=150;lev.21=100;
i=1
while(i<=21)
lev=lev.i
'set lev 'lev''
'd hgtprs'
i=i+1
endwhile
'disable fwrite'
'c'
;

2.针对提取出来的数据的ctl文件
dset d:/barnes/out/hgtprs.grd
undef 9.999E+20
title 20120719_00_00 hgt700

ydef 181 linear -90.000000 1.00000
xdef 360 linear 0.000000 1.000000
zdef 21 levels 1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100
tdef 1 linear 00Z19jul2012 1mo
vars 1
HGTprs    26,100  0,3,5 ** Geopotential Height [gpm]

endvars


3.用idx文件直接画图的gs，画的是700hpa的图
'reinit'
'open D:\Barnes\fnl_20120719_00_00.ctl'
'enable print d:/grads/data/700.gmf'
'set lat -90 90.000000'
'set lon  0 360.000000'
'set lev 700'
'd HGTprs'
'print'
'disable print'
;


4.用提取出来的资料画图的gs
'reinit'
'open D:\Barnes\program\hgt700.ctl'
'enable print d:/grads/data/700out.gmf'
'set lat -90.00000 90.000000'
'set lon  0 360.000000'
lev.1=1000;lev.2=975;lev.3=950;lev.4=925;lev.5=900;lev.6=850;lev.7=800;lev.8=750;lev.9=700;lev.10=650;lev.11=600;
lev.12=550;lev.13=500;lev.14=450;lev.15=400;lev.16=350;lev.17=300;lev.18=250;lev.19=200;lev.20=150;lev.21=100;
i=1
while(i<=21)
lev=lev.i
'set lev 'lev''
'd hgtprs'
'print'
'c'
i=i+1
endwhile
'disable print'
;



恳请各位大神解答！！！

沙颖凯

目测经度范围定义错了，把set lon  0 360.000000改成-180 180试试？

Mr.Yes

沙颖凯 发表于 2014-5-9 19:48
目测经度范围定义错了，把set lon  0 360.000000改成-180 180试试？

好，我是看idx文件的ctl文件里面这么写的，那中午回去试试看！You are so sweet!

Mr.Yes

沙颖凯 发表于 2014-5-9 19:48
目测经度范围定义错了，把set lon  0 360.000000改成-180 180试试？

还是不对- -   感觉就是经度错位了一样

Mr.Yes

沙颖凯 发表于 2014-5-9 19:48
目测经度范围定义错了，把set lon  0 360.000000改成-180 180试试？

还是不对- -   感觉就是经度错位了一样   改完以后错位的更离谱了

栤唥de嶶笑

描述性文件写错了吧？HGTprs&#160; &#160; 26,100&#160;&#160;0,3,5 ** Geopotential Height [gpm]   不是只有21层么？你怎么写26层？？？

green

请问楼主问题解决了没？

Mr.Yes

green 发表于 2014-5-12 18:25
请问楼主问题解决了没？

没解决呢= =

Mr.Yes

green 发表于 2014-5-12 18:25
请问楼主问题解决了没？

没解决呢= =

Mr.Yes

栤唥de嶶笑 发表于 2014-5-10 15:10
描述性文件写错了吧？HGTprs&#160; &#160; 26,100&#160;&#160;0,3,5 ** Geopotential Height [gpm]   不是 ...

21也不对- -