24小时热门版块排行榜    

查看: 676  |  回复: 3
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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)

[ Last edited by sunxiao on 2009-4-30 at 11:36 ]
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tjpm

金虫 (正式写手)

是需要设置初值的,不然就会随机给出一个值为初值。
LZ可能每次运行都出来一样的结果,所以觉得是每次没有赋初值的浮点数的值是一样的。
实际上只是因为没有随机数种子。如果先call random_seed()则可能每次的情况不一样。
不过我用的intel的编译器,初值都是默认设置为0了的
4楼2009-03-21 10:54:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

spin

运行这个程序。
其中的可变空数组DC , C 都是很大的一个数  -6.277438562204192E+066
2楼2009-03-19 11:01:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)


sunxiao(金币+1,VIP+0):欢迎参与讨论,鼓励 3-21 09:37
这个问题我也经常遇到,但为什么我不知道,呵呵,所以定义的数组,我一般都设置初值为0.
3楼2009-03-20 19:52:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见