薄膜表面粗糙度fortran程序.docVIP

c phase model(2012)explicit finite difference method 2012.2.29. c this programm debuged through on 2012.2.29. c cc !!!!!!!!!!!!!!!!!!!!!!! PROGRAM MAIN!!!!!!!!!!!!!!!!!!!!!!! IMPLICIT DOUBLE PRECISION (A-H,O-Z) INTEGER I,J,N,k,h DIMENSION FAI(351,351),CATA(351,351) 1 ,FAI1(301,301),CATA1(301,301) 1 ,dfai(301,301),dcata(301,301),dfai1(301,301) 1 ,dcata1(301,301),A(300),A1(300) c c OPEN (1,FILE=D:\TEST\OUTFL1) OPEN (6,FILE=D:\TEMP4\OF6.DAT) OPEN (4,FILE=D:\TEMP4\OF4.DAT) OPEN (5,FILE=D:\TEMP4\OF5.DAT) OPEN (7,FILE=D:\TEMP4\OF7.DAT) C C !!!!!!!!!!!!!!!!!!!!!!!!!!!!parameters!!!!!!!!!!!!!!!!!!!!!!!!!!!!! C LMD=0 LMDN=358.70 PI=3.14 D=800 W=2 TAO=1 n=300 h=300 C c !!!!!!!!!!!!!!!!!!!!!initial value!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! C DO I=1,300 DO J=1,300 FAI(i,j)=0.0 CATA(i,j)=0.0 END DO END DO OPEN(100,FILE=D:\TEMP4\OF24.DAT,STATUS=OLD) DO k=1,h READ(100,(F10.0))a(k) END DO OPEN(101,FILE=D:\TEMP4\OF240.DAT,STATUS=OLD) DO k=1,h READ(101,(F10.0))A1(k) END DO DO k=1,h do i=1,300 do j=1,300 q1=(i-A(k))**2+(j-A1(k))**2 if(q1.LT.28)then FAI(i,j)=1.0 CATA(i,j)=1.0 end if end do end do end do WRITE(6,20)((cata(I,j),j=1,300),i=1,300) C C !!!!!!!!!!!!!!!!!!!!!!!纽曼边界条件!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! C do i=1,n do j=1,n dfai(-1,j)=dfai(1,j) !!!!!!!!!!!!!!!!!!!!!一阶导数为零!!!!!!!!!!!!!!!!!!! dfai(N+1,j)=dfai(N-1,j) dfai(I,-1)=dfai(I,1) dfai(I,N-1)=dfai(I,N+1) dcata(-1,j)=dcata(1,j) dcata(N-1,j)=dcata(N+1,j) dcata(I,-1)=dcata(I,1) dcata(I,N-1)=dcata(I,N+1) end do end do c C C BEGAIN CACULATION C DTIME=0.000001 c dtime=0.000005 dx=2 C DO KK=1,1000 C C CACLATE THE VALUE OF DFAI AND DCATA OF EACH POINT C C !!!!!!!!!!!!!!!差分格式!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! C do i=2,n-1 do J=2,n-1 dfai(i,j)=FAI(I-1,J)+FAI(