爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3241|回复: 8

[脚本编辑] 求助,grd有数据但grads出图显示值为0

[复制链接]

新浪微博达人勋

发表于 2016-5-13 19:48:04 | 显示全部楼层 |阅读模式

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

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

x
fortran第一个程序处理逐日四时次风场资料后得出风切变点的数据(有切变1,无切变0),得出结果再用第二个程序将切变点按月/年叠加,和grd一起生成的txt有数据,但是grads画出图显示value=0,两个程序生成的结果都这样,我认为不是grd的问题,因为如果第一个程序生成的grd有问题,第二个程序就转不起来啊,求助/(ㄒoㄒ)/~~我尝试画的是一年的叠加图
ctl
DSET C:\lunwen\700\data\2012.year.h.grd
TITLE 700hPa.h
undef 1111
xdef 15 linear 95 1
ydef 8 linear 35 1
zdef 1 levels 700
tdef 1 linear 01jan2012 1yr
vars 1
hqbd    0   9999  heng qie bian dian
endvars

gs就是最简单的
'reinit'
'open c:\lunwen\700\qiebiandian\h.y.ctl'
'set lon 95 110'
'set lat 35 43'
'set grads off'
'set lev 700'
'set t 1'
'set gxout contour'
'd hqbd'
'enable print c:\lunwen\700\qiebiandian\x.gmf'
'print'
'disable print'
;

x.png






密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-5-13 20:40:10 | 显示全部楼层
求指点= =告诉我一点可能的原因啊
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-5-13 21:23:37 | 显示全部楼层
我不是很懂,你把ctl里的9999改为999试试?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-5-13 21:37:37 | 显示全部楼层
小小书童 发表于 2016-5-13 21:23
我不是很懂,你把ctl里的9999改为999试试?

改过了 然后缺测值各种也都试过了 都没用。。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-5-14 08:31:46 | 显示全部楼层
hqbd,这个不是坐标吗?怎么用CONTOUR
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-5-14 08:49:21 | 显示全部楼层
yuyumage 发表于 2016-5-14 08:31
hqbd,这个不是坐标吗?怎么用CONTOUR

我fortran判断的时候存的是h(i,j,t)
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-5-14 09:28:43 | 显示全部楼层
天坑啊 grads2.0画出的图 文件什么的全没改 但昨天和今天还不一样 最后用grads1.9就好了!!!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-5-17 19:10:56 | 显示全部楼层
fortran生成同时生成grd和txt,txt有数据,但是grd无法画出图,用fortran将txt读入grd,逐年数据可以,逐月不可以,最后是用matlab解决的。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-5-17 22:16:09 | 显示全部楼层
得要看下fortran程序~
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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