fortran循环调用子程序，运行出结果20多次后就出现error,求助大神

ruilier

我的fortran程序为调用子程序（6S模式-一个小模式）计算光学厚度的
数据从文件中读取，然后计算，但每次call 6s至第25次调用时都会出错，不显示结果，但是如果从文件中第25行执行，则至第50次时再次出错。。。。我不明白，为什么会是这样？如果程序本身有错，那这种错误是什么呢？
求指教
将程序运行的结果（前24次）和25次出错提示附下面，求各位fortran老师们指点一二，不甚感激！
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.020               user-defined    *
*                2        0.5             0.020               user-defined    *
*                3        2.5             0.060               user-defined    *
*                4        8.0             0.001               user-defined    *
*                                                                             *
第           1 次运算
  0.1000000      0.1561989
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.040               user-defined    *
*                2        0.5             0.040               user-defined    *
*                3        2.5             0.120               user-defined    *
*                4        8.0             0.002               user-defined    *
*                                                                             *
第           2 次运算
  0.2000000      0.1654775
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.060               user-defined    *
*                2        0.5             0.060               user-defined    *
*                3        2.5             0.180               user-defined    *
*                4        8.0             0.002               user-defined    *
*                                                                             *
第           3 次运算
  0.3000000      0.1749900
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.080               user-defined    *
*                2        0.5             0.080               user-defined    *
*                3        2.5             0.240               user-defined    *
*                4        8.0             0.003               user-defined    *
*                                                                             *
第           4 次运算
  0.4000000      0.1843241
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.100               user-defined    *
*                2        0.5             0.100               user-defined    *
*                3        2.5             0.300               user-defined    *
*                4        8.0             0.004               user-defined    *
*                                                                             *
第           5 次运算
  0.5000000      0.1931361
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.120               user-defined    *
*                2        0.5             0.120               user-defined    *
*                3        2.5             0.360               user-defined    *
*                4        8.0             0.005               user-defined    *
*                                                                             *
第           6 次运算
  0.6000000      0.2012857
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.140               user-defined    *
*                2        0.5             0.140               user-defined    *
*                3        2.5             0.420               user-defined    *
*                4        8.0             0.006               user-defined    *
*                                                                             *
第           7 次运算
  0.7000000      0.2087223
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.160               user-defined    *
*                2        0.5             0.160               user-defined    *
*                3        2.5             0.480               user-defined    *
*                4        8.0             0.006               user-defined    *
*                                                                             *
第           8 次运算
  0.8000000      0.2154519
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.180               user-defined    *
*                2        0.5             0.180               user-defined    *
*                3        2.5             0.540               user-defined    *
*                4        8.0             0.007               user-defined    *
*                                                                             *
第           9 次运算
  0.9000000      0.2215103
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.200               user-defined    *
*                2        0.5             0.200               user-defined    *
*                3        2.5             0.600               user-defined    *
*                4        8.0             0.008               user-defined    *
*                                                                             *
第          10 次运算
   1.000000      0.2269316
           1
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.020               user-defined    *
*                2        0.5             0.020               user-defined    *
*                3        2.5             0.060               user-defined    *
*                4        8.0             0.001               user-defined    *
*                                                                             *
第           1 次运算
  0.1000000      0.1559616
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.040               user-defined    *
*                2        0.5             0.040               user-defined    *
*                3        2.5             0.120               user-defined    *
*                4        8.0             0.002               user-defined    *
*                                                                             *
第           2 次运算
  0.2000000      0.1653229
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.060               user-defined    *
*                2        0.5             0.060               user-defined    *
*                3        2.5             0.180               user-defined    *
*                4        8.0             0.002               user-defined    *
*                                                                             *
第           3 次运算
  0.3000000      0.1749104
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.080               user-defined    *
*                2        0.5             0.080               user-defined    *
*                3        2.5             0.240               user-defined    *
*                4        8.0             0.003               user-defined    *
*                                                                             *
第           4 次运算
  0.4000000      0.1843147
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.100               user-defined    *
*                2        0.5             0.100               user-defined    *
*                3        2.5             0.300               user-defined    *
*                4        8.0             0.004               user-defined    *
*                                                                             *
第           5 次运算
  0.5000000      0.1931854
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.120               user-defined    *
*                2        0.5             0.120               user-defined    *
*                3        2.5             0.360               user-defined    *
*                4        8.0             0.005               user-defined    *
*                                                                             *
第           6 次运算
  0.6000000      0.2013879
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.140               user-defined    *
*                2        0.5             0.140               user-defined    *
*                3        2.5             0.420               user-defined    *
*                4        8.0             0.006               user-defined    *
*                                                                             *
第           7 次运算
  0.7000000      0.2088694
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.160               user-defined    *
*                2        0.5             0.160               user-defined    *
*                3        2.5             0.480               user-defined    *
*                4        8.0             0.006               user-defined    *
*                                                                             *
第           8 次运算
  0.8000000      0.2156360
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.180               user-defined    *
*                2        0.5             0.180               user-defined    *
*                3        2.5             0.540               user-defined    *
*                4        8.0             0.007               user-defined    *
*                                                                             *
第           9 次运算
  0.9000000      0.2217261
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.200               user-defined    *
*                2        0.5             0.200               user-defined    *
*                3        2.5             0.600               user-defined    *
*                4        8.0             0.008               user-defined    *
*                                                                             *
第          10 次运算
   1.000000      0.2271739
           1
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.020               user-defined    *
*                2        0.5             0.020               user-defined    *
*                3        2.5             0.060               user-defined    *
*                4        8.0             0.001               user-defined    *
*                                                                             *
第           1 次运算
  0.1000000      0.1564872
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.040               user-defined    *
*                2        0.5             0.040               user-defined    *
*                3        2.5             0.120               user-defined    *
*                4        8.0             0.002               user-defined    *
*                                                                             *
第           2 次运算
  0.2000000      0.1657971
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.060               user-defined    *
*                2        0.5             0.060               user-defined    *
*                3        2.5             0.180               user-defined    *
*                4        8.0             0.002               user-defined    *
*                                                                             *
第           3 次运算
  0.3000000      0.1753326
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.080               user-defined    *
*                2        0.5             0.080               user-defined    *
*                3        2.5             0.240               user-defined    *
*                4        8.0             0.003               user-defined    *
*                                                                             *
第           4 次运算
  0.4000000      0.1846837
