Fortran Coder

查看: 44|回复: 5

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

[复制链接]

6

帖子

2

主题

0

精华

入门

F 币
40 元
贡献
23 点
发表于 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] 纯文本查看 复制代码
program main
    implicit none
    character(*),parameter :: infile='Depth.01',outfile='Depth2.01'
    integer :: i,j
    real,dimension(521,1440):: m
    real,dimension(521,1440):: b

    OPEN(UNIT=20,FILE=infile)
    OPEN(UNIT=30,FILE=outfile,status='replace',action='write')
!读取文件
    do i=1,521
        read(20,'(1440(I6))',end=99) (m(i,j),j=1,1440)
    enddo

!赋值
do i=1,521
       do j=4,1437
        if(m(i,j)==-999) then 
            b(i,j)=1
       
           else if(m(i,j)/=-999.and.m(i,j)==-999) then         !999陆地
               b(i,j)=2                                                        !2 海岸线
               b(i,j-1)=2                                                      !1 海洋
               b(i,j-2)=2                     
               b(i,j-3)=2
           else if(m(i,j)==-999 .and. m(i,j+1)/=-999) then
               b(i,j)=2
               b(i,j+1)=2
               b(i,j+2)=2
               b(i,j+3)=2
           else if (m(i,j)<200) then
               b(i,j)=2
        else 
         b(i,j)=m(i,j)
        endif
       enddo
   enddo

!写入新文件
  do i=1,521
        write(30,'(1440(I6))') (b(i,j),j=1,1440)
        end do

99 stop
end program main
回复

使用道具 举报

1292

帖子

12

主题

5

精华

论坛跑堂

Fcode跑堂伙计

F 币
56 元
贡献
159 点

新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

发表于 2018-10-7 22:47:48 | 显示全部楼层
是格式问题呀。
b 是浮点数,不能用 I 编辑符。(因为 I 编辑符只针对整型)
其他意见,见下面代码的注释内容。

[Fortran] 纯文本查看 复制代码
program main
  implicit none
  Integer , parameter :: NX = 521 , NY = 1440 , SLAND = -999!//定义常量,以便容易修改
  character(*),parameter :: infile='Depth.01',outfile='Depth2.01'
  integer :: i,j
  real,dimension(NX,NY):: m, b

  OPEN(UNIT=20,FILE=infile)
  OPEN(UNIT=30,FILE=outfile,status='replace',action='write')
  !读取文件
  do i=1,NX
    read(20,*,end=99) m(i,:)!//数组片段,简化代码
  enddo
  !赋值
  do i=1,NX
    do j=1+3,NY-3
      if(m(i,j)==SLAND) then
        b(i,j)=1
      else if(m(i,j)/=SLAND.and.m(i,j+1)==SLAND) then!//此处是否m(i,j+1)??
        b(i,j-3:j)=2!//数组片段,简化代码
      else if(m(i,j)==SLAND .and. m(i,j+1)/=SLAND) then
        b(i,j:j+3)=2!//数组片段,简化代码
      else if (m(i,j)<200) then
        b(i,j)=2
      else
        b(i,j)=m(i,j)
      endif
    enddo
  enddo
  !写入新文件
  do i=1,NX
    write(30,'(*(f13.4))') b(i,:)!// *(f13.4) 表示有多少个算多少个。
  end do
99 close(20)!//良好的关闭文件习惯
   close(30)
end program main

6

帖子

2

主题

0

精华

入门

F 币
40 元
贡献
23 点
 楼主| 发表于 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。这个要重新再写两个循环吗,还是有更优化的代码呢?

1292

帖子

12

主题

5

精华

论坛跑堂

Fcode跑堂伙计

F 币
56 元
贡献
159 点

新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

发表于 2018-10-8 14:45:43 | 显示全部楼层
write(30,'(*(f13.4))') b(i,1+3:NY-3)

6

帖子

2

主题

0

精华

入门

