作者RuinGgg (坦克)
看板Fortran
标题[问题] 新手请教~"~
时间Fri Oct 31 00:16:16 2008
首先 我先写了一个副程式
SUBROUTINE IDMAT(NFIX,IDND,NDN,NNOD,NEQ)
.
.
.
DO I=1,NNOD
DO J=1,NDN
.
.
.
IDND(J,I)=xxx
.
.
.
END DO
END DO
WRITE(*,*) IDND
^^^^这个IDND矩阵的值有跑出来，没问题
接下来我要再写一个副程式
SUBROUTINE MEMDOF(IDBC,NDN,NDE,NBC,NNOD,IDND,LM)
IMPLICIT REAL*8 (A-H,O-Z)
IMPLICIT INTEGER*4 (I-N)
INTEGER*4 IDBC(5,NBC),IDND(NDN,NNOD),LM(NDE,NBC),NDN,NNOD,NDE,NBC
WRITE(*,*) IDND
^^^^我想去测试有没有办法读到上面副程式所写的IDND
但结果显示读不出东西，请问为什麽呢???
---------------------------------------------------------------
Sorry 附上部份主程式...
PROGRAM
IMPLICIT REAL*8 (A-H,O-Z)
IMPLICIT INTEGER*4 (I-N)
REAL*8 A(200000)
INTEGER*4 MAXMEM,IEND,IRI
INTEGER*4 NNOD,NBC,NMAT,NSEC,ITP,NNE,IFORCE
INTEGER*4 ICOOR,IVECTY,INFIX,IEXLD,IPROP,ISECT,IIDBC,IIDND,IFEF
INTEGER*4 ILM,IGLOAD,IGLK,IDISP,IELFOR
.
.
.
C
C Compute pointers for the input arrays
C
ICOOR = 1
INFIX = ICOOR + NCO*NNOD
IEXLD = INFIX +(NDN*NNOD+IRI-1)/IRI
IIDBC = IEXLD +NDN*NNOD
IVECTY = IIDBC+(5*NBC+IRI-1)/IRI
IF (ITP .EQ. 6) THEN
IFEF = IVECTY+ 3*NBC
ELSE
IFEF = IVECTY
ENDIF
IF (IFORCE .EQ. 1) IPROP=IFEF
IF (IFORCE .EQ. 2) IPROP=IFEF+NDE*NBC
ISECT =IPROP+5*NMAT
IEND =ISECT+5*NSEC
C
C Memory check
C
PLACE='INPUT'
CALL MEMCHECK(IEND,MAXMEM,PLACE)
C
C Read the remaining data
C
CALL INPUT(A(ICOOR),A(INFIX),A(IEXLD),A(IIDBC),A(IVECTY),A(IFEF),
1 A(IPROP),A(ISECT),NNOD,NBC,NMAT,NSEC,ITP,NCO,NDN,NDE,
2 IFORCE)
CLOSE(UNIT=1)
C***** UP TO HERE --- PROG 1 *****
IIDND = IEND
IEND = IIDND +(NDN*NNOD+IRI-1)/IRI
C
PLACE='IDMAT'
CALL MEMCHECK(IEND,MAXMEM,PLACE)
CALL IDMAT(A(INFIX),A(IIDND),NDN,NNOD,NEQ)
C Compute LM(NDE,NBC)
ILM = IEND
IEND = ILM+(NDE*NBC+IRI-1)/IRI
C
PLACE='MEMDOF'
CALL MEMCHECK(IEND,MAXMEM,PLACE)
C
CALL MEMDOF(A(IIDBC),A(IIDND),A(ILM),NBC,NDE,NDN,NNOD)
.
.
.
end
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.112.217.19
1F：→ DMFC:没主程式嬷 @_@? 线索只有这样~谁解的出来 ...... =_= 10/31 06:17
※ 编辑: RuinGgg 来自: 125.233.1.90 (11/01 01:40)
2F：→ blc:错误讯息是？ 11/02 22:13
3F：→ DMFC:我的猜测...IDND宣告问题乎?!?! 11/03 01:10
4F：→ DMFC:检查看看主程式跟sub对IDND的宣告是否一致 11/03 01:10
5F：→ RuinGgg:嗯嗯~谢谢^^ 我跑出来了 我发现我回传顺序和宣告都有点误 11/03 16:18