У меня проблема с моей простой программой на Фортране. Я работаю в Fortran 77, используя Compaq Visual Fortran. Структура программы должна быть в форме основной и подпрограммы, потому что она является частью большой программы, связанной с методом конечных элементов.

Моя проблема в том, что я хотел бы установить значения 10000 и 10000 дляNHELE а такжеNVELE соответственно, но когда я запускаю код, программа останавливается и выдает следующую ошибку:

forrt1: server <170>: program Exception - stack overflow

Я пытался итеративно уменьшать требуемые значения, пока не достиг 507 и 507. На этом этапе код работает без ошибок.

Однако увеличение значений до 508 и 508 приводит к повторному появлению той же ошибки.

Я думаю, что проблема связана с подпрограммойNIGTEE, потому что, когда я переставить программу без него, все работает нормально.

Я попытался увеличить размер стека до максимума с помощью менюproject>>settings>>link>>output>>reserve & commit но это не имело значения.

Как я могу решить эту проблему?

Вот моя программа:

PARAMETER(NHELE=508,NVELE=508)
PARAMETER(NHNODE=NHELE+1,NVNODE=NVELE+1)
PARAMETER(NTOTALELE=NHELE*NVELE)

DIMENSION MELE(NTOTALELE,4)

    CALL NIGTEE(NHELE,NVELE,NHNODE,NVNODE,NTOTALELE,MELE)

OPEN(UNIT=7,FILE='MeshNO For Rectangular.TXT',STATUS='UNKNOWN')
WRITE(7,500) ((MELE(I,J),J=1,4),I=1,NTOTALELE)
500 FORMAT(4I20)

    STOP
END

  SUBROUTINE NIGTEE(NHELE,NVELE,NHNODE,NVNODE,NTOTALELE,MELE)
DIMENSION NM(NVNODE,NHNODE),NODE(4)
DIMENSION MELE(NTOTALELE,4)

KK=0
DO 20 I=1,NVNODE
DO 20 J=1,NHNODE
KK=KK+1
NM(I,J)=KK
20  CONTINUE
  KK=0
DO 30 I=1,NVELE
DO 30 J=1,NHELE
NODE(1)=NM(I,J)
NODE(2)=NM(I,J+1)
NODE(3)=NM(I+1,J+1)
NODE(4)=NM(I+1,J)
KK=KK+1
DO 50 II=1,4
50  MELE(KK,II)=NODE(II)

30  CONTINUE
  RETURN
END

Благодарю.