| 查看: 119 | 回复: 0 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】为什莫fortran的空数组自动赋值-6.277438562204192E+066
|
|||
|
IMPLICIT NONE INTEGER(KIND=8)::li=10**65 REAL*8,PARAMETER:: PI=3.1415926535 COMPLEX*16,ALLOCATABLE:: speq1( ,AC(:, ,speqd( ,DC(:,![]() REAL*8,ALLOCATABLE:: A(:, ,AA(:, ,C(:, ,D(:,![]() INTEGER i,j,k,NX,NY,NXX,NYY,M,N,FNN,FNM REAL*8 dkx,dky,L,W,kx,ky,dx,dy REAL*8 cx,cy,ms,dum,NyqL,NyqW,NLW,krad,corner,KCx,KCy,UpperK INTEGER ci,cj,lL,rL,bW,tW,dw OPEN(101,FILE='slipgen.out') OPEN(105,FILE='specx.txt') OPEN(106,FILE='specy.txt') OPEN(130,FILE='slipgen.in') OPEN(120,FILE='slipx.txt') OPEN(121,FILE='slipy.txt') !Startup settings CALL RANDOM_SEED() read(130,*)L,W read(130,*)M,N read(130,*)UpperK FNN=N/2+1;FNM=M/2+1 ALLOCATE(speq1(N),AC(M/2,N),speqd(N),DC(M/2,N),A(M,N),AA(M,N),C(M,N),D(M,N)) dkx=1./L;dky=1./W dx=1./dkx/real(M) dy=1./dky/real(N) KCx=UpperK/L !Corner wave-number for along-strike direction KCy=UpperK/W !Corner wave-number for down-dip direction print*,DC(1,1)IMPLICIT NONE INTEGER(KIND=8)::li=10**65 REAL*8,PARAMETER:: PI=3.1415926535 COMPLEX*16,ALLOCATABLE:: speq1( ,AC(:, ,speqd( ,DC(:,![]() REAL*8,ALLOCATABLE:: A(:, ,AA(:, ,C(:, ,D(:,![]() INTEGER i,j,k,NX,NY,NXX,NYY,M,N,FNN,FNM REAL*8 dkx,dky,L,W,kx,ky,dx,dy REAL*8 cx,cy,ms,dum,NyqL,NyqW,NLW,krad,corner,KCx,KCy,UpperK INTEGER ci,cj,lL,rL,bW,tW,dw OPEN(101,FILE='slipgen.out') OPEN(105,FILE='specx.txt') OPEN(106,FILE='specy.txt') OPEN(130,FILE='slipgen.in') OPEN(120,FILE='slipx.txt') OPEN(121,FILE='slipy.txt') !Startup settings CALL RANDOM_SEED() read(130,*)L,W read(130,*)M,N read(130,*)UpperK FNN=N/2+1;FNM=M/2+1 ALLOCATE(speq1(N),AC(M/2,N),speqd(N),DC(M/2,N),A(M,N),AA(M,N),C(M,N),D(M,N)) dkx=1./L;dky=1./W dx=1./dkx/real(M) dy=1./dky/real(N) KCx=UpperK/L !Corner wave-number for along-strike direction KCy=UpperK/W !Corner wave-number for down-dip direction print*,DC(1,1) [ Last edited by sunxiao on 2009-4-30 at 11:36 ] |
» 猜你喜欢
球磨粉体时遇到了大的问题,请指教!
已经有15人回复
过年走亲戚时感受到了所开私家车的鄙视链
已经有5人回复
情人节自我反思:在爱情中有过遗憾吗?
已经有5人回复
江汉大学解明教授课题组招博士研究生/博士后
已经有3人回复













,AC(:,
回复此楼