gkd.for

有限元求解柏松方程。本文采用FORTRAN语言编制程序。程序中大部分变量采用有名公共区存储方式存储

	SUBROUTINE GKD !(NG,NP,ND,LEE,II,KD,NS)
	INTEGER N,M,NG,NP,ND,II(3,2*80*80),LEE,NS,
	+        KD(81*81+1)
	DOUBLE PRECISION XY(2,81*81)
	COMMON /C1/N,M,NP,ND/C2/NS,KD,NG,LEE
	COMMON /C6/XY,II
	KD(1)=0
	PRINT*,'NG=',NG
	DO 10, L=1,NG
	   MIN=NG
	   DO 20, LE=1,LEE
	     DO 30, J=1,ND
	       IF(L.EQ.II(J,LE)) GOTO 40
   30	     CONTINUE
           GOTO 70
   40      DO 50, K=1,ND
             IF(II(K,LE).LT.MIN) MIN=II(K,LE)
   50      CONTINUE
   70		 CONTINUE
   20   CONTINUE
        LI=NP*(L-MIN)
	  I=NP*(L-1)
	  DO 60, N=1,NP
	    KD(I+N+1)=KD(I+N)+LI+N
   60   CONTINUE
   10 CONTINUE
      NS=KD(NG*NP+1)
	END