1 |
subroutine Grav(X, M, A)
|
2 | !HPF ALIGN Xdelta1, Vdelta1 . . . WITH X |
3 | . . . |
4 | end subroutine |
5 | ! main program |
6 | program nbody |
7 | use nbodyvars |
8 | . . . |
9 | call setup ( ) |
10 | !HPF DISTRIBUTE X ( :, BLOCK) |
11 | !HPF ALIGN V, M WITH X |
12 | . . . |
13 | do k = 1, np |
14 | call runge-kutta(timestep, ns) |
15 | call print_state () |
16 | end do |
17 | end program |