Changeset 4022 for palm/trunk/SOURCE/parin.f90

Timestamp:

Jun 12, 2019 11:52:39 AM (5 years ago)

Author:

suehring

Message:

Synthetic turbulence generator: Revise bias correction of imposed perturbations (correction via volume flow can create instabilities in case the mean volume flow is close to zero); Introduce lower limits in calculation of coefficient matrix, else the calculation may become numerically unstable; Impose perturbations every timestep, even though no new set of perturbations is generated in case dt_stg_call /= dt_3d; Implement a gradual decrease of Reynolds stress and length scales above ABL height (within 1 length scale above ABL depth to 1/10) rather than an abrupt decline; Bugfix in non-nested case: use ABL height for parametrized turbulence; Offline nesting: Rename subroutine for ABL height calculation and make it public; Change top boundary condition for pressure back to Neumann

File:

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

palm/trunk/SOURCE/parin.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Change default top boundary condition for pressure to Neumann in offline
+! nesting case
+! 
+! 4017 2019-06-06 12:16:46Z schwenkel
 ! Introduce alternative switch for debug output during timestepping
 ! 
@@ -983 +987 @@
              bc_pt_t  = 'nesting_offline'
              bc_q_t   = 'nesting_offline'
-           !  bc_s_t   = 'nesting_offline'  ! scalar boundary condition is not clear
+           !  bc_s_t   = 'nesting_offline'  ! scalar boundary condition is not clear yet
            !  bc_cs_t  = 'nesting_offline'  ! same for chemical species
-!
-!--          For the pressure set Dirichlet conditions, in contrast to the 
-!--          self nesting. This gives less oscilations within the 
-!--          free atmosphere since the pressure solver has more degrees of 
-!--          freedom. In constrast to the self nesting, this might be 
-!--          justified since the top boundary is far away from the domain
-!--          of interest. 
-             bc_p_t   = 'dirichlet' !'neumann'
+             bc_p_t   = 'neumann'
           ENDIF