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