Mathematical Physics - Volume II - Numerical Methods

3.6 Program for solving of elliptical problems

95

& maxband=25) COMMON /ccon / nnode , nelem , nmat , nbce , nbcn , npoint , & nband COMMON /cmatl / prop(3,maxm) WRITE (6,99001) READ (5,*) READ (5,*) nmat

IF ( nmat.LT.1 .OR. nmat.GT.maxm ) THEN PRINT * , ’error in rmat, nmat=’ , nmat STOP ENDIF DO 100 i = 1 , nmat READ (5,*) (prop(j,i),j=1,3) WRITE (6,99002) i , (prop(j,i),j=1,3)

100 CONTINUE RETURN 99001 FORMAT (/’ mat no 99002 FORMAT (I4,5X,3F15.4) END

k

b

f’/)

SUBROUTINE rbc PARAMETER (maxn=400,maxe=100,maxm=5,maxbce=70,maxbcn=40,maxpt=10, & maxband=25) COMMON /ccon / nnode , nelem , nmat , nbce , nbcn , npoint , & nband COMMON /cbc / nodbc(maxbce) , vbce(maxbce) , nelbc(maxbcn) , & nside(maxbcn) , vbcn(2,maxbcn) , npt(maxpt) , & vpt(maxpt) c read point loads ir = 0 READ (5,*) READ (5,*) nrec npoint = nrec IF ( nrec.NE.0 ) THEN ir = ir + 1 DO 50 i = 1 , nrec READ (5,*) n , v

npt(i) = n vpt(i) = v

50

CONTINUE

ENDIF c read essential boundary condition data READ (5,*) READ (5,*) nrec nbce = 0 IF ( nrec.NE.0 ) THEN ir = ir + 1 DO 100 j = 1 , nrec READ (5,*) n1 , num , inc , v