DO I = 1, N_STEP CALL PREDICT(DT) CALL CALCULATE_LJ() CALL CORRECT(DT) CALL BOUNDARY_CONDITIONS() CALL CORRECT_LJ() K = 0.0 DO J = 1, N K = K + VX(I) ** 2 + VY(I) ** 2 + VZ(I) ** 2 END DO K = 0.5 * M * K T = T + 2.0 / 3.0 * K / N WRITE (*, '(''ENERGIES AND TEMPERATURES '')', ADVANCE = 'NO') WRITE (*, *) I, V, K, VC, (T / I) END DO