24小时热门版块排行榜    

查看: 112  |  回复: 0
当前主题已经存档。

spin

[交流] 【求助】为什莫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 ]
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 spin 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见