【已解决】Cannot draw vectors/barbs -- U field all undefined

say苹果有虫

问题：grads 报错：Cannot draw vectors/barbs -- U field all undefined
描述：
我有一个 .nc文件，一个wrf输出文件.dat配有.ctl，两个文件含有相同的变量。想在一个page里画两张图以便对比，然而总是报错：Cannot draw vectors/barbs -- U field all undefined，只能画出来第一张图；但是两个文件分开画就可以正常输出……这是为什么？求各位大侠指导！

附gs脚本及配图：
1.合在一起画的：

'reinit'
'open F:\ex\haze_0901_cv5\plot_type_press\haze_0901_press\case1\press.ctl'
'sdfopen F:\data\ERA-interim\north2009_0101__0228.nc'
'enable print F:\ex\haze_0901_cv5\plot_type_press\haze_0901_press\case1\compare.gmf'
*'define_colors'

'set vpage 0 11 0 8.5'
'set parea 1 5.5 1 7'
'set lat 18 33'
'set lon 108 123'
'set z 1'
'set time 00z22jan2009'
'set grads off'
'set map 9 1 6'
'set gxout vector'
'set arrscl 0.5 10'
'set arrowhead 0.08'
'd skip(U10.1,15,15);skip(V10.1,15,15)'
'draw title 00UTC22Jan \\ WRF wind at 10m'

'set vpage 0 11 0 8.5'
'set parea 6 10.5 1 7'
'set lat 18 33'
'set lon 108 123'
'set z 1'
'set time 00z22jan2009'
'set grads off'
'set map 9 1 6'
'set gxout vector'
'set arrscl 0.5 10'
'set arrowhead 0.08'
'd u10.2;v10.2'
'draw title 00UTC22Jan \\ ERA-interim wind at 10m'
'print'

'reinit'

2.分开画：
2.1  .nc文件
'reinit'
'sdfopen F:\data\ERA-interim\north2009_0101__0228.nc'
'enable print F:\ex\haze_0901_cv5\plot_type_press\haze_0901_press\case1\compare_era.gmf'
*'define_colors'

'set vpage 0 11 0 8.5'
'set parea 1 10 1 7'
'set lat 18 33'
'set lon 108 123'
'set z 1'
'set time 00z22jan2009'
'set grads off'
'set map 9 1 6'
'set gxout vector'
'set arrscl 0.5 10'
'set arrowhead 0.08'
'd u10;v10'
'draw title 00UTC22Jan \\ ERA-interim wind at 10m'
'print'
'reinit'


2.2  .dat
'reinit'
'open F:\ex\haze_0901_cv5\plot_type_press\haze_0901_press\case1\press.ctl'
'enable print F:\ex\haze_0901_cv5\plot_type_press\haze_0901_press\case1\compare_wrf.gmf'
*'define_colors'

'set vpage 0 11 0 8.5'
'set parea 1 10 1 7'
'set lat 18 33'
'set lon 108 123'
'set z 1'
'set time 00z22jan2009'
'set grads off'
'set map 9 1 6'
'set gxout vector'
'set arrscl 0.5 10'
'set arrowhead 0.08'
'd skip(U10.1,15,15);skip(V10.1,15,15)'
'draw title 00UTC22Jan \\ WRF wind at 10m'
'print'
'reinit'


#------------------------------------------------------------------------------------------------------------------------------------------------
#------------------------------------------       已解决   -------------------------------------------------------------------------------------
#------------------------------------------------------------------------------------------------------------------------------------------------
'reinit'
'open F:\ex\haze_0901_cv5\plot_type_press\haze_0901_press\case1\press.ctl'
'sdfopen F:\data\ERA-interim\north2009_0101__0228.nc'
'enable print F:\ex\haze_0901_cv5\plot_type_press\haze_0901_press\case1\compare.gmf'
*'define_colors'

'set vpage 0 11 0 8.5'
'set parea 1 5.5 1 7'
'set dfile 1'
'set lat 18 33'
'set lon 108 123'
'set z 1'
'set time 00z22jan2009'
'set grads off'
'set map 9 1 6'
'set gxout vector'
'set arrscl 0.5 10'
'set arrowhead 0.08'
'd skip(U10.1,15,15);skip(V10.1,15,15)'
'draw title 00UTC22Jan \\ WRF wind at 10m'


