| 查看: 112 | 回复: 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 ] |
» 猜你喜欢
实验室接单子
已经有7人回复
假如你的研究生提出不合理要求
已经有11人回复
全日制(定向)博士
已经有5人回复
萌生出自己或许不适合搞科研的想法,现在跑or等等看?
已经有4人回复
Materials Today Chemistry审稿周期
已经有4人回复
参与限项
已经有3人回复
对氯苯硼酸纯化
已经有3人回复
求助:我三月中下旬出站,青基依托单位怎么办?
已经有12人回复
所感
已经有4人回复
要不要辞职读博?
已经有7人回复












,AC(:,
回复此楼