*             user-defined aerosol profile: 4 layers                          *
*              Layer   Height(km)   Opt. thick.(at 0.55 mkm)     Model        *
*                1        0.2             0.100               user-defined    *
*                2        0.5             0.100               user-defined    *
*                3        2.5             0.300               user-defined    *
*                4        8.0             0.004               user-defined    *
*                                                                             *

run-time error M6201: MATH
- atan2: DOMAIN error
Image              PC        Routine            Line        Source
main.exe           0049A119  Unknown               Unknown  Unknown
Incrementally linked image--PC correlation disabled.

topmad

看看是不是有什么内存没有释放 以至于溢出了

言深深

看不到子程序，难以判断。单单从出错来看，猜测很可能是在调用过程中，反复运算一个变量，导致在经过25次迭代之后，超出阈值范围了。

ruilier

言深深 发表于 2014-12-24 11:07
看不到子程序，难以判断。单单从出错来看，猜测很可能是在调用过程中，反复运算一个变量，导致在经过25次迭 ...

恩恩，好的，我找找atan2的变量看看

ruilier

topmad 发表于 2014-12-24 11:01
看看是不是有什么内存没有释放 以至于溢出了

好的，我再仔细找找，谢谢

ruilier

ruilier 发表于 2014-12-24 11:34
恩恩，好的，我找找atan2的变量看看

谢谢啦

ruilier

ruilier 发表于 2014-12-24 11:34
恩恩，好的，我找找atan2的变量看看

果然是这样，按照这个方向改了下就对了。。。。。。。。。