nc直接出图和nc转成grd 出图 不一样吗

nickbsb

这是nc直接出图的gs文件 1980年12月日平均资料 我画12月1号的风场图
'reinit'
'sdfopen f:\data\uwind\uwnd.1980.nc'
'set lat 20 90'
'set lon 0 360'
'set lev 500'
'set t 336'
'd uwnd'
*'set t 336 366'
*'d ave(uwnd,t=336,t=366)'
'enable print f:\198012\198012nc.gmf'
'print'
'disable print'
;
这是nc转成的对应的grd的gs文件

'reinit'
'sdfopen f:\data\uwind\uwnd.1980.nc'
'set fwrite f:\198012\198012.grd'
'set gxout fwrite'
'set x 1 145'
'set y 44 73'
'set lev 500'
day=336
while(day<=366)
'set t 'day''
'd uwnd'
day=day+1

endwhile
'close 1'

'disable fwrite'
;
俩个文件画出的图不一样 怎么回事 有错误吗

mofangbao

不要写  set lon 0 360   这样就多了一个点，用 set x start end  来表示你就理解了

叩容

比较一下你的图，看看你的ctl设置上经纬度有没有问题。
不一定照着NC文件给的信息，那个有可能有错。

冰枫

mofangbao 发表于 2012-2-24 14:17

                               
登录/注册后可看大图

不要写  set lon 0 360   这样就多了一个点，用 set x start end  来表示你就理解了

那在用转的grd画图的时候也不能设set lon嘛？

mofangbao

set lon 可以用，写对起止值就行，0 和360是同一条线，比如2.5的分辨率写0 357.5就可以