1      subroutine foo(mxgtot,mxsh)
2      logical b
3      dimension ex(mxgtot),cs(mxgtot)
4         do 500 jg = k1,ig
5            u = ex(ig)+ex(jg)
6            z = u*sqrt(u)
7            x = cs(ig)*cs(jg)/z
8            if (ig .eq. jg) go to 480
9               x = x+x
10  480       continue
11            y = y+x
12  500    continue
13      if(y.gt.t) z=1/sqrt(y)
14      if (b) then
15         write(9) z
16      endif
17      end
18