forrrt1：severe<408>数组超限怎么整呀

xiao. · 发表于 2021-1-13 10:50:35

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

01program co2h20n2
02implicit none
03integer,parameter::M=2966297,N=2,L=1227853
04real ::cob(3000001,2)
05real::co2(L,N)
06real::h2o(M,N)
07real ::k
08integer::i,j,q,p,count
09 
10 
11open(10,file='T_1800_P_25.0_CO2Mol_0.1_N2mol_0.9.data')
12 do i=1,1227852
13 
14    read(10,*)co2(i,:)
15 
16 end do
17 close(10)
18 
19 
20 open(11,file='T_1800_P_25.0_H2OMol_0.1_N2mol_0.9.data')
21 do i=1,2966296
22 
23    read(11,*)h2o(i,:)
24    
25 end do
26 close(11)
27 
28 
29 q=1
30 count=1
31 k=1000
32 do k=0,30000,0.01
33     do p=count,1227852
34      if (k-co2(p,1)<0)then
35             cycle
36      else if (k-co2(p,1)>0.01)then
37      
38     else
39        cob(q,2)=cob(q,2)+co2(p,2)  
40        count=count+1
41      end if
42     end do
43     cob(q,1)=k
44    q=q+1
45     
46  end do 
47 
48 
49 
50 
51  q=1
52   count=1
53 do k=0,30000,0.01
54  do p=count,2966296
55     if (k-co2(p,1)<0)then
56    cycle
57       else if (k-co2(p,1)>0.01)then
58       else
59        cob(q,2)=cob(q,2)+co2(p,2)
60             count=count+1
61 
62       end if
63  end do
64 
65  q=q+1
66 end do
67 
68 
69 open(12,file='hecheng.txt')
70  do i=1,3000000
71 
72    write (12,*)cob(i,:)
73 
74  end do
75  close(12)
76  end 

fcode · 发表于 2021-1-13 10:59:54

错误提示仔细看看，英文翻译一下。
第57行，数组越界

co2 大小是 1227853，但p循环到 2966296

此外，既然你定义了 L 和 M，为什么代码里还要写一大堆 2966296 呢？这样你改这些数字多累啊？？为何不写成 do p = count , M-1 呢？

do k=0,30000,0.01
  do p=count,2966296
   if (k-co2(p,1)<0)then
cycle
   else if (k-co2(p,1)>0.01)then
   else
      cob(q,2)=cob(q,2)+co2(p,2)
         count=count+1
   end if
  end do

[求助] forrrt1：severe<408>数组超限怎么整呀

美女勋章

热心勋章

星光勋章

新人勋章

贡献勋章

管理勋章

帅哥勋章

爱心勋章

规矩勋章

元老勋章

水王勋章