'set vpage 0 11 0 8.5'
'set parea 6 10.5 1 7'
'set dfile 2'
'set lat 18 33'
'set lon 108 123'
'set z 1'
'set time 00z22jan2009'
'set grads off'
'set map 9 1 6'
'set gxout vector'
'set arrscl 0.5 10'
'set arrowhead 0.08'
'd u10.2;v10.2'
'draw title 00UTC22Jan \\ ERA-interim wind at 10m'
'print'
'c'

'reinit'

say苹果有虫

给力的小伙伴儿指出，多文件绘图，之前最好加上“'set dfile file#”，file#=1,2,3...问题解决~


修改后：
'reinit'
'open F:\ex\haze_0901_cv5\plot_type_press\haze_0901_press\case1\press.ctl'
'sdfopen F:\data\ERA-interim\north2009_0101__0228.nc'
'enable print F:\ex\haze_0901_cv5\plot_type_press\haze_0901_press\case1\compare.gmf'
*'define_colors'

'set vpage 0 11 0 8.5'
'set parea 1 5.5 1 7'
'set dfile 1'
'set lat 18 33'
'set lon 108 123'
'set z 1'
'set time 00z22jan2009'
'set grads off'
'set map 9 1 6'
'set gxout vector'
'set arrscl 0.5 10'
'set arrowhead 0.08'
'd skip(U10.1,15,15);skip(V10.1,15,15)'
'draw title 00UTC22Jan \\ WRF wind at 10m'


'set vpage 0 11 0 8.5'
'set parea 6 10.5 1 7'
'set dfile 2'
'set lat 18 33'
'set lon 108 123'
'set z 1'
'set time 00z22jan2009'
'set grads off'
'set map 9 1 6'
'set gxout vector'
'set arrscl 0.5 10'
'set arrowhead 0.08'
'd u10.2;v10.2'
'draw title 00UTC22Jan \\ ERA-interim wind at 10m'
'print'
'c'

'reinit'
C:\Users\Administrator\Desktop\QQ截图20160714154153.png

WZZ809525579

为什么我只打开了一个文件，也报这个错误

versatile

WZZ809525579 发表于 2017-3-29 16:52
为什么我只打开了一个文件，也报这个错误

我也出现这种错误，你解决了吗？

orangez

WZZ809525579 发表于 2017-3-29 16:52
为什么我只打开了一个文件，也报这个错误

我也是诶，有人知道吗？

Kyrie_x

单文件报错有人解决了吗 解决方案为把'set t'改为'set time ...'……）不知道有无参考价值。
'reinit'
'sdfopen D:\GrADS\test8\uwnd.mon.mean.nc'
'sdfopen D:\GrADS\test8\vwnd.mon.mean.nc'
'set vpage 5.5 11.0 0 8.5'
'set parea 1.0 5.0 0 8.5'
'set lat 10 50'
'set lon 100 135'
'set lev 1000'
'set t 00z01JAN2008'

'set gxout vector'
'set xlint 20'
'set ylint 10'
'set xlopts 1 5 0.2'
'set ylopts 1 5 0.2'
'set map 1 1 5'
'set arrscl 0.5 3'
'd uwnd.1; vwnd.2'
'draw title (b) JAN2008 Anomalies'
;

Houckin

Kyrie_x 发表于 2022-05-20 09:32
本帖最后由 Kyrie_x 于 2022-5-20 09:40 编辑

单文件报错有人解决了吗 解决方案为把'set t'改为'set time ...'……）不知道有无参考价值。
'reinit'
'sdfopen D:\GrADS\test8\uwnd.mon.mean.nc'
'sdfopen D:\GrADS\test8\vwnd.mon.mean.nc'
'set vpage 5.5 11.0 0 8.5'
'set parea 1.0 5.0 0 8.5'
'set lat 10 50'
'set lon 100 135'
'set lev 1000'
'set t 00z01JAN2008'

'set gxout vector'
'set xlint 20'
'set ylint 10'
'set xlopts 1 5 0.2'
'set ylopts 1 5 0.2'
'set map 1 1 5'
'set arrscl 0.5 3'
'd uwnd.1; vwnd.2'
'draw title (b) JAN2008 Anomalies'
;

好像是啊，x,y,z,t对应，lon,lat,lev,time对应？

来自安卓APP客户端