[포트란 초보입니다. "array bounds exceeded" 관련 질문드립니다.]

[alstom]

Y bus matrix 변환 프로그램인데 "array bounds exceeded" 에러가 발생해서 며칠째 고생하고 있습니다.

여기에 올라와 있는 글들을 읽었는데 잘 이해가 안되네요.

소스하고 입력파일 첨부합니다.

고수님들의 조언 부탁드립니다.











C MAXIMUM NUMBER OF BUSES

PARAMETER (NMAX=2000)

C MAXIMUM NUMBER OF BRANCHES

C PARAMETER (MXBRCH=2500)
PARAMETER (MXBRCH=2500)
C MAXIMUM NUMBER OF TYPE 2, 3, AND 4 TRANSFORMERS

PARAMETER (MXFMR=1000)

C MAXIMUM NUMBER OF GENERATORS TO MONITER DURING CONTINUATION

C PARAMETER (MONMX=300)

C MAXIMUM NUMBER OF WEAK BUSES TO MONITER DURING CONTINUATION

C PARAMETER (NWKMX=100)

C NUMBER OF BRANCH PARTICIPATIONS TO CALCULATE

C PARAMETER (NBPRT=20)
C PARAMETER (NBPRT=6)

C MAXIMUM NUMBER OF ADMITTANCE MATRIX COMPONENTS

PARAMETER (MAXY=NMAX*8)



C IMPLICIT REAL*8 (A-H,O-Z)

C COMPLEX*16 YSHUNT(NMAX),YVAL(MAXY),ZBRAN(MXBRCH),B,ADMIT,TCMPLX


IMPLICIT REAL*8 (A-H,O-Z)

COMPLEX*16 YSHUNT(NMAX),YVAL(MAXY),ZBRAN(MXBRCH),B,ADMIT,TCMPLX


REAL*8 CHARGE(MXBRCH),RATE(MXBRCH),TAP(MXBRCH),PHASE(MXBRCH),
*RATES(3),RATMIN(MXFMR),RATMAX(MXFMR),STEP(MXFMR),VMIN(MXFMR),
*VMAX(MXFMR),DMY(MXFMR)

INTEGER NPOINT(2000),JCOLY(MAXY),NEXTY(MAXY),IFROM(MXBRCH),
*ITO(MXBRCH),NCKT(MXBRCH),IBRTP(MXBRCH),LINE(MXFMR),NREG(MXFMR),
*NSIDE(MXFMR)




C DEFINITIONS:
C
C N - NUMBER OF BUSES
C NMAX - MAXIMUM NUMBER OF BUSES
C YSHUNT - SHUNT ADMITTANCE AT EACH BUS (FROM BUS DATA)
C NPOINT - MAPS BUS NUMBERS TO INTERNAL BUS NUMBERS
C MAXY - MAXIMUM NUMBER OF ELEMENTS IN Y BUS
C YVAL - COMPLEX ADMITTANCE VALUES IN Y BUS
C JCOLY - COLUMN NUMBER OF Y BUS ENTRIES
C NEXTY - LOCATION OF NEXT VALUE IN SAME ROW OF Y BUS
C NSTORY - ACTUAL NUMBER OF ELEMENTS IN Y BUS
C MXBRCH - MAXIMUN NUMBER OF BRANCHES ALLOWED
C IFROM - BUS NUMBER WHERE BRANCH ORIGINATES
C ITO - BUS NUMBER WHERE BRANCH ENDS
C NCKT - CIRCUIT NUMBER OF BRANCH
C IBRTP - BRANCH TYPE
C 0 - TRANSMISSION LINE
C 1 TO 4 - TRANSFORMER
C ZBRAN - BRANCH IMPEDENCE
C CHARGE - BRANCH CHARGING
C RATE - MVA RATING OF BRANCH
C NRATE - WHICH OF 3 MVA RATINGS IS TO BE USED
C TAP - VOLTAGE RATIO OF TRANSFORMER
C PHASE - PHASE SHIFT OF TRANSFORMER
C NBRCH - ACTUAL NUMBER OF BRANCHES
C LIUNIT - NUMBER OF INPUT FILE



OPEN(UNIT=5,FILE='DATA.IN')
OPEN(UNIT=1000,FILE='DATA.OUT')


READ(5, *)


DO 100 I=1,N
YVAL(I)=YSHUNT(I)
JCOLY(I)=I
NEXTY(I)=0
100 CONTINUE

C SET Y BUS STORAGE POINTER

ISTORE=N+1

C INITIALIZE REGULATING TRANSFORMER COUNT

