求助：将6个变量的grd合并成一个grd时发生漂移

zxholystar

我有6个变量abcdef,想画qnet=d+e-a-b-c-f，但为什么把qnet生成一个grd文件画出的图形比不写grd直接画出的图形要明显向东偏移很多？因为模式比较多，所以我还是想生成一个grd文件。实在不知道为什么会这样第二张图为生成grd文件后画出的，明显漂移，第一张图是不生成grd直接画出
生成qnet.grd的gs
'reinit'
'set gxout fwrite'
'set fwrite E:\MIROC5-net.grd'
'open E:\25\MIROC5\1\hfls25.ctl'
'define A=ave(hfls.1,t=1,t=360)'
'open E:\25\MIROC5\1\hfss25.ctl'
'define B=ave(hfss.2,t=1,t=360)'
'open E:\25\MIROC5\1\rlus25.ctl'
'define C=ave(rlus.3,t=1,t=360)'
'open E:\25\MIROC5\1\rlds25.ctl'
'define D=ave(rlds.4,t=1,t=360)'
'open E:\25\MIROC5\1\rsds25.ctl'
'define E=ave(rsds.5,t=1,t=360)'
'open E:\25\MIROC5\1\rsus25.ctl'
'define F=ave(rsus.6,t=1,t=360)'
'define qnet=d+e-a-b-c-f'
'd qnet'
'disable fwrite'

生成的qnet.grd的ctl
dset E:\MIROC5-net.grd
undef 32766
title  MIROC5-net
xdef 144 linear 0 2.5
ydef 73 linear  -90  2.5
zdef 1 levels 0 0
tdef 1 linear 00z01JUL1979  1mon
vars 1
qnet 0 -999  net heat flux 1979-2000
endvars



直接画出的gs
'reinit'
'open E:\25\MIROC5\1\hfls25.ctl'
'define A=ave(hfls.1,t=1,t=360)'
'open E:\25\MIROC5\1\hfss25.ctl'
'define B=ave(hfss.2,t=1,t=360)'
'open E:\25\MIROC5\1\rlus25.ctl'
'define C=ave(rlus.3,t=1,t=360)'
'open E:\25\MIROC5\1\rlds25.ctl'
'define D=ave(rlds.4,t=1,t=360)'
'open E:\25\MIROC5\1\rsds25.ctl'
'define E=ave(rsds.5,t=1,t=360)'
'open E:\25\MIROC5\1\rsus25.ctl'
'define F=ave(rsus.6,t=1,t=360)'
'define qnet=d+e-a-b-c-f'
'set gxout shaded'
'set clevs -150 -120 -90 -60 -30 0 30 60 90 120 150 180'
'set ccols 9 14 4 11 5 13 3 10 7 12 8 2'
'd qnet'

lqouc

原始数据ctl呢？尽量保持一致。

zxholystar

lqouc 发表于 2014-10-21 13:00
原始数据ctl呢？尽量保持一致。

dset E:\25\ACCESS1-0\1\hfls25.grd
undef 1e+20
title  hfls
xdef 144 linear 0 2.5
ydef 73 linear  -90  2.5
zdef 1 levels 0 0
tdef 360 linear 12z14jan1979  1mon
vars 1
hfls 0 -999  latent heat flux
endvars
这是原始数据的ctl之一，其余5个都差不多，只改了变量名

river

什么都没设置啊，在你的d前面加上set x 1 144   set y 1 73  set t 1 试试

zxholystar

river 发表于 2014-10-21 16:30
什么都没设置啊，在你的d前面加上set x 1 144   set y 1 73  set t 1 试试

Ahhhhhhh! That works. river, u are really life saver! Thank you!