F 币
40 元
贡献
23 点
 楼主| 发表于 2018-10-8 15:27:07 | 显示全部楼层
fcode 发表于 2018-10-8 14:45
write(30,'(*(f13.4))') b(i,1+3:NY-3)

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


[Fortran] 纯文本查看 复制代码
program main
  implicit none
  Integer , parameter :: NX = 521 , NY = 1440 , SLAND = -999
  character(*),parameter :: infile='Depth.01',outfile='Depth2.01'
  integer :: i,j
  integer,dimension(NX,NY):: m, b

  OPEN(UNIT=20,FILE=infile)
  OPEN(UNIT=30,FILE=outfile,status='replace',action='write')
  !读取文件
  do i=1,NX
    read(20,'(*(I6))',end=99) m(i,:) !//数组片段,简化代码
  enddo

  !赋值
  do i=1,NX
    do j=1,4
    I if(m(i,j)==SLAND) then
        b(i,j)=m(i,j)
      else if(m(i,1)/=SLAND.and.m(i,2)==SLAND) then 
        b(i,1)=1 
      else if(m(i,2)/=SLAND.and.m(i,3)==SLAND) then 
        b(i,1:2)=1
      else if(m(i,3)/=SLAND.and.m(i,4)==SLAND) then 
        b(i,1:3)=1
      else if(m(i,1)==SLAND.and. m(i,2)/=SLAND) then
        b(i,2:3)=1 
      else if(m(i,2)==SLAND.and. m(i,3)/=SLAND) then 
        b(i,3)=1
      else if (m(i,j)<200) then
        b(i,j)=1
      else
        b(i,j)=555
      endif 
    enddo
   enddo

  do i=1,NX
    do j=1+3,NY-3
      if(m(i,j)==SLAND) then
        b(i,j)=m(i,j)
      else if(m(i,j)/=SLAND.and.m(i,j+1)==SLAND) then !//此处是否m(i,j+1)??
        b(i,j-3:j)=1 !//数组片段,简化代码
      else if(m(i,j)==SLAND .and. m(i,j+1)/=SLAND) then
        b(i,j:j+3)=1 !//数组片段,简化代码
      else if (m(i,j)<200) then
        b(i,j)=1
      else
        b(i,j)=555
      endif
    enddo
  enddo
  
    do i=1,NX   
     do j=NY-3,NY
    I if(m(i,j)==SLAND) then
        b(i,j)=m(i,j)
      else if(m(i,NY-2)/=SLAND.and.m(i,NY-1)==SLAND) then 
        b(i,NY-2)=1 
      else if(m(i,NY-1)/=SLAND.and.m(i,NY)==SLAND) then 
        b(i,NY-2:NY-1)=1
      else if(m(i,NY-3)==SLAND.and. m(i,NY-3)/=SLAND) then
        b(i,NY-2:NY)=1 
      else if(m(i,NY-2)==SLAND.and. m(i,NY-1)/=SLAND) then 
        b(i,NY-1:NY)=1
      else if(m(i,NY-1)==SLAND.and. m(i,NY)/=SLAND) then 
        b(i,NY)=1
      else if (m(i,j)<200) then
        b(i,j)=1
      else
        b(i,j)=555
      endif
     enddo
  enddo

 

  !写入新文件
  do i=1,NX
    write(30,'(*(I6))') b(i,:) 
  end do
99 close(20)
    close(30)
end program main

1292

帖子

12

主题

5

精华

论坛跑堂

Fcode跑堂伙计

F 币
56 元
贡献
159 点

新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

发表于 2018-10-8 16:36:49 | 显示全部楼层
1. 出错请给错误提示。
2. 让别人调试,请提供输入数据文件。
3. 你说“能不能保留呢?”,我不知道保留成什么?涉及到程序逻辑问题,请描述清楚。
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

QQ|捐赠本站|Archiver|关于我们 About Us|群聊|Fcode

GMT+8, 2018-10-23 23:57

Powered by Discuz! X3.2

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表