| 1 | INTEGER, PARAMETER :: LDX = 129, NX1 = 128, NX2 = 256 |
|---|
| 2 | INTEGER, PARAMETER :: LDY = 64, NY1 = 64, NY2 = 256 |
|---|
| 3 | REAL(KIND=8), DIMENSION(LDX, NX2) :: X |
|---|
| 4 | COMPLEX(KIND=8), DIMENSION(LDY, NY2) :: Y |
|---|
| 5 | REAL(KIND=8), DIMENSION(512*NX2) :: WORK |
|---|
| 6 | REAL(KIND=8), DIMENSION(100+2*(NX1+NX2)) :: TABLE |
|---|
| 7 | INTEGER :: IRC |
|---|
| 8 | CHARACTER(LEN=80) :: MESSAGE |
|---|
| 9 | |
|---|
| 10 | ! ... (Definition de X) |
|---|
| 11 | |
|---|
| 12 | CALL JMSETERREUR(.false.) |
|---|
| 13 | CALL SCFFT2D(0, 128, 256, 1.d0, X, 129, Y, 64, TABLE, WORK, 0) |
|---|
| 14 | CALL JMGETCODE(IRC) |
|---|
| 15 | IF (IRC .NE. 0) THEN |
|---|
| 16 | CALL JMGETMESSAGE(IRC,MESSAGE) |
|---|
| 17 | PRINT *,'Code',IRC,', ',TRIM(MESSAGE) |
|---|
| 18 | CALL HANDLER() |
|---|
| 19 | END IF |
|---|
| 20 | CALL SCFFT2D(1, 128, 256, 1.d0, X, 129, Y, 64, TABLE, WORK, 0) |
|---|
| 21 | |
|---|
| 22 | ! ... (Impression de Y) |
|---|
| 23 | |
|---|
| 24 | END |
|---|
| 25 | |
|---|
| 26 | SUBROUTINE HANDLER() |
|---|
| 27 | STOP |
|---|
| 28 | END SUBROUTINE HANDLER |
|---|
