| 查看: 147 | 回复: 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 ] |
» 猜你喜欢
售一区SCI文章T0P,我:8O.551.O54,科目全,可十急
已经有7人回复
售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急
已经有10人回复
售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
售SCI文章,我:8O.5.5.1O.54,科目全,可十急
已经有8人回复
售SCI文章,我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急
已经有10人回复
售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急
已经有8人回复
国社科又开始会评了,不知道这次命运如何
已经有14人回复
售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
麻烦专家们看看评委们的意见(F口面上)
已经有13人回复










,AC(:,
回复此楼
30