Changeset 4855 for palm/trunk/SOURCE/pres.f90

Timestamp:

Jan 25, 2021 12:30:54 PM (4 years ago)

Author:

raasch

Message:

bugfix: mean w removal not applied to ghost points of the total domain in case of non-cyclic setups (pres), bugfix for correct identification of indices of extreme values in case of non-cyclic boundary conditions

File:

• palm/trunk/SOURCE/pres.f90 (modified) (3 diffs)

palm/trunk/SOURCE/pres.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! bugfix: mean w removal not applied to ghost points of the total domain in case of non-cyclic
+! setups
+!
+! 4828 2021-01-05 11:21:41Z Giersch
 ! Bugfix in OpenMP directives - intel compiler do not allow reduction operations on array elements
 !
@@ -357 +361 @@
           w_l(k) = w_l(k) / ngp_2dh_outer(k,0)
        ENDDO
-       DO  i = nxlg, nxrg
-          DO  j = nysg, nyng
+       DO  i = nxl, nxr
+          DO  j = nys, nyn
              DO  k = nzb+1, nzt
                 w(k,j,i) = w(k,j,i) - w_l(k)                                                       &
@@ -365 +369 @@
           ENDDO
        ENDDO
+!
+!--    Instead of running the above loop over ghost points, they are set via exchange_horiz,
+!--    in order to correctly consider non-cyclic boundary conditions, where ghost boundaries
+!--    of the total domain must not be set. Otherwise w may continuously increase/decrease
+!--    at these points.
+       CALL exchange_horiz( w, nbgp )
     ENDIF