已解决 二维数组读取时，外循环不运行

wutongyu

在读取一个（71,73）的二维数组过程中，程序只执行内循环，外循环不执行，自己写的代码如图，不知道哪儿错了，求指导

wc619lf

a，b没赋初值，读和写的循环最好分开写，我是这样想的，你试试

言深深

楼主如何判断外层循环没有运作？
write的时候format和定义需要相匹配，a默认是实型的，写成I型的不妥，再者下面的3X···建议修改，另外，status语句行去掉···

wutongyu

言深深 发表于 2012-12-29 11:02

                               
登录/注册后可看大图

楼主如何判断外层循环没有运作？
write的时候format和定义需要相匹配，a默认是实型的，写成I型的不妥，再者 ...

额，对，从新看了下数据，数据格式是这样的
33    32   37    33    33    33  .....
33    33   33    33    33    38   ....
20    50   33    33    33    74 ......
总共71行，73列，我自己写的程序只能读出来第一列，比如上面的33,33,20.。。。
这是写的完整代码，不知哪儿出的问题，求指导一下
program ex19951
implicit none
integer,parameter::inputfileid=10,outputfileid=11
integer m,n,t,i,j,k
character(len=80)::inputfile,outputfile
parameter(m=71,n=73)
integer x(m,n)
integer a
real b
a=-180
b=90.0
inputfile="0010.txt"
outputfile="1.txt"

open(10,file=inputfile,status='old')
open(11,file=outputfile,status="replace")

do t=1,71
do j=1,73
read(10,*)(x(t,j))
b=b-2.5
a=a+5
write(11,"(I8,3XF5.1,I6,I6,I6)")a,b,(x(1,j)),j,t
end do
end do
stop
end program  

言深深

wutongyu 发表于 2012-12-29 11:24

                               
登录/注册后可看大图

额，对，从新看了下数据，数据格式是这样的
33    32   37    33    33    33  .....
33    33   33  ...

如果是这样就比较容易理解了···
read的时候该用隐循环，具体如下
do i=1,71
read(1,*)(x(i,j),j=1,73)
enddo
语句的意思就是读取71行，每行读取73个···

你原来的程序应该还会报错的，错误应该是end of read

wutongyu

言深深 发表于 2012-12-29 11:53

                               
登录/注册后可看大图

如果是这样就比较容易理解了···
read的时候该用隐循环，具体如下
do i=1,71

嗯，你猜的太对了，就是运行时DOS窗口弹出这个错误，你说的改正方法可以，不过不能逐个输出啊，我想在每个输出的数前面加上经纬度，a,b分别表示的就是经纬度。

言深深

wutongyu 发表于 2012-12-29 14:20

                               
登录/注册后可看大图

嗯，你猜的太对了，就是运行时DOS窗口弹出这个错误，你说的改正方法可以，不过不能逐个输出啊，我想在每个 ...

可以逐个输出的。可以用i,j*间隔来输出的···输出还是用
do i=1,m
do j=1,n
write(2,"(2f12.2,f12.4)") i*2.5,j*2.5,x(i,j)
enddo
enddo

wutongyu

言深深 发表于 2012-12-29 17:23

                               
登录/注册后可看大图

可以逐个输出的。可以用i,j*间隔来输出的···输出还是用
do i=1,m
do j=1,n

好滴，十分感谢，问题解决啦，哈哈

言深深

wutongyu 发表于 2012-12-29 21:07

                               
登录/注册后可看大图

好滴，十分感谢，问题解决啦，哈哈

建议你查看一下，输出的经纬度有可能对应不上，需要做相应的修改，比如90-j*2.5等等···具体原因就是和之前数据存放的方式有关，这些细节想一下就明白了

wutongyu

言深深 发表于 2012-12-29 21:30

                               
登录/注册后可看大图

建议你查看一下，输出的经纬度有可能对应不上，需要做相应的修改，比如90-j*2.5等等···具体原因就是和 ...

嗯，经纬度的问题已经解决啦