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

Timestamp:

Jan 18, 2011 4:19:48 PM (13 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/prognostic_equations.f90 (modified) (11 diffs)

palm/trunk/SOURCE/prognostic_equations.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-!
+! Consideration of the pressure gradient (steered by tsc(4)) during the time 
+! integration removed.
 ! Former revisions:
 ! -----------------
@@ -212 +213 @@
                           dt_3d * (                                            &
                                    tsc(2) * tend(k,j,i) + tsc(3) * tu_m(k,j,i) &
-                                 - tsc(4) * ( p(k,j,i) - p(k,j,i-1) ) * ddx    &
                                   ) -                                          &
                            tsc(5) * rdf(k) * ( u(k,j,i) - ug(k) )
@@ -297 +297 @@
                           dt_3d * (                                            &
                                    tsc(2) * tend(k,j,i) + tsc(3) * tv_m(k,j,i) &
-                                 - tsc(4) * ( p(k,j,i) - p(k,j-1,i) ) * ddy    &
                                   ) -                                          &
                           tsc(5) * rdf(k) * ( v(k,j,i) - vg(k) )
@@ -383 +382 @@
                           dt_3d * (                                            &
                                    tsc(2) * tend(k,j,i) + tsc(3) * tw_m(k,j,i) &
-                              - tsc(4) * ( p(k+1,j,i) - p(k,j,i) ) * ddzu(k+1) &
                                   ) -                                          &
                           tsc(5) * rdf(k) * w(k,j,i)
@@ -899 +897 @@
 !------------------------------------------------------------------------------!
 ! Version with one optimized loop over all equations. It is only allowed to
-! be called for the standard Piascek-Williams advection scheme.
+! be called for the Wicker and Skamarock or Piascek-Williams advection scheme.
 !
 ! Here the calls of most subroutines are embedded in two DO loops over i and j,
@@ -983 +981 @@
                              dt_3d * (                                         &
                                    tsc(2) * tend(k,j,i) + tsc(3) * tu_m(k,j,i) &
-                                 - tsc(4) * ( p(k,j,i) - p(k,j,i-1) ) * ddx    &
                                      ) -                                       &
                              tsc(5) * rdf(k) * ( u(k,j,i) - ug(k) )
@@ -1050 +1047 @@
                              dt_3d * (                                         &
                                    tsc(2) * tend(k,j,i) + tsc(3) * tv_m(k,j,i) &
-                                 - tsc(4) * ( p(k,j,i) - p(k,j-1,i) ) * ddy    &
                                      ) -                                       &
                              tsc(5) * rdf(k) * ( v(k,j,i) - vg(k) )
@@ -1116 +1112 @@
                           dt_3d * (                                         &
                                 tsc(2) * tend(k,j,i) + tsc(3) * tw_m(k,j,i) &
-                           - tsc(4) * ( p(k+1,j,i) - p(k,j,i) ) * ddzu(k+1) &
                                   ) -                                       &
                           tsc(5) * rdf(k) * w(k,j,i)
@@ -1525 +1520 @@
                           dt_3d * (                                            &
                                    tsc(2) * tend(k,j,i) + tsc(3) * tu_m(k,j,i) &
-                                 - tsc(4) * ( p(k,j,i) - p(k,j,i-1) ) * ddx    &
                                   ) -                                          &
                           tsc(5) * rdf(k) * ( u(k,j,i) - ug(k) )
@@ -1617 +1611 @@
                           dt_3d * (                                            &
                                    tsc(2) * tend(k,j,i) + tsc(3) * tv_m(k,j,i) &
-                                 - tsc(4) * ( p(k,j,i) - p(k,j-1,i) ) * ddy    &
                                   ) -                                          &
                           tsc(5) * rdf(k) * ( v(k,j,i) - vg(k) )
@@ -1706 +1699 @@
                           dt_3d * (                                            &
                                    tsc(2) * tend(k,j,i) + tsc(3) * tw_m(k,j,i) &
-                              - tsc(4) * ( p(k+1,j,i) - p(k,j,i) ) * ddzu(k+1) &
                                   ) -                                          &
                           tsc(5) * rdf(k) * w(k,j,i)