Changeset 2617 for palm


Ignore:
Timestamp:
Nov 16, 2017 12:47:24 PM (7 years ago)
Author:
suehring
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/time_integration.f90

    r2563 r2617  
    2525! -----------------
    2626! $Id$
     27! Bugfix, assure that the reference state does not become zero.
     28!
     29! 2563 2017-10-19 15:36:10Z Giersch
    2730! Variable wind_turbine moved to module control_parameters
    2831!
     
    549552                ref_state(:)  = hom(:,1,44,0)
    550553             ENDIF
    551 
     554!
     555!--          Assure that ref_state does not become zero at any level
     556!--          ( might be the case if a vertical level is completely occupied
     557!--            with topography ).
     558             ref_state = MERGE( MAXVAL(ref_state), ref_state,                  &
     559                                ref_state == 0.0_wp )
    552560          ENDIF
    553561
Note: See TracChangeset for help on using the changeset viewer.