14Fermer16
WarptenLe 19/04/2017 à 10:36
    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

Mon tour, ce bout de code a un bug, K devient 0 ou infinity au pif, et ca sigsegv au pif++ ...