13       numat = 0
14       call valdis(var,ipntr,numat,iop,iasel,istat)
15       if (istat.eq.0) call eemcalc(var,ipntr,numat)
22       subroutine valdid(var,ipntr,numat,iop,iasel,istat,  argument
255       numat = icnt
260       subroutine eemcald(var,ipntr,numat,ianz,coo,qat)  argument
272       do i=1,(numat-1)
274          x(i,numat+1) = -1
275          x(numat+1,i) = 1
277          do j=i+1,numat
286       x(numat,numat)     = 2*var(numat,2)
287       x(numat,numat+1)   = -1
288       x(numat+1,numat)   = 1
289       x(numat+1,numat+1) = 0
291       y(numat+1)       = dble(itot)
293       do i=1,numat
297       call dgefa(x,mxeat+1,numat+1,ipvt,info)
298       call dgedi(x,mxeat+1,(numat+1),ipvt,det,work,01)
300       call mtmul(x,y,q,numat)
306       do i=1,(numat)
320       subroutine mtmul(a,y,b,numat)  argument
324       integer i,j,numat  local
327       do i=1,numat+1
331          do j=1,numat+1