Ignore:
Timestamp:
Jan 18, 2011 4:19:48 PM (10 years ago)
Author:
suehring
Message:

Right computation of the pressure using Runge-Kutta weighting coefficients. Consideration of the pressure gradient during the time integration removed. Removed bugfix concerning velocity variances.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/timestep_scheme_steering.f90

    r484 r673  
    44! Current revisions:
    55! -----------------
    6 !
     6! No pressure term during time integration (tsc(4)=0.0).
    77!
    88! Former revisions:
     
    6060             IF ( timestep_scheme == 'leapfrog+euler'  .OR. &
    6161                  timestep_scheme == 'euler' .OR. simulated_time == 0.0 )  THEN
    62                 tsc(1:5) = (/ 1.0, 1.0, 0.0, 1.0, 1.0 /)
     62                tsc(1:5) = (/ 1.0, 1.0, 0.0, 0.0, 1.0 /)
    6363             ELSE
    64                 tsc(1:5) = (/ 0.0, 2.0, 0.0, 1.0, 2.0 /)
     64                tsc(1:5) = (/ 0.0, 2.0, 0.0, 0.0, 2.0 /)
    6565             ENDIF
    6666          ELSE
     
    6969!--          user.
    7070             IF ( timestep_scheme == 'euler' )  THEN
    71                 tsc(1:5) = (/ 1.0, 1.0, 0.0, 1.0, 1.0 /)
     71                tsc(1:5) = (/ 1.0, 1.0, 0.0, 0.0, 1.0 /)
    7272             ELSE
    73                 tsc(1:5) = (/ 0.0, 2.0, 0.0, 1.0, 2.0 /)
     73                tsc(1:5) = (/ 0.0, 2.0, 0.0, 0.0, 2.0 /)
    7474             ENDIF
    7575          ENDIF
     
    8484          IF ( simulated_time == 0.0 )  THEN
    8585             dt_changed = .TRUE.
    86              tsc(1:5) = (/ 1.0, 1.0, 0.0, 1.0, 1.0 /)
     86             tsc(1:5) = (/ 1.0, 1.0, 0.0, 0.0, 1.0 /)
    8787          ELSE
    8888             dt_changed = .FALSE.
    8989             IF ( timestep_scheme == 'euler' )  THEN
    90                 tsc(1:5) = (/ 1.0, 1.0, 0.0, 1.0, 1.0 /)
     90                tsc(1:5) = (/ 1.0, 1.0, 0.0, 0.0, 1.0 /)
    9191             ELSE
    92                 tsc(1:5) = (/ 0.0, 2.0, 0.0, 1.0, 2.0 /)
     92                tsc(1:5) = (/ 0.0, 2.0, 0.0, 0.0, 2.0 /)
    9393             ENDIF
    9494          ENDIF
Note: See TracChangeset for help on using the changeset viewer.