爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 3144|回复: 3

[脚本编辑] 为什么使用set xaxis命令后出图与实际不符

[复制链接]
发表于 2016-5-8 21:43:15 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 wty118 于 2016-5-8 21:47 编辑

'reinit'
'open e:\grib\201604\fnl_201604.ctl'
'enable print e:\picture\thetase(huanghua).gmf'
tt=57
while(tt<=61)
'set t 'tt
'set lat 28.18'
'set lon 113.22'
'set lev 1000 200'
'set xaxis 315 355'   

'define prs=lev'
'set grads off'
'set grid off'
'set mpdset cnworld'
'set font 1'
'define es=(6.112*exp(17.67*(Tmpprs-273.15)/(Tmpprs-29.65)))'
'define q=RHprs*(0.62197*es/(prs-es))/100.'
'define e=prs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(Tmpprs)-log(e)-4.805)'
'define theta=Tmpprs*pow((1000./prs),(0.2854*(1.0-0.28*q)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'
'set gxout contour'
'set ccolor 1'
'd eqt'

'q time'date=subwrd(result,3)
'draw title thetase(huanghua)at 'date''
'print'
'c'
tt=tt+1
endwhile
'disable print'
'reinit'

用set xaxis命令设置有坐标数据范围以后出图和实际数据不一致,只是坐标改变了,图并没有改变,这是为什么?  
QQ截图20160508213453.png
QQ截图20160508213553.png
QQ截图20160508213528.png
密码修改失败请联系微信:mofangbao
发表于 2016-5-9 08:53:37 | 显示全部楼层
楼主能发过来你的数据及描述文件吗,我想调试一下
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2016-5-9 10:18:59 | 显示全部楼层
风之精灵 发表于 2016-5-9 08:53
楼主能发过来你的数据及描述文件吗,我想调试一下

怎么发给你?我看回复发不了附件。
数据就是FNL1*1
CTL也是对应的,论坛里好多
密码修改失败请联系微信:mofangbao
发表于 2016-10-17 14:04:25 | 显示全部楼层
请问楼主解决问题了吗?我也遇到了同样的问题,求解答
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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