爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 2373|回复: 1

[脚本编辑] 求助大神,grads求水汽通量出图结果为什么都是一样的?

[复制链接]

新浪微博达人勋

发表于 2016-4-17 17:15:27 | 显示全部楼层 |阅读模式

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

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

x
求的是10年内季节平均水汽通量代码如下:
'reinit'
'set grid off'
'set grads off'
'set poli on'
'sdfopen D:\Grads\shum.mon.mean.nc'
'sdfopen D:\Grads\uwnd.mon.mean.nc'
'sdfopen D:\Grads\vwnd.mon.mean.nc'
'sdfopen D:\Grads\pres.mon.mean.nc'
'set fwrite d:\aut_shuiqi.grd'
'set gxout fwrite'
'set lat -10 55'
'set lon 10 140'
'set xlpos -20'                                                                           
'set ylpos -20'
'set mpdset cnworld'
'set map 1 1 1'
'set cthick 5'

'set t 601 721'
'set z 1 8'
'define uq=shum.1*uwnd.2'
'define vq=shum.1*vwnd.3'
'define v=vwnd.3'
'define u=uwnd.2'

'set z 1'
'define uqs=vint(pres.4(z=1),uq,300)/1000'
'define vqs=vint(pres.4(z=1),vq,300)/1000'
'define sq=mag(uqs,vqs)'

'set t 1 12'
'define w1=ave(sq,t=601,t=721,12)'
'define w2=ave(uqs,t=601,t=721,12)'
'define w3=ave(vqs,t=601,t=721,12)'

'set t 1'
'define w1ave=sum(w1,t=9,t=11)'
'define w2ave=sum(w2,t=9,t=11)'
'define w3ave=sum(w3,t=9,t=11)'

*发现修改t=9,t=11、t=6,t=8或者t=3,t=5,出来的图都一样,求助如何修改,谢谢  

'set grads off'
'set gxout shaded'
'set ccolor 1'
'set gxout vector'
'set ccolor 1'
'd w2ave;w3ave'
'cbarn 1 0 5.5 7.8'
'printim e:\aut.png x1000 y800 white'
'disable fwrite'
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-4-17 20:29:46 | 显示全部楼层
你这程序和我之前用的求水汽通量不一样啊,论坛搜一搜吧
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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