Fortran Coder

查看: 12760|回复: 1

[求助] 写入冲突问题

[复制链接]

6

帖子

4

主题

0

精华

入门

F 币
48 元
贡献
26 点
发表于 2019-10-16 18:47:25 | 显示全部楼层 |阅读模式
异常.jpg
异常如上图所示,编译没有错误和警告,百度上实在没找到解决办法。。故来此请教各位这个问题是由什么引起的,如何解决。
代码如下:
[Fortran] 纯文本查看 复制代码
!本程序主要用于网格模型的转化,2019年10月16日
program main
implicit none
integer, parameter :: FN=200,MN=200
integer counter1,counter2,counter3,dn,i
integer,allocatable :: NE(:)
integer,allocatable :: NSEQ(:)
integer,allocatable :: NADD(:)
integer,allocatable :: MA2(:)
character(len=3),allocatable :: EL(:),MA1(:)
character(len=10) title1,title2,variables,zone
real(kind=8) d,dm
real(kind=8),allocatable :: VOLX(:),AHTX(:),PMX(:),X(:),Y(:),Z(:),Xf(:),Yf(:),Zf(:),f(:)
real(kind=4),allocatable :: M(:,:)
!打开地质模型生成的岩相数据文件
open(1,file='Zone.txt',status='unknown')
!打开正演模型网格文件
open(2,file='MESH.txt',status='unknown')
read(1,*) title1
write(*,*) title1
read(1,*) variables
read(1,*) zone
read(2,*) title2
write(*,*) title2
dm=10000000.0
i=0
!开始循环读取地质模型岩相数据文件
Do counter1=1,FN
  write(*,*) "good luck"
  i=i+1
  read(1,10) Xf(i),Yf(i),Zf(i),f(i)
10 format(E10.4,E10.4,E10.4,I2)
  write(*,*) Xf(i),Yf(i),Zf(i),f(i)
end do
!开始循环读取网格模型数据文件并修改
Do counter2=1,MN
  read(2,20) EL(counter2),NE(counter2),NSEQ(counter2),NADD(counter2),MA1(counter2),MA2(counter2),VOLX(counter2),AHTX(counter2),PMX(counter2),X(counter2),Y(counter2),Z(counter2)
20 format(A3,I2,I5,I5,A3,I2,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4)  
!开始计算最短距离
   Do counter1=1,MN
     d=(X(counter2)-Xf(counter1))**2+(X(counter2)-Yf(counter1))**2+(Z(counter2)-Zf(counter1))**2
  if(dm>=d)then
    dm=d
    dn=counter1
  end if 
   end do
   MA2(counter2)=f(dn)
end do
!替换正演网格模型
write(2,*) title2
Do counter2=1,MN,1
  write(2,30) EL(counter2),NE(counter2),NSEQ(counter2),NADD(counter2),MA1(counter2),MA2(counter2),VOLX(counter2),AHTX(counter2),PMX(counter2),X(counter2),Y(counter2),Z(counter2)
  30 format(A3,I2,I5,I5,A3,I2,E10.4,E10.4,E10.4,E10.4,E10.4,E10.4)  
end do
stop
end



Zone.txt

13.76 KB, 下载次数: 1

MESH.txt

30.04 KB, 下载次数: 1

1948

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
1298 元
贡献
547 点

美女勋章热心勋章星光勋章新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

发表于 2019-10-17 13:54:46 | 显示全部楼层
Xf(i),Yf(i),Zf(i),f(i)
这些可分配数组,需要先分配再使用
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-3-29 01:10

Powered by Tencent X3.4

© 2013-2024 Tencent

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