二维数组赋值出错

吃饭叫我 · 发表于 2018-10-7 21:45:56

求助！！！！！！！！！！请大家帮忙看看，很简单的赋值，不知道哪里出错，卡在这里了。

背景：一个规则矩阵的数据文件，矩阵大小为【521,1440】。存储了某地区海洋和陆地的高度信息。其中-999表示陆地，其余数值表示水深，水深用正值表示。如：26，即表示此处水深26米;  -999则表示该地区为陆地。
需求：1.将文件信息读入数组（已完成）
      2。提取海岸线信息，
                  （1）将海岸线的点标注为2，
                  （2）同时，扩大范围，将与海岸线相邻的连续三个点赋值为2，
                  （3）且若水水深<200米，赋值为2
                  （4）其余表示海洋的点赋值为1

      3.将更新后的数组存入到另一个文件（已完成）

问题：只要开始赋值，就出错变成  **********，是格式问题还是赋值出错？难道是编译器问题？？？

代码如下：

[Fortran] 纯文本查看 复制代码

01program main
02    implicit none
03    character(*),parameter :: infile='Depth.01',outfile='Depth2.01'
04    integer :: i,j
05    real,dimension(521,1440):: m
06    real,dimension(521,1440):: b
07 
08    OPEN(UNIT=20,FILE=infile)
09    OPEN(UNIT=30,FILE=outfile,status='replace',action='write')
10！读取文件
11    do i=1,521
12        read(20,'(1440(I6))',end=99) (m(i,j),j=1,1440)
13    enddo
14 
15！赋值
16do i=1,521
17       do j=4,1437
18        if(m(i,j)==-999) then
19            b(i,j)=1
20        
21           else if(m(i,j)/=-999.and.m(i,j)==-999) then         !999陆地
22               b(i,j)=2                                                        !2 海岸线
23               b(i,j-1)=2                                                      !1 海洋
24               b(i,j-2)=2                    
25               b(i,j-3)=2
26           else if(m(i,j)==-999 .and. m(i,j+1)/=-999) then
27               b(i,j)=2
28               b(i,j+1)=2
29               b(i,j+2)=2
30               b(i,j+3)=2
31           else if (m(i,j)<200) then
32               b(i,j)=2
33        else
34         b(i,j)=m(i,j)
35        endif
36       enddo
37   enddo
38 
39！写入新文件
40  do i=1,521
41        write(30,'(1440(I6))') (b(i,j),j=1,1440)
42        end do
43 
4499 stop
45end program main

fcode · 发表于 2018-10-7 22:47:48

是格式问题呀。
b 是浮点数，不能用 I 编辑符。（因为 I 编辑符只针对整型）
其他意见，见下面代码的注释内容。

[Fortran] 纯文本查看 复制代码

01program main
02  implicit none
03  Integer , parameter :: NX = 521 , NY = 1440 , SLAND = -999!//定义常量，以便容易修改
04  character(*),parameter :: infile='Depth.01',outfile='Depth2.01'
05  integer :: i,j
06  real,dimension(NX,NY):: m, b
07 
08  OPEN(UNIT=20,FILE=infile)
09  OPEN(UNIT=30,FILE=outfile,status='replace',action='write')
10  !读取文件
11  do i=1,NX
12    read(20,*,end=99) m(i,:)!//数组片段，简化代码
13  enddo
14  !赋值
15  do i=1,NX
16    do j=1+3,NY-3
17      if(m(i,j)==SLAND) then
18        b(i,j)=1
19      else if(m(i,j)/=SLAND.and.m(i,j+1)==SLAND) then!//此处是否m(i,j+1)??
20        b(i,j-3:j)=2!//数组片段，简化代码
21      else if(m(i,j)==SLAND .and. m(i,j+1)/=SLAND) then
22        b(i,j:j+3)=2!//数组片段，简化代码
23      else if (m(i,j)<200) then
24        b(i,j)=2
25      else
26        b(i,j)=m(i,j)
27      endif
28    enddo
29  enddo
30  !写入新文件
31  do i=1,NX
32    write(30,'(*(f13.4))') b(i,:)!// *(f13.4) 表示有多少个算多少个。
33  end do
3499 close(20)!//良好的关闭文件习惯
35   close(30)
36end program main

吃饭叫我 · 发表于 2018-10-8 14:26:48

fcode 发表于 2018-10-7 22:47
是格式问题呀。
b 是浮点数，不能用 I 编辑符。（因为 I 编辑符只针对整型）
其他意见，见下面代码的注释内 ...

