原来读入的变量值怎么自动改变了啊？？？？

pingnuaa

在第一个do 循环后面和第二个do循环后面设置断点，运行完第一个do前面读入的变量值都是对的，
为什么运行完第二个do之后，前面所有变量值都自动变了呢？太奇怪了，求助高手啊，非常感谢！
源程序截取的一部分如下，输入文件见附件，要读入的数据有点多，麻烦高手了，谢谢

PROGRAM FEM2D
      DIMENSION IJK_ELE(500,3),X(500),Y(500),IJK_U(100,14),
     &          AK(2000,2000),
     &THTA(10),ZKD(10),AL(10),CK(10),TX(10),TY(10),TZ(10)
OPEN(4,FILE='BASIC1.IN')
OPEN(5,FILE='NODE_ANSYS1.IN')
OPEN(6,FILE='ELEMENT_ANSYS1.IN')
READ(4,*) N_ELE,N_NODE,N_BC,N_LAYER,N_GAS,N_GAS2
CALL READ_IN(N_ELE,N_NODE,N_BC,N_BAND,E1,E2,G12,V12,G23,
     &            G31,IJI_ELE,X,Y,IJK_U,N_LAYER,THTA,ZKD,HD,XL,
     &            N_GAS,N_GAS2,AL,CK,TX,TY,TZ)
      CLOSE(4)
CLOSE(5)
CLOSE(6)
  END
c   
c to get the original data in order to model the problem
SUBROUTINE READ_IN(N_ELE,N_NODE,N_BC,N_BAND,E1,E2,G12,V12,G23,
     &           G31,IJK_ELE,X,Y,IJK_U,N_LAYER,THTA,ZKD,HD,XL,
     &           N_GAS,N_GAS2,AL,CK,TX,TY,TZ)
      DIMENSION IJK_ELE(N_ELE,3),X(N_NODE),Y(N_NODE),IJK_U(N_BC,14),
     & NE_ANSYS(N_ELE,14),THTA(N_LAYER),ZKD(N_LAYER+1),
     &    AL(N_GAS2),CK(N_GAS),TX(N_GAS),TY(N_GAS),TZ(N_GAS2)  
REAL ND_ANSYS(N_NODE,3)
READ(4,*)E1,E2,G12,G23,V12,G31,HD,XL
READ(4,*)(THTA(I),I=1,N_LAYER)
READ(4,*)(ZKD(I),I=1,N_LAYER+1)
      READ(4,*)((IJK_U(I,J),J=1,14),I=1,N_BC)
READ(4,*)(AL(I),I=1,N_GAS2)
READ(4,*)(CK(J),J=1,N_GAS)
READ(4,*)(TX(K),K=1,N_GAS)
READ(4,*)(TY(K),K=1,N_GAS)
READ(4,*)(TZ(K),K=1,N_GAS2)
READ(5,*)((ND_ANSYS(I,J),J=1,3),I=1,N_NODE)
READ(6,*)((NE_ANSYS(I,J),J=1,14),I=1,N_ELE)
DO L=1,N_NODE
X(L)=ND_ANSYS(L,2)
Y(L)=ND_ANSYS(L,3)
      ENDDO
      DO MM=1,N_ELE
  DO NN=1,3
       IJK_ELE(MM,NN)=NE_ANSYS(MM,NN)
       ENDDO
ENDDO  
c   执行完上面这句前面所有读入的值就都变了，怎么回事？？？？？？？？？？                                                                     
  N_BAND=0
DO  IE=1,N_ELE
   DO  I=1,3
    DO  J=1,3
       IW=IABS(IJK_ELE(IE,I)-IJK_ELE(IE,J))
       IF(N_BAND.LT.IW)N_BAND=IW
         ENDDO
   ENDDO
ENDDO
      N_BAND=(N_BAND+1)*13
  RETURN
  END

铁道科学

1.     PROGRAM FEM2D
   
2.     DIMENSION IJK_ELE(500,3),X(500),Y(500),IJK_U(100,14),
   
3.   & AK(2000,2000),
   
4.   & THTA(10),ZKD(10),AL(10),CK(10),TX(10),TY(10),TZ(10)
   
5.     OPEN(4,FILE='BASIC1.IN')
   
6.     OPEN(5,FILE='NODE_ANSYS1.IN')
   
7.     OPEN(6,FILE='ELEMENT_ANSYS1.IN')
   
8.     READ(4,*) N_ELE,N_NODE,N_BC,N_LAYER,N_GAS,N_GAS2
   
9.     CALL READ_IN(N_ELE,N_NODE,N_BC,N_BAND,E1,E2,G12,V12,G23,
   
10.   & G31,IJI_ELE,X,Y,IJK_U,N_LAYER,THTA,ZKD,HD,XL,
    
11.   & N_GAS,N_GAS2,AL,CK,TX,TY,TZ)
    

复制代码

你程序似乎有点问题，在call那个里面，第十行第二个参数IJI_ELE，这个应该是IJk_ELE吧，你仔细查一下。

pingnuaa

怎么没人解答啊，大侠们快快现身啊，拜托了

铁道科学

1. 
   
2. DO MM=1,N_ELE
   
3. DO NN=1,3
   
4.    IJK_ELE(MM,NN)=NE_ANSYS(MM,NN)
   
5. ENDDO
   
6. ENDDO
   

复制代码

你的意思是执行完这一句后，本来这句没有赋值的变量的值也发生了变化？

pingnuaa

3# 铁道科学

就是执行完这句，之前赋的所有变量值都变得乱七八糟，执行以下就会发现了，附件里有输入文件。就这么看着没什么问题吧？

pingnuaa

3# 铁道科学
麻烦帮忙看看怎么回事好吗，纠结啊

pingnuaa

有没有人帮帮忙啊！！！！！！！！！！！！！！！！！！

pingnuaa

实在是找不出问题，很诡异，大侠快现身啊

铁道科学

稍等片刻，我先运行一下完后给你消息

pingnuaa

9# 铁道科学
真的是啊，太马虎了，看了好久都找不到问题，整个程序很大的，这只是一部分，这么粗心真愁人……。非常非常感谢！！！！！！

铁道科学

10# pingnuaa
建议你在编程序的时候加入implicit none，这样可以避免此类问题的。

pingnuaa

11# 铁道科学
非常感谢！

pingnuaa

2# 铁道科学
大侠再麻烦问一下，如果在主程序里CALL  READ_IN的后面再随便加一句其他的语句比如输出IJK_ELE，为什么IJK_ELE数组里的值还是会变呢？谢谢啊！

铁道科学

建议在程序开头加上implicint none，我怀疑你程序还有一些其他错误。