画水汽通量为什么直接画图和写成数据后再画图不相同呢？（解决）

洗雨

运用NCEP2.5*2.5的资料计算水汽通量，直接用资料画出来的图和把水汽通量资料写为dat数据文件后再画图，图形不匹配，弄了好久，不知道怎么办啊
*下面是把水汽通量写到一个文件中
'reinit'
'open G:\NCEP\shum.1992.nc'
'open G:\NCEP\uwnd.1992.nc'
'open G:\NCEP\vwnd.1992.nc'
'set gxout fwrite'
'set fwrite g:\data\qu.dat'
mt=1;while(mt<=31)
'set t 'mt
z=1;while(z<=8)
'set x 1 144'
'set y 1 73'
'set z 'z
'define qu=(shum.1*uwnd.2)'
*'define qv=(vwnd.3*shum.1)'
'd qu'
*'d qv'
z=z+1;endwhile
mt=mt+1;endwhile
'disable fwrite'
'reinit'
'open G:\NCEP\shum.1992.nc'
'open G:\NCEP\uwnd.1992.nc'
'open G:\NCEP\vwnd.1992.nc'
'set gxout fwrite'
'set fwrite g:\data\qv.dat'
mt=1;while(mt<=31)
'set t 'mt
z=1;while(z<=8)
'set x 1 144'
'set y 1 73'
'set z 'z
'define qv=(shum.1*vwnd.3)'
'd qv'
z=z+1;endwhile
mt=mt+1;endwhile
'disable fwrite'
描述文件写错，所以后面怎么画图都不对，当时应该将图的范围画大一点，这样就容易知道是什么问题了（我出错的原因是y的起始纬度设错了）
dset g:\data\qu.dat
title daily qu at every lev(1992)
undef -9.99e+33
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 8 levels 1000 925 850 700 600 500 400 300
tdef 31 linear 00z01jan1992 1dy
vars 1
qu 8 99 vapor flux at jingxiang
endvars

DSET m:\data\qv.dat
title daily qv at every lev(1992)
undef -9.99e+33
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 8 levels 1000 925 850 700 600 500 400 300
tdef 31 linear 00z01jan1992 1dy
vars 1
qv 8 99 vapor flux at weixiang
endvars
上图为成生文件画图，下图为NCEP资料直接画图
‘reinit'
'open g:\gs\1992\qu.ctl'
'set lon 30 160'
'set lat 20 60'
'set t 1'
'set z 1'
'd qu'
'printim g:\picture\z1t11992qu.jpg white'


'reinit'
'sdfopen g:\NECP2_DAILY\shum.1992.nc'
'sdfopen g:\NECP2_DAILY\uwnd.1992.nc'
'sdfopen g:\NECP2_DAILY\vwnd.1992.nc'
'set lon 30 160'
'set lat 20 60'
'set mpdset cnworld'
'set z 1'
'set t 1'
'define qu=(uwnd.2*shum.1*1e3)'
'define qv=(vwnd.3*shum.1*1e3)'
*'set gxout shaded'
*'d mag(qu,qv)'
*'d qu;qv'
'd qu'
*'d qv'
'printim g:\picture\qu23.jpg white'
图放置的顺序乱了，但是验证是对的

wheelar

你只给出了生成qu.dat,qv.dat的gs，画出两个图形的gs呢？

洗雨

wheelar 发表于 2015-5-11 21:14
你只给出了生成qu.dat,qv.dat的gs，画出两个图形的gs呢？

你好，谢谢你的回复，出图的gs我只设置了经纬度，时间 层次均为1，直接画图。当时画图还试过ctl的options修改，画出的图都不对