谢谢回复！
用您的代码输出成功了。但是我的数都是整形的，所以我改成了integer，dimention（NX,NY）:: m,b。
格式还是用的 I 编辑符
read(20,'(*(I6))',end=99) m(i,:)
write(30,'(*(I6))') b(i,:)
这样输出就没问题了。
但是现在还有一个问题，那就是每列开头三个和结尾三个数都是0，因为循环的范围是 j= 1+3，NY-3。这个要重新再写两个循环吗，还是有更优化的代码呢？

fcode · 发表于 2018-10-8 14:45:43

write(30,'(*(f13.4))') b(i,1+3:NY-3)

吃饭叫我 · 发表于 2018-10-8 15:27:07

fcode 发表于 2018-10-8 14:45
write(30,'(*(f13.4))') b(i,1+3:NY-3)

谢谢您，这样就少输出了一些数据，能不能保留呢，而且这些点的值也要按要求修改的。我改了一下，又出错了，能不能帮忙看下

[Fortran] 纯文本查看 复制代码

01program main
02  implicit none
03  Integer , parameter :: NX = 521 , NY = 1440 , SLAND = -999
04  character(*),parameter :: infile='Depth.01',outfile='Depth2.01'
05  integer :: i,j
06  integer,dimension(NX,NY):: m, b
07 
08  OPEN(UNIT=20,FILE=infile)
09  OPEN(UNIT=30,FILE=outfile,status='replace',action='write')
10  !读取文件
11  do i=1,NX
12    read(20,'(*(I6))',end=99) m(i,:) !//数组片段，简化代码
13  enddo
14 
15  !赋值
16  do i=1,NX
17    do j=1,4
18    I if(m(i,j)==SLAND) then
19        b(i,j)=m(i,j)
20      else if(m(i,1)/=SLAND.and.m(i,2)==SLAND) then
21        b(i,1)=1
22      else if(m(i,2)/=SLAND.and.m(i,3)==SLAND) then
23        b(i,1:2)=1
24      else if(m(i,3)/=SLAND.and.m(i,4)==SLAND) then
25        b(i,1:3)=1
26      else if(m(i,1)==SLAND.and. m(i,2)/=SLAND) then
27        b(i,2:3)=1
28      else if(m(i,2)==SLAND.and. m(i,3)/=SLAND) then
29        b(i,3)=1
30      else if (m(i,j)<200) then
31        b(i,j)=1
32      else
33        b(i,j)=555
34      endif
35    enddo
36   enddo
37 
38  do i=1,NX
39    do j=1+3,NY-3
40      if(m(i,j)==SLAND) then
41        b(i,j)=m(i,j)
42      else if(m(i,j)/=SLAND.and.m(i,j+1)==SLAND) then !//此处是否m(i,j+1)??
43        b(i,j-3:j)=1 !//数组片段，简化代码
44      else if(m(i,j)==SLAND .and. m(i,j+1)/=SLAND) then
45        b(i,j:j+3)=1 !//数组片段，简化代码
46      else if (m(i,j)<200) then
47        b(i,j)=1
48      else
49        b(i,j)=555
50      endif
51    enddo
52  enddo
53   
54    do i=1,NX   
55     do j=NY-3,NY
56    I if(m(i,j)==SLAND) then
57        b(i,j)=m(i,j)
58      else if(m(i,NY-2)/=SLAND.and.m(i,NY-1)==SLAND) then
59        b(i,NY-2)=1
60      else if(m(i,NY-1)/=SLAND.and.m(i,NY)==SLAND) then
61        b(i,NY-2:NY-1)=1
62      else if(m(i,NY-3)==SLAND.and. m(i,NY-3)/=SLAND) then
63        b(i,NY-2:NY)=1
64      else if(m(i,NY-2)==SLAND.and. m(i,NY-1)/=SLAND) then
65        b(i,NY-1:NY)=1
66      else if(m(i,NY-1)==SLAND.and. m(i,NY)/=SLAND) then
67        b(i,NY)=1
68      else if (m(i,j)<200) then
69        b(i,j)=1
70      else
71        b(i,j)=555
72      endif
73     enddo
74  enddo
75 
76  
77 
78  !写入新文件
79  do i=1,NX
80    write(30,'(*(I6))') b(i,:)
81  end do
8299 close(20)
83    close(30)
84end program main

fcode · 发表于 2018-10-8 16:36:49

1. 出错请给错误提示。
2. 让别人调试，请提供输入数据文件。
3. 你说“能不能保留呢？”，我不知道保留成什么？涉及到程序逻辑问题，请描述清楚。

[输入输出] 二维数组赋值出错

美女勋章

热心勋章

星光勋章

新人勋章

贡献勋章

管理勋章

帅哥勋章

爱心勋章

规矩勋章

元老勋章

水王勋章