官方的tfm案例为什么算出来有错误

MFiX
1

当我跟着官网的TFM案例一步步过来的时候，到最后计算的时候会报错。我试着重新去操作了一边，依旧会报错。

image
image1256×473 57.3 KB

image
image486×153 33.7 KB

2

Hi. Sorry, I can’t read Chinese. If you want help from the development team, please attach your input files, as well as all log files.
The code at set_icbc_flags.f:211 looks like this:

! If an error is detected, have each rank open a log file and write
! it's own message. Otherwise, we need to send all the data back to
! PE_IO and that's too much work!
      IF(ERROR) THEN

         CALL OPEN_PE_LOG(IER)

         WRITE(ERR_MSG, 1100) trim(iVal(myPE))
         CALL LOG_MESSAGE(__FILE__, __LINE__, LOGLEVEL_INFO, FOOTER=.FALSE.)

         DO K = kStart2, kEnd2
         DO J = jStart2, jEnd2
         DO I = iStart2, iEnd2

            IF(ICBC_FLAG(FUNIJK(I,J,K)) == '   ') THEN
               WRITE(ERR_MSG,1101) I, J, K
               CALL LOG_MESSAGE(__FILE__, __LINE__, LOGLEVEL_INFO, HEADER=.FALSE., FOOTER=.FALSE.)
            ENDIF

         ENDDO
         ENDDO
         ENDDO

         CALL LOG_ERROR()

which shows that the actual error message is in a different file OPEN_PE_LOG(IER)