用include来读取一个数组的时候碰到的问题
在debug一个程序的时候发现了这样的一个bug：有个头文件output_step.h的内容大致如下：
!====2D output
nt2d(1)=1 ; nt2d(2)=243
! nt2d(1)=3286 ; nt2d(2)=3650
!====time-series output(nts分記入,stepは重複しないようにする)
nnts(1,1)= 1 ; nnts(1,2)=243
nnts(2,1)=366 ; nnts(2,2)=730
nnts(3,1)=731 ; nnts(3,2)= 1095
nnts(4,1)= 1096 ; nnts(4,2)= 1460
nnts(5,1)= 1461 ; nnts(5,2)= 1825
nnts(6,1)= 1826 ; nnts(6,2)= 2190
nnts(7,1)= 2191 ; nnts(7,2)= 2555
nnts(8,1)= 2556 ; nnts(8,2)= 2920
nnts(9,1)= 2921 ; nnts(9,2)= 3285
nnts(10,1)= 3286 ; nnts(10,2)= 3650
主程序中用include来读取一个头文件，大概内容如下
implicit none
(省略部分)
include "output_step.h"
!====debug==
write(*,*)nnts(1,1)
write(*,*)nnts(1,2)
write(*,*)nnts(2,1)
write(*,*)nnts(2,2)
（省略部分）
end
电脑是64位，windows7，编译器是Microsoft Visual Studio Fortran2008
用ifort 编译，结果是：
1
366
366
731
也就是说读入了数据以后马上就输出，没有经过任何数值计算，结果竟然错了。一开始以为Fortran数据储存的方式是先循环前面的，然后出现了这个bug，可是看这组数据貌似不是那么简单了。
我瞬间就凌乱了，难道我电脑也是醉了吗？
这让我以后还敢用include吗？
求各位大神指点到底是什么原因``
传所有代码吧 由于write 后面的部分没有参与所以省略不计
程序如下
program xx
implicit none
include "param.h"
include "TECIO_mod.f90"
real(8) :: X(ngi,ngj),Y(ngi,ngj),T(ni,nj),SS(ni,nj),TURB(ni,nj),U(ni,nj),V(ni,nj),yvd(ngj),dx(ni)
real(8),allocatable :: XT(:,:),YT(:,:),TT(:,:),SST(:,:),TURBT(:,:),UT(:,:),VT(:,:)
real(8) :: dys,ytmp1,ytmp2
integer :: i,j,ii,jj,tloop,io
integer :: dts,dto,tts1,tts2
integer :: iyr,imon,iday,ihour,niu,njs
integer :: js(ni)
!====Variables for Tecplot-Binary Output
integer(4) :: ti,tiii,imax,jmax,kmax,debug,disdouble,visdouble
integer(4) :: strandid,parentzn
integer(4) :: pvl(7),valloc(7),sharingzone(7)
real(8) :: soltime1,soltime2
character(len=1) :: nulchar
character(len=15) :: timchar1,timchar2
include "output_step.h"
!====debug==
write(*,*)nnts(1,1)
write(*,*)nnts(1,2)
write(*,*)nnts(2,1)
write(*,*)nnts(2,2)
end
output_step.h 文件的内容已经完整传上来了。可以看到前面的代码都只是定义变量的部分，跟这部分应该没有关系。
本帖最后由 楚香饭 于 2014-10-17 17:44 编辑
根据你已有的代码来看，不会出错，实际上在我这里也是正常的。
建议你无省略的上传所有代码（可打包，上传附件）
页:
[1]