grads计算水汽通量求助！！！！（已解决，并更新）

slc2005011020

你好 能帮帮忙看看那错了吗？
'reinit'
'sdfopen f:\shuju2\hgt.mon.mean.nc'
'sdfopen f:\shuju2\uwnd.mon.mean.nc'
'sdfopen f:\shuju2\vwnd.mon.mean.nc'
'sdfopen f:\shuju2\shum.mon.mean.nc'
'sdfopen f:\shuju2\pres.mon.mean.nc'

'enable print F:\tu\sqtl.gmf'
'set lat 10 70'
'set lon 80 160'
'set z 1 8'
'set t 7 800'
'set grid off'
'set xlint 20'
'set ylint 10'
'set grads off'

'set mpdset cnworld'
'define qu=uwnd.2*shum.4*1000'
'define qv=vwnd.3*shum.4*1000'
'define qus=vint(pres.5(lev=0)/100,qu,300)/9.8'
'define qvs=vint(pres.5(lev=0)/100,qv,300)/9.8'

'set z 1'
'set t 7'

'define qum1=ave(qus,t+0,t=799,12)'
'define qvm1=ave(qvs,t+0,t=799,12)'

'set t 8'
'define qum2=ave(qus,t+0,t=801,12)'
'define qvm2=ave(qvs,t+0,t=801,12)'

'define qum=(qum1+qum2)/2'
'define qvm=(qvm1+qvm2)/2'


'd hdivg(qum,qvm)/100'

'print'
'disable print'

slc2005011020

错误提示是：Cannot contour grid -all undefined values

任我行

谢谢分享，很不错的资料哦！

柚子皮

好帖，几个帖子结合起来看，收获很大

a1s1

太阳雨 发表于 2012-10-20 20:06
十分感谢楼主的指教，我现在把修改过后的gs脚本和图贴出来，还希望楼主帮我看看是否正确，还有就是右下角 ...

请问一下，为什么我只是替换了一下nc文件，就出错呢
‘reinit'
'sdfopen e:\uwnd.mon.mean.nc'
'sdfopen e:\vwnd.mon.mean.nc'
'sdfopen e:\shum.mon.mean.nc'
'sdfopen e:\pres.sfc.mon.mean.nc'
'set lat -10 70'
'set lon 40 160'
'set z 1 8'
'set t 163 178'
'define qu=uwnd*shum.3*1000'
'define qv=vwnd.2*shum.3*1000'
'define qus=vint(pres.4(lev=0)/100,qu,300)/9.8'
'define qvs=vint(pres.4(lev=0)/100,qv,300)/9.8'
'define qum=ave(qus,t=163,t=178,1)'
'define qvm=ave(qvs,t=163,t=178,1)'
'set z 1'
'set t 163'
'd qum;qvm'
'print'

却出不了图，提示 Cannot draw vectors/barbs--U field all undefined

limoxi

wet510 发表于 2012-9-24 08:19
其实如果算整层的话不用循环液可以，如果你要求确定时间的一个区域的整层水汽通量，可以如下设置维度，比 ...

感谢楼主分享经验，最近正琢磨这个问题，想请问楼主按上面的程序
'define uq=shum.2*uwnd.3*1000'是直接得到了一个数组其中包括了每一层从6月1号到8月31号的比湿乘风分量的值吗？如果是这样，那'set z 1 8'和
'set time 12z01jun1998 12z31aug1998'
如果都是求和的意思还需要后面用vint做垂直积分吗?

limoxi

太阳雨 发表于 2012-10-20 20:06
十分感谢楼主的指教，我现在把修改过后的gs脚本和图贴出来，还希望楼主帮我看看是否正确，还有就是右下角 ...

您好，我最近正需要这个水汽垂直积分，对这个算法不是很明白，求指点！
我的资料只是1年的6,7,8三个月的月平均资料，所以t从1到3，我的gs文件如下：'reinit'
'sdfopen f:/data/u2014sum.nc'
'sdfopen f:/data/v2014sum.nc'
'sdfopen f:/data/shum2014sum.nc'
'sdfopen f:/data/pressur2014sum.nc'
'set lat 0 60'
'set lon 20 180'
'set z 1 8'
'set t 1 3'
'define uq=shum.3*uwnd.1*1000'
'define vq=shum.3*vwnd.2*1000'
'define uqs=vint(pres.4(lev=0)/100,uq,300)/9.8'
'define vqs=vint(pres.4(lev=0)/100,vq,300)/9.8'
'define qum=ave(uqs,t=1,t=3,1)'
'define qvm=ave(vqs,t=1,t=3,1)'
'set z 1'
'set t 1'
'set gxout shaded'
'define a=mag(qum,qvm)'
'd a'
'set gxout vector'
'd qum;qvm'
'gxprint'
;


得到却是这样的warning


出的图也就不是很可靠：


下面的单位也不是很明白
看了很多程序都是用循环的，这个程序没有用循环，也没有另写ctl文件，不明白定义的uq,vq是一个数组还是一个对z和t维求和之后的数字？

limoxi

a1s1 发表于 2015-4-14 09:13
请问一下，为什么我只是替换了一下nc文件，就出错呢
‘reinit'
'sdfopen e:\uwnd.mon.mean.nc'

看一下是不是数据的范围不合适？

wet510

a1s1 发表于 2015-4-14 09:13
请问一下，为什么我只是替换了一下nc文件，就出错呢
‘reinit'
'sdfopen e:\uwnd.mon.mean.nc'

在输出前，要设定你的图像形式，也就是要加一句
'set gxout vector' 这个语句放到'd ups;vps'之前

wet510

limoxi 发表于 2015-5-1 20:44
感谢楼主分享经验，最近正琢磨这个问题，想请问楼主按上面的程序
'define uq=shum.2*uwnd.3*1000'是直接 ...

垂直积分还是需要的，四维数据，如果你不固定任一维，那么计算得到的就是四维数组，你固定了一维，得出的就是3维数组，依次类推