Changeset 2617

Timestamp:

Nov 16, 2017 12:47:24 PM (6 years ago)

Author:

suehring

Message:

Avoid that the mean reference state in buoyancy term does not become zero

File:

• palm/trunk/SOURCE/time_integration.f90 (modified) (2 diffs)

palm/trunk/SOURCE/time_integration.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix, assure that the reference state does not become zero.
+! 
+! 2563 2017-10-19 15:36:10Z Giersch
 ! Variable wind_turbine moved to module control_parameters
 ! 
@@ -549 +552 @@
                 ref_state(:)  = hom(:,1,44,0)
              ENDIF
-
+!
+!--          Assure that ref_state does not become zero at any level
+!--          ( might be the case if a vertical level is completely occupied 
+!--            with topography ).
+             ref_state = MERGE( MAXVAL(ref_state), ref_state,                  &
+                                ref_state == 0.0_wp )
           ENDIF