NXFMR=0


DO 200 ITER=1,MXBRCH+1


READ(5,10)IFROM(ITER),ITO(ITER),NCKT(ITER),IBRTP(ITER),
*ZBRAN(ITER),CHARGE(ITER),RATES(1),RATES(2),RATES(3),NBREG,
*ISIDE,TAP(ITER),PHASE(ITER),TAPMIN,TAPMAX,TSTEP,VLTMIN,VLTMAX,IDX

10 FORMAT(2(I4,1X),6X,2(I1,1X),2F10.6,F9.5,3(1X,F5.0),1X,I4,1X,
*I1,2X,F6.4,1X,F7.2,3F7.4,1X,2F7.4,I6)




C CHECK FOR END OF DATA

IF(IFROM(ITER) .EQ. -999)THEN
NBRCH=ITER-1
NSTORY=ISTORE-1
C RETURN
END IF

C FIND INTERNAL BUS NUMBERS

I=NPOINT(IFROM(ITER))
J=NPOINT(ITO(ITER))


IF(IBRTP(ITER) .EQ. 0)THEN

ADMIT=(1.0D0,0.0D0)/ZBRAN(ITER)
YVAL(I)=YVAL(I)+ADMIT
YVAL(J)=YVAL(J)+ADMIT
CALL YSTORE(I,J,ADMIT,ISTORE,MAXY,YVAL,JCOLY,NEXTY,N)
CALL YSTORE(J,I,ADMIT,ISTORE,MAXY,YVAL,JCOLY,NEXTY,N)
END IF

200 CONTINUE

WRITE(1000,40)NBRCH
40 FORMAT(1X,I5,' BRANCHES READ IN')
WRITE(6,*)

WRITE(1000,60)NSTORY
60 FORMAT(1X,I5,' ELEMENTS STORED IN ADMITTANCE MATRIX')
WRITE(6,*)

END



C*************************************************************************
C---------------- YSTORE -------------------------------------------------
C*************************************************************************
SUBROUTINE YSTORE(I,J,VALUE,ISTORE,MAXY,YVAL,JCOLY,NEXTY,N)

IMPLICIT REAL*8 (A-H,O-Z)

COMPLEX*16 YVAL(MAXY),VALUE

INTEGER JCOLY(MAXY),NEXTY(MAXY)

C DEFINITIONS:
C
C I - ROW NUMBER OF ELEMENT TO BE STORED
C J - COLUMN NUMBER OF ELEMENT TO BE STORED
C VALUE - ADMITTANCE VALUE TO BE STORED
C ISTORE - POINTER FOR LAST ELEMENT STORED
C MAXY - MAXIMUM NUMBER OF ELEMENTS ALLOWED IN Y BUS
C YVAL - ELEMENTS OF Y BUS
C JCOLY - COLUMN NUMBERS OF Y BUS ELEMENTS
C NEXTY - NUMBER OF NEXT Y BUS VALUE IN SAME ROW

C STORAGE IF ONLY DIAGONAL ELEMENT EXISTS

IF(NEXTY(I) .EQ. 0)THEN

NEXTY(I)=ISTORE
NEXTY(ISTORE)=0
YVAL(ISTORE)=-VALUE
JCOLY(ISTORE)=J
ISTORE=ISTORE+1
RETURN

C STORAGE IF MORE THAN DIAGONAL EXISTS

ELSE

IPOS=NEXTY(I)
DO 20 ITER=1,N

C IF THIS ELEMENT CORRRESPONDS TO THE ONE TO BE STORED

IF(JCOLY(IPOS) .EQ. J)THEN

YVAL(IPOS)=YVAL(IPOS)-VALUE
RETURN

C IF THIS ELEMENT IS NOT THE RIGHT ONE AND THERE ARE NO MORE

ELSE IF(NEXTY(IPOS) .EQ. 0)THEN

NEXTY(IPOS)=ISTORE
NEXTY(ISTORE)=0
YVAL(ISTORE)=-VALUE
JCOLY(ISTORE)=J
ISTORE=ISTORE+1
RETURN

C IF FIRST OPTIONS FAIL, GO TO NEXT ELEMENT IN ROW

ELSE

IPOS=NEXTY(IPOS)

END IF

20 CONTINUE

END IF

END

[경제운]

array bounds exceeded => 포트란FAQ에 있습니다.

[꿈꾸는가인]

배열 변수 정의할 때 잡아 놓았던 것보다 클 경우 발생합니다. 즉 a(10)이라고 정의를 하였는데 a(11)을 사용하려면 문제가 되겠지요? ^^