Changeset 678 for palm/trunk/SOURCE

Timestamp:

Feb 2, 2011 2:31:56 PM (14 years ago)

Author:

raasch

Message:

New:
---

further adjustments on Tsubame and concerning openMP usage
(mrun, mbuild, subjob)

Changed:

---

Errors:

---

Bugfix in calculation of divergence of vertical flux of resolved scale
energy, pressure fluctuations, and flux of pressure fluctuation itself
(flow statistics)

Bugfix: module pegrid was missing. (user_statistics)

Location:

palm/trunk/SOURCE

Files:

• flow_statistics.f90 (modified) (3 diffs)
• user_statistics.f90 (modified) (2 diffs)

palm/trunk/SOURCE/flow_statistics.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
+! Bugfix in calculation of divergence of vertical flux of resolved scale
+! energy, pressure fluctuations, and flux of pressure fluctuation itself
 !
 ! Former revisions:
@@ -802 +804 @@
                    sums_ll(k,1) = sums_ll(k,1) + 0.5 * w(k,j,i) * (        &
                   ( 0.25 * ( u(k,j,i)+u(k+1,j,i)+u(k,j,i+1)+u(k+1,j,i+1)   &
-                           - 2.0 * ( hom(k,1,1,sr) + hom(k+1,1,1,sr) )     &
+                           - 0.5 * ( hom(k,1,1,sr) + hom(k+1,1,1,sr) )     &
                            ) )**2                                          &
                 + ( 0.25 * ( v(k,j,i)+v(k+1,j,i)+v(k,j+1,i)+v(k+1,j+1,i)   &
-                           - 2.0 * ( hom(k,1,2,sr) + hom(k+1,1,2,sr) )     &
+                           - 0.5 * ( hom(k,1,2,sr) + hom(k+1,1,2,sr) )     &
                            ) )**2                                          &
                    + w(k,j,i)**2                                  )
@@ -820 +822 @@
           sums_ll(nzt+1,2) = 0.0
 
-          DO  k = nzb_s_inner(j,i)+1, nzt
+          DO  k = nzb+1, nzt
              sums_l(k,55,tn) = ( sums_ll(k,1) - sums_ll(k-1,1) ) * ddzw(k)
              sums_l(k,56,tn) = ( sums_ll(k,2) - sums_ll(k-1,2) ) * ddzw(k)

palm/trunk/SOURCE/user_statistics.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix: module pegrid was missing
 !
 ! Former revisions:
@@ -35 +35 @@
     USE indices
     USE netcdf_control
+    USE pegrid
     USE statistics
     USE user