fortran 读取micaps数值预报结果出错？

zdf2075

欧洲中心数值预报结果：
diamond 4 2015012820-0102_EC_2米气温预报
2015 01 28 20 78 0
0.125000 -0.125000 103.000000 105.000000 32.000000 30.000000 17 17 1.000000 -100.000000 88.000000 1.000000 0.000000
-8.61 -5.83 -3.86 -4.61 -6.98 -7.33 -8.17 -8.83 -7.20 -5.39 -4.89 -4.51 -3.58 -2.58 -1.58 -0.67  0.17
-11.80 -8.83 -6.05 -5.45 -6.17 -5.73 -6.33 -7.36 -6.98 -5.05 -3.70 -2.92 -2.42 -1.23  0.11  1.11  1.61
-12.17 -9.95 -7.73 -7.76 -7.95 -6.58 -5.86 -6.01 -6.58 -5.64 -3.67 -1.76 -0.39  0.77  1.74  2.05  1.83
-9.92 -8.42 -6.80 -6.92 -6.80 -5.55 -5.64 -6.58 -6.51 -5.11 -2.42  0.05  1.42  2.02  2.33  2.27  1.89
-7.11 -6.17 -5.51 -5.55 -5.01 -4.67 -6.20 -7.70 -5.70 -3.17 -0.39  1.58  2.27  2.45  2.45  2.36  2.17
-7.36 -6.39 -6.17 -5.95 -5.45 -5.98 -6.95 -6.86 -3.51 -1.01  0.86  1.99  2.36  2.64  2.67  2.52  2.27
-11.48 -9.17 -7.36 -5.89 -5.64 -5.86 -5.23 -3.55 -1.23  0.39  1.49  2.11  2.36  2.77  2.74  2.42  1.99
-11.33 -8.55 -6.95 -6.42 -4.86 -3.36 -2.11 -1.01  0.24  1.27  2.20  2.55  2.20  2.27  2.45  2.30  1.83
-8.55 -6.11 -5.80 -6.20 -2.58 -0.08  0.95  0.95  0.92  1.89  2.20  2.02  1.61  1.49  2.11  2.20  2.08
-9.67 -7.01 -5.11 -3.58  0.02  1.33  1.49  1.27  1.61  1.67  1.39  0.86  0.49  1.80  2.30  2.27  2.14
-7.51 -5.83 -4.05 -1.76  0.42  1.36  1.74  1.70  1.36  1.05  0.58  0.36  0.64  1.80  2.30  2.45  2.55
-6.26 -4.11 -1.86  0.42  1.49  1.95  2.17  2.17  1.74  1.14  0.70  0.77  1.36  2.20  2.45  2.74  3.08
-5.08 -2.33  0.33  2.20  2.45  2.58  2.74  2.80  2.52  1.70  1.36  1.61  2.30  2.86  2.92  3.11  3.36
-3.01 -0.48  1.86  3.11  2.83  2.95  3.20  3.30  2.99  2.27  2.17  2.55  3.11  3.49  3.55  3.55  3.55
-0.67  0.99  2.49  3.27  2.89  3.24  3.49  3.49  3.05  2.67  2.89  3.30  3.61  3.83  3.89  3.83  3.74
0.77  1.58  2.42  3.24  3.05  3.52  3.77  3.64  3.17  3.08  3.42  3.80  3.92  3.95  3.92  3.83  3.80
0.83  1.77  2.86  3.58  3.36  3.80  4.05  3.83  3.24  3.30  3.70  4.08  4.20  4.27  4.08  3.95  3.86

fortran程序如下，为了显示哪个地方错误我把程序是一个一个数据读取的，显示只读了第一列的数据 而没有读取第二列，所以出错，请大神看看怎么了？
parameter hs=17,ls=17
character*5000  a
real zl(hs,ls)
open(1,file='1.txt')
read(1,'(a)')a
print*,trim(adjustl(a))
read(1,'(a)')a
print*,trim(adjustl(a))
read(1,'(a)')a
print*,trim(adjustl(a))
read(1,*)b
print*,b
read(1,*)b
print*,b
read(1,*)b
print*,b
read(1,*)b
print*,b
read(1,*)b
print*,b
read(1,*)b
print*,b
read(1,*)b
print*,b
read(1,*)b
print*,b
read(1,*)b
print*,b
read(1,*)b
print*,b
read(1,*)b
print*,b
read(1,*)b
print*,b
read(1,*)b
print*,b
read(1,*)b
print*,b
read(1,*)b
print*,b
read(1,*)b
print*,b
read(1,*)b
print*,b
read(1,*)b
print*,b
close(1)
end

龍_彳亍_天_下

b又没定义数组，当然读取第一列了。。。

zdf2075

龍_彳亍_天_下 发表于 2015-2-2 08:28
b又没定义数组，当然读取第一列了。。。

即便没有定义数组不是也应该往后面读吗？

lqouc

zdf2075 发表于 2015-2-2 10:31
即便没有定义数组不是也应该往后面读吗？

事实不是已经告诉你不会往后面读了么？
而且已经声明了数组real zl(hs,ls)，干嘛还用个b。

zdf2075

lqouc 发表于 2015-2-2 11:06
事实不是已经告诉你不会往后面读了么？
而且已经声明了数组real zl(hs,ls)，干嘛还用个b。

那是因为事实是 我先用的 zl(hs,ls)，不行，我才用b看一下到底是哪里出错了

srflyh

是要声明数组的哦

zdf2075

srflyh 发表于 2015-2-3 08:31
是要声明数组的哦

定义了 还不行吗？  应该怎么声明呢

wysyhd

zdf2075 发表于 2015-2-3 09:59
定义了 还不行吗？  应该怎么声明呢

貌似你得写一个循环读取数据，你这样一行一行read和print只读第一行第一个

1530718856

那个b是啥东西，