- 积分
- 210
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-6-26
- 最后登录
- 1970-1-1
![未绑定新浪微博用户 新浪微博达人勋](source/plugin/sina_login/img/gray.png)
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
使用re函数出现问题。我想把1°*1°海温资料使用re函数转化成2°*2°的资料海温原始数据资料的ctl描述文件
dset g:\sst_1.dat
undef -999.0
title hsst
xdef 360 linear -180.0 1.0
ydef 180 linear -90. 1.0
zdef 1 levels 1000
tdef 612 linear 00z01jan1950 1mon
vars 1
sst 0 99 globle ocean sst
endvars
我读取时写的re的gs文件
'reinit'
'open G:\sst1.ctl'
'set gxout fwrite'
'set fwrite G:\sst2.grd'
'set t 1'
'set lev 1000'
'set x 1 360'
'set y 1 180'
*'set lat -90 89'
*'set lon 0 360'
*'define s1=re(sst,2,2,bl)'
'define s1=re(sst,180,linear,0,2,91,linear,-90,2,bl)'
'd s1'
'disable fwrite'
;
数据读出后搭配的ctl描述文件
dset g:\sst2.grd
undef -999.0
title sst
xdef 180 linear 0 2.0
ydef 91 linear -90. 2.0
zdef 1 levels 1000
tdef 1 linear 00z01jan1950 1mon
vars 1
sst 0 99 globle ocean sst
endvars
同一时次画出的图
原数据图
读出图
请教,从图中看感觉数据应该是读入进去了,就是经纬度偏移太厉害,海洋上的数据都到陆地上了。不管是re函数的简写还是所有参数都设置的写法都用了,不管是设置经纬度set lat,lon还是设置格点set x,y都试了,结果都是这种。
请问这种re函数应该怎么使用?我的gs和ctl有什么问题?谢谢!!
|
|