- 积分
- 2114
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-9-6
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 tvfantasy 于 2017-5-26 12:46 编辑
已经解决
- 'reinit'
- 'sdfopen D:\data\uwnd.2016.nc'
- 'sdfopen D:\data\vwnd.2016.nc'
- *set map*
- 'set mpdset cnworld'
- t=737
- while(t<741)
- <font color="red">'set lon 70 170'
- 'set lat 0 60'</font>
- 'set lev 700'
- 'set t 't
- 'define vort=hcurl(uwnd.1,vwnd.2)*10000'
- <font color="red">'set lon 80 160'
- 'set lat 10 50'
- 'd vort'</font>
- *set title with elements and time*
- 'q time'
- rc=sublin(result,1)
- tt=subwrd(rc,3)
- 'draw title 700hpa_Vorticity_'%tt%''
- 'printim D:\by_pic\700hpa_Vorticity_'%tt%'.png white x1000 y800'
- 'c'
- t=t+1
- endwhile
- 'reinit'
复制代码
求助!绘制的700hpa水平涡度分布图,四周有一圈空白是什么情况?
脚本
- 'reinit'
- 'sdfopen D:\data\uwnd.2016.nc'
- 'sdfopen D:\data\vwnd.2016.nc'
- *set map*
- 'set mpdset cnworld'
- 'set lon 80 160'
- 'set lat 10 50'
- t=737
- while(t<741)
- 'set lev 700'
- 'set t 't
- 'define vort=hcurl(uwnd.1,vwnd.2)*100000'
- 'd vort'
- *set title with elements and time*
- 'q time'
- rc=sublin(result,1)
- tt=subwrd(rc,3)
- 'draw title 700hpa_Vorticity_'%tt%''
- 'printim D:\by_pic\700hpa_Vorticity_'%tt%'.png white x1000 y800'
- 'c'
- t=t+1
- endwhile
- 'reinit'
复制代码
|
|