1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
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