1*59599516SKenneth E. Jansen module spebc 2*59599516SKenneth E. Jansen 3*59599516SKenneth E. Jansen integer nfin, nelint, npin, npint, nfint 4*59599516SKenneth E. Jansen real*8 sang 5*59599516SKenneth E. Jansen real*8 xnrml, ynrml, znrml, aR, aI 6*59599516SKenneth E. Jansen 7*59599516SKenneth E. Jansen 8*59599516SKenneth E. Jansen 9*59599516SKenneth E. Jansen real*8, allocatable :: xyn(:), xynin(:) 10*59599516SKenneth E. Jansen real*8, allocatable :: xcyl(:,:),xintl(:,:,:),xsinfin(:,:,:) 11*59599516SKenneth E. Jansen 12*59599516SKenneth E. Jansen 13*59599516SKenneth E. Jansen integer, allocatable :: ien2D(:,:), nen1(:), elcnfin(:,:) 14*59599516SKenneth E. Jansen integer, allocatable :: nrint(:), imax(:) 15*59599516SKenneth E. Jansen 16*59599516SKenneth E. Jansen end module 17*59599516SKenneth E. Jansen 18*59599516SKenneth E. Jansenc----------------------------------------------------------------------- 19*59599516SKenneth E. Jansenc allocate the spebc arrays 20*59599516SKenneth E. Jansenc----------------------------------------------------------------------- 21*59599516SKenneth E. Jansen subroutine setSPEBC(numnp,nsd) 22*59599516SKenneth E. Jansen 23*59599516SKenneth E. Jansen use spebc 24*59599516SKenneth E. Jansen 25*59599516SKenneth E. Jansen allocate (xyn(numnp)) 26*59599516SKenneth E. Jansen allocate (xynin(numnp)) 27*59599516SKenneth E. Jansen allocate (xcyl(numnp,nsd)) 28*59599516SKenneth E. Jansen allocate (nen1(numnp)) 29*59599516SKenneth E. Jansen 30*59599516SKenneth E. Jansenc allocate (elcnpin(numnp)) 31*59599516SKenneth E. Jansenc allocate (xsi(numnp,nsd)) 32*59599516SKenneth E. Jansen 33*59599516SKenneth E. Jansen 34*59599516SKenneth E. Jansen 35*59599516SKenneth E. Jansen return 36*59599516SKenneth E. Jansen end 37*59599516SKenneth E. Jansen 38