Ignore:
Timestamp:
Mar 5, 2018 2:44:20 PM (4 years ago)
Author:
suehring
Message:

Bugfix in init log-law correction

File:
1 edited

Legend:

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

    r2841 r2853  
    2525! -----------------
    2626! $Id$
     27! Bugfix in init log-law correction.
     28!
     29! 2841 2018-02-27 15:02:57Z knoop
    2730! Bugfix: wrong placement of include 'mpif.h' corrected
    2831!
     
    15781581       IMPLICIT NONE
    15791582
    1580        INTEGER(iwp) ::  direction      !<  Wall normal index: 1=k, 2=j, 3=i.
     1583       INTEGER(iwp) ::  direction      !< Wall normal index: 1=k, 2=j, 3=i.
     1584       INTEGER(iwp) ::  dum            !< dummy value for reduce operation
    15811585       INTEGER(iwp) ::  i              !<
    15821586       INTEGER(iwp) ::  icorr          !<
    1583        INTEGER(iwp) ::  inc            !<  Wall outward-normal index increment -1
     1587       INTEGER(iwp) ::  ierr           !< MPI status
     1588       INTEGER(iwp) ::  inc            !< Wall outward-normal index increment -1
    15841589                                       !< or 1, for direction=1, inc=1 always
    15851590       INTEGER(iwp) ::  iw             !<
     
    17161721          nzt_topo_nestbc_n = nzt_topo_nestbc_n + 1
    17171722       ENDIF
     1723
     1724#if defined( __parallel )
     1725!
     1726!--       Determine global topography-top index along child boundary.
     1727          dum = nzb
     1728          CALL MPI_ALLREDUCE( nzt_topo_nestbc_l, dum, 1, MPI_INTEGER,          &
     1729                              MPI_MAX, comm1dy, ierr )
     1730          nzt_topo_nestbc_l = dum
     1731
     1732          dum = nzb
     1733          CALL MPI_ALLREDUCE( nzt_topo_nestbc_r, dum, 1, MPI_INTEGER,          &
     1734                              MPI_MAX, comm1dy, ierr )
     1735          nzt_topo_nestbc_r = dum
     1736
     1737          dum = nzb
     1738          CALL MPI_ALLREDUCE( nzt_topo_nestbc_n, dum, 1, MPI_INTEGER,          &
     1739                              MPI_MAX, comm1dx, ierr )
     1740          nzt_topo_nestbc_n = dum
     1741
     1742          dum = nzb
     1743          CALL MPI_ALLREDUCE( nzt_topo_nestbc_s, dum, 1, MPI_INTEGER,          &
     1744                              MPI_MAX, comm1dx, ierr )
     1745          nzt_topo_nestbc_s = dum
     1746#endif
    17181747!
    17191748!--    Then determine the maximum number of near-wall nodes per wall point based
Note: See TracChangeset for help on using the changeset viewer.