求大神帮忙看看程序哪儿出错了~~用fortran编程将单列数据按要求变为矩阵

竹外桃

打算将50年的列数据变为日对齐的矩阵。（为百分法求阈值做准备）
fortran程序如下：
program main
real a(366,50),b(18300)
integer i,j,k,n
open(1,file = 'F:\in.dat')
open(2,file = 'F:\out.dat')
read(1,*) (b(i),i=1,18262)
k=0
j=0
do n=1961,2010
j=j+1
if(mod(n,4)>0) then
do i=1,365
  k=k+1
a(i,j)=b(k)
end do
else
do i=1,366
k=k+1
a(i,j)=b(k)
end do
end if
end do
do i=1,366
write(2,'(1x,366i6)'）(a(i,j),j=1,50)
end do
end
程序能够运行，不过结果不是预想的。我想将年作为列，日作为行，把单列数据变成366*50的行列式，运行结果却是366*50的星号。
谢谢大家~~~

lqouc

结果怎么预想不一样啊？好歹描述一下啊。花10个金币把数据下载下来还要替你算一遍，一般人是不会这么做的吧。

竹外桃

哦，其实数据不看也行。
50年的温度数据，我想把以年作为列，日为行，就应该是366*50或者365*50的行列式。但是结果是3660*5的行列式。
谢谢啊~

竹外桃

lqouc 发表于 2013-4-1 21:51

                               
登录/注册后可看大图

结果怎么预想不一样啊？好歹描述一下啊。花10个金币把数据下载下来还要替你算一遍，一般人是不会这么做的吧 ...

谢谢提醒~~，我还不知道下这个附件要金币呢~~~

lqouc

竹外桃 发表于 2013-4-1 21:58

                               
登录/注册后可看大图

哦，其实数据不看也行。
50年的温度数据，我想把以年作为列，日为行，就应该是366*50或者365*50的行列式。 ...

那个是因为文本文件的自动换行问题，其实我看你的矩阵部分是没有问题的，得到的数组也是对的。两个办法解决，一个是在write的时候加个换行控制，另一个是输出为二进制。

竹外桃

lqouc 发表于 2013-4-1 22:12

                               
登录/注册后可看大图

那个是因为文本文件的自动换行问题，其实我看你的矩阵部分是没有问题的，得到的数组也是对的。两个办法 ...

好的，我试试看~~谢了哦

kongfeng0824

请问楼主，你的百分法程序是如何处理的呢？很想学习一下

qxtlyf

用matlab导入数据，在进行数组转换，就ok了

mofangbao

竹外桃 发表于 2013-4-1 21:58

                               
登录/注册后可看大图

哦，其实数据不看也行。
50年的温度数据，我想把以年作为列，日为行，就应该是366*50或者365*50的行列式。 ...

下次请在提问专区发帖

尽头的尽头

楼主可以直接把代码贴上来，这样大家方便给你看，不用耗费金币