【求助】站点插值数据后经过计算想用set fwrite写入数据，出现如图所示，求解！

Sylvia

站点已经处理成格点数据，grid.grd和station.grd都已生成，经过插值图也成功出来了。但是我现在想用插值后的数据，即把图不显示出来，而将数据写到文件中，再画图，结果，发现图总是错位。并且set fwrite后，再run的话，没有反应，请问什么原因呢，怎么解决，是不是站点插值格点的数据用set fwrite保存不下来？.如图所示：


ga-> reinit
No hardcopy metafile open
All files closed; all defined objects released;
All GrADS attributes have been reinitialized
ga-> open f:\air\grid2.ctl
Scanning description file:  f:\air\grid2.ctl
Data file f:\air\grid.grd is open as file 1
LON set to 70 149
LAT set to 15 59
LEV set to 1000 1000
Time values set: 1951:1:1:0 1951:1:1:0
E set to 1 1
ga-> open f:\air\t.ctl
Scanning description file:  f:\air\t.ctl
Data file f:\air\65ave.grd is open as file 2
ga-> run f:\air\tt.gs
ga->


很奇怪，run之后就立即出现ga->，无其他提示信息，也出来了grd，但很小，然后再画图就是错位的，请问怎么解决。
多谢各位了

mofangbao

楼主下次提问前请先搜索：
http://bbs.06climate.com/forum.php?mod=viewthread&tid=7357

Sylvia

mofangbao 发表于 2012-9-8 18:23

                               
登录/注册后可看大图

楼主下次提问前请先搜索：
http://bbs.06climate.com/forum.php?mod=viewthread&tid=7357

我已经看过这个帖子了，只是它输出的是txt，我想要的不是这个格式的数据，我想要保存为grd，本来打算用set fwrite的，但是似乎用不了

mofangbao

你的gs贴出来

Sylvia

mofangbao 发表于 2012-9-8 20:38

                               
登录/注册后可看大图

你的gs贴出来

gs:将set fwrite改成enable print可以出图，但这样没有结果
'open f:\air\grid2.ctl'
'open f:\air\t.ctl'
'set fwrite f:\air\czj.grd'
'set gxout fwrite'
'set grads off'
'set grid off'
'set cthick 5'
'set lat 15 60'
'set lon 70 150'
'set mpdset cnworld'
'define a=oacres(g,t.2,1.5)'
'define a1=maskout(a,g-0.5)'
'define aa=smth9(a1)'
'set lat 16 56'
'set lon 75 135'
'set cint 2'
*'d aa*0.1'
'cnbasemap aa*0.1'
'draw title 2012 China Spring Anomaly'
'run china -p'
'disable fwrite'
;


站点转成格点的文件的ctl：
dset f:\air\65anos.grd
dtype station
stnmap f:\air\ano.map
undef -9.99e33
title rain
tdef 62 linear JAN1951 1yr
vars 1
t 0 99 rainfall data
endvars


格点文件的ctl
dset f:\air\grid.grd
undef -9.99e33
title Sample Grid Data
xdef 80 linear 70 1
ydef 45 linear 15 1
zdef 1 linear 1000 1
tdef 1 linear JAN1951 1yr
vars 1
g 0 99 grid data
endvars

Sylvia

Sylvia 发表于 2012-9-8 22:11

                               
登录/注册后可看大图

gs:将set fwrite改成enable print可以出图，但这样没有结果
'open f:\air\grid2.ctl'
'open f:\air\t.c ...

我现在的问题就是set fwrite出的结果grd不对或者出不来

Sylvia

mofangbao 发表于 2012-9-8 20:38

                               
登录/注册后可看大图

你的gs贴出来

我看过你写的那个站点数据处理的pdf，很全很好，我也一一对应了下，没有发现问题，我之前也试图做过很多尝试的。就出在结果了，图是出来了，和你的结果是一样的，可惜你没有写怎么处理出来的grd，哎，我整了半天也没出来，本来还有点数据的，后来整的都没有了，图总是错位。。。。感谢大神啊，能解决我万分感谢啊，好好学学怎么彻底的处理站点数据。

mofangbao

把你的数据传上看看，如果不是保密级别的数据的话，还有你的ctl什么的，所有的打个包上来，数据如果很大的话，就提取一个时次的上来就行了

Sylvia

mofangbao 发表于 2012-9-9 15:45

                               
登录/注册后可看大图

把你的数据传上看看，如果不是保密级别的数据的话，还有你的ctl什么的，所有的打个包上来，数据如果很大的话 ...

不好意思啊，这会才看到回复。那个我弄出来了，准备闲下来的时候写个帖子注明一下。后面的set fwrite 用法是一样用，只是个别地方需要修改下就可以了。麻烦了。资料一点都不保密，呵呵，就是国家160站的资料。

mofangbao

Sylvia 发表于 2012-9-13 19:31

                               
登录/注册后可看大图

不好意思啊，这会才看到回复。那个我弄出来了，准备闲下来的时候写个帖子注明一下。后面的set fwrite 用法 ...

好的，别忘了哈~