Changeset 2853 for palm

Timestamp:

Mar 5, 2018 2:44:20 PM (7 years ago)

Author:

suehring

Message:

Bugfix in init log-law correction

File:

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

palm/trunk/SOURCE/pmc_interface_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix in init log-law correction.
+! 
+! 2841 2018-02-27 15:02:57Z knoop
 ! Bugfix: wrong placement of include 'mpif.h' corrected
 ! 
@@ -1578 +1581 @@
        IMPLICIT NONE
 
-       INTEGER(iwp) ::  direction      !<  Wall normal index: 1=k, 2=j, 3=i.
+       INTEGER(iwp) ::  direction      !< Wall normal index: 1=k, 2=j, 3=i.
+       INTEGER(iwp) ::  dum            !< dummy value for reduce operation
        INTEGER(iwp) ::  i              !<
        INTEGER(iwp) ::  icorr          !<
-       INTEGER(iwp) ::  inc            !<  Wall outward-normal index increment -1
+       INTEGER(iwp) ::  ierr           !< MPI status
+       INTEGER(iwp) ::  inc            !< Wall outward-normal index increment -1
                                        !< or 1, for direction=1, inc=1 always
        INTEGER(iwp) ::  iw             !<
@@ -1716 +1721 @@
           nzt_topo_nestbc_n = nzt_topo_nestbc_n + 1
        ENDIF
+
+#if defined( __parallel )
+!
+!--       Determine global topography-top index along child boundary. 
+          dum = nzb
+          CALL MPI_ALLREDUCE( nzt_topo_nestbc_l, dum, 1, MPI_INTEGER,          &
+                              MPI_MAX, comm1dy, ierr )
+          nzt_topo_nestbc_l = dum
+
+          dum = nzb
+          CALL MPI_ALLREDUCE( nzt_topo_nestbc_r, dum, 1, MPI_INTEGER,          &
+                              MPI_MAX, comm1dy, ierr )
+          nzt_topo_nestbc_r = dum
+
+          dum = nzb
+          CALL MPI_ALLREDUCE( nzt_topo_nestbc_n, dum, 1, MPI_INTEGER,          &
+                              MPI_MAX, comm1dx, ierr )
+          nzt_topo_nestbc_n = dum
+
+          dum = nzb
+          CALL MPI_ALLREDUCE( nzt_topo_nestbc_s, dum, 1, MPI_INTEGER,          &
+                              MPI_MAX, comm1dx, ierr )
+          nzt_topo_nestbc_s = dum
+#endif
 !
 !--    Then determine the maximum number of near-wall nodes per wall point based