Changeset 1652 for palm/trunk/SOURCE

Timestamp:

Sep 17, 2015 8:12:24 AM (9 years ago)

Author:

raasch

Message:

bugfix in calculation of vertical profile of energy production by turbulent transport of TKE
bugfix for output of configuration settings (mbuild)
output of login init command (mbuild)

File:

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

palm/trunk/SOURCE/flow_statistics.f90

@@ -21 +21 @@
 ! Current revisions:
 ! -----------------
-! 
+! bugfix in calculation of energy production by turbulent transport of TKE
 ! 
 ! Former revisions:
@@ -993 +993 @@
 
                    sums_ll(k,1) = sums_ll(k,1) + 0.5_wp * w(k,j,i) * (         &
-                  ( 0.25_wp * ( u(k,j,i)+u(k+1,j,i)+u(k,j,i+1)+u(k+1,j,i+1)    &
-                            - 0.5_wp * ( hom(k,1,1,sr) + hom(k+1,1,1,sr) )     &
-                              ) )**2                                           &
-                + ( 0.25_wp * ( v(k,j,i)+v(k+1,j,i)+v(k,j+1,i)+v(k+1,j+1,i)    &
-                            - 0.5_wp * ( hom(k,1,2,sr) + hom(k+1,1,2,sr) )     &
-                              ) )**2                                           &
+                  ( 0.25_wp * ( u(k,j,i)+u(k+1,j,i)+u(k,j,i+1)+u(k+1,j,i+1) )  &
+                            - 0.5_wp * ( hom(k,1,1,sr) + hom(k+1,1,1,sr) ) )**2&
+                + ( 0.25_wp * ( v(k,j,i)+v(k+1,j,i)+v(k,j+1,i)+v(k+1,j+1,i) )  &
+                            - 0.5_wp * ( hom(k,1,2,sr) + hom(k+1,1,2,sr) ) )**2&                                           &
                 + w(k,j,i)**2                                        )
 
@@ -2867 +2865 @@
 
                    sums_ll(k,1) = sums_ll(k,1) + 0.5_wp * w(k,j,i) * (         &
-                  ( 0.25_wp * ( u(k,j,i)+u(k+1,j,i)+u(k,j,i+1)+u(k+1,j,i+1)    &
-                              - 0.5_wp * ( hom(k,1,1,sr) + hom(k+1,1,1,sr) )   &
-                              ) )**2                                           &
-                + ( 0.25_wp * ( v(k,j,i)+v(k+1,j,i)+v(k,j+1,i)+v(k+1,j+1,i)    &
-                              - 0.5_wp * ( hom(k,1,2,sr) + hom(k+1,1,2,sr) )   &
-                              ) )**2                                           &
+                  ( 0.25_wp * ( u(k,j,i)+u(k+1,j,i)+u(k,j,i+1)+u(k+1,j,i+1) )  &
+                            - 0.5_wp * ( hom(k,1,1,sr) + hom(k+1,1,1,sr) ) )**2&                                           &
+                + ( 0.25_wp * ( v(k,j,i)+v(k+1,j,i)+v(k,j+1,i)+v(k+1,j+1,i) )  &
+                            - 0.5_wp * ( hom(k,1,2,sr) + hom(k+1,1,2,sr) ) )**2&                                          &
                 + w(k,j,i)**2                                        )