Changeset 2478 for palm/trunk/SOURCE/production_e.f90

Timestamp:

Sep 18, 2017 1:37:24 PM (7 years ago)

Author:

suehring

Message:

Bugfixes concerning top fluxes and TKE production

File:

• palm/trunk/SOURCE/production_e.f90 (modified) (6 diffs)

palm/trunk/SOURCE/production_e.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix, consider case where no constant-flux layer and no surfaces fluxes 
+! are used
+! 
+! 2329 2017-08-03 14:24:56Z knoop
 ! Bugfix: added division by density as kinematic fluxes are needed
 ! 
@@ -432 +436 @@
              ENDDO
 
-          ELSEIF ( use_surface_fluxes )  THEN
+          ELSE
 
              DO  j = nys, nyn
@@ -1182 +1186 @@
           ENDDO
 
-       ELSEIF ( use_surface_fluxes )  THEN
+       ELSE
 !
 !--       Calculate TKE production by shear. Here, no additional
@@ -1235 +1239 @@
 !--                the bottom and top surface layer
                    DO  k = nzb+1, nzt
+ 
                       tend(k,j,i) = tend(k,j,i) +                              &
                                     kh(k,j,i) * g / rho_reference *            &
@@ -1241 +1246 @@
                                       MERGE( 1.0_wp, 0.0_wp,                   &
                                              BTEST( wall_flags_0(k,j,i), 0 )   &
-                                           )                                   
+                                           )                   
                    ENDDO
 
@@ -1323 +1328 @@
                                       MERGE( 1.0_wp, 0.0_wp,                   &
                                              BTEST( wall_flags_0(k,j,i), 0 )   &
-                                           )
+                                           )           
                    ENDDO