|
subroutine foo(mxgtot,mxsh)
|
|
logical b
|
|
dimension ex(mxgtot),cs(mxgtot)
|
|
do 500 jg = k1,ig
|
|
u = ex(ig)+ex(jg)
|
|
z = u*sqrt(u)
|
|
x = cs(ig)*cs(jg)/z
|
|
if (ig .eq. jg) go to 480
|
|
x = x+x
|
|
480 continue
|
|
y = y+x
|
|
500 continue
|
|
if(y.gt.t) z=1/sqrt(y)
|
|
if (b) then
|
|
write(9) z
|
|
endif
|
|
end
|