继承上次的ctl文件不能打开的问题之后，现在出现了新的问题~~~

wk1021

继承上次的ctl文件不能打开的问题之后，现在出现了新的问题~~~gs文件编写如下，照着论坛里的台风路径的gs自己修改的一个程序：
'open e:\2.ctl' ;
i=1
while (i<13)
'set t ' i
t=subwrd(result,4)
tc=substr(t,8,5)
'd sp'
'define locx=minloc(min(sp,lat=15,lat=35),lon=110,lon=135)'
'define locy=minloc(min(sp,lon=110,lon=135),lat=15,lat=35)'
'd locx'
locxx=subwrd(result,4)
'd locy'
locyy=subwrd(result,4)
'q gr2w 'locxx' 'locyy ;
x=subwrd(result,3)
y=subwrd(result,6)
'set x 'locxx
'set y 'locyy
minslp=subwrd(result,4)
rc=write('E:\position.txt',x ' 'y ''minslp' ' tc,append) ;
*say rc
i=i+1 ;
endwhile
close('E:\position.txt') ;
'!cp position.txt position1.txt'

'reinit' ;
'open e:\2.ctl'
'set mproj latlon'
'set mpdset hires' ;
'set lon 110 135' ;
'set lat 15 35'
'set cmax -200' ;
'set xlint 2' ;
'set ylint 2' ;
'set xlopts 1 4 0.22' ;
'set ylopts 1 4 0.20' ;
'set clopts -1 -1 0.20' ;

'd sp' ;
read_file=read('E:\position1.txt') ;
read_line=sublin(read_file,2) ;
lon0=subwrd(read_line,1) ;
lat0=subwrd(read_line,2) ;
date0=subwrd(read_line,3) ;
'q w2xy ' lon0 ' ' lat0 ;
x0=subwrd(result,3) ;
y0=subwrd(result,6) ;
'draw mark 3 'x0' 'y0' 0.12' ;
'set string 2 '
'set strsiz 0.080 0.080' ;
'draw string 'x0-0.1' 'y0+0.2' 'date0''

i=1
while(i<=2000) ;
xbefore=x0
ybefore=y0
read_file=read('E:\position.txt')
read_code=sublin(read_file,1) ;
if(read_code>0);break;endif;
read_line=sublin(read_file,2) ;
lon0=subwrd(read_line,1) ;
lat0=subwrd(read_line,2) ;
date0=subwrd(read_line,3) ;
'q w2xy ' lon0 ' ' lat0 ;
x0=subwrd(result,3) ;
y0=subwrd(result,6) ;
'draw title 0808 "TIANCHENG" track by 20121508018 zx'
'draw mark 3 'x0' 'y0' 0.12'
'draw line 'xbefore' 'ybefore' 'x0' 'y0''
'draw string 'x0-0.1' 'y0+0.2' 'date0'';
i=i+1
endwhile

'printim E:\typhoon.gif  x400 y300' ;
'print'
'disable print

现在出现了以下的情况，sp是不合理的变量描述，“position.txt” can not stat

但是E盘里面生成了position.TXT文件，是这样的

，请问如何解决？谢了！

wk1021

是gs文件那里编写错误了吗？是照着论坛里的程序一步步修改调试的。。。

wk1021

谢谢上次帮忙的@夏天很美，@惠子两位小伙伴，谢谢你们了~

wk1021

希望这次也能有好心人麻烦帮我看一下~

lanlaniris

我用过这个gs画路径，不过我用的不是nc资料不知道~你后面的问题主要就是文件没读进去导致sp变量不能用吧~你要不试一下，专门写一个gs来读读资料，等成功了再带到别的语句里用~

wk1021

lanlaniris 发表于 2015-4-25 10:18
我用过这个gs画路径，不过我用的不是nc资料不知道~你后面的问题主要就是文件没读进去导致sp变量不能用吧~你 ...

恩，那应该如何操作呢？先写一个什么样的gs呢？然后代入这个gs文件么？文件格式有问题么？nc，grib原始资料文件我都有啊，感觉应该和文件类型关系不大啊。。。

lanlaniris

wk1021 发表于 2015-4-25 11:51
恩，那应该如何操作呢？先写一个什么样的gs呢？然后代入这个gs文件么？文件格式有问题么？nc，grib原始资 ...

后面的问题都是你的文件打开失败,没有成功读取变量造成的.我的意思是你单纯写一个gs文件，里面只有打开NC数据和它的描述文件，然后看能不能成功运行。只是单纯的测试。因为我看你在GRADS里面直接打开时是成功的，但是用批处理文件就有打开错误。
另外，我今天也运行了这个gs，也有出现语法错误的提示(图上的query error：Syntax is QUERY W2XY……)这个好像是文件读到最尾了，没有成功读入经纬度的数据。原本这个脚本里面是有筛选的（就是if语句那一行），但不知道为什么没用，你得自己设定一下循环的次数。

wk1021

lanlaniris 发表于 2015-4-25 19:59
后面的问题都是你的文件打开失败,没有成功读取变量造成的.我的意思是你单纯写一个gs文件，里面只有打开NC ...

ok,我来试试，thanks~

wk1021

lanlaniris 发表于 2015-4-25 19:59
后面的问题都是你的文件打开失败,没有成功读取变量造成的.我的意思是你单纯写一个gs文件，里面只有打开NC ...

编写了测试gs如下：’sdfopen e:\2010tc.nc' 'open e:\2010tc.ctl'。运行正确，循环次数也改了，图出来是错误的，应该还是文件没有正确打开，变量没有正确读取导致的，苦恼啊，怎么修改呢~

wk1021

lanlaniris 发表于 2015-4-25 19:59
后面的问题都是你的文件打开失败,没有成功读取变量造成的.我的意思是你单纯写一个gs文件，里面只有打开NC ...