想请问下关于fortran提取数据的问题

Daisy°L

各位大神 大家好 毕业论文有一部分是要求副高指数的
老师给的资料是74项环流特征量资料
老师给了一个程序
parameter(iva=74,imon=12,iyr=100)
dimension hc068(iva,imon,iyr)
open(10,file='hc068',err=200)
     read(10,100) (((hc068(k,i,j),k=1,iva),i=1,imon),j=1,iyr)
100 format(37I5)
200  continue
close(10)
是这样的
然后file那里我加了文件路径
程序最后一行应该加一个end对吧
然后运行结果是这样的

想知道哪里出了问题
（嗯 是这样的 我的fortran基础也不是很好 是后来学grads的时候才知道了点提取数据的程序 如果提取数据的话 是不是应该把提取的数据存入到新的文件中啊？还有这个程序里面有两个地方不太懂 一个是err=200 一个是format（37I5）这个查书也没有查到 嗯 麻烦大家看下了哈）

Daisy°L

哦哦哦 那个err在百度里面查到了 还有就是百度到的是formatted的用法 想问下和format一样么 后面的括号内的数字是什么意思啊

半个瞎子

不一样的，formatted是对于文件属性的描述，意思是十进制的，format（）是读或是写的格式，你这里的format（37I5）是指连续37次读取都是以整数型读取的，这要和你读取的文件的写入格式一致了。

Daisy°L

半个瞎子 发表于 2014-5-21 11:22
不一样的，formatted是对于文件属性的描述，意思是十进制的，format（）是读或是写的格式，你这里的format ...

老师给的资料说明里面确实是说都是整型数来着 那这个程序有什么问题么

半个瞎子

没什么问题啊，你这段程序就是把文件里的数据读到你的数组了，没有任何输出语句，你可以 加入一句： print *，hc068（：，：，1） ，看读取是否正确。

lqouc

上面的同学解释的差不多了，报错是因为你的数据没有你认为的那么长，需要进行修改，可以把你的数据截个图上来。

普鲁斯特

这样改：
read(10,100) (hc068(k,I,j),k=1,37)
read(10,100) (hc068(k,I,j),k=38,74)
自己加I,j的循环，其他不变

Daisy°L

lqouc 发表于 2014-5-21 12:34
上面的同学解释的差不多了，报错是因为你的数据没有你认为的那么长，需要进行修改，可以把你的数据截个图上 ...

                               
登录/注册后可看大图

                               
登录/注册后可看大图

额 回复的时候不知道怎么添加图片 嗯 这个数据只知道是从1951年1月开始的 程序是老师给的 说可以直接提取

Daisy°L

lqouc

Daisy°L 发表于 2014-5-21 15:14

7楼的方法是正确的，修改下有问题可以再发上来，而且后面的年份是没有数据的，你可以不用设置100年那么长。