时间数组在建模中超出范围？

这些天我正在尝试使用新的气象数据（以 netcdf 格式而不是旧的 cray 格式给出）运行我的气候模型。模型编译顺利，但是当进行模拟时，模型在第一天运行良好，但在模拟的第二天就停止了，总是在同一时间步长，无论我使用哪个开始日期。错误是：

forrtl: strict (408): fort: (2): 数组 TIMEVALS 的下标 #1 的值为 141，该值大于上限 140。

所以我做了一些研究，浏览了我的源代码，看看 timevals 指的是什么样的数组，我发现它指的是这个具有时间维度的新气象数据，而这样的时间数组由 140 个元素组成。每个元素都是模型应该用于模拟的气象数据的特定日期和时间......所以我开始相信这是我的代码的问题，但我的同事已经能够在没有任何情况下运行模型问题，这对我来说很奇怪。他在 Makefile 中使用一些不同的设置编译了模型，我不知道这是否重要，我对 fortran 等还不是很熟悉。但是下面是使用这个 TIMEVALS 数组的代码部分：

           CASE(2)        ! nudging data is in netcdf-format

           cfile = str_filter(ndg_file_nc,yr,mo,dy,hr,mi,se,ndgblock)
           CALL message('   Adjust date using file: ',TRIM(cfile))

           IF (p_parallel_io) THEN
             INQUIRE(file=cfile,exist=found)
             IF (.NOT.found) &
                CALL finish('NudgingInit','Nudging data file not found.')
             ndgfile%format = NETCDF
             CALL IO_open (cfile, ndgfile, IO_READ)
             CALL IO_INQ_DIMID(ndgfile%file_id, 'time', ndimid)
             CALL IO_INQ_DIMLEN(ndgfile%file_id, ndimid, nts)
             CALL IO_INQ_VARID(ndgfile%file_id, 'time', nvarid)
             ALLOCATE (timevals(nts))
             CALL IO_GET_VAR_DOUBLE (ndgfile%file_id, nvarid, timevals)
             ihead_nc(1) = FLOOR(timevals(1))                          ! ihead_nc(1) is YYYYMMDD
             ihead_nc(2) = INT((timevals(1)-ihead_nc(1))*24._dp)       ! ihead_nc(2) is HH
             DEALLOCATE (timevals)
           ENDIF

           IF (p_parallel) CALL p_bcast(ihead_nc, p_io)
           CALL inp_convert_date(ihead_nc(1),ihead_nc(2)*10000, ndg_date0)

             IF (p_parallel_io) THEN
             ! skip first record and read second header
             ALLOCATE (timevals(nts))
             CALL IO_GET_VAR_DOUBLE (ndgfile%file_id, nvarid, timevals)
             ihead_nc(1) = FLOOR(timevals(2))                          ! ihead_nc(1) is YYYYMMDD
             ihead_nc(2) = INT((timevals(2)-ihead_nc(1))*24._dp)       ! ihead_nc(2) is HH
             DEALLOCATE (timevals)
             CALL IO_close(ndgfile)
           ENDIF

           IF (p_parallel) CALL p_bcast(ihead_nc, p_io)
           CALL inp_convert_date(ihead_nc(1),ihead_nc(2)*10000, ndg_date1)

ndg_file and ndg_date参考微推（气象数据） 你们知道什么可能会导致这个错误吗？

我现在有一些时间来详细阐述我之前的评论。 （请注意，我使用斜体表示您可能想阅读的术语。）

您报告的错误是运行时错误，而不是编译器在编译时能够看到的错误。如果您不明白两者之间的区别run-time（即当代码执行时）和编译时（即当编译器将您的源代码转换为可执行代码时）进行一些研究。此外，很明显，您（或某人）已指示编译器创建一个代码版本，用于检查数组元素访问是否在数组范围内。在测试新软件时，这是一个非常重要的安全功能，但是当代码执行时，一旦通过测试，就在没有进行此检查的情况下编译了如此多的代码，从而造成了性能损失。

我不知道您使用的是什么编译器，但请查看其文档以找到打开的选项数组边界检查在运行时。

错误消息非常明确——在代码中的某个时刻，它尝试访问只有 140 个元素的数组的第 141 个元素。我们无法告诉您这是如何发生的，即使我们看到了您的整个代码，也可能无法告诉您。当加载的数据不符合程序员的期望时，这种事情经常发生。当程序员制作时也经常发生这种情况差一编写循环时出错。我们可能会通过查看您的整个代码来发现这一点，但您比我们处于更有利的位置。

你写

但我的同事能够毫无问题地运行该模型，这对我来说很奇怪。他在 Makefile 中使用一些不同的设置编译了模型，我不知道这是否重要，

嗯，是的，这很重要，非常重要。如果您编写的代码访问具有 140 个元素的数组的元素 141，则 Fortran 与许多其他编译语言一样，将很乐意访问内存中元素 140 之后的下一个位置。通常，您不知道程序正在干扰哪些数据。如果幸运的话，内存中的下一个位置位于操作系统分配给程序的地址空间之外，操作系统会立即停止程序并报告分段错误。

如果你运气不好，程序会继续愉快地读取和写入元素 141，无论它是什么。

我推测您的同事尚未对其代码版本实现数组边界检查。是否告诉他他的代码（很可能）被破坏了，这取决于你。

那么你该怎么办呢？你debug该程序。您可以通过多种方式执行此操作，其中最简单的方法是（在我看来）插入一些write语句在代码中的关键点打印出变量值，以测试您对它们可能、可以或实际采用的值的假设。更困难，但在未来解决问题方面值得付出最初的努力，那就是在一个人的控制下运行代码debugger。有几个很好的调试器可用于所有主要平台上的 Fortran 程序。