[½»Á÷] ¡¾ÇóÖú¡¿ÎªÊ²Äª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 ]

»Ø¸´´ËÂ¥