1 SUBROUTINE garc(tp,xp,zp,csx,csz,npmax,arc,npc)
5 REAL(rprec) :: tp(*),xp(*),zp(*),arc(*)
6 REAL(rprec) :: csx(*),csz(*)
8 REAL(rprec) :: yg(4),wg(4)
9 REAL(rprec) :: el,summ,a,b
11 REAL(rprec) :: ws1,ws2,ws,delta
15 data yg/.069431844202974,.330009478207572,.669990521792428,
16 1 .930568155797027/,wg/.173927422568727,
17 2 2*.326072577431273,.173927422568727/,ngaus/4/
32 ws1=(xp(j)-xp(is))/delta-(3.*a**2-1.)/6.*delta*csx(is)
33 $ +(3.*b**2-1.)/6.*delta*csx(j)
34 ws2=(zp(j)-zp(is))/delta-(3.*a**2-1.)/6.*delta*csz(is)
35 $ +(3.*b**2-1.)/6.*delta*csz(j)
39 ws=sqrt(ws1*ws1+ws2*ws2)