grads画风速距平场，报错问题

momo_2013

'reinit'

'sdfopen E:\method\\uwnd.mon.mean.nc'

'sdfopen E:\method\\vwnd.mon.mean.nc'

'set grads off'

'set grid off'

'set lon 40 180'

'set lat -10 70'

'set lev 850'

'set t 1'

*计算1981-2010年各月风速U常年值----------------------------

'set t 1 12'

'haveu=ave(uwnd.1,t+396,t=756,12)'

*计算1961-2014年12月各月风速U距平值

'modify have seasonal'

'set t 157 804'

'u=uwnd.1-haveu'

'set t 1'

*计算2014年夏季风速距平值

'define us=ave(u,t=798,t=800)'

*计算1981-2010年各月风速V常年值----------------------------

'set t 1 12'

'havev=ave(vwnd.2,t+396,t=756,12)'

*计算1961-2014年12月各月风速V距平值

'modify have seasonal'

'set t 157 804'

'v=vwnd.2-havev'

'set t 1'

*计算2014年夏季风速距平值

'define vs=ave(v,t=798,t=800)'

'set gxout vector'

'd  us;vs'

'cbar 1  1'

'draw title UV anomaly at 850hPa in summer of 2014'

原始风速场可出图，用上述方法画水汽输送距平场也能出图，看着也正常，但计算风速气候场画距平场怎么不出图啊？报错如下：“U field all undefiled”

风之精灵

你的程序中前后两次使用'modify have seasonal'语句，根据你的意思，第一次应该是'modify haveu seasonal',第二次应该是'modify havev seasonal'吧

mofangbao

如果不能熟练的用隐循环，请改为while循环，代码便也会清晰多了，自己循环计算平均值，再计算距平，最后计算的时候两个变量的t要统一，你set t只是指定了第一个变量的，第二个是自己计算出来的，你指定的那些自然是缺测的

momo_2013

mofangbao 发表于 2015-12-31 13:54
如果不能熟练的用隐循环，请改为while循环，代码便也会清晰多了，自己循环计算平均值，再计算距平，最后计 ...

没学过grads，只是照搬了实习范例和参考论坛的一个同学的程序，可是画别的要素都出图了，原样不动换成风速怎么就不出图呢，我疑惑点在这里。画水汽输送距平场可以出图

momo_2013

画水汽输送距平场可以出图，也是矢量的

momo_2013

风之精灵 发表于 2015-12-31 14:19
你的程序中前后两次使用'modify have seasonal'语句，根据你的意思，第一次应该是'modify haveu seasonal', ...

太感激了，果然是这个问题！万分感谢！

风之精灵

momo_2013 发表于 2015-12-31 14:26
太感激了，果然是这个问题！万分感谢！

不客气，

夏天夏夏夏

求问楼主，你的水汽输送距平图是怎么画的呢