| 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 |