Changeset 3022 for palm/trunk/SOURCE/prognostic_equations.f90

Timestamp:

May 18, 2018 11:12:35 AM (6 years ago)

Author:

suehring

Message:

Revise recent bugfix in nested runs at left and south boundary; bugfix in advection of u in case of OpenMP parallelization; bugfix in plant transpiration

File:

• palm/trunk/SOURCE/prognostic_equations.f90 (modified) (4 diffs)

palm/trunk/SOURCE/prognostic_equations.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Revise recent bugfix for nesting
+! 
+! 3021 2018-05-16 08:14:20Z maronga
 ! Bugfix in IF clause for nesting
 ! 
@@ -298 +301 @@
                intermediate_timestep_count_max, large_scale_forcing,           &
                large_scale_subsidence, microphysics_morrison,                  &
-               microphysics_seifert, microphysics_sat_adjust, nest_bound_s,    &
-               neutral, nudging,                                               &
+               microphysics_seifert, microphysics_sat_adjust, neutral, nudging,&
                ocean, outflow_l, outflow_s, passive_scalar, plant_canopy,      &
                prho_reference, prho_reference,                                 &
@@ -527 +529 @@
        DO  j = nys, nyn
 !
-!--       Tendency terms for u-velocity component
-          IF ( .NOT. outflow_l  .OR.  i > nxl )  THEN
+!--       Tendency terms for u-velocity component. Please note, in case of 
+!--       non-cyclic boundary conditions the grid point i=0 is excluded from
+!--       the prognostic equations for the u-component.    
+          IF ( i >= nxlu )  THEN
 
              tend(:,j,i) = 0.0_wp
@@ -598 +602 @@
 
           ENDIF
-
-!
-!--       Tendency terms for v-velocity component
-          IF ( (.NOT. outflow_s .AND. .NOT. nest_bound_s ) .OR.  j > nys )  THEN
+!
+!--       Tendency terms for v-velocity component. Please note, in case of 
+!--       non-cyclic boundary conditions the grid point j=0 is excluded from
+!--       the prognostic equations for the v-component. !--       
+          IF ( j >= nysv )  THEN
 
              tend(:,j,i) = 0.0_wp