Changeset 106 for palm/trunk/SOURCE/flow_statistics.f90

Timestamp:

Aug 16, 2007 2:30:26 PM (17 years ago)

Author:

raasch

Message:

preliminary update of bugfixes and extensions for non-cyclic BCs

File:

• palm/trunk/SOURCE/flow_statistics.f90 (modified) (7 diffs)

palm/trunk/SOURCE/flow_statistics.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! Prescribed momentum fluxes at the top surface are used
+! Prescribed momentum fluxes at the top surface are used,
+! profiles for w*p* and w"e are calculated
 !
 ! Former revisions:
@@ -620 +621 @@
 !
 !--    Divergence of vertical flux of resolved scale energy and pressure
-!--    fluctuations. First calculate the products, then the divergence.
+!--    fluctuations as well as flux of pressure fluctuation itself (68).
+!--    First calculate the products, then the divergence.
 !--    Calculation is time consuming. Do it only, if profiles shall be plotted.
-       IF ( hom(nzb+1,2,55,0) /= 0.0 )  THEN
+       IF ( hom(nzb+1,2,55,0) /= 0.0  .OR.  hom(nzb+1,2,68,0) /= 0.0 )  THEN
 
           sums_ll = 0.0  ! local array
@@ -654 +656 @@
              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)
+             sums_l(k,68,tn) = sums_ll(k,2)
           ENDDO
           sums_l(nzb,55,tn) = sums_l(nzb+1,55,tn)
           sums_l(nzb,56,tn) = sums_l(nzb+1,56,tn)
-
-       ENDIF
-
-!
-!--    Divergence of vertical flux of SGS TKE
-       IF ( hom(nzb+1,2,57,0) /= 0.0 )  THEN
+          sums_l(nzb,68,tn) = 0.0    ! because w* = 0 at nzb
+
+       ENDIF
+
+!
+!--    Divergence of vertical flux of SGS TKE and the flux itself (69)
+       IF ( hom(nzb+1,2,57,0) /= 0.0  .OR.  hom(nzb+1,2,69,0) /= 0.0 )  THEN
 
           !$OMP DO
@@ -669 +673 @@
                 DO  k = nzb_s_outer(j,i)+1, nzt
 
-                   sums_l(k,57,tn) = sums_l(k,57,tn) + (                       &
+                   sums_l(k,57,tn) = sums_l(k,57,tn) - 0.5 * (                 &
                    (km(k,j,i)+km(k+1,j,i)) * (e(k+1,j,i)-e(k,j,i)) * ddzu(k+1) &
                  - (km(k-1,j,i)+km(k,j,i)) * (e(k,j,i)-e(k-1,j,i)) * ddzu(k)   &
-                                                 ) * ddzw(k)
+                                                             ) * ddzw(k)
+
+                   sums_l(k,69,tn) = sums_l(k,69,tn) - 0.5 * (                 &
+                   (km(k,j,i)+km(k+1,j,i)) * (e(k+1,j,i)-e(k,j,i)) * ddzu(k+1) &
+                                                              )
 
                 ENDDO
@@ -678 +686 @@
           ENDDO
           sums_l(nzb,57,tn) = sums_l(nzb+1,57,tn)
+          sums_l(nzb,69,tn) = sums_l(nzb+1,69,tn)
 
        ENDIF
@@ -832 +841 @@
        hom(:,1,55,sr) = sums(:,55)     ! w*u*u*/dz
        hom(:,1,56,sr) = sums(:,56)     ! w*p*/dz
-       hom(:,1,57,sr) = sums(:,57)     ! w"e/dz
+       hom(:,1,57,sr) = sums(:,57)     ! ( w"e + w"p"/rho )/dz
        hom(:,1,58,sr) = sums(:,58)     ! u"pt"
        hom(:,1,59,sr) = sums(:,59)     ! u*pt*
@@ -843 +852 @@
        hom(:,1,66,sr) = sums(:,66)     ! w*sa*
        hom(:,1,67,sr) = sums(:,65) + sums(:,66)    ! wsa
+       hom(:,1,68,sr) = sums(:,68)     ! w*p*
+       hom(:,1,69,sr) = sums(:,69)     ! w"e + w"p"/rho
 
        hom(:,1,pr_palm-1,sr) = sums(:,pr_palm-1)