绘制直方图时总提示无效变量，跪求大神帮助~~~~~

怪力金鱼姬

绘制直方图总是提示定义好的变量是无效变量，Grads小白试了好多办法都没有解决，跪求大神帮助~
以下是fortran程序：
program main
implicit none
integer i,j
integer,parameter::it=10,tt=12
integer yy(it),mm(tt)
open(1,file='f:\lunwen\first\zl\year.txt')
open(2,file='f:\lunwen\first\zl\mon.txt')
open(11,file='f:\lunwen\first\cx\cishu\y.grd',form='binary')
open(12,file='f:\lunwen\first\cx\cishu\m.grd',form='binary')
read(1,*)(yy(i),i=1,it)
do i=1,it
write(11)yy(i)
print *,yy(i)
enddo
read(2,*)(mm(j),j=1,tt)
do j=1,tt
write(12)mm(j)
print *,mm(j)
enddo
close(1)
close(2)
close(11)
close(12)
end


ctl文件：
dset F:\lunwen\first\cx\cishu\y.grd
undef -9.99E+33
xdef 1 linear 1 1
ydef 1 linear 1 1
zdef 1 levels 1
tdef 10 linear JAN2010 1yr
vars 1
yy 1 99
endvars
；


gs文件：
'reinit'
'open f:\lunwen\first\cx\cishu\y.ctl'
'set parea 1 9 1 5'
'set x 1'
'set y 1'
'set z 1'
'set grads off'
'set gxout bar'
'set vrange 0 8'
'set xlint 1'
'set ylint 1'
'set bargap 60'
'set baropts filled'
'set barbase bottom'
'set grid off'
'set t 0 10'
'set ccolor 3'
'd yy'
'printim F:/lunwen/first/cx/cishu/y.png white'
;

十分感谢！！

伽蓝鸟

set t 1 10？

怪力金鱼姬

伽蓝鸟 发表于 2021-5-30 17:11
set t 1 10？

谢谢您的回复，不过好像不是这个原因。时间t从0开始，画图的时候X轴的起点是0；时间t从1开始，X轴的起点是2010年。

river

上GrADS报错的信息，截图

怪力金鱼姬

river 发表于 2021-5-30 21:46
上GrADS报错的信息，截图

您好，现在没有提示错误了，就是只能画出来坐标轴，变量画不出来。

                               
登录/注册后可看大图

怪力金鱼姬

river 发表于 2021-5-30 21:46
上GrADS报错的信息，截图

C:\Users\dell\Desktop\1.png

苏拉苏拉

有txt文件直接拿excel画多方便，干嘛搞成二进制啊

怪力金鱼姬

苏拉苏拉 发表于 2021-5-31 14:41
有txt文件直接拿excel画多方便，干嘛搞成二进制啊

感谢您的回复，个人需求。

付亚男

你写出来的二进制文件不对，这个文件里的值应该全部集中于0附近，就是比如1e-10这种极小的值，你把set vrange去掉，就能看见错误的结果了。
在fortran代码里，你把变量声明为real就可以解决这个问题了。

怪力金鱼姬

付亚男 发表于 2021-5-31 18:56
你写出来的二进制文件不对，这个文件里的值应该全部集中于0附近，就是比如1e-10这种极小的值，你把set vran ...

哇！是这个原因，终于画出来正确的图了~真的非常感谢您的帮助