GRADS 画格点数据的地形场和站点数据的风场 求问

yanglei_nuist

最近从论坛里找到大家提供的地形数据，学习GRADS画风场和地形场，但是画完的图里只有地形数据，没有风场，请大家给予帮助。
'reinit'
*分别打开地形和风场的描述文件
'sdfopen D:\hourt\china.nc'
'open d:\hourt\HOURT.ctl'
m=1
while(m<25)
   'set grads off'
   'set grid off'  
  'set lon 118 126'
  'set lat 38  44'
  'set mpdset cnworld'
   'set mpdset liaoning liaoning_q cnworld'
   'set csmooth on'
*因为地形数据里时间只有一个时次，所以这里设置时间t=1
   'set t 1'
  'set lon 118 126'
  'set lat 38  44'
   'set gxout shaded'
   'run D:\hourt\define_colors.gs'
   'set clevs  200 300 400 500 600 700  800'
   'set ccols 39 37  33   21  23  25  27  29'
   'd china.1'
   'cbarn'

*现在开始画风场
   'set t ' m
   'set gxout barb'
   'set cthick 5'
   'set digsiz 0.04'
   'D UUU.2*2.5;VVV.2*2.5'
  'printim D:\hourt\wind_height\hour_'m'.png x1280 y1024 white'
   m=m+1
   'c'
endwhile
'print'  
'disable print'
;


china.nc (1.79 MB, 下载次数: 144)

2014-2-24 16:40 上传

点击文件名下载附件 HOURT.ctl (221 Bytes, 下载次数: 53)

2014-2-24 16:40 上传

点击文件名下载附件 define_colors.gs (1.53 KB, 下载次数: 56)

2014-2-24 16:40 上传

点击文件名下载附件 grid.ctl (200 Bytes, 下载次数: 34)

2014-2-24 16:40 上传

点击文件名下载附件 wind height.gs (790 Bytes, 下载次数: 32)

2014-2-24 16:40 上传

点击文件名下载附件





在大家的帮助下弄好了！
'reinit'
m=1
while(m<25)
   'sdfopen D:\hourt\china.nc'
   'set grads off'
   'set grid off'   
   'set lon 118 126'
   'set lat 38  44'
   'set mpdset cnworld'
   'set mpdset liaoning liaoning_q cnworld'
   'set csmooth on'
   'set t 1'
   'set lon 118 126'
   'set lat 38  44'
   'set gxout shaded'
   'run D:\hourt\define_colors.gs'
   'set clevs  200 300 400 500 600 700  800'
   'set ccols 39 37  33   21  23  25  27  29'
   'd china.1'
   'cbarn'
   'close 1'

   'open d:\hourt\HOURT.ctl'
   'set lon 118 126'
   'set lat 38  44'
   'set t ' m
   'set gxout barb'
   'set cthick 5'
   'set digsiz 0.04'
   'D uuu;vvv'
   'close 1'

   'printim D:\hourt\wind_height\hour_'m'.emf x1280 y1024 white'
   m=m+1
   'c'

endwhile
'print'  
'disable print'
;

river

是不是这句'set t ' m 多了一个空格啊，应该是'set t 'm 吧，试一下看看。顺便问一下emf是什么格式的图片，可以用什么软件打开啊？

容与

才学grads 的飘过帮顶一个，弱弱的问一下那个'D UUU.2*2.5;VVV.2*2.5' 中，不应该是小写的D吗

mofangbao

把各种辅助设置都去掉 只保留最必须的绘制语句看看行不行 然后再贴一下grads的提示和你现在出来的图看看

北京小脚

关注一下，希望楼主按着管理员给的建议给个提示和图

yanglei_nuist

river 发表于 2014-2-24 19:10

                               
登录/注册后可看大图

是不是这句'set t ' m 多了一个空格啊，应该是'set t 'm 吧，试一下看看。顺便问一下emf是什么格式的图片， ...

不是空格的问题啊，我弄了下还是只有地形数据

yanglei_nuist

容与 发表于 2014-2-25 06:18

                               
登录/注册后可看大图

才学grads 的飘过帮顶一个，弱弱的问一下那个'D UUU.2*2.5;VVV.2*2.5' 中，不应该是小写的D吗

这里大写也可以

yanglei_nuist

mofangbao 发表于 2014-2-25 08:33

                               
登录/注册后可看大图

把各种辅助设置都去掉 只保留最必须的绘制语句看看行不行 然后再贴一下grads的提示和你现在出来的图看看

都去掉了，还是不行啊，附件里有我按照原来程序画的图，只有地形数据，没有风场。我的这个风场是站点数据，是不是站点数据和NC文件一起画有什么要求啊？

yanglei_nuist

river 发表于 2014-2-24 19:10

                               
登录/注册后可看大图

是不是这句'set t ' m 多了一个空格啊，应该是'set t 'm 吧，试一下看看。顺便问一下emf是什么格式的图片， ...

EMF也是图片的一种格式，电脑都可以打开，放大后更清晰些。

qixiangren34

"我的这个风场是站点数据，是不是站点数据和NC文件一起画有什么要求啊？”这句是关键，站点数据经过差值转换为格点数据后，就可以画等值线了