爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 11311|回复: 19

[秀图] 求助水汽通量散度的画图问题

[复制链接]
发表于 2012-10-18 19:08:36 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
画出的水汽通量散度图,如何将周边没有画到的部分与图边的经纬度连接起来呢?file:///C:/Users/iap-nju-qian/AppData/Roaming/Tencent/Users/12270645/QQ/WinTemp/RichOle/FKGIMKUMAO%25(2RW%60J%60EJO(O.jpg

密码修改失败请联系微信:mofangbao
发表于 2012-10-18 20:04:26 | 显示全部楼层
图挂了,请楼主上图,没明白你的意思
密码修改失败请联系微信:mofangbao
发表于 2012-10-18 20:25:36 | 显示全部楼层
我也没明白,楼主上图,把问题说明白点儿啊
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2012-10-19 09:29:17 | 显示全部楼层
图是这样的 850hPa-divergence-of-moisture-flux00Z21JUL2012.gif ,有什么办法能把水汽通量散度覆盖到整个要画的图面?
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2012-10-19 09:29:42 | 显示全部楼层
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2012-10-19 09:30:01 | 显示全部楼层
天目神眉 发表于 2012-10-18 20:04
图挂了,请楼主上图,没明白你的意思

图示已发,望高手解答
密码修改失败请联系微信:mofangbao
发表于 2012-10-19 20:53:39 | 显示全部楼层
风暴侦察兵 发表于 2012-10-19 09:29
图是这样的,有什么办法能把水汽通量散度覆盖到整个要画的图面?

楼主这个是什么资料画的啊,如果是通过FORTRAN计算得来的话,应该是差分之后边界就少掉一圈。具体的我也不太清楚
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2012-10-20 16:54:27 | 显示全部楼层
river 发表于 2012-10-19 20:53
楼主这个是什么资料画的啊,如果是通过FORTRAN计算得来的话,应该是差分之后边界就少掉一圈。具体的我也不 ...

就是NCEP资料画的水汽通量散度,只是后续处理不太会,想把散度的分布覆盖到整个图面
密码修改失败请联系微信:mofangbao
发表于 2012-10-20 18:34:39 | 显示全部楼层
风暴侦察兵 发表于 2012-10-20 16:54
就是NCEP资料画的水汽通量散度,只是后续处理不太会,想把散度的分布覆盖到整个图面

如果是NC资料的话是不应该出现这种问题的。你把gs贴上来吧。我怀疑是你计算的时候设置的经纬度比你要画的经纬度范围小了
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2012-10-20 18:49:22 | 显示全部楼层
river 发表于 2012-10-20 18:34
如果是NC资料的话是不应该出现这种问题的。你把gs贴上来吧。我怀疑是你计算的时候设置的经纬度比你要画的 ...

'reinit'
'sdfopen e:\data\ncep\uwnd.2012.nc'
'sdfopen e:\data\ncep\vwnd.2012.nc'
'sdfopen e:\data\ncep\shum.2012.nc'

'set grads off'
'set grid off'
'set mpdset cnworld cnriver'
'set map 15 1 2'
'set xlopts 1 4 0.2'
'set ylopts 1 4 0.2'
'set parea 1.5 10 1 8'

nt=809
while(nt<=812)
'set t 'nt''
'set lat 20 45'
'set lon 100 135'
'set lev 850'
'set ylint 5'
'set xlint 5'
'set csmooth on'

*'run define_color.gs'
*'set ccols 4 11 5 13 3 10 12 8 2 27 68'
*'set clevs 2 4 6 8 10 12 14 16 18 20'
'set gxout shaded'
'd hdivg(shum.3*uwnd.1/9.8,shum.3*vwnd.2/9.8)'
'cbarn 0.85 0 5.75 0.35'

'q time'
rc=sublin(result,1)
tt=subwrd(rc,3)

'draw title 850hPa Divergence of Moisture Flux '
'printim e:\IAP\721\850hPa-divergence-of-moisture-flux'tt'.gif white'
'enable print e:\IAP\721\850hPa-divergence-of-moisture-flux'tt'.gmf'

'print'
'disable print'
'c'
'set grads off'
nt=nt+1
endwhile
;
经纬度就设置了一次,不太理解你所说的。。。。请指教
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表