Changeset 673 for palm/trunk/SOURCE/timestep_scheme_steering.f90

Timestamp:

Jan 18, 2011 4:19:48 PM (14 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:

• palm/trunk/SOURCE/timestep_scheme_steering.f90 (modified) (4 diffs)

palm/trunk/SOURCE/timestep_scheme_steering.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-! 
+! No pressure term during time integration (tsc(4)=0.0).
 !
 ! Former revisions:
@@ -60 +60 @@
              IF ( timestep_scheme == 'leapfrog+euler'  .OR. &
                   timestep_scheme == 'euler' .OR. simulated_time == 0.0 )  THEN
-                tsc(1:5) = (/ 1.0, 1.0, 0.0, 1.0, 1.0 /)
+                tsc(1:5) = (/ 1.0, 1.0, 0.0, 0.0, 1.0 /)
              ELSE
-                tsc(1:5) = (/ 0.0, 2.0, 0.0, 1.0, 2.0 /)
+                tsc(1:5) = (/ 0.0, 2.0, 0.0, 0.0, 2.0 /)
              ENDIF
           ELSE
@@ -69 +69 @@
 !--          user.
              IF ( timestep_scheme == 'euler' )  THEN
-                tsc(1:5) = (/ 1.0, 1.0, 0.0, 1.0, 1.0 /)
+                tsc(1:5) = (/ 1.0, 1.0, 0.0, 0.0, 1.0 /)
              ELSE
-                tsc(1:5) = (/ 0.0, 2.0, 0.0, 1.0, 2.0 /)
+                tsc(1:5) = (/ 0.0, 2.0, 0.0, 0.0, 2.0 /)
              ENDIF
           ENDIF
@@ -84 +84 @@
           IF ( simulated_time == 0.0 )  THEN
              dt_changed = .TRUE.
-             tsc(1:5) = (/ 1.0, 1.0, 0.0, 1.0, 1.0 /)
+             tsc(1:5) = (/ 1.0, 1.0, 0.0, 0.0, 1.0 /)
           ELSE
              dt_changed = .FALSE.
              IF ( timestep_scheme == 'euler' )  THEN
-                tsc(1:5) = (/ 1.0, 1.0, 0.0, 1.0, 1.0 /)
+                tsc(1:5) = (/ 1.0, 1.0, 0.0, 0.0, 1.0 /)
              ELSE
-                tsc(1:5) = (/ 0.0, 2.0, 0.0, 1.0, 2.0 /)
+                tsc(1:5) = (/ 0.0, 2.0, 0.0, 0.0, 2.0 /)
              ENDIF
           ENDIF