@@ -91,16 +91,19 @@ SUBROUTINE velocity_verlet(pos, acc, dt, time, vel)
9191 ! Verlet
9292 CALL verlet(pos, acc, dt, time)
9393
94- WRITE (* ,* ) ' [#] Time for Verlet integration:'
94+ WRITE (* ,* ) ' [#]'
95+ WRITE (* ,* ) ' Time for Verlet integration:'
9596 WRITE (* ,* ) time
9697
9798 ! stormer Verlet
9899 pos = 5d0
99100 CALL stormer_verlet(pos, acc, dt, time, vel)
100101
101- WRITE (* ,* ) ' [#] Time for Stormer Verlet integration:'
102+ WRITE (* ,* ) ' [#]'
103+ WRITE (* ,* ) ' Time for Stormer Verlet integration:'
102104 WRITE (* ,* ) time
103- WRITE (* ,* ) ' [#] Velocity for Stormer Verlet integration:'
105+ WRITE (* ,* ) ' [#]'
106+ WRITE (* ,* ) ' Velocity for Stormer Verlet integration:'
104107 WRITE (* ,* ) vel
105108
106109
@@ -109,9 +112,11 @@ SUBROUTINE velocity_verlet(pos, acc, dt, time, vel)
109112 pos = 5d0
110113 CALL velocity_verlet(pos, acc, dt, time, vel)
111114
112- WRITE (* ,* ) ' [#] Time for velocity Verlet integration:'
115+ WRITE (* ,* ) ' [#]'
116+ WRITE (* ,* ) ' Time for velocity Verlet integration:'
113117 WRITE (* ,* ) time
114- WRITE (* ,* ) ' [#] Velocity for velocity Verlet integration:'
118+ WRITE (* ,* ) ' [#]'
119+ WRITE (* ,* ) ' Velocity for velocity Verlet integration:'
115120 WRITE (* ,* ) vel
116121
117122END PROGRAM verlet_integration
0